Homogeneous Isotropic Turbulence Research Articles

Turbulent vortex pair at equilibrium and its interaction with the ground at

A turbulent two-vortex system (T-2VS) is obtained by inserting analytical model wake vortices into very weak homogeneous isotropic turbulence (HIT) and by evolving them in time using large-eddy simulation until a turbulent state at statistical equilibrium is reached. The T-2VS is characterised as follows: circulation distribution of the vortices; energy of the mean and fluctuating fields; energy dissipation rate. It is also verified that essentially the same T-2VS is obtained when varying the initial model or initial HIT perturbation. A wall-resolved simulation of the T-2VS further interacting with a smooth ground is then performed at $Re_\varGamma = 2 \times 10^5$ ; this is $10 \times$ higher than in previous works, which allows us to better capture the high Reynolds number behaviour. The high release height of the T-2VS also ensures a physically correct approach to the ground. The results are compared with the literature and also to what is obtained for the case of non-turbulent vortices interacting with the same ground at the same Reynolds number. The flow topologies are discussed, and significant differences are highlighted regarding the separation of the boundary layer generated at the ground, and the way this secondary vorticity interacts with the primary vortices and makes them decay. The vortex trajectories are also measured, together with their circulation distribution and global circulation evolution, and the differences are discussed.

Modeling the Unsteady Wake of an Impulsively Started Circular Cylinder Using Refined Potential Flow Theory

Abstract Cylindrical structures find usage in many engineering applications including tethered oil drums and engine canisters slung beneath helicopters in flight. The motion of air around such circular cylindrical structures and the helicopter presents interesting phenomena including flow separation, wakes and turbulence. The physics of these are enshrined in the continuity equation and the Navier-Stokes equations. Therefore, their studies are not only important in mathematics and physics, but they are also required for efficient helicopter operations. In practice, reduced-order models of these operations that take in aerodynamics models of the tethered loads are utilized for stability analysis, flight certification and pilot training because of the prohibitive cost of experimentation and computational analyses of these configurations. However, there is a dearth of realistic analytical models of finite cylinder flows because of the Navier-Stokes problem. Classical potential flow theory provides an avenue to develop such models, but the extant gaps in its predictions significantly preclude its usage for engineering applications. Attempting to bridge these gaps, this article introduces refined potential flow theory in which the governing equations and boundary conditions are satisfied. Viscous effects, fluctuations of the mean flow and three-dimensional effects are also incorporated. For characterization, refined potential flow theory is employed on an incompressible flow over an impulsively started circular cylinder for Reynolds numbers and non-dimensional times in the range 30 < Re < 10^4 and 0.2 ≤ T≤ 77, 047 respectively. There is an excellent prediction of 0.209 for the Strouhal number at Re = 3, 900. At this transitional Reynolds number, the harmonics of the Strouhal frequency are also captured, and the characteristic irregular fluctuations at sub-Strouhal frequencies are discernible in the velocity spectra. As the flow becomes more turbulent, these become more pronounced at Re = 9, 500 when the predicted Strouhal number is within 10% of experimental result. In the fully developed stage, spectra analyses of the wake velocity components at some downstream locations also display Kolmogorov’s Five-Thirds law of homogeneous isotropic turbulence. The present model can thus aid the development of reduced-order models of helicopter operations that feature tethered cylindrical loads.

Clustering of settling microswimmers in turbulence

Abstract. Clustering of plankton plays a vital role in several biological activities, including feeding, predation, and mating. Gyrotaxis is one of the mechanisms that induces clustering. A recent study (Candelier et al., 2022) reported a fluid inertial torque acting on a spherical microswimmer, which has the same effect as a gyrotactic torque. In this study, we model plankton cells as microswimmers that are subject to gravitational sedimentation as well as a fluid inertial torque. We use direct numerical simulations to obtain the trajectories of swimmers in homogeneous isotropic turbulence. We also investigate swimmers' clustering using Voronoï analysis. Our findings indicate that fluid inertial torque leads to notable clustering, with its intensity depending on the swimming and settling speeds of swimmers. Using Voronoï analysis, we demonstrate that swimmers preferentially sample downwelling regions where clustering is more prevalent.

Experimental study on modulation of homogeneous isotropic turbulence by bubbles of different sizes

The mechanism of turbulence modulation by bubbles is crucial for understanding and predicting turbulent bubbly flow. In this study, we conducted an experimental investigation of turbulence modulation by bubbles of different sizes in homogeneous isotropic turbulence using two-phase stereo-particle image velocimetry measurement techniques. Two bubble generation methods, electrolysis and porous medium, were employed to generate bubbles in micrometer and millimeter sizes, respectively. The oscillating grid system was utilized to generate homogeneous isotropic turbulence, allowing precise control of turbulent boundary conditions. The ratio of the fluctuating velocities and the comparison between turbulent kinetic energy and average kinetic energy indicated that the generated turbulence was nearly homogeneous and isotropic. With increasing turbulence intensity, micron-sized bubbles transition from suppressing turbulence to enhancing it, while millimeter-sized bubbles exhibit the opposite behavior. Turbulence modulation by millimeter-sized bubbles appears to be nearly isotropic, whereas micrometer-sized bubbles do not exhibit isotropy.

An Implicit Factorized Transformer with Applications to Fast Prediction of Three-dimensional Turbulence

Transformer has achieved remarkable results in various fields, including its application in modeling dynamic systems governed by partial differential equations. However, transformer still face challenges in achieving long-term stable predictions for three-dimensional turbulence. In this paper, we propose an implicit factorized transformer (IFactFormer) model, which enables stable training at greater depths through implicit iteration over factorized attention. IFactFormer is applied to large eddy simulation of three-dimensional homogeneous isotropic turbulence (HIT), and is shown to be more accurate than the FactFormer, Fourier neural operator (FNO), and dynamic Smagorinsky model (DSM) in the prediction of the velocity spectra, probability density functions of velocity increments and vorticity, temporal evolutions of velocity and vorticity root-mean-square value and isosurface of the normalized vorticity. IFactFormer can achieve long-term stable predictions of a series of turbulence statistics in HIT. Furthermore, IFactFormer showcases superior computational efficiency compared to the conventional DSM in large eddy simulation.

Temperature statistics of settling particles in homogeneous isotropic turbulence

Heat transfer of settling particles in homogeneous isotropic turbulence is investigated in one-way coupled Eulerian–Lagrangian direct numerical simulations. In the absence of gravity, our observations highlight that there was a correlation between particles concentration and temperature fronts, which are characterized by high-temperature gradients. The particle clustering is mainly influenced by the preferential concentration and the non-local effect, the latter being the memory of their interaction with the flow along their trajectories. Nevertheless, gravity significantly affects the clustering of particles, which further influence the heat transfer of particles with the fluid. We find that the heat transfer of particles is intricately related to particle dynamic inertia and thermal inertia, which can be quantified by the particle Stokes number St and particle thermal Stokes number Stθ, respectively. In this study, we observe that the variance of particle temperature rate of change is independent of Stθ at small thermal inertia but inversely proportional to Stθ2 for large thermal inertia. In the presence of gravity, the variance converges to Stθ−1 for particles with negligible thermal inertia. Our findings reveal that the second-order structure function of the particle temperature, indicating the enhancement of Lagrangian scalar intermittency, adheres to a power-law behavior with limited thermal inertial particles in the dissipation region, denoted as r2. However, as Stθ increases, the power law of the structure function of the particle temperature is disrupted leading to ‘thermal caustics’. The present findings of the heat transfer behavior of settling particles advance the understanding of gravity effect, and are valuable for the modeling of heat transfer in particle-laden turbulent flows.

Navier-Stokes characteristic boundary conditions for high-enthalpy compressible flows in thermochemical non-equilibrium

This fundamental study presents Navier-Stokes characteristic boundary conditions (NSCBCs) for high-enthalpy hypersonic flows in thermochemical non-equilibrium. In particular, the relevant locally one-dimensional inviscid (LODI) relations are derived within a two-temperature framework for high-enthalpy hypersonic flows undergoing finite-rate thermochemical processes, including air dissociation and vibrational relaxation. Using these LODI relations, a set of NSCBCs are proposed and later demonstrated in canonical test cases, including the interaction of homogeneous isotropic turbulence with a shock wave subject to high-enthalpy thermochemical non-equilibrium effects.

A Method to Simulate an Outwardly Propagating Turbulent Premixed Flame at Constant Pressure

AbstractAn outwardly propagating premixed flame in homogeneous isotropic turbulence at constant pressure is considered one of canonical configurations to study turbulent premixed flames. In this paper, a surface forcing method to prevent the undesirable influence of the boundary-condition-induced backflow on the flame evolution, while maintaining the constant pressure, in the simulation of the outwardly propagating flame is presented. The method is validated for laminar and turbulent flames. The results show that the present method well preserves the characteristics of turbulence and of an outwardly propagating flame, without the undesirable influence of the boundary condition, by feeding the homogeneous turbulence relative to the velocity field induced by the volume expansion due to heat release to the domain in which the flame develops.

Departure from the statistical equilibrium of large scales in forced three-dimensional homogeneous isotropic turbulence

We study the statistically steady states of the forced dissipative three-dimensional homogeneous isotropic turbulence at scales larger than the forcing scale in real separation space. The probability density functions (p.d.f.s) of longitudinal velocity difference at large separations are close to, but deviate from, Gaussian, measured by their non-zero odd parts. The analytical expressions of the third-order longitudinal structure functions derived from the Kármán–Howarth–Monin equation prove that the odd-part p.d.f.s of velocity differences at large separations are small but non-zero. Specifically, when the forcing effect in the displacement space decays exponentially as the displacement tends to infinity, the odd-order longitudinal structure functions have a power-law decay with an exponent of $-$ 2, implying a significant coupling between large and small scales. Under the assumption that forcing controls the large-scale dynamics, we propose a conjugate regime to Kolmogorov's inertial range, independent of the forcing scale, to capture the odd parts of p.d.f.s. Thus, dynamics of large scales departs from the absolute equilibrium, and we can partially recover small-scale information without explicitly resolving small-scale dynamics. The departure from the statistical equilibrium is quantified and found to be viscosity-independent. Even though this departure is small, it is significant and should be considered when studying the large scales of the forced three-dimensional homogeneous isotropic turbulence.

Universal relative scaling of longitudinal structure functions in shear-dominated turbulence

Shear significantly influences turbulence in the energy-containing range of shear-dominated flows, and the longitudinal structure functions do not have a universal form as they do in homogeneous isotropic turbulence. Despite this, the relative scaling of structure functions exhibits universal sub-Gaussian behaviour in shear-dominated flows, in particular for turbulent boundary layers, channels and Taylor–Couette flows. Our investigation of a turbulent vertical buoyancy layer at $Pr = 0.71$ using direct numerical simulation shows this universality even in moderate-Reynolds-number buoyancy-driven but shear-dominated boundary layers. It is demonstrated that the universality is related to the energy density of the eddies, which attains a hierarchical equilibrium in the energy-containing range of shear-dominated turbulence. We conjecture that the universal sub-Gaussian behaviour of the energy density of the energy-containing range, which was considered to be non-trivial in prior studies, is related to the universal anomalous scaling exponents of the inertial subrange turbulence. Based on this conjecture, we propose a hypothesis that relates large-scale eddies and the intermittent dissipation field in shear-dominated turbulence, highlighting a relationship between large and small scales. A phenomenological model is also developed to predict the scaling, which is verified using data from a turbulent boundary layer, half-channel and vertical buoyancy layer at friction Reynolds numbers spanning four orders of magnitude. Excellent agreement is observed.

Temporal vortex tracking in isotropic turbulence

An algorithm to allow the temporal vortex tracking of the intense vorticity structures is introduced and applied to direct numerical simulations (DNS) of homogeneous isotropic turbulence (HIT). The algorithm allows the identification of the individual life events of the worms, including events of splits into several structures and of merging between different structures. Results obtained from a relatively small DNS, with Taylor based Reynolds number equal to Reλ = 54, for a simulation lasting more than 6 integral time scales, shows that the complement of the CDF of the life time of the IVS is well approximated by an exponential distribution function. The mean life time of the IVS which is approximately 3.6τη , where τη is the Kolmogorov time, and the birth rate of these structures is close to 4.5/τη (number of new structures per Kolmogorov time). A population model is built to predict the number of existing structures at any time, which gives a number of existing IVS equal to ≈ 16, which agrees well with the results from the DNS.

A robust phase-field method for two-phase flows on unstructured grids

A phase-field method for unstructured grids that is accurate, conservative, and robust is proposed in this work. The proposed method also results in bounded transport of volume fraction, and the interface thickness adapts automatically to local grid size. In addition to this, we present a novel formulation for two-phase flows on collocated grids that is provably energy stable, which is a critical feature for robust simulations of two-phase turbulent flows. The proposed method is first evaluated on canonical test cases, including scenarios like drop advection and drop in a shear flow. The accuracy of the proposed method in the presence of grid transitions, which is important for simulations in complex geometries using unstructured grids, is also evaluated. We assess the robustness of our method by performing simulations of a drop in homogeneous isotropic turbulence at infinite Reynolds number with varying density ratios. Furthermore, as a step toward verification and validation, simulations of test cases in complex geometries and conditions, such as, a damped surface wave, infinitesimal interface distortion on a liquid jet, and the atomization of a liquid jet from the engine combustion network's Spray A nozzle are presented. These simulations illustrate the accuracy, robustness, and applicability of the proposed method in various kinds of complex two-phase flow applications of engineering interest.

Detailed characterization of extreme clustering at near-contact scales in isotropic turbulence

Recent measurements of inertial particles in isotropic turbulence (Hammond & Meng, J. Fluid Mech., vol. 921, 2021, A16) revealed surprising extreme clustering of particles at near-contact separations $(r)$ , whereby the radial distribution function, $g(r)$ , grows from $O(10)$ to $O(10^3)$ with a $(r/a)^{-6}$ scaling (where $a$ is the particle radius), and a surprising upturn of the mean inward particle-pair relative velocity (MIRV). Hydrodynamic interactions (HIs) were proposed to explain the extreme clustering, but despite predicting the correct scaling $(r/a)^{-6}$ , the HI theory underpredicted $g(r)$ by at least two orders of magnitude (Bragg et al., J. Fluid Mech., vol. 933, 2022, A31). To further understand the extreme clustering phenomenon and the relevance of HI, we characterize $g(r)$ and particle-pair kinematics for Stokes numbers $0.07 \leq St \leq 3.68$ in a homogeneous isotropic turbulence chamber using three-dimensional (3-D) particle tracking resolved to near–contact. A drift–diffusion equation governing $g(r)$ is presented to investigate the kinematic mechanisms of particle pairs. Measurements in all 24 conditions show that when $r/a\lessapprox 20$ , extreme clustering consistently occurs, scaling as $g(r) \sim (r/a)^{-k}$ with $4.5 \leq k \leq 7.6$ , which increases with $St$ . Here $g(r)$ varies with $St$ , particle size, density and polydispersity in ways that HI cannot explain. The extreme clustering region features an inward drift contributed by particle-pair turbophoresis and an inward radial relative acceleration. The latter indicates an interparticle attractive force at these separations that HI also cannot explain. The MIRV turns upward when approaching the extreme clustering region, opposite to direct numerical simulation predictions. These observations further support our previous assessment that extreme clustering arises from particle–particle interactions, but HI is not the main mechanism.

Generation of nearly homogeneous isotropic turbulence using a novel oscillating grid system

The creation and control of homogeneous isotropic turbulence in laboratory settings is of practical significance. In this paper, we propose a novel design for a turbulence generation apparatus driven by two pairs of oscillating grids, with the object of striking a balance between the stochastics of momentum injection and the flexibility of control of turbulent kinematic energy (TKE). To validate the performance of this system, tomographic particle image velocimetry measurements are conducted. Based on Reynolds decomposition, we characterize the secondary flow pattern and examine the impact of the mean velocity on the TKE transport. Analysis of one- and two-point statistics indicates that the generated turbulence exhibits reasonable approximations of homogeneity and isotropy. The TKE and the dissipation rate are linearly proportional to the grid oscillation frequency, with the Taylor-scale Reynolds number Reλ ranging from 47 to 154, suggesting a convenient control of turbulence with low-to-moderate intensity.

A novel vector potential random flow generation method for synthesizing divergence-free homogeneous isotropic turbulence with arbitrary spectra

A novel method, known as the vector potential random flow generation (VPRFG) method, is introduced for synthesizing divergence-free homogeneous isotropic turbulence with arbitrary spectra. First, the proposed approach employs the random-flow-generation-based method to create a vector potential field. Subsequent application of the curl operation to this field produces a turbulent flow that inherently satisfies the divergence-free condition. In the formulas of the proposed method, we explicitly impose arbitrary homogeneous isotropic three-dimensional spatial cross-spectral density (CSD) and Taylor's frozen hypothesis. This ensures that the generated turbulence conforms to prescribed statistical characteristics, including energy spectra, one-dimensional spatial power spectral density (PSD), temporal PSD, spatial coherence function, turbulent kinetic energy, and Reynolds stress. Additionally, the theoretical accuracy of the proposed method is validated through numerical examples, employing the von Kármán energy spectrum as the target value. Finally, large eddy simulations of homogeneous isotropic temporal-decaying box turbulence generated by the VPRFG method demonstrate a close alignment with the experimental results.

Clustering, rotation, and swirl of inertial particles in turbulent channel flow

Clustering dynamics of inertial particles in turbulent channel flow are studied via tessellation-based analysis of high-fidelity simulation data at Reτ≈230 with various values of mass loading (10%−100%) and the Stokes number (St+=[1−60]). We then characterise the solenoidal, rotational, and swirling motions of clusters by computing the probability density functions (PDFs) of the divergence, curl, and helicity of the particle velocity, as well as their dependence on wall-normal distance, using the methods of Oujia et al. (2020); Maurel–Oujia et al. (2023). Particle inertia gives heavier tails to the PDFs of divergence and curl, suggesting enhanced intermittency in the convergence/divergence of clusters, and in their rotational motions. The fluctuations of the divergence and curl are most intense in the buffer layer, due to the stronger fluctuations of fluid velocity there. Similarities are identified between the cluster dynamics in the logarithmic region and those in homogeneous isotropic turbulence, including the dependence of divergence, curl, and helicity on Stokes number. The effect of increasing mass loading on cluster dynamics is relatively small except in the viscous sublayer, where attenuation of clustering, rotation, and swirling motions are observed. The effect of increasing Stokes number on the viscous sublayer is different, resulting in more intense convergence/divergence and rotation of particle clusters, as the particles become more independent of the carrier fluid.

Turbulence ingestion noise generation in rotating blades

The interaction between ingested turbulence and rotating blades is a key source of broadband noise in engineering applications. In this paper, a far-field noise model accounting for source correlation across the span of the blade and between blades is developed and applied to the study of homogeneous isotropic turbulence ingestion by a model cooling fan, wind turbine and aircraft propeller. The widely used theory of Amiet is revisited and it is shown that previous works produce conflicting results when attempting to account for blade-to-blade correlation. Central to the model is the calculation of the time between blade chops of the same turbulent eddy as heard by the observer. In this paper it is shown that, when derived correctly, Amiet's theory accounts for correlated sources between blades and, thus, can predict haystacking tones. Comparisons with the new rotational formulation and with experimental data enable us to show that Amiet's theory can be used to accurately predict turbulence ingestion noise from open rotors. In particular, it is found that the infinite-span assumption in strip theory and the neglect of correlation effects across the span do not undermine the accuracy of this theory. This is of great importance because, unlike Amiet's theory, models which treat rotational effects and source correlation exactly are expensive to evaluate routinely at high frequencies due to the slow convergence of infinite series with Bessel functions.

Effects of electrostatic interaction on clustering and collision of bidispersed inertial particles in homogeneous and isotropic turbulence

In sandstorms and thunderclouds, turbulence-induced collisions between solid particles and ice crystals lead to inevitable triboelectrification. The charge segregation is usually size dependent, with small particles charged negatively and large particles charged positively. In this work, we perform numerical simulations to study the influence of charge segregation on the dynamics of bidispersed inertial particles in turbulence. Direct numerical simulations of homogeneous isotropic turbulence are performed with the Taylor Reynolds number ${Re}_{\lambda }=147.5$ , while particles are subjected to both electrostatic interactions and fluid drag, with Stokes numbers of 1 and 10 for small and large particles, respectively. Coulomb repulsion/attraction is shown to effectively inhibit/enhance particle clustering within a short range. Besides, the mean relative velocity between same-size particles is found to rise as the particle charge increases because of the exclusion of low-velocity pairs, while the relative velocity between different-size particles is almost unaffected, emphasizing the dominant roles of differential inertia. The mean Coulomb-turbulence parameter, ${Ct}_0$ , is then defined to characterize the competition between the Coulomb potential energy and the mean relative kinetic energy. In addition, a model is proposed to quantify the rate at which charged particles approach each other and to capture the transition of the particle relative motion from the turbulence-dominated regime to the electrostatic-dominated regime. Finally, the probability distribution function of the approach rate between particle pairs is examined, and its dependence on the Coulomb force is further discussed using the extended Coulomb-turbulence parameter.

Machine learning-based vorticity evolution and super-resolution of homogeneous isotropic turbulence using wavelet projection

A wavelet-based machine learning method is proposed for predicting the time evolution of homogeneous isotropic turbulence where vortex tubes are preserved. Three-dimensional convolutional neural networks and long short-term memory are trained with a time series of direct numerical simulation (DNS) data of homogeneous isotropic turbulence at the Taylor microscale Reynolds number 92. The predicted results are assessed by using the flow visualization of vorticity and statistics, e.g., probability density functions of vorticity and enstrophy spectra. It is found that the predicted results are in good agreement with DNS results. The small-scale flow topology considering the second and the third invariants of the velocity gradient tensor likewise shows an approximate match. Furthermore, we apply the pre-trained neural networks to coarse-grained vorticity data using super-resolution. It is shown that the super-resolved flow field well agrees with the reference DNS field, and thus small-scale information and vortex tubes are well regenerated.

Estimating the turbulent kinetic energy dissipation rate from one-dimensional velocity measurements in time

Abstract. The turbulent kinetic energy dissipation rate is one of the most important quantities characterizing turbulence. Experimental studies of a turbulent flow in terms of the energy dissipation rate often rely on one-dimensional measurements of the flow velocity fluctuations in time. In this work, we first use direct numerical simulation of stationary homogeneous isotropic turbulence at Taylor-scale Reynolds numbers 74≤Rλ≤321 to evaluate different methods for inferring the energy dissipation rate from one-dimensional velocity time records. We systematically investigate the influence of the finite turbulence intensity and the misalignment between the mean flow direction and the measurement probe, and we derive analytical expressions for the errors associated with these parameters. We further investigate how statistical averaging for different time windows affects the results as a function of Rλ. The results are then combined with Max Planck Variable Density Turbulence Tunnel hot-wire measurements at 147≤Rλ≤5864 to investigate flow conditions similar to those in the atmospheric boundary layer. Finally, practical guidelines for estimating the energy dissipation rate from one-dimensional atmospheric velocity records are given.