Turbulent Mixing Layer Research Articles

Carrier-phase direct numerical simulation and flamelet modeling of NOx formation in a pulverized coal/ammonia co-firing flame

A carrier-phase direct numerical simulation (CP-DNS) of a coal/ammonia co-firing flame in a turbulent mixing layer is performed with a detailed chemical reaction mechanism including NOx formation. The combustion characteristics and NOx formation mechanism of the coal/ammonia co-firing flame are investigated in detail through a conditional reaction pathway analysis. A new four-fuel-stream flamelet model is proposed to adapt to the piloted pulverized coal/ammonia co-combustion system, in which the complex mixing among the volatiles, char off-gases, ammonia, and the pilot stream is characterized with four mixture fractions. The flamelet solutions in the fourfold mixture fraction space are transformed into a uniform space to facilitate the access to the flamelet table. The performance of the flamelet model is evaluated through an a priori analysis, a combustion mode analysis, and a chemical timescale analysis. For the turbulent coal/ammonia co-firing flame studied, two distinct flame fronts within the shear layer are observed with the upper part being dominated by ammonia combustion and the lower part being governed by coal combustion. The conditional analyses show that NO species is mainly generated by the lower layer governed by pulverized coal combustion where the nitrogen-containing species from volatile combustion are highly involved in the NO formation. The NO concentration first increases with decreasing the fraction of ammonia in the local computational cell, and then decreases towards pure pulverized coal combustion, which is characterized with a newly introduced manifold coordinate. While the temperature and the major and intermediate species mass fractions in the coal/ammonia co-firing flame can be reasonably predicted by the proposed flamelet model, the NOx species mass fractions are not well predicted in the region where pulverized coal combustion dominates. It is clarified that the inaccurate prediction of the NOx species is not caused by interpolation errors, multiple combustion modes or the kinetics of NOx formation, but mainly attributed to the definition of the progress variable.Novelty and significance statementThe novelty of this research is the first DNS of a coal/ammonia co-firing flame in a turbulent mixing layer, and the proposal of a new flamelet tabulation method for predicting NOx formation in the coal/ammonia co-firing flame. It is significant because• Coal/ammonia co-firing promises to provide continuous, secure power supply at a much reduced carbon footprint compared to pure coal combustion;• The underlying physics governing the coal/ammonia co-firing are analyzed using carrier-phase DNS;• The performance of the newly proposed flamelet model in predicting the NOx species is evaluated based on the DNS dataset.

Artificial neural network aided vapor–liquid equilibrium model for multi-component high-pressure transcritical flows with phase change

In this work, an artificial neural network (ANN) aided vapor–liquid equilibrium (VLE) model is developed and coupled with a fully compressible computational fluid dynamics (CFD) solver to simulate the transcritical processes occurring in high-pressure liquid-fueled propulsion systems. The ANN is trained in Python using TensorFlow, optimized for inference using Open Neural Network Exchange Runtime, and coupled with a C++ based CFD solver. This plug-and-play model/methodology can be used to convert any multi-component CFD solver to simulate transcritical processes using only open-source packages, without the need of in-house VLE model development. The solver is then used to study high-pressure transcritical shock-droplet interaction in both two- and four-component systems and a turbulent temporal mixing layer (TML), where both qualitative and quantitative agreement (maximum relative error less than 5%) is shown with respect to results based on both direct evaluation and the state-of-the-art in situ adaptive tabulation (ISAT) method. The ANN method showed a 6 times speed-up over the direct evaluation and a 2.2-time speed-up over the ISAT method for the two-component shock-droplet interaction case. The ANN method is faster than the ISAT method by 12 times for the four-component shock-droplet interaction. A 7 times speed-up is observed for the TML case for the ANN method compared to the ISAT method while achieving a data compression factor of 2881. The ANN method also shows intrinsic load balancing, unlike traditional VLE solvers. A strong parallel scalability of this ANN method with the number of processors was observed for all the three test cases. Code repository for 0D VLE solvers, and C++ ANN interface—https://github.com/UMN-CRFEL/ANN_VLE.git.

Mean-flow structures of the turbulent boundary layers bounding a two-dimensional separation bubble

Understanding the mean-flow structures of a separated turbulent boundary layer (TBL) is crucial for turbulence modeling. This study investigates the spatial scaling properties of the total shear stress and mixing length in the TBLs bounding a two-dimensional (2D) separation bubble, aiming to derive analytical descriptions for the entire mean-velocity profiles of the TBLs. For the adverse pressure gradient (APG) TBL upstream of the separation bubble, the total shear stress possesses a two-layer structure with an inner layer adhering to a linear law and an outer layer conforming to a defect power law. In contrast, the mixing length profile consists of four layers, namely the viscous sublayer, the buffer layer, the overlap layer, and the wake region. Each of the layers exhibits a power law or a defect power law relationship with the spatial coordinate normal to the wall. In the four-layer structure, three parameters are sensitive to the variation of the APG: the buffer-layer thickness, the relative magnitude of the mixing length at the boundary layer edge, and a defect power law exponent quantifying the extent of the wake region. For the reattached TBL downstream of the separation bubble, the total shear stress consists of two parts. One part is induced by the pressure gradient and retains the two-layer structure, while the other, engendered by the intense turbulence advected from the separated shear layer, exhibits a dual-power-law distribution. The advected turbulence also significantly alters the four-layer structure of the mixing length, resulting in an augmented buffer layer, a diminished overlap layer, and a wake region that mimics a turbulent mixing layer. Via a dilation ansatz to describe the scaling transition between adjacent layers, the study formulates the complete profiles of the total shear stress and mixing length. The formulation leads to the derivation of novel analytical expressions for the entire mean-velocity profiles of the TBLs. The expressions are in precise accord with the direct numerical simulations of an incompressible 2D separation-bubble flow and a 2D impinging shock wave/TBL interaction. The elucidation of the mean-flow structures through this study is anticipated to facilitate the analysis of turbulence models, thereby enhancing their performance in simulating separated TBLs. The construction of the mean-flow descriptions by inspecting the spatial scaling properties of turbulence paves a promising way for theoretical exploration of complex nonequilibrium TBLs.

Rapid Restratification Processes Control Mixed Layer Turbulence and Phytoplankton Growth in a Deep Convection Region

AbstractThe Gulf of Lion, Northwestern Mediterranean Sea, is one of few oceanic regions where deep convection occurs. We investigate the restratification following a convection event using measurements from an ocean glider equipped with turbulence microstructure sensors. This unique combination of instruments provides a high‐resolution description of the mixed layer with regard to turbulence, stratification and chlorophyll. We observe a rapid restratification process that proceeds over a timescale of days to one week. We find that restratification exerts a leading order control on surface mixed layer turbulence variability, as abrupt changes in turbulence dissipation rates are associated with the formation of near‐surface stratification. The near‐surface formation of stratification occurs through both the diurnal variability in surface buoyancy fluxes and through lateral advective processes. We conclude that daily near‐surface processes that influence stratification control mixed layer turbulence levels, and thus the phytoplankton response in the critical transition period to spring bloom.

Taming the TuRMoiL: The Temperature Dependence of Turbulence in Cloud–Wind Interactions

Turbulent radiative mixing layers play an important role in many astrophysical contexts where cool (≲104 K) clouds interact with hot flows (e.g., galactic winds, high-velocity clouds, infalling satellites in halos and clusters). The fate of these clouds (as well as many of their observable properties) is dictated by the competition between turbulence and radiative cooling; however, turbulence in these multiphase flows remains poorly understood. We have investigated the emergent turbulence arising in the interaction between clouds and supersonic winds in hydrodynamic enzo-e simulations. In order to obtain robust results, we employed multiple metrics to characterize the turbulent velocity, v turb. We find four primary results when cooling is sufficient for cloud survival. First, v turb manifests clear temperature dependence. Initially, v turb roughly matches the scaling of sound speed on temperature. In gas hotter than the temperature where cooling peaks, this dependence weakens with time until v turb is constant. Second, the relative velocity between the cloud and wind initially drives rapid growth of v turb. As it drops (from entrainment), v turb starts to decay before it stabilizes at roughly half its maximum. At late times, cooling flows appear to support turbulence. Third, the magnitude of v turb scales with the ratio between the hot phase sound-crossing time and the minimum cooling time. Finally, we find tentative evidence for a length scale associated with resolving turbulence. Underresolving this scale may cause violent shattering and affect the cloud’s large-scale morphological properties.

LES of autoignition process in a turbulent droplet-laden mixing layer – impact of an evaporation model and discretisation method

This paper focuses on the numerical modelling of the autoignition process in an ethanol-air mixture within a turbulent droplet-laden temporally evolving mixing layer, utilising the Large Eddy Simulation (LES) method. We compare results obtained by implementing two notably distinct discretisation methods (second-order Total Variation Diminishing – TVD, and fifth-order Weighted Essentially Non-Oscillatory – WENO) in the species and energy transport equations. Additionally, we examine three widely recognised evaporation models (the ‘ D 2 ’ model, the Abramzon–Sirignano model, and a non-equilibrium model based on the Langmuir-Knudsen law). For modelling the combustion process, we employ the Eulerian Stochastic Field method and a laminar chemistry model. The simulations encompass two droplet diameter distributions, one with low and the other with significant diameter non-uniformity. Furthermore, nine distinct initial conditions are considered. Our findings indicate that adjusting the discretisation scheme or the evaporation model can yield simulation outcomes differing to a substantial degree, similar to the impact of modifying the actual properties of the spray droplets or their initial distribution. This is particularly pronounced in the case of sprays featuring a narrowed droplet diameter distribution. We found that the Abramzon–Sirignano model results in the most rapid autoignition. The TVD scheme contributes to a more even spatial distribution of the evaporated fuel, consequently facilitating quicker evaporation and faster autoignition. Using the WENO scheme generates higher fuel vapor gradients near the droplets, leading to less intense evaporation and slower fuel transport in space.

Intermediate Gas Phases within Turbulent Radiative Mixing Layers

Turbulent radiative mixing layers (TRMLs) occur ubiquitously in astrophysical environments; e.g., TRMLs are prevalent within galactic outflows at the intersections between hot supernovae ejecta and cold molecular clouds. A velocity shear between the rapidly outflowing hot gas and cold clouds drives the Kelvin–Helmholtz instability producing TRMLs, with radiative cooling dominating the heat transfer between the gas phases. Using hydrodynamic simulations, we have modeled TRMLs for a range of overdensities (100, 1000, 3000) applied to cold phase temperatures of 400, 103, and 104 K. The production of an intermediate gas phase at the interface between the hot and cold phases is consistently observed at ∼104 K for molecular clouds, in agreement with larger-scale wind-tunnel simulations.

MÔ PHỎNG XOÁY TÁCH RỜI TRONG NGHIÊN CỨU CÁC ĐẶC TRƯNG KHÍ ĐỘNG CỦA DÒNG CHẢY RỐI TRÊN VÙNG TƯƠNG TÁC

The turbulent mixing layer is a phenomenon that occurs and has many applications in aeronautical engineering. The study of mixing layer turbulence by the numerical method is suitable under domestic conditions. Two numerical simulation methods are used in most studies including, direct numerical simulation DNS and large eddy simulation LES, while there are very few studies using detached eddy simulation DES. In this article, the authors use numerical simulation methods evolving Reynolds average Navier-Stokes simulation RANS and detached eddy simulation DES to propose a numerical simulation model that allows to study of the aerodynamic characteristics of the turbulent flow interaction. RANS is used to evaluate grid independence. With the selected mesh, two models counting the DES-2D model and the DES-3D model are conducted. Processing velocity field data with Matlab code, the results show that the average flow behavior has closely followed the experimental results, while the turbulent kinetic energy still has a large error. POD and Q-criterion analysis has shown some large coherent structures in the interaction region.

Investigating Ionization in the Intergalactic Medium

The intergalactic medium (IGM) contains >50% of the baryonic mass of the Universe, yet the mechanisms responsible for keeping the IGM ionized have not been fully explained. Hence, we investigate ion abundances from the largest blind QSO absorption catalog for clouds that show C iv, N v, and O vi simultaneously. The wavelength range of present UV spectrographs, however, makes it possible to probe C iv and O vi only over a small range of redshift (z ≈ 0.12–0.15). As a result, we only have five IGM absorbing clouds, yet these provide a powerful and representative tool to probe the IGM ionization state. We found one cloud to be in collisional ionization equilibrium while three of the five showed signs of being produced by nonequilibrium processes, specifically conductive interfaces and turbulent mixing layers. None of the models we explore here were able to reproduce the ionization state of the remaining system. Energetic processes, such as galactic feedback from star formation and active galactic nucleus winds, would be excellent candidates that can cause such widespread ionization.

Opposed jet burner and Tsuji burner for representative laminar flamelet with differential molecular diffusion for non-premixed combustion modeling

Modeling differential molecular diffusion in turbulent non-premixed combustion remains challenging for flamelet models. Laminar flamelet is a key building block of flamelet models for turbulent combustion. One significant challenge that has not been adequately addressed is the representativity of laminar flamelet for the characteristics of differential molecular diffusion in turbulent combustion. Laminar flamelet is typically generated based on two conceptual burner configurations, the opposed jet burner and the Tsuji burner. The two burners are commonly viewed to be equivalent for the description of laminar flamelet structures. In this work, a major difference between them is revealed for the first time when they are used to represent differential molecular diffusion. The traditional opposed jet burner yields almost fixed equal diffusion locations defined as the points with equal values of mixture fractions based on different elements in the mixture fraction space. The Tsuji burner, with a slight extension, can produce a continuous variation of the equal diffusion locations in the mixture fraction space. This variation of the equal diffusion locations is shown to be an important feature of non-premixed combustion, as demonstrated in a laminar jet mixing layer, a turbulent jet mixing layer, and a turbulent jet non-premixed flame. Theoretical analysis is conducted based on an idealized opposed jet mixing layer to identify the major cause of the difference between the two burners. It is found that the inlet diffusion caused the major difference in differential molecular diffusion in the two burners. The curvature difference, another main difference between the two burners, however, is not found to cause any major difference in differential molecular diffusion in the two burners. Based on the finding, a modified opposed jet burner with disabled inlet diffusion is introduced and is shown to be able to reproduce the feature of differential molecular diffusion observed in the Tsuji burner. Both burners, after properly generalized with suitable boundary conditions and additional parameters (like the velocity ratio of fuel and oxidizer), are expected to be suitable and equivalent choices for the generation of representative laminar flamelets that can be used eventually in flamelet modeling of differential molecular diffusion in turbulent non-premixed combustion.Novelty and significance statement: (1) A significant difference between the opposed jet burner and the Tsuji burner is revealed for generating representative laminar flamelet with differential molecular diffusion (under specific conditions); (2) The flamelet generated from the Tsuji burner is shown to be representative of differential molecular diffusion found in practically relevant mixing and combustion problems; (3) The treatment of inlet diffusion is identified to be the major cause of the difference observed between the two burners; (4) Modifications are introduced to the opposed jet burner to reconcile the two different burners for producing representative flamelet that can be used in the modeling of differential molecular diffusion in turbulent non-premixed combustion.

Interfaces of high- and low-speed large-scale structures in compressible turbulent mixing layers: compressibility effects and structures

Direct numerical simulations of temporally developing mixing layers have been performed to investigate the effects of compressibility on statistics and structures near the interfaces of high- and low-speed large-scale structures (LSSs), covering a range of convective Mach numbers from $M_c=0.2$ to $1.8$ and Taylor Reynolds numbers up to 290. The interfaces of LSSs are directly defined by the isosurface of zero fluctuating streamwise velocity. The characteristic velocity jump at the interfaces grows rapidly in the transition stage and then decreases until reaching a gradual self-similar stage. The evolution of velocity jump is evidently delayed as the convective Mach number increases. The interface layer is formed by the shearing motion of neighbouring LSSs. A small vortical motion characterized by the Kolmogorov scale is induced in the interface layer by shear-dominated outer regions. The strengths of outer shearing motion and central vortical motion are greater at the leading edge, but smaller at the trailing edge, which is also reflected in a larger velocity jump at the leading edge. As the convective Mach number increases, the small-scale vortical structure is obviously suppressed by compressibility. At high convective Mach number $M_c=1.8$ , the compressive shear-dominated outer regions are linked with a sheet-like structure passing through the centre of the expansion region corresponding to a small-scale vortical structure. The compressibility and shearing strength near the interface are highly dependent on the interface orientation.

Carrier-Phase DNS of Ignition and Combustion of Iron Particles in a Turbulent Mixing Layer

Three-dimensional carrier-phase direct numerical simulations (CP-DNS) of reacting iron particle dust clouds in a turbulent mixing layer are conducted. The simulation approach considers the Eulerian transport equations for the reacting gas phase and resolves all scales of turbulence, whereas the particle boundary layers are modelled employing the Lagrangian point-particle framework for the dispersed phase. The CP-DNS employs an existing sub-model for iron particle combustion that considers the oxidation of iron to FeO and that accounts for both diffusion- and kinetically-limited combustion. At first, the particle sub-model is validated against experimental results for single iron particle combustion considering various particle diameters and ambient oxygen concentrations. Subsequently, the CP-DNS approach is employed to predict iron particle cloud ignition and combustion in a turbulent mixing layer. The upper stream of the mixing layer is initialised with cold particles in air, while the lower stream consists of hot air flowing in the opposite direction. Simulation results show that turbulent mixing induces heating, ignition and combustion of the iron particles. Significant increases in gas temperature and oxygen consumption occur mainly in regions where clusters of iron particles are formed. Over the course of the oxidation, the particles are subjected to different rate-limiting processes. While initially particle oxidation is kinetically-limited it becomes diffusion-limited for higher particle temperatures and peak particle temperatures are observed near the fully-oxidised particle state. Comparing the present non-volatile iron dust flames to general trends in volatile-containing solid fuel flames, non-vanishing particles at late simulation times and a stronger limiting effect of the local oxygen concentration on particle conversion is found for the present iron dust flames in shear-driven turbulence.

Numerical Investigation of Diffusion Flame in Transonic Flow with Large Pressure Gradient

A finite-volume method for the steady, compressible, reacting, turbulent Navier–Stokes equations is developed by using a steady-state-preserving splitting scheme for the stiff source terms in chemical reactions. Laminar and turbulent reacting flows in a mixing layer with a large streamwise pressure gradient and acceleration are studied and compared to boundary-layer solutions. It reveals that chemical reactions strongly enhance turbulent transport due to the intensive production of turbulence by the increased velocity gradients and thus produce large turbulent viscosities in the reaction region. The influence of vitiated air on the combustion process and aerodynamic performance is also investigated in the cases of turbulent mixing layers and highly loaded transonic turbine cascades. Both cases indicate the viability of the turbine-burner concept.

Characterization of the mixing layer self-similarity with multiple parameters

In the paper, results from direct numerical simulation of a planar incompressible mixing layer spatially developing incompressible between two co-flowing laminar boundary layers are used to analyze a possibility for multiple flow parameters to achieve self-similarity within the same flow region. The Reynolds numbers for the boundary layers are 3930 and 2412 based on the free-stream velocities far above and below the splitter plate and the boundary layer thicknesses at the splitter plate trailing edge. The three-dimensional computational domain is sufficiently large for the mixing layer transition to fully turbulent far upstream the domain exit. The mixing layer growth is characterized using various definitions of the mixing layer thickness. It is shown that the proposed mixing layer thickness based on the gradient of the streamwise mean velocity in the transverse direction defines more accurately the area of turbulent mixing. Three regions of the flow linear growth are discovered using a rigorous approach, with only one of them being located within the fully turbulent mixing layer. Other parameters included in the flow self-similarity analysis are the streamwise and transverse mean velocities along with the Reynolds stresses. New normalization is proposed to observe self-similarity of the transverse mean velocity. The flow region where all considered parameters exhibit self-similarity is determined. It is shown that this region is limited by the “pulsating” streamwise distribution of the transverse mean velocity. The computational domain dimension along with the boundary conditions in the transverse direction for all considered parameters are suggested for the Reynolds-averaged Navier–Stokes simulations.

Carrier-phase DNS study of particle size distribution effects on iron particle ignition in a turbulent mixing layer

The ignition and combustion of iron particles in a turbulent mixing layer is studied by means of three-dimensional carrier-phase direct numerical simulations (CP-DNS). A particular focus is set on particle size distribution (PSD) effects on the ignition behaviour by comparing CP-DNS results from using a realistic experimental PSD to DNS data based on a monodisperse (MD) particle cloud with the same equivalence ratio. The CP-DNS solves the Eulerian transport equations of the reacting gas phase and resolves all turbulent scales, while the particle boundary layers are modelled in the Lagrangian point-particle framework. A previously validated sub-model for the oxidation of iron to Wüstite (FeO) that accounts for both diffusion- and kinetically-limited combustion is employed. The mixing layer is initialised with an upper stream of air carrying cold iron particles and an opposed lower stream of hot air. Simulation results show distinct differences in the ignition behaviour between the MD and PSD cases. The ignition of the PSD case is delayed compared to the MD case and does not show any significant particle clustering prior to ignition. Further investigations indicate that the particle size has a crucial effect on the mixing process and ignition time. Small particles start their oxidation process early and already consume some of the available oxygen, while not crucially affecting the gas temperature due to their limited iron mass contribution. Conversely, a slower entrainment into the lower stream combined with higher thermal inertia and the prior oxygen depletion by the small particles leads to a delayed oxidation of the larger particles. As a net result, the PSD case shows a wide spread of individual particle ignition delay times and overall delayed bulk ignition compared to the MD case, where the majority of the particles ignites over a shorter period of time.

High-speed 1-D and 2-D Raman scattering measurement for quantitative characterization of transient hydrogen jets

This paper reports the first quantitative characterization of transient hydrogen (H2) jets issued from a single-hole injector into nitrogen (N2), using high-speed 1-D and 2-D Raman scattering techniques. By employing the pulse-burst laser as the light source together with the low-noise EMCCD cameras running in subframe burst gating mode for signal collection, single-shot 1-D mole fraction of H2 (xH2) deduced from the ratio of H2 and N2 Raman scattering intensities, is obtained at repetition rate of 50 kHz. By changing the injection conditions, transient jet behaviors of both subsonic and highly under-expanded jets are well captured in single shots of 1-D, and the dynamic of the oblique shocks is visualized by local peaks in relative H2 number density. Statistical information regarding the jet behaviors is also obtained by repeating the measurement at same injection conditions, and the mixing between H2 and N2 can be directly evaluated by the average xH2. A coefficient of variance (COV) of 1.6 to 6 % is measured for all the conditions near the core region of the jets at late injection time, which demonstrates a good repeatability of the 1-D line measurements during the quasi-steady injection period. Using high-speed CMOS cameras, 2-D Raman imaging of transient H2 jets is performed at repetition rate of 10 kHz. Results in single shots demonstrate the vortices formation at the start of injection leading to turbulent mixing layers at later stages. Due to the low SNR, the accuracy and precision of the 2-D measurements is limited, and a one-to-one comparison of statistical results shows a maximum deviation of approximately 15 % from the more accurate 1-D measurements.

Cloud atlas: navigating the multiphase landscape of tempestuous galactic winds

ABSTRACT Galaxies comprise intricate networks of interdependent processes which together govern their evolution. Central among these are the multiplicity of feedback channels, which remain incompletely understood. One outstanding problem is the understanding and modelling of the multiphase nature of galactic winds, which play a crucial role in galaxy formation and evolution. We present the results of three-dimensional magnetohydrodynamical simulations of tall–box interstellar medium (ISM) patches with clustered supernova-driven outflows. Dynamical fragmentation of the ISM during superbubble breakout seeds the resulting hot outflow with a population of cool clouds. We focus on analyzing and modelling the origin and properties of these clouds. Their presence induces large-scale turbulence, which, in turn, leads to complex cloud morphologies. Cloud sizes are well described by a power-law distribution and mass growth rates can be modelled using turbulent radiative mixing layer theory. Turbulence provides significant pressure support in the clouds, while magnetic fields only play a minor role. We conclude that many of the physical insights and analytic scalings derived from idealized small-scale simulations of turbulent radiative mixing layers and cloud–wind interactions are directly translatable and applicable to these larger scale cloud populations. This opens the door to developing effective subgrid recipes for their inclusion in global-scale galaxy models where they are unresolved.

Magnetic fields in multiphase turbulence: impact on dynamics and structure

ABSTRACT Both multiphase gas and magnetic fields are ubiquitous in astrophysics. However, the influence of magnetic fields on mixing of the different phases is still largely unexplored. In this study, we use both turbulent radiative mixing layer (TRML) and turbulent box simulations to examine the effects of magnetic fields on cold gas growth rates, survival, and the morphology of the multiphase gas. Our findings indicate that, in general, magnetic fields suppress mixing in TRMLs, while turbulent box simulations show comparatively marginal differences in growth rates and survival of the cold gas. We reconcile these two seemingly contrasting results by demonstrating that similar turbulent properties result in comparable mixing, regardless of the presence or absence of magnetic fields. We, furthermore, find the cold gas clump size distribution to be independent of the magnetic fields, but the clumps are more filamentary in the MHD case. Synthetic Mg ii absorption lines support this picture being marginally different with and without magnetic fields; both cases align well with observations. We also examine the magnetic field strength and structure in turbulent boxes. We generally observe a higher mean magnetic field in the cold gas phase due to flux freezing and reveal fractal-like magnetic field lines in a turbulent environment.

Simulations of weakly magnetized turbulent mixing layers

ABSTRACT Radiative turbulent mixing layers (TMLs) are expected to form pervasively at the phase boundaries in multiphase astrophysical systems. This inherently small-scale structure is dynamically crucial because it directly regulates the mass, momentum, and energy exchanges between adjacent phases. Previous studies on hydrodynamic TMLs have revealed the interactions between cold and hot phases in the context of the circumgalactic medium, offering important insight into the fate of cold clouds traveling through hot galactic winds. However, the role of magnetic field has only been sparsely investigated. We perform a series of 3D magnetohydrodynamics simulations of such mixing layers in the presence of weak to modest background magnetic field. We find that due to field amplification, even relatively weak background magnetic fields can significantly reduce the surface brightness and inflow velocity of the hot gas in the mixing layer. This reduction is attributed to a combination of magnetic pressure support and direct suppression of turbulent mixing, both of which alter the phase structures. Our results are largely independent of thermal conduction and converged with resolution, offering insights on the survival of cold gas in multiphase systems.

Interface instability of the thermal plasma jet

This work is first focused to experimentally study the interface instability and expansion mechanism of thermal plasma jet and provide a better understanding of the complex fluid-dynamic interactions occurring on the surface of the plasma bubble due to the Kelvin–Helmholtz effect. The experimental techniques used include a plasma generator, a pulse-forming network based on the capacitive energy storage, pressure measurement system along the capillary tube, and high-speed camera system to trace the development processes of the plasma interface. Results indicate that the plasmas jet has a better advantage of radial expansion with a high light at the beginning. However, the axial expansion velocity is larger than the radial one with time going on; thus, a torch-shaped jet body occurs under the Rayleigh–Taylor effect and can be divided into two parts including a plasma head and tail. With a dissipation of the initial energy and turbulent mixing between the plasmas and the gas, the jet boundary is broken and even the local rupture phenomena occur on the plasma jet surface. The turbulent dissipation is also very violent when the discharge voltage increases to 3000 V; thus, the turbulent mixing layer between the plasma jet and the gas is quite thicker and the plasma jet boundary is also fuzzy resulting in that the fold surfaces with much larger degree exist earlier. These experimental phenomena are also explained further from the mechanism by deriving the momentum equations of the interface of the plasma jet into the gas. Finally, a fitting formula of the surface area as an important factor in the expansion process of the plasma is obtained to analyze the interface characteristic of the plasma jet.