Kelvin-Helmholtz Instability Research Articles

Effects of Kelvin–Helmholtz instability on the material mixing in the double-cone ignition

The Kelvin–Helmholtz instability (KHI) occurs on the interface of gold cones and embedded fuels for fusion schemes with gold cones. The development of KHI on the inner surface of gold cones in the double-cone ignition scheme is investigated with two-dimensional radiation hydrodynamic simulations. It has been found that the colliding high-density fuel plasma between the tips of the two cones is spatiotemporally separated from the mixed gold ions from the inner surface of the gold cones due to the KHI. Furthermore, it is found that fuel layers coated on the inner surface of the cones can effectively mitigate the energy loss in the compression process. These results could provide a reference for fast ignition schemes with gold cones.

Study on the Performance of a Surface with Coupled Wettability Difference and Convex-Stripe Array for Improved Air Layer Stability.

The existence of an air layer reduces friction drag on superhydrophobic surfaces. Therefore, improving the air layer stability of superhydrophobic surfaces holds immense significance in reducing both energy consumption and environmental pollution caused by friction drag. Based on the properties of mathematical discretization and the contact angle hysteresis generated by the wettability difference, a surface coupled with a wettability difference treatment and a convex-stripe array is developed by laser engraving and fluorine modification, and its performance in improving the air layer stability is experimentally studied in a von Kármán swirling flow field. The results show that the destabilization of the air layer is mainly caused by the Kelvin-Helmholtz instability, which is triggered by the density difference between gas and liquid, as well as the tangential velocity difference between gas and liquid. When the air layer is relatively thin, tangential wave destabilization occurs, whereas for larger thicknesses, the destabilization mode is coupled wave destabilization. The maximum Reynolds number that keeps the air layer fully covering the surface of the rotating disk (with drag reduction performance) during the disk rotation process is defined as the critical Reynolds number (Rec), which is 1.62 × 105 for the uniform superhydrophobic surface and 3.24 × 105 for the superhydrophobic surface with a convex stripe on the outermost ring (SCSSP). Individual treatments of wettability difference and a convex-stripe array on the SCSSP further improve the air layer stability, but Rec remains at 3.24 × 105. Finally, the coupling of the wettability difference treatment with a convex-stripe array significantly improves the air layer stability, resulting in an increase of Rec to 4.05 × 105, and the drag reduction rate stably maintained around 30%.

Growth and Breakdown of Kelvin–Helmholtz Billows in the Stable Atmospheric Boundary Layer

The development and breakdown of Kelvin–Helmholtz (KH) waves (billows) in the stable atmospheric boundary layer (SABL) and their impact on vertical transport of momentum and scalars have been examined utilizing large eddy simulations. These simulations are initialized with a vertically uniform geostrophic wind and a constant potential temperature lapse rate. An Ekman type of boundary layer develops, and an inflection point forms in the SABL, which triggers the KH instability (KHI). KHI develops with the kinetic energy (KE) in the KH billows growing exponentially with time. The subsequent onset of secondary shear instability along S-shaped braids leads to the turbulent breakdown of the KH billow cores and braids. The frictional ground surface tends to slow down the growth of KE near the surface, reduce the KH billow core depth, and likely suppress other types of secondary instability. KH billows induce substantial down-gradient transport of momentum and sensible heat, which can be further enhanced by the onset of secondary shear instability. Although the KHI-induced strong transport only lasts for around 10–20 min, it reduces vertical shear and stratification in the SABL, enhances surface winds, and results in a 2–3-fold increase in the SABL depth.

Vertical shear instability in two-moment radiation-hydrodynamical simulations of irradiated protoplanetary disks

Context. The vertical shear instability (VSI) is a hydrodynamical instability predicted to produce turbulence in magnetically inactive regions of protoplanetary disks. The regions in which this instability can occur and the physical phenomena leading to its saturation are a current matter of research. Aims. We explore the secondary instabilities triggered by the nonlinear evolution of the VSI and their role in its saturation. We also expand on previous investigations on stability regions by considering temperature stratifications enforced by stellar irradiation and radiative cooling, and including the effects of dust-gas collisions and molecular line emission. Methods. We modeled the gas-dust mixture in a circumstellar disk around a T Tauri star by means of high-resolution axisymmetric radiation-hydrodynamical simulations including stellar irradiation with frequency-dependent opacities, considering different degrees of depletion of small dust grains. Results. The flow pattern produced by the interplay of the axisymmetric VSI modes and the baroclinic torque forms bands of nearly uniform specific angular momentum. In the high-shear regions in between these bands, the Kelvin–Helmholtz instability (KHI) is triggered. A third instability mechanism, consisting of an amplification of eddies by baroclinic torques, forms meridional vortices with Mach numbers up to ∼0.4. Our stability analysis suggests that protoplanetary disks can be VSI-unstable in surface layers up to tens of au for reasonably high gas emissivities. Conclusions. The significant transfer of kinetic energy to small-scale eddies produced by the KHI and possibly even the baroclinic acceleration of eddies limit the maximum energy of the VSI modes, likely leading to the saturation of the VSI. Depending on the gas molecular composition, the VSI can operate at surface layers even in regions where the midplane is stable. This picture is consistent with current observations of disks showing thin midplane millimeter-sized dust layers while appearing vertically extended in optical and near-infrared wavelengths.

Evolution of a shock-impacted reactive liquid fuel droplet with evaporation effects: A numerical study

We describe detailed numerical simulations of a liquid fuel droplet impacted by a Mach 5 shock wave, considering the effects of chemical reactions and phase change due to evaporation. In our baseline case, a 5μm, n-Dodecane fuel droplet is preheated to 460K, and surrounded by preheated O2 gas at 700K. The simulations investigate the low Ohnesorge number regime (Oh<0.1), where viscous effects are insignificant and therefore not considered in this work. The fuel droplet undergoes significant deformation and morphological changes following shock impingement, as the droplet surface becomes unstable to the Kelvin-Helmholtz instability. The droplet core is also observed to eject a thin sheet near the equatorial plane, which is then stretched by the high-speed post-shock gas flow, affecting the late-time behavior. We find the observed dominant modes associated with the Kelvin–Helmholtz instability are in reasonable agreement with inviscid linear theory (Chandrasekhar, 1981; Jepsen et al., 2006), when the local conditions at the droplet surface are considered. Furthermore, an evaporation-induced Stefan flow is established which blows off the hot post-shock gases surrounding the droplet, leading to droplet cooling. As the fuel vapors react, a diffusion flame is formed on the droplet-windward side, leading to intense droplet heating and enhanced vapor production in this region. We investigated the effect of the Damkohler number on droplet evolution by varying the fuel reactivity, and found the flame thickness decreased with increasing reactivity in agreement with trends predicted by laminar diffusion flame theory. At the highest reactivity, secondary burning of fuel vapors was observed in the droplet wake, while the resulting flame eventually reattached to the droplet surface. Our results show significant spatial inhomogeneities are present in the droplet flowfield in all the cases investigated, which must be considered in the development of reduced order point-particle models for system-level simulations of detonation engines.

Analysis of Atmospheric Elements in Near Space Based on Meteorological-Rocket Soundings over the East China Sea

As an important means of in situ detection in near space, meteorological rockets can provide a high-precision distribution analysis of atmospheric elements. However, there are currently few studies on the principles of meteorological-rocket detection and the application of rocket-sounding data. The purpose of this paper is to fill this gap by providing a detailed introduction to the detection principle of a meteorological rocket launched in the East China Sea in November 2022. Moreover, empirical models, satellite data, and reanalysis data were selected for comparison and verification with the rocket-sounding data. Furthermore, the accuracy of these widely used datasets was studied based on the rocket-sounding data in the near space over the East China Sea. Additionally, gravity-wave power–frequency spectra were extracted using the maximum entropy method from both the rocket-sounding data and the remote-sensing data. Furthermore, the relationship between gravity waves and Kelvin–Helmholtz instability (KHI) was investigated by analyzing the gravity-wave energy and the Richardson number. The research findings indicate that among the remote-sensing data describing the atmospheric environment over the launch site, the COSMIC occultation data is more accurate compared with the SABER data. The wind-field distribution derived from rocket detection is consistent with the Modern-Era Retrospective analysis for Research and Applications (MERRA) reanalysis data, while also providing a more detailed description of the wind field. The main wavelengths of gravity waves extracted from rocket-sounding data are consistently smaller than those obtained from satellite remote-sensing data, indicating that rocket sounding is capable of capturing more intricate structures of gravity waves. The good correspondence between the peaks of gravity-wave energy and the regions where KHI occurs indicates that there is a strong interaction between gravity waves and KHI in the middle atmosphere.

Identification of slow waves in the evolution of KHI near the Venusian ionopause

The MHD slow waves in the evolution of Kelvin–Helmholtz instability (KHI) near the Venusian ionopause has been identified by using the magnetohydrodynamics equations. The diagnosis of slow waves utilizes two criteria: the phase speed and the ratio of density oscillation to the velocity oscillation. With the speed of sound greater than Alfvén speed, the phase velocity of slow waves along the initial magnetic field is approximately equal to Alfvén speed. For slow waves, the ratio of the density oscillation to the velocity oscillation along the initial magnetic field is independent of the angle between wave vector and the magnetic field. The data from the side of low-density were analyzed around the KH vortex. The numerical simulation results show that there are slow waves in the induced magnetosphere near the magnetopause of Venus. We also note that during the nonlinear growth stage, the parallel oscillations contribute approximately 76% ∼ 93% to the wave energy. This work will provide more clues that the KHI is a possible source of slow waves observed near Venus-like planets.

Heat transfer across an air curtain for heat confinement in road tunnels: Influence of the heat source

Experiments were carried out to study the influence of the heat source using double-flow double-jet curtains (DS-TJ) to confine heat in fire safety problems. For this purpose, a small-scale road tunnel facility in 1:34 was built. The heat source used is an electric resistance. Simultaneous two-dimensional temperature and velocity measurements were made across the DS-TJ air curtain and in its immediate vicinity using a fine thermocouple and a 2D laser Doppler anemometer (LDA), respectively. Four cases were studied, by modifying the air curtain velocities at the nozzle exits and the heat source power. It is noted that hot-cold jets and heat source play a role in the curtain heat confinement by modifying the development of the two jets. Both the entrainment of the upward flow rising from the heat source as well as the Kelvin-Helmholtz instabilities contribute to this. On the other hand, the nozzle exit velocities determine the behavior of the Reynolds stress profiles as well as the production of the velocities fluctuations across the curtain, showing a non-correlating effect on the hot side, so that the profiles take values around zero, suggesting a heat source influence. In contrast in the cold side negative values indicate energy-draining behavior. The intensity of thermal turbulence reveals that the hot jet has higher values than the cold jet, constituting a thermal barrier to turbulent heat flow, and all activity is relegated to the zone of interaction between the two jets. In fact, the production derived from temperature fluctuations occurs mainly in the interaction of hot-cold jets, and negative values are due to large turbulent eddies, such as Kelvin-Helmholtz instabilities and Görtler-type structures. The Richardson number shows that buoyancy forces prevail on both sides of the curtain (Ri > 0.1), especially on the hot side, due to the heat source. Lower Ri values within the curtain indicate the presence of a thermal barrier, attenuating heat transfer through the curtain and favoring thermal sealing.

PRIMARY BREAKUP INSTABILITY OF LIQUID JET IN CROSSFLOW

The liquid jet in crossflows (LJICF) has been analyzed using the compressible volume of fluid-Lagrangian particle tracking (VOF-LPT) coupled solver for the instabilities that result in the primary breakup. It is understood that the dominant force driving the instabilities changes with the Weber number and momentum flux ratio. The Kelvin- Helmholtz (KH) instability is found to be prevalent at low momentum flux ratio, whereas the Rayleigh-Taylor (RT) instability is dominant at higher values. In the present work, the instability causing the primary breakup is analyzed for a range of Weber numbers and momentum flux ratios, where the breakup is predominantly caused by either KH or RT instability. It is observed that the transition from KH waves to RT waves happens for the momentum flux ratio values ranging from 20 to 50. Also, the lower Weber number cases show the domination of long KH waves on the liquid jet column with negligible turbulence.

CHARACTERIZATION OF SEPARATED SHEAR LAYER OF ELEVATED JET IN CROSSFLOW AT LOW-VELOCITY RATIO

The present study reports an experimental investigation of the flow dynamics of separated shear layers of an elevated jet-in-crossflow (EJICF). The current work aims to enhance the understanding of the downwash of jet fluid into the stack-wake region, akin to plume downwash in chimney flows, particularly at a low-velocity ratio. One end of the square stack with the aspect ratio of 9.0 is kept on a flat plate while the jet emerges from the free-end open to crossflow. Flow visualization and particle image velocimetry (PIV) experiments are conducted in a recirculating, low-speed water tunnel for three Reynolds numbers, namely Re &amp;#61; 1000, 2000, and 4000, based on the side of the square stack (&lt;i&gt;D&lt;/i&gt;) and the freestream velocity (&lt;i&gt;U&lt;/i&gt;&lt;sub&gt;&amp;#8734;&lt;/sub&gt;). The study mainly focuses on the jet shear layer characteristics and the downstream wake near the free-end formed by the interaction of crossflow and the jet shear layer. The instantaneous flow behavior includes the analysis of the shedding of the jet shear layer in a symmetric plane (&lt;i&gt;z&lt;/i&gt; &amp;#61; 0) at a constant velocity ratio of &lt;i&gt;R&lt;/i&gt; &amp;#61; 0.5. On the windward side, the jet shear layer reveals the formation of small-scale vortices due to the Kelvin-Helmholtz (KH) instability for Re &amp;ge; 2000. Also, it is found that at a high Reynolds number (Re &amp;#61; 4000), an unstable vortex appears above the jet exit on the upstream side, mentioned as a hovering vortex. The mean and fluctuating flow behavior are investigated using streamline and contour plots of mean velocity and fluctuating components.

Magnetic field suppression characteristics in interaction process between shock wave and light gas cylinder

Based on ideal compressible magnetohydrodynamics (MHD) equations, the interface instabilities induced by the interaction between planar shock wave and the light gas (Helium) cylinder under the influence of the magnetic fields with different directions are investigated numerically by using the CTU(corner transport upwind)+CT (constrained transport) algorithm. The numerical results elucidate the evolution of flow field characteristics and wave structures with and without magnetic field. Moreover, we examine the influence of the magnetic field direction on a characteristic scales (including the length, height and width of the central axis of gas cylinder), as well as the volume compressibility. Then, the mechanism of the magnetic field direction affecting the interface instability is studied in depth by integrating the analyses of the circulation, energy, velocity and magnetic force distribution within the flow field. The core of this study, is to explore the suppression mechanism of interface instability by magnetic field force. The results show that the magnetic pressure plays a crucial role in driving vorticity away from the interface, thereby reducing its deposition on the density interface. Simultaneously, it adheres to the divided vortex layer, thereby effectively isolating the influence of Richtmyer-Meshkov instability on the interface. On the other hand, the magnetic tension adheres to the separated vortex layer, and its direction is opposite to that of the vorticity generated by the shear of interface velocity. This action effectively suppresses the Kelvin-Helmholtz instability and the rolling-up of vortices on the density interface. Additionally, under the control of a longitudinal magnetic field, the direction of magnetic tension is opposite to the direction of the central jet, effectively suppressing the development of Rayleigh-Taylor instability.

Numerical Simulation Study on the Dynamics of Bluff-Body Flames under Oxygen-Lean Conditions

As modern aeroengine combustors advance towards high temperatures, afterburners are inevitably affected by diminished oxygen content in incoming flows, thus affecting combustion efficiency, instability, and flammability limits. In this study, the dynamic combustion characteristics of V-shaped bluff body-stabilised diffusion flames were investigated using a large eddy simulation method with an oxygen mass fraction range of 14–23% and temperatures ranging from 900 to 1100 K. The results show the significant effects of oxygen content and inflow temperature on the flame/flow behaviours downstream of the bluff-body flame holder. In a separated shear layer, two distinct modes of flow/flame shedding are observed when varying the oxygen content and inflow temperature. The results show that BVK instability governs the far-field wake flow/flame features, whereas the oxygen concentration and temperature significantly affect their oscillation amplitudes. In addition, variations in the incoming oxygen content and temperature shift the axial position of the transition from KH instability to BVK instability. Finally, a spectral analysis is conducted to investigate the characteristics of pressure and heat release pulsations under different scenarios. This study highlights the importance of oxygen content on the combustion dynamics of bluff body-stabilised diffusion flames at various temperatures, which is essential for optimising combustion efficiency and stability in practical applications.

Demonstrating the Kelvin-Helmholtz Instability Using a Low-Cost Experimental Apparatus and Computational Fluid Dynamics Simulations

A Kelvin-Helmholtz instability is formed when two fluids of different densities exert a shear on one another at their interface when flowing in opposite directions. This paper presents a step-by-step guide for the design of a low-cost, small-scale, experimental tilt tube apparatus and a corresponding computational fluid dynamics (CFD) model that can be used to introduce the Kelvin-Helmholtz instability to undergraduate mechanical engineering students in several courses. A thermal-fluids laboratory course is taken by our fourth-year mechanical engineering students, and the overall variety of experiments has been limited by the cost of commercial teaching equipment. The tilt tube apparatus allows students to induce and record the Kelvin-Helmholtz instability, and no ongoing costs are involved in incorporating this experiment into the course. In our introductory CFD course, students perform CFD simulations as part of the design and analysis process. Developing a two-dimensional (2D) CFD model with two different fluids is well within their capabilities after completing initial software and simulation tutorial exercises and homework. Representative experiments were conducted with fresh water and salt water of different densities, and results showed that both the amplitude of the waves and the amount of time the instability was visible decreased with increasing salt water salinity. Results from a 2D CFD model developed in Ansys Fluent exhibited the same trends as the experimental data.

Swirling instability of coaxial liquid jet in gas surroundings

Linear instability analysis of an inviscid coaxial swirling jet is carried out by deriving an analytical dispersion relation of perturbation growth. The azimuthal Rankine vortex and the axial discontinuous velocity distribution are utilized as the jet basic flow. Due to the existence of double interfaces, the instability mechanisms of the coaxial swirling jet are much more complex than those of the single-layered swirling jet. The effects of control parameters (including the swirling ratio, the Weber number, the jet radius ratio, the velocity ratios between different fluids, and the azimuthal velocity jump at the inner interface) on the temporal instability of coaxial swirling jet with different azimuthal modes are studied. By comparing the growth rate of different azimuthal modes, the predominant mode that determines the jet breakup is identified. The results indicate that an increase in the swirling ratio, the Weber number, and the radius ratio can lead to predominant mode transition to larger azimuthal wavenumbers. The velocity ratio between the inner jet and the annular jet and that between the surrounding fluid and the annular jet mainly affect the axial Kelvin–Helmholtz (KH) instability. An enhancement of the KH instability leads to the jet breakup with smaller azimuthal wavenumbers. The azimuthal velocity jump affects the azimuthal KH instability, the centrifugal instability, and the Coriolis instability simultaneously, thus leading to a multiple influence on modes transition. The phase-diagram of the predominant modes is further given, showing that the relative importance between the centrifugal force and the interfacial tension plays a significant role on the transition of predominant modes.

Turbulence Transitions in Kelvin–Helmholtz Instability “Tube” and “Knot” Dynamics: Vorticity, Helicity, and Twist Waves

We address the sources and dynamics of vorticity and helicity and their relations in transitions to turbulence arising due to Kelvin–Helmholtz instability (KHI) “Tube” and “Knot” (T&amp;K) events. Such events are common in the atmosphere and oceans, and initial numerical simulations reveal that T&amp;K dynamics significantly accelerate turbulence transitions and enhance KHI peak and mean energy dissipation rates. KHI T&amp;K events arise where emerging KH billows exhibit varying wavelengths, phases, amplitudes, and/or discontinuities along their axes. As the KH billows intensify, these regions evolve roughly orthogonal billow cores and induced vortex tubes in close proximity. Their mutual advection as they intensify induces large-amplitude Kelvin vortex waves, or “twist waves”, that arise where locally uniform vortices are distorted by axial or radial advection. The twist waves propagate along, and fragment, the vortex tubes and billow cores, thus accounting for the emergence of helicity and the down-scale energy, enstrophy, and helicity fluxes within the turbulence inertial range. We describe the results of four direct numerical simulations (DNS) addressing KHI T&amp;K dynamics in large and idealized small domains. The large-domain vorticity fields reveal the character and diversity of KHI T&amp;K dynamics, the emergence of twist waves at larger and smaller scales, and their driving of turbulence transitions. Two small-domain DNS exhibit idealized KHI T&amp;K events arising from KH billows that are mis-aligned and that exhibit phase variability along their axes. A third examines the interactions of two vortex tubes in close proximity. These reveal that twist waves drive the character and evolutions of the vorticity and helicity fields.

Analysis of the E-Region Low-Altitude Quasi-Periodic Event in Low-Latitudes of China

The low-altitude quasi-periodic (LQP) E-region echoes observed with the upgraded Hainan COherent scatter Phased Array Radar (HCOPAR) located in Hainan Island of China (19.5°N, 109.1°E) are presented. With the interferometry technique, interferometry equations were developed to investigate the structures of the irregularities. The results show that the plasma structures of the LQP echoes with a negative slope mainly drifted southwestward almost horizontally, but they also drifted southeastward in some striations. According to the movement of the plasma structures, we found strong wind shear in the region in which the LQP echoes were produced. Meanwhile, the periodicity of the LQP echoes may have been related to the strength of the wind shear. It is proposed that the plasma structures were most likely produced by the Kelvin–Helmholtz instability (KHI).

Shock induced variable density flows in the vacuum microchannel: I. medium laser fluence

Laser-matter interactions with metal target cause plasma explosion and shock accelerated variable density flow instabilities in the Semiconfined Configuration (SCC). Their study gives deeper insight into the flow instabilities present in all microchannel devices. Blast wave motion along the SCC microchanel causes the Kelvin–Helmholtz (KH) instability and formation of vortex filaments for the critical Reynolds number. Appearing in all shear layers—it affects the fluid transport efficiency. Shear layer acceleration causes a Raleigh-Taylor instability (RTI). Oriented bubble growth by discrete merging indicates anisotropic RTI mixing. Similar RTI flame instability appears in the conversion of chemical energy into electricity affecting microcombustion efficiency. Another case of anisotropic RTI is the flow boiling for cooling of chips and microelectronic devices. The RTI boiling which appears for the critical heat flux is based on rising surface vapor columns (oriented bubble growth) with liquid counterflow (spike prominences) for the critical wavelength at density interface. The RT bubble merging graph trees determine turbulent mixing which affects the heat transfer rates. Bottom-wall turbulent flow in the SCC microchannel causes streaks of the low momentum fluid and formation of hairpin vortex packets with lattice organization. This makes possible to quantify parameters responsible for the evolution of hairpin vortex packets in the microchannel devices. Appearing from the low to the high Reynolds numbers they affect the transport properties, control of the fluid motion, enhancement of mixing, or the separation of fluids. Fluid particle ejecta—thin supersonic jets - evolve into long needle-like jets which start spiraling, helical pairing and swirling in the field of thermal gradients. Such instabilities appear in the microcombustion flame instability and in the space micropropulsion systems. Oscillating and spiral flames appear in the presence of thermal gradient in the microchannel, due to the combined effects of thermal gradient fields and the mixture flow rates.

Formation of Billowed Structures Along Margins of Basaltic Intrusions: Insights from Big Bend National Park, Texas and 71 Gulch, Idaho

Billowed margins of basaltic intrusions provide opportunity to study the conditions of non-explosive magma-pore water interactions and how magma propagates and interacts with unconsolidated sediment in the subsurface. Basaltic intrusions exhibiting billowed structures that formed at different paleodepths in Big Bend National Park, Texas (200–500 m) and 71 Gulch, Idaho (≤24 m) have been examined and reconstructed in three dimensions (3D) using structure for motion (SfM) techniques. To characterize billowed structure diversity, a classification system was developed utilizing 21 high-resolution 3D models of outcrops from both locations. Seven structure types are recognized: linear, sinuous, teardrop, bulbous, circular, transitional, and overprinted. Relative abundances of each structure type are similar between both localities. All structures identified were recognized along dikes, apart from one sill. Each site hosts ≥3 structure types, no one structure type is dependent on another, and the dimensions of structure types are similar at both localities with insignificant variation in mean wavelength and amplitude of crests and troughs (≤6 cm), suggesting that depth of intrusion is not a major influence in the occurrence of different structure types. Aligned plagioclase laths and stretched vesicles along margins are evidence that billow formation occurred while margins were fluid, deformable, and flowing. These observations indicate billowed structures likely form during the initial stages of intrusion while the margins are experiencing Kelvin-Helmholtz instabilities during cooling. As billowed intrusions from varying paleodepths at the two localities both occur near phreatomagmatic diatreme structures, changes in the rate or flux of magma intrusion or local heterogeneities in the host sediment and/or pore fluid content may be more responsible for explosive vs. non-explosive interactions rather than simple arrival of magma.

Instability of mixing layers: Momentum and thermal transport in the continuum breakdown regime.

We examine the momentum and thermal transport in the continuum breakdown regime of a mixing layer flow, which exhibits Kelvin-Helmholtz instability under ideal continuum conditions. The Grad 13 moment model is used as it provides an adequate description of the flow physics (second-order accurate in Knudsen number) in the transition regime. Analytical solutions are developed under breakdown conditions for two-dimensional, compressible, parallel shear flows. It is shown that the deviation of viscous stress and heat flux from the Navier-Stokes-Fourier system follows two different scaling regimes depending upon the Mach number. At low Mach numbers, the departure of all stress and heat-flux components depends only upon the Knudsen number. At high Mach number, the scaling of shear stress and transverse heat flux depends on the product of the Knudsen and Mach numbers. The normal stresses depend individually on the Knudsen and Mach number. The scaling results are verified against numerical simulations of compressible mixing layers performed using the unified gas kinetic scheme for various degrees of rarefaction.

Heat transfer correlation of smooth inverted annular film boiling in a circular tube

Experiments on steady film boiling heat transfer are conducted in a circular tube with test conditions covering the range of 10–50 K inlet subcooling, 0.2–2.5 MPa pressure, and 100–1000 kg/m2s mass flux. Through analysis of experimental data and comparison with corresponding predictions by the existing correlations, they are demonstrated to be insufficient for accurately capturing the varying trends of the experimental data within inverted annular flow boiling (IAFB) regime. It is discovered that IAFB is divided into three regions, namely, smooth IAFB, rippled IAFB and turbulent IAFB, by the distribution of wall temperature and its rate of change along axial direction, which is consistent with relevant observations from previous study. Only the vapor film in the smooth IAFB can be attributed to laminar vapor film flow, which leads to the discrepancy in the overall trend of existing correlations and experimental data. As fluid flows downstream, interface interaction due to velocity difference between vapor film and liquid core accounts for the fluctuations. The assumption that transition between smooth and rippled IAFB is result from Kelvin-Helmholtz (K-H) instability accounts for the flow characteristics of vapor film in smooth IAFB. Heat transfer at interface is obtained based on vapor–liquid velocity. Additionally, a new heat transfer correlation for smooth IAFB is proposed. The MAPE of the new correlation is 16.4%, and the RMSPE is 18.5%. The new correlation could more accurately describe the heat transfer characteristic in smooth IAFB.