Homogeneous Turbulence Research Articles

Assessment of a multiphase formulation of one-dimensional turbulence using direct numerical simulation of a decaying turbulent interfacial flow

The interaction between turbulence and surface tension is studied numerically using the one-dimensional-turbulence (ODT) model. ODT is a stochastic model simulating turbulent flow evolution along a notional one-dimensional line of sight by applying instantaneous maps that represent the effects of individual turbulent eddies on property fields. It provides affordable high resolution of interface creation and property gradients within each phase, which are key for capturing the local behavior as well as overall trends, and has been shown to reproduce the main features of an experimentally determined regime diagram for primary jet breakup. Here ODT is used to investigate the interaction of turbulence with an initially planar interface. The notional flat interface is inserted into a periodic box of decaying homogeneous isotropic turbulence, simulated for a variety of turbulent Reynolds and Weber numbers. Unity density and viscosity ratios serve to focus solely on the interaction between fluid inertia and the surface-tension force. Statistical measures of interface surface density and spatial structure along the direction normal to the initial surface are compared to corresponding direct-numerical-simulation (DNS) data. Allowing the origin of the lateral coordinate system to follow the location of the median interface element improves the agreement between ODT and DNS, reflecting the absence of lateral nonvortical displacements in ODT. Beyond the DNS-accessible regime, ODT is shown to obey the predicted parameter dependencies of the Kolmogorov critical scale in both the inertial and dissipative turbulent-cascade subranges. Notably, the probability density function of local fluctuations of the critical scale is found to collapse to a universal curve across both subranges. Published by the American Physical Society 2024

Two-point statistics in non-homogeneous rotating turbulence

Time-resolved two-dimensional two-component particle image velocimetry measurements with high spatial resolution are carried out in a water tank agitated by four blades rotating at constant speed. Different blade geometries and rotation speeds are tested for the purpose of modifying turbulent flow conditions. In all cases where no baffles are used to break the rotation, the Zeman length is an order of magnitude smaller than the Taylor length. Compared with the cases with baffles which break the rotation, in the unbaffled cases turbulence production and/or mean advection are significant and the turbulence nonlinearity is dramatically reduced for the horizontal (i.e. normal to the axis of rotation) two-point turbulence fluctuating velocities. This nonlinearity reduction is manifest not only in the interscale turbulent energy transfer but also in the interspace turbulent energy transfer, which nearly vanishes. However, the nonlinearity is not reduced for the vertical two-point turbulence fluctuating velocities: the corresponding interscale turbulent transfer rate is in fact intensified, and its dependence on the two-point separation distance, as well as that of the corresponding interspace turbulent transfer rate, which does not vanish, is significantly modified. Even though non-homogeneities are very different for different blades and rotation speeds in the unbaffled cases, the horizontal fluctuating velocity's second-order structure function collapses with scalings which resemble predictions for homogeneous turbulence subject to strong rotation. The vertical fluctuating velocity's second-order structure function does not collapse for different blade geometries by neither these nor the Kolmogorov predictions.

Scale dependence of local shearing motion in decaying turbulence generated by multiple-jet interaction

Scale dependence of local shearing motion is investigated experimentally in decaying homogeneous isotropic turbulence generated through multiple-jet interaction. The turbulent Reynolds number, based on the Taylor microscale, is between approximately 900 and 400. Velocity fields, measured using particle image velocimetry, are analysed through the triple decomposition of a low-pass filtered velocity gradient tensor, which quantifies the intensities of shear and rigid-body rotation at a given scale. These motions manifest predominantly as layer and tubular vortical structures, respectively. The scale dependence of the moments of velocity increments, associated with shear and rigid-body rotation, exhibits power-law behaviours. The scaling exponents for shear are in quantitative alignment with the anomalous scaling of the velocity structure functions, suggesting that velocity increments are influenced predominantly by shearing motion. In contrast, the exponents for rigid-body rotation are markedly smaller than those predicted by Kolmogorov scaling, reflecting the high intermittency of rigid-body rotation. The mean flow structure associated with shear at intermediate scales is investigated with conditional averages around locally intense shear regions in the filtered velocity field. The averaged flow field exhibits a shear layer structure with aspect ratio approximately 4.5, surrounded by rotating motion. The analysis at different scales reveals the existence of self-similar structures of shearing motion across various scales. The mean velocity jump across the shear layer increases with the layer thickness. This relationship is well predicted by Kolmogorov's second similarity hypothesis, which is useful in predicting the mean characteristics of shear layers across a wide range of scales.

Posterior comparison of model dynamics in several hybrid turbulence model forms

Hybrid turbulence models that can accurately reproduce unsteady three-dimensional flow physics across the entire range of grid scales and turbulence dynamics from Reynolds-averaged Navier–Stokes (RANS), through large-eddy simulation (LES), down to direct numerical simulations (DNS) are of increasing interest to the turbulence modeling community. However, despite decades of research and development, the basic tasks of eliminating poor-performing hybrid RANS-LES models and accelerating adoption of superior models through well-designed validation and verification have yet to occur. As a step in this direction, in this work we evaluate thirteen different hybrid RANS-LES models via systematic grid refinement of decaying homogeneous isotropic turbulence. We further derive a novel mathematical framework for assessing the energy partitioning dynamics of each Hybrid RANS-LES model, wherein model-to-model variations in energy partitioning can be interpreted as different feedback mechanisms operating on a low-dimensional nonlinear dynamical system. We found that model forms similar to the flow simulation methodology—also often termed very-large eddy simulation—are dynamically inconsistent with DNS at all resolutions. Additionally, we found a strong dynamical similarity in the feedback mechanisms of all models related to detached eddy simulation and partially averaged Navier–Stokes that is inherent to their general model forms.

Experimental Study of Turbulent Premixed Flames of Gas-to-Liquids (GTL) Fuel in a Fan-Stirred Combustion Bomb

This study experimentally investigates the turbulent flame speeds (St) of Gas-to-Liquids (GTL) fuel and a 50/50 blend with diesel across turbulence intensities (u') ranging from 0.5 to 3.0 m/s and equivalence ratios (Φ) from 0.7 to 1.3. Experiments were conducted using a cylindrical fan-stirred combustion bomb at an initial temperature (Ti) of 463 K and atmospheric pressure under near homogeneous and isotropic turbulence (HIT) conditions. A pressure transducer measured the St of the propagating GTL flame. High-speed imaging reveals that changes in flame brightness are associated with variations in flame temperature and soot incandescence. Compared to diesel, stoichiometric GTL and the 50/50 diesel-GTL blend showed peak combustion pressure reductions of 8.9% and 4.9%, respectively. Rich diesel fuel and lean GTL exhibited higher pressure rise rates, flame propagation, and St due to their Lewis numbers (Le) being less than one, enhancing flame-turbulence interaction. At u` of 0.5, 1.5, and 3.0 m/s, St for GTL increased by approximately 3.6%, 5.3%, and 2.8%, respectively, when compared to diesel, with these increases observed at Φ < 1.1. Transition from wrinkled to corrugated flamelets with increasing u` supports models predicting flame stability and potential industrial combustion efficiency improvements. Comparisons with numerical St results using the Zimont Turbulent Flame Speed Closure (Zimont TFC) model showed good agreement at lower (u' = 0.5 m/s) and mid-range (u' = 2.0 m/s) turbulence intensities but discrepancies at higher intensities (u' = 3.0 m/s), highlighting the need to refine numerical models for more accurate predictions across all turbulence levels.

Strong evidence for universality in homogeneous compressible turbulence

Quantifying the degree of universality in compressible turbulence is challenging due to the existence of different modes and their complex interactions. For a restricted family of flows, Donzis and John [Phys. Rev. Fluids 5, 084609 (2020)] showed that universal behavior is indeed observed in compressible turbulence if the ratio of dilatational to solenoidal root mean square (rms) velocities (δ=u′d/u′s) is incorporated as a scaling parameter along with the traditional turbulent Mach number (Mt=u′/〈c〉, where u′ is the rms velocity and 〈c〉 is the mean speed of sound). In this paper, we argue for the generality of those results by analyzing a wide range of compressible turbulent flows spanning a variety of flow configurations and setups to assess the degree of universal behavior. These include, among others, reacting flows, flows with solenoidal, thermal, and dilatational forcing, and flows with mean shear and bulk viscosity. We also performed new direct numerical simulations, which include turbulence in situations where vibrational modes of constitutive molecules are not in thermal equilibrium. Collectively, we offer the largest comparison across studies in the literature to date. We find that despite the wide range of forcing conditions and physical processes, universality holds across all these turbulent flows to a very satisfactory degree when both δ and Mt are considered as intrinsic compressibility parameters. The statistics investigated here—single-point statistics up to order four—are chosen such that they represent different ranges across the spectrum of dynamically relevant turbulence scales. We discuss the applicability of the purposed universal behavior for other key statistics in these turbulent flows, including two-point statistics and inhomogeneity effects, and the perspective it opens for modeling them.

Super-resolution reconstruction of turbulence for Newtonian and viscoelastic fluids with a physical constraint

Super-resolution reconstruction (SR) of turbulent flow fields with high physical fidelity from low-resolution turbulence data is a novel and cost-effective way in a turbulence study. However, some naive image-to-image machine learning methods often produce nonphysical features inconsistent with the physical characteristics of turbulence. The present work proposes, respectively, convolutional neural network and generative adversarial network-based turbulence SR models using the kinetic energy spectra of turbulence flow as a physical constraint. The models have been validated in turbulence SR reconstruction for a Newtonian fluid under the flow condition of homogeneous isotropic turbulence at Reynolds number, Re=3140 and 4710 and viscoelastic fluid at the same Re numbers and Weisenberg number, Wi=0.796 and 1.194 (the elasticity number El=Wi/Re=0.000 254), respectively. The results show that with the energy spectra constraint (ESC), not only the nonphysical features occurred in the energy spectra of velocity field could be eliminated by the SR models, the errors of their reconstructed vorticity fields in comparison with the results of direct numerical simulation are also significantly smaller than those of the SR models without ESC. Therefore, incorporation of physical constraints is vital in preserving physical characteristics in SR of turbulent flow.

Reynolds stress decay modeling informed by anisotropically forced homogeneous turbulence

Reynolds stress decay modeling informed by anisotropically forced homogeneous turbulence

Stationary and nonstationary energy cascades in homogeneous ferrofluid turbulence

Stationary and nonstationary energy cascades in homogeneous ferrofluid turbulence

A Comparative Study of the Hydrogen Auto-Ignition Process in Oxygen–Nitrogen and Oxygen–Water Vapor Oxidizer: Numerical Investigations in Mixture Fraction Space and 3D Forced Homogeneous Isotropic Turbulent Flow Field

In this paper, we analyze the auto-ignition process of hydrogen in a hot oxidizer stream composed of oxygen–nitrogen and oxygen–water vapor with nitrogen/water vapor mass fractions in a range of 0.1–0.9. The temperature of the oxidizer varies from 1100 K to 1500 K and the temperature of hydrogen is assumed to be 300 K. The research is performed in 1D mixture fraction space and in a forced homogeneous isotropic turbulent (HIT) flow field. In the latter case, the Large Eddy Simulation (LES) method combined with the Eulerian Stochastic Field (ESF) combustion model is applied. The results obtained in mixture fraction space aim to determine the most reactive mixture fraction, maximum flame temperature, and dependence on the scalar dissipation rate. Among others, we found that the ignition in H2-O2-H2O mixtures occurs later than in H2-O2-N2 mixtures, especially at low oxidizer temperatures. On the other hand, for a high oxidizer temperature, the ignitability of H2-O2-H2O mixtures is extended, i.e., the ignition occurs for a larger content of H2O and takes place faster. The 3D LES-ESF results show that the ignition time is virtually independent of initial conditions, e.g., randomness of an initial flow field and turbulence intensity. The latter parameter, however, strongly affects the flame evolution. It is shown that the presence of water vapor decreases ignitability and makes flames more prone to extinction.

A theoretical study of a radiofrequency wave propagation through a realistic turbulent marine atmospheric boundary layer based on large eddy simulation

This study aims at modeling and investigating the impact of realistic turbulence inhomogeneity on radiowave propagation by utilizing large eddy simulations (LES). An up-to-date version of the X-LES method for generating realistic turbulent phase screens is first introduced. It combines atmospheric simulations with the classical Tatarski statistical modeling. This method naturally incorporates vertical turbulence inhomogeneity into phase screens at scales resolved by LES. It is then extended to sub-grid scales by weighting statistically generated phase variations with the vertical profile of the turbulent structure constant extracted from atmospheric data. This method is applied to replicate turbulence of a classical tropical marine atmospheric boundary layer. The impact of the generated medium on the propagation of a 10 GHz spherical wave emanating from a Gaussian aperture is analyzed through a statistical study of log-amplitude profiles performed for three different source altitudes. Results first show that contrary to a classical homogeneous turbulence modeling, log-amplitude profiles resulting from the propagation into inhomogeneous turbulence exhibit a statistical heterogeneity strongly dependent of source altitude. Furthermore, classical stochastic phase screen generation from a homogeneous Von-Kàrmàn Kolmogorov spectrum seems to give a statistically significant underestimation of the actual impact of turbulence compared to the X-LES method.

Effects of Turbulence on Upper-Tropospheric Ice Supersaturation

Abstract Effects of turbulence on ice supersaturation at cirrus heights (&gt;8 km) remain unexplored. Small-scale mixing processes become important for high Reynolds number flows, which may develop below the buoyancy length scale (10–100 m). The current study couples a stochastic turbulent mixing model with reduced dimensionality to an entraining parcel model to investigate, in large-ensemble simulations, how supersaturation evolves due to homogeneous turbulence in the stably stratified, cloud-free upper troposphere. The rising parcel is forced by a mesoscale updraft. The perturbation of an initially homogeneous vertical distribution of supersaturation is studied after a 36-m ascent in a baseline case and several sensitivity scenarios. Turbulent mixing and associated temperature fluctuations alone lead to changes in ensemble-mean distributions with standard deviations in the range 0.001–0.006, while mean values are hardly affected. Large case-to-case variability in the supersaturation field is predicted with fluctuation amplitudes of up to 0.03, although such large values are rare. A vertical gradient of supersaturation (≈10−3 m−1) is generated for high turbulence intensities due to the development of a dry-adiabatic lapse rate. Entrainment of slightly warmer (less than 0.1 K) environmental air into the parcel decreases the mean supersaturation by less than 0.01. Supersaturation fluctuations are substantially larger after entrainment events with an additional small offset in absolute humidity (by ±3.5%) between the parcel and environmental air. The predicted perturbations of ice supersaturation are significant enough to motivate studies of turbulence–ice nucleation interactions during cirrus formation that abandon the assumption of instantaneous mixing inherent to traditional parcel models. Significance Statement The purpose of our study is to investigate the effects of microscale turbulence on ice supersaturation in the upper troposphere. The associated variability in temperature and moisture fields is not resolved in cloud models and cannot easily be represented in terms of large-scale flow variables. We specify the conditions in which turbulent mixing and entrainment cause substantial variations in distributions of supersaturation. These include high turbulence intensity, strong atmospheric stability, and large moisture gradients. Our results suggest that turbulence may affect the strongly supersaturation-dependent ice formation processes in high-altitude clouds, pointing to the need to investigate cirrus formation in the presence of turbulence.

Theory and Modelling of Isotropic Turbulence: From Incompressible through Increasingly Compressible Flows

Homogeneous isotropic turbulence (HIT) has been a useful theoretical concept for more than fifty years of theory, modelling, and calculations. Some exact results are revisited in incompressible HIT, with special emphasis on the 4/5 Kolmogorov law. The finite Reynolds number effect (FRN), which yields corrections to that law, is investigated, using both Kármán–Howarth-type equations and a statistical spectral closure of the Eddy-Damped Quasi-Normal Markovian (EDQNM)-type. This discussion offers an opportunity to give an extended review of such spectral closures, from weak turbulence, as in wave turbulence theory, to a strong one. Extensions of the 4/5 or 4/3 Kolmogorov/Monin laws to anisotropic cases, such as stably stratified and MHD turbulence, are briefly touched on. Before addressing more recent work on compressible isotropic turbulence, the simplest case of quasi-incompressible turbulence subjected to externally imposed isotropic compression or dilatation is presented. Rapid distortion theory is found to be a poor model in this isotropic case, in contrast with its relevance in strongly anisotropic flow cases. Accordingly, a fully nonlinear approach based on a rescaling of all fluctuating variables is used, in order to show its interplay with the linear operator. This opens the discussion on the cases of homogeneous incompressible turbulence, where RDT and nonlinear models are relevant, provided that anisotropy is accounted for. Finally, isotropic compressible flows of increasing complexity are considered. Recent studies using weak turbulence theory, modelling, and DNS are discussed. A final unpublished study involves interactions between the solenoidal mode, inherited from incompressible turbulence, and the acoustic and entropic modes, which are specific to the compressible problem. An approach to acoustic wave turbulence, with resonant triads, is revisited on this occasion.

Effects of spatiotemporal plasma power distribution on the modeling of ignition kernel evolution in quiescent and turbulent methane/air mixtures

Abstract The present work improves a phenomenological plasma-assisted combustion model by integrating the spatiotemporal distribution of plasma power density, thereby considering the evolution of plasma streamers in the modeling, and subsequently, better predicting the ignition kernel evolution. The improved phenomenological model is validated against experiments representing the plasma discharge and post-discharge ignition kernel evolution. Specifically, the new model demonstrates a more accurate prediction of ultrafast gas heating and O2 dissociation during the plasma discharge, compared to the original model. In addition, the new model is found to closely match the experimental pressure wave and heated channel profiles post-discharge without the need for tuning the energy deposition (unlike the original model), highlighting its accuracy of post-discharge ignition kernel dynamics. The improved phenomenological model is then employed to investigate ignition kernel evolution for a stoichiometric methane-air discharge across various discharge gap configurations. Simulations reveal a non-uniform temperature and streamer distribution progressing from the electrode tips toward the center, contrasting uniform cylindrical discharges previously described in the original model. Streamer propagation is observed to be faster for larger gaps when maintained at the same average electric field for different discharge gaps. The tendency of smaller gaps to produce detached toroidal ignition kernels is observed, while larger gaps promote cylindrical and attached ignition kernels. Interactions between successive ignition kernels from consecutive discharges varied significantly, with the smallest gap (1 mm) promoting the quenching of the preceding ignition kernel due to the initial kernel–kernel separation. The intermediate gap (2 mm) promotes detached kernel growth. In contrast, in the largest gap (4 mm), kernels consistently combine and expand attached to electrodes. The impact of homogeneous isotropic turbulence is also explored, showing the persistence of ignition kernels early on but eventually quenching due to enhanced radical and heat losses with pronounced turbulence intensity.

Turbulence effect on disk settling dynamics

Turbulence can have a strong effect on the fall speed of snowflakes and ice crystals. In this experimental study, the behaviour of thin disks falling in homogeneous turbulence is investigated, in a range of parameters relevant to plate crystals. Disks ranging in diameter from 0.3 to 3 mm, and in Reynolds number $Re = 10\unicode{x2013}435$ , are dispersed in two air turbulence levels, with velocity fluctuations comparable to the terminal velocity. For each case, thousands of trajectories are captured and reconstructed by high-speed laser imaging, allowing for statistical analysis of the translational and rotational dynamics. Air turbulence reduces the disk terminal velocities by up to 35 %, with the largest diameters influenced most significantly, which is primarily attributed to drag nonlinearity. This is evidenced by large lateral excursions of the trajectories, which correlate with cross-flow-induced drag enhancement as reported previously for falling spheres and rising bubbles. As the turbulence intensity is increased, flat-falling behaviour is progressively eliminated and tumbling becomes prevalent. The rotation rates of the tumbling disks, however, remain similar to those displayed in still air. This is due to their large moment of inertia compared to the surrounding fluid, in stark contrast with studies conducted in water. In fact, the observed reduction of settling velocity is opposite to previous findings on disks falling in turbulent water. This emphasizes the importance of the solid-to-fluid density ratio in analogous experiments that aim to mimic the behaviour of frozen hydrometeors.

Robust control of charged particles position in a homogeneous isotropic turbulence using an electric field

This research investigates the position control of spherical charged particles suspended in a directly simulated homogeneous isotropic turbulence, which has the potential to improve the efficiency of the electrostatic precipitators (ESPs) in removing gaseous pollutants. To do this, the electric field used in ESPs is controlled by the sliding mode control method to steer the particles toward their desired positions. The particles are charged using the field charging method, which leads to electrical interactions between them. Employing the linked-list method reduces the computational cost related to these interaction forces. The dynamics of particles is simulated using the Maxey–Riley equation including the control force. Coupling the direct numerical simulation solver with the proposed control method and the linked-list scheme, along with solving the particle dynamics, complicates the problem from a control theory point of view. The robustness of the proposed model-free control method leads to overcoming disturbances and uncertainties in the simulation due to turbulence and repulsive electric forces. Hence, steering a large number of charged particles to places with more pollution concentration or increasing the duration of exposure of the particles to gaseous pollutants becomes possible. The findings of the present research are pertinent to improving the performance of air purification equipment.

The Relative Prevalence of Wave Packets and Coherent Structures in the Inertial and Kinetic Ranges of Turbulence as Seen by Solar Orbiter

The Solar Orbiter (SO) mission provides the opportunity to study the evolution of solar wind turbulence. We use SO observations of nine extended intervals of homogeneous turbulence to determine when turbulent magnetic field fluctuations may be characterized as: (i) wave packets and (ii) coherent structures (CSs). We perform the first systematic scale-by-scale decomposition of the magnetic field using two wavelets known to resolve wave packets and discontinuities, the Daubechies 10 (Db10) and Haar, respectively. The probability distribution functions (PDFs) of turbulent fluctuations on small scales exhibit stretched tails, becoming Gaussian at the outer scale of the cascade. Using quantile–quantile plots, we directly compare the wavelet fluctuations PDFs, revealing three distinct regimes of behavior. Deep within the inertial range (IR) both decompositions give essentially the same fluctuation PDFs. Deep within the kinetic range (KR) the PDFs are distinct as the Haar decompositions have larger variance and more extended tails. On intermediate scales, spanning the IR–KR break, the PDF is composed of two populations: a core of common functional form containing ∼97% of fluctuations, and tails that are more extended for the Haar decompositions than the Db10 decompositions. This establishes a crossover between wave-packet (core) and CS (tail) phenomenology in the IR and KR, respectively. The range of scales where the PDFs are two-component is narrow at 0.9 au (4–16 s) and broader (0.5–8 s) at 0.4 au. As CS and wave–wave interactions are both candidates to mediate the turbulent cascade, these results offer new insights into the distinct physics of the IR and KR.

Settling velocity characteristics of inertial particles in turbulent and wave-induced environments

In the present study, a theoretical model is proposed to calculate the settling velocity of solid particles in turbulence generated by both oscillatory and nearly horizontal flow motions. In contrast to other previous models that typically compartmentalize two flow conditions in their studies, this study takes a more comprehensive approach by considering the combined effects of both conditions. Taking into account the influence of drag nonlinearity, virtual mass and Basset history forces, the new theoretical model is formulated. The proposed model is obtained by solving the particle motion in fluid flow under some reasonable assumptions. Accordingly, we obtain a new dimensionless term to better take into account the effect of turbulence anisotropy on the settling velocity and the role of the sediment damping coefficient. Application of this term for other conditions is discussed in the paper. The present model shows satisfactory agreement with a wide range of experimental and numerical data and with different flow conditions found in the existing literature. These data include homogeneous isotropic turbulence with high resolution direct numerical simulations (DNS), turbulent open channel and vertical oscillation.

Lifting relations for a generalized total-energy double-distribution-function kinetic model and their impact on compressible turbulence simulation

Recently, Qi et al. (2022) and Guo et al. (2023) proposed two alternative designs of an efficient mesoscopic method using the total-energy double-distribution-function (DDF) formulation, hereafter referred to as the Qi model and the Guo model. The two models share the same advantage of using only 40 discrete particle velocities to fully reproduce the Navier–Stokes-Fourier (NSF) system. However, the Guo model is based on a more rigorous kinetic consideration, while the Qi model relies on a more general design of the source term to allow for adjustable bulk-to-shear viscosity ratio. In this paper, we derive lifting relations for the Qi model based on two alternative approaches, namely, the Hermite expansion and the Chapman–Enskog expansion, which can be used to construct the boundary and initial conditions for the mesoscopic method. For three-dimensional compressible turbulence simulations, including compressible decaying homogeneous isotropic turbulence and Taylor–Green vortex flows, the derived two sets of lifting relations are applied to the initialization distribution function to study their impacts. Interestingly, for the Qi model, the two sets of lifting relations yield the same results without numerical artifacts, whereas for the Guo model, an appropriate lifting relation must be specified to avoid numerical artifacts resulting from the flow initialization (Qi et al., 2023).

Mixing Of A Starting Turbulent Jet In Homogeneous And Isotropic Turbulence

The mixing of a starting turbulent jet, injected with an injection velocity of U = 18.8 m/s from a tube with a diameter of D = 4mm (Reynolds number Re = 4800, characteristic time scale t* = D/U = 0.21 ms at the nozzle exit), is experimentally investigated in the presence of homogenous and isotropic turbulence (HIT) using Planar Laser-induced Fluorescence (PLIF). The root-mean-square (RMS) velocity of the ambient HIT is 0.38 and 0.63 m/s, corresponding to turbulent Reynolds numbers of Reλ = 190 and 250, respectively. Spatial distributions of mixture fraction and scalar dissipate rate (SDR) were measured within an axial range of y = 0 – 10 D and at different times after the start of injection (ASIs) up to 21.4 t*. Ambient HIT did not alter the location or size of the leading vortex ring in the initial formation stage. However, it was able to enhance the entrainment of ambient fluids in this stage and deform/transport both the leading and the trailing vortex rings at later stages, despite that the ASIs examined were only comparable to its Kolmogorov timescale. Ambient HIT also introduced a significant increase to the fluctuations of the jet penetrations length but only a marginal reduction (less than 10%) to the averaged jet penetration speed. PDFs of the ensemble averaged mixture fraction within the whole jet region showed increased probability of low mixture fraction values – and therefore indicating enhanced entrainment and mixing – in the presence of ambient HIT. Ensemble averaged distributions of the SDR showed no significant alterations by the ambient HIT in the initial formation stage of the leading vortex ring. At later stages, while the high SDRs in the leading vortex ring remained concentrated in a layer along the edge of the vortex ring for Reλ = 0, it became increasingly uniformly distributed with a lower maximum value as the ambient HIT increased. Similar, but less significant, impact of the ambient HIT on the SDR fields was observed in the trailing vortex rings, while no noticeable impact was found in the shear layer close to the nozzle exit.