Fractal modelling of turbulent mixing

The aim of this work is to propose a new model for turbulent flows, called the fractal model (FM), applicable both in a Reynolds averaged Navier–Stokes (RANS) and a large-eddy simulation (LES) formulation, with the ultimate goal of applying it to simulate turbulent combustion irrelevant of its mode (premixed or non-premixed). The model is able to turn itself off in the laminar zones of the flow, and in particular near walls. It is based on the fractal theory. It describes the physics of the smaller spatial scales and therefore represents a small-scales model.

FM describes the physics of the small scales of turbulence based on the phenomenological concept of vortex cascade and on the self-similar behaviour of turbulence in the inertial range. Such a model is used in each cell of a numerical calculation. A characteristic length Δ is associated to each cell, and the local energy u ³ _Δ/Δ is distributed over a certain number of eddies, which depends on the local Reynolds number Re _Δ. Each vortex of the cascade generates N _c vortices; the recursive process of vortex generation terminates at the dissipative scale level, i.e. when the eddy Reynolds number is equal to one. FM is also able to estimate the volume fraction occupied by the dissipative fine structures of turbulence; this quantity is critical in reactive turbulent flows.

The physics of small scales is summarized by a turbulent ‘viscosity’ μ_t, to be added to the molecular one. μ_t is zero where the flow is laminar and, in particular, goes to zero at solid walls. Assuming μ_t to be isotropic, FM is applicable in a RANS formulation (IFM, isotropic fractal model). The model can be extended to the anisotropic case (AFM, anisotropic fractal model) and therefore used to close the transport equations in an LES approach. In the present paper, the model (IFM) is used in a RANS approach and is validated through a test case studied experimentally by Johnson and Bennett, and numerically (with LES) by Akselvoll and Moin. The results obtained are in good agreement both with the experimental and the numerical ones. Other tests are being performed.     

Vortex rings moving into turbulence

We report on an experimental study of turbulent vortex rings injected with velocity U _v0 into a grid-generated turbulent flow (with RMS streamwise velocity u _*) and followed relative to the mean flow. The initial Reynolds number of the vortices varies from 4500 to 11,500. The turbulence was characterised by an intensity I_t =u _*/U _v0, which varied over the range 0<I_t <0.03. A mathematical model based on a stochastic model of the vortex core is developed to explain and interpret the results. The vortex radius grows diffusively in time with the rate of increase of the square of the vortex radius increasing linearly with I_t . As the vortices grow, they slow down sufficiently rapidly in a manner that they penetrate a finite distance into the turbulence. The vortex velocity, averaged over many experiments, showed an initial t ^?1 decay, consistent with Maxworthy’s experiments. The analysis and experiments show that such vortices ultimately only move a finite distance from their point of generation and this distance varies inversely with I_t .     

The self-similarity of wall-bounded temporally accelerating turbulent flows

From the study of viscous flow it is known that certain time-dependent laminar problems, such as the impulsively started flat plate and the diffusion of a vortex sheet, possess self-similar solutions. Previous studies of turbulent channel and pipe flows accelerating between two steady states have shown that the flow field evolves in three distinct stages. Furthermore, recent direct numerical simulations have shown that the perturbation velocity, i.e. the surplus velocity from the initial value, in an impulsively accelerating turbulent channel and pipe flow also possesses a self-similar distribution during the initial stage. In here, these results are developed analytically and it is shown that accelerating flows in which the centreline velocity develops as U^∧_c(t) = U₀(t/t₀)^m will possess a self-similar velocity distribution during the initial stage. The displacement thickness of the perturbation velocity is shown to be dependent only on the type of acceleration, and not on the initial Reynolds number, the acceleration rate or the change in Reynolds number. The derived formulas are verified with good agreement against measurements performed in a linearly accelerating turbulent pipe flow and with data from channel flow simulations.     

Stabilization mechanism of turbulent premixed flames in strongly swirled flows

The stabilization of turbulent premixed flames in strongly swirled flows undergoing vortex breakdown is studied in the case of the ALSTOM En-Vironmental (EV) double cone burner using a simple one-dimensional boundary layer type model and computational fluid dynamics, mainly at the level of large-eddy simulation. The analysis shows that, due to flame curvature effects, the flame speed on the combustor axis is 2 D _t/R _F lower than the turbulent burning rate, where D _t is a characteristic turbulent diffusion coefficient and R _F the flame radius of curvature. Flame propagation with negative speed observed in the experiments, i.e. the flame completely embedded in the central recirculation zone on the symmetry axis, is explained with the one-dimensional model as caused by the factor 2 D _t/R _F being larger than the characteristic turbulent burning rate. A peculiar sudden displacement of the flame anchoring location deep into the burner, which takes place experimentally at a critical value of the equivalence ratio, cannot however be explained with the present one-dimensional approach due to the modelling assumptions. The mathematical analysis is supported in this case with large-eddy simulation which can accurately reproduce the flame behaviour across the full operating range. It is finally shown that steady RANS methods cannot cope with the problem due to their inability to correctly predict the velocity flowfield in this burner.     

Diffusion and transport of a magnetic field in a turbulent medium with helicity

An analysis is made of the relation between accurate formulas for the coefficients of turbulent diffusion D _T and the alpha effect α _T for a magnetic field in the Lagrange and Euler representations. It is shown that the quadratic term with respect to α _T in the diffusion coefficient derived by Moffatt and Kraichnan is incorrect and should be dropped. First, a numerical solution of the nonlinear equation (DIA equation) for the Green function is presented, describing the transport of a magnetic field for the case of incompressible, uniform, isotropic, steady-state turbulence possessing helicity. These solutions are used to calculate the steady-state coefficients D _T and α _T for various values of the parameters ξ ₀=u ₀ σ ₀/R ₀, a=H ₀/u ₀ ² p ₀, σ ₀/σ ₁, and R ₀/R ₁, where u ₀, σ ₀, and R ₀ are the characteristic velocity, lifetime, and scale of the turbulent pulsations, and H ₀, σ ₁, and R ₁ are similar values describing the helicity of the medium h(1,2)=〈u(1)· (∇×u(2))〉, and the parameter α characterizes the degree of helicity. The DIA values of D _T and α _T and the self-consistent values of these quantities calculated using the Green tensor in the diffusion approximation are in qualitative agreement. It is shown that the coefficient of turbulent diffusion is always positive for all the types of turbulence studied. Nonsteady-state values of D _T(t) and α _T(t) calculated by a self-consistent method are given. Zh. éksp. Teor. Fiz. 112, 1312–1331 (October 1997)     

Comparison of methods for calculating turbulent diffusion coefficients

Different methods for calculating the turbulent diffusion coefficient D _T of a passive scalar impurity in an infinite homogeneous isotropic stationary turbulent medium are examined. The values of D _T calculated by these methods are compared for two limiting types of turbulence, viz., turbulence with a δ-function spectrum and turbulence with a Kolmogorov-type spectrum. The temporal dependence of the velocity correlators is assumed to be exponential. It is shown that the most accurate method is based on the use of the solution of the nonlinear equation for the averaged Green’s function with consideration of the contribution from the four-point turbulent velocity correlators. A comparison with the results of other methods that are simpler from the mathematical standpoint shows that some of them also permit the calculation of D _T with relatively good accuracy. Zh. éksp. Teor. Fiz. 111, 871–881 (March 1997)     

Fractal modelling of turbulent combustion

In a previous paper we proposed a new model for turbulent flows, called the fractal model (FM), which is applicable both to RANS and LES formulations. Here, the model is extended to the reactive case with the goal of simulating turbulent flames, both premixed and non-premixed.

FM is a subgrid model that describes the physics of the small scales of turbulence building on the phenomenological concept of vortex cascade and on fractal theory. The physics of the small scales is summarized by a turbulent ‘viscosity’ μ_t, to be added to the molecular one. μ_t is zero where the flow is laminar and, in particular, goes to zero at solid walls.

The fundamental assumption in treating combustion in this work is that chemical reactions take place only at the dissipative scales of turbulence, i.e. near the so-called ‘fine structures’ (the eddy dissipation concept). FM predicts the growth of dissipative scales due to heat release; therefore, it enables a local DNS in the hot regions of the flow where the dissipative scale may grow up to the cell dimension. FM can also estimate the volume fraction γ* occupied by the ‘fine structures’; this quantity is critical for modelling the reaction rate, and therefore the source terms in the energy and species equations. FM can also estimate the local surface of the reactive ‘fine structures’, that is, the local flame front area. It also takes into account, although in approximate manner, the formation of islands of unburnt mixture. In this paper, the model (in the isotropic formulation (IFM)) is used in conjunction with a time-dependent LES (but with the limitations of an isotropic model) approach and is validated through a three-dimensional axisymmetric diffusion flame studied experimentally by Correa and Gulati and numerically by many researchers. The time-dependent results obtained are in good agreement with the experiments. Moreover, the IFM solution offers a possible explanation for the stabilization process of this flame, which undergoes local stretching of the order of 46 000 s^?1.     

Vortex dynamics of a trapezoidal bluff body placed inside a circular pipe

The influence of Reynolds number and blockage ratio on the vortex dynamics of a trapezoidal bluff body placed inside a circular pipe is studied experimentally and numerically. Low aspect ratio, high blockage ratio, curved end conditions (junction of pipe and bluff body), axisymmetric upstream flow with shear and turbulence are some of the intrinsic features of this class of bluff body flows which have been scarcely addressed in the literature. A large range (200:200,000) of Reynolds number (Re_D) is covered in this study, encompassing all the three pipe flow regimes (laminar, transition and turbulent). Four different flow regimes are defined based on the distinct features of Strouhal number (St)–Re_D relation: steady, laminar irregular, transition and turbulent. The wake in the steady regime is stationary with no oscillations in the shear layer. The laminar regime is termed as irregular owing to irregular vortex shedding. The vortex shedding in this regime is observed to be symmetric. The emergence of separation bubble downstream of the bluff body on either side is another interesting feature of this regime, which is further observed to be symmetric. Two pairs of mean streamwise vortices are noticed in the near-wake regime, which are termed as reverse dipole-type wake topology. Beyond the irregular laminar regime, the Strouhal number falls gradually and vortex shedding becomes more periodic. This regime is named transition and occurs close to the Reynolds number at which transition to turbulence takes place in a fully developed pipe. The turbulent regime is characterised by a nearly constant Strouhal number. Typical Karman-type vortex shedding is noticed in this regime. The convection velocity, wake width formation length and irrecoverable pressure loss are quantified to highlight the influence of blockage ratio. These results will be useful to develop basic understanding of vortex dynamics of confined bluff body flow for several practical applications.     

Level set simulations of turbulent thermonuclear deflagration in degenerate carbon and oxygen

We study the dynamics of thermonuclear flames propagating in fuel stirred by stochastic forcing. The fuel consists of carbon and oxygen in a state which is encountered in white dwarfs close to the Chandrasekhar limit. The level set method is applied to represent the flame fronts numerically. The computational domain for the numerical simulations is cubic, and periodic boundary conditions are imposed. The goal is the development of a suitable flame speed model for the small-scale dynamics of turbulent deflagration in thermonuclear supernovae. Because the burning process in a supernova explosion is transient and spatially inhomogeneous, the localized determination of subgrid scale closure parameters is essential. We formulate a semi-localized model based on the dynamical equation for the subgrid scale turbulence energy k _sgs. The turbulent flame speed s _t is of the order √2k _sgs. In particular, the subgrid scale model features a dynamic procedure for the calculation of the turbulent energy transfer from resolved toward subgrid scales, which has been successfully applied to combustion problems in engineering. The options of either including or suppressing inverse energy transfer in the turbulence production term are compared. In combination with the piece-wise parabolic method for the hydrodynamics, our results favour the latter option. Moreover, different choices for the constant of proportionality in the asymptotic flame speed relation, s _t∝√2k _sgs, are investigated.     

Propagation of orifice- and nozzle-generated vortex rings in air

The effect of the exit geometry of a vortex ring generator was studied experimentally. Two types of exit geometries were chosen: an orifice and a nozzle. Vortex rings were generated by pushing a solenoid-valve-controlled, pressurized-air jet through the circular opening of the orifice or nozzle. Experiments were performed over a wide range of initial Reynolds number (450≤Re≤4580) and length-to-diameter ratio (0.7≤L/D≤7.0) of the air jet. The exit geometry was found to significantly influence the entire course of propagation of the vortex ring. The orifice-generated vortex ring had superior characteristics to that produced by the nozzle under the same conditions. The vorticity generated along the wall in the orifice exit plane had a negligible effect on the circulation of the vortex ring within the specified range of Reynolds number. Compared to the nozzle-generated vortex ring, the orifice-generated ring showed reduced initial vorticity losses and less diffusive entrainment of ambient fluid. The vortex rings produced by the orifice attained more circulation, less entrainment of ambient fluid and hence rapidly propagated through longer distances in comparison to the nozzle-generated rings.     

Flow and mixing fields of turbulent bluff-body jets and flames

The mean structure of turbulent bluff-body jets and flames is presented. Measurements of the flow and mixing fields are compared with predictions made using standard turbulence models. It is found that two vortices exist in the recirculation zone; an outer vortex close to the air coflow and an inner vortex between the outer vortex and the jet. The inner vortex is found to shift downstream with increasing jet momentum flux relative to the coflow momentum flux and gradually loses its circulation pattern. The momentum flux ratio of the jet to the coflow in isothermal flows is found to be the only scaling parameter for the flow field structure. Three mixing layers are identified in the recirculation zone. Numerical simulations using the standard k-? and Reynolds stress turbulence models underpredict the length of the recirculation zone. A simple modification to the C₁ constant in the dissipation transport equation fixes this deficiency and gives better predictions of the flow and mixing fields. The mixed-is-burnt combustion model is found to be adequate for simulating the temperature and mixing field in the recirculation zone of the bluff-body flames.     

Observation of vortex packets in direct numerical simulation of fully turbulent channel flow

A number of experimental studies have inferred the existence of packets of inclined, hairpinlike vortices in wall turbulence on the basis of observations made in two-dimensional x−y planes using visualization and particle image velocimetry (PIV). However, there are very few observations of hairpins in existing three-dimensional studies made using direct numerical simulation (DNS), and no such study claims to have revealed packets. We demonstrate, for the first time, the existence of hairpin vortex packets in DNS of turbulent flow. The vortex packet structure found in the present study at low Reynolds number,Re _t=300, is consistent with and substantiates the observations and the results from twodimensional PIV measurements at higher Reynolds numbers in channel, pipe and boundary layer flows. Thus, the evidence supports the view that vortex packets are a universal feature of wall turbulence, independent of effects due to boundary layer trips or critical conditions in the aforementioned numerical studies. Visualization of the DNS velocity field and vortices also shows the close association of hairpin packets with long low-momentum streaks and the regions of high Reynolds shear stress.     

Homogeneous Isotropic Turbulence: Geometric and Isometry Properties of the 2-point Velocity Correlation Tensor

The emphasis of this review is both the geometric realization of the 2-point velocity correlation tensor field B_ij (x,x′,t) and isometries of the correlation space K³ equipped with a (pseudo-) Riemannian metrics ds²(t) generated by the tensor field. The special form of this tensor field for homogeneous isotropic turbulence specifies ds²(t) as the semi-reducible pseudo-Riemannian metric. This construction presents the template for the application of methods of Riemannian geometry in turbulence to observe, in particular, the deformation of length scales of turbulent motion localized within a singled out fluid volume of the flow in time. This also allows to use common concepts and technics of Lagrangian mechanics for a Lagrangian system (M^t, ds²(t)), M^t ? K³. Here the metric ds²(t), whose components are the correlation functions, evolves due to the von Kármán-Howarth equation. We review the explicit geometric realization of ds²(t) in K³ and present symmetries (or isometric motions in K³) of the metric ds²(t) which coincide with the sliding deformation of a surface arising under the geometric realization of ds²(t). We expose the fine structure of a Lie algebra associated with this symmetry transformation and construct the basis of differential invariants. Minimal generating set of differential invariants is derived. We demonstrate that the well-known Taylor microscale λ_g is a second-order differential invariant and show how λ_g can be obtained by the minimal generating set of differential invariants and the operators of invariant differentiation. Finally, we establish that there exists a nontrivial central extension of the infinite-dimensional Lie algebra constructed wherein the central charge is defined by the same bilinear skew-symmetric form c as for the Witt algebra which measures the number of internal degrees of freedom of the system. For turbulence, we give the asymptotic expansion of the transversal correlation function for the geometry generated by a quadratic form.     

Identification of lagrangian coherent structures in the turbulent boundary layer

Using Finite-Time Lyapunov Exponents (FTLE) method, Lagrangian coherent structures (LCSs) in a fully developed flat-plate turbulent boundary layer are successfully identified from a two-dimensional (2D) velocity field obtained by time-resolved 2D PIV measurement. The typical LCSs in the turbulent boundary layer are hairpin-like structures, which are characterized as legs of quasi-streamwise vortices extending deep into the near wall region with an inclination angle θ to the wall, and heads of the transverse vortex tube located in the outer region. Statistical analysis on the characteristic shape of typical LCS reveals that the probability density distribution of θ accords well with t-distribution in the near wall region, but presents a bimodal distribution with two peaks in the outer region, corresponding to the hairpin head and the hairpin neck, respectively. Spatial correlation analysis of FTLE field is implemented to get the ensemble-averaged inclination angle θ _R of typical LCS. θ _R first increases and then decreases along the wall-normal direction, similar to that of the mean value of θ. Moreover, the most probable value of θ saturates at y ⁺=100 with the maximum value of about 24°, suggesting that the most likely position where hairpins transit from the neck to the head is located around y ⁺=100. The ensemble- averaged convection velocity U _c of typical LCS is finally calculated from temporal-spatial correlation analysis of FTLE field. It is found that the wall-normal profile of the convection velocity U _c(y) accords well with the local mean velocity profile U(y) beyond the buffer layer, evidencing that the downstream convection of hairpins determines the transportation properties of the turbulent boundary layer in the log-region and beyond. Supported by the National Natural Science Foundation of China (Grant Nos. 10425207 and 10832001)     

Dynamic magnetic permeability of a thin, high-T c superconducting wafer

The field dependence of the vibrational contribution to the dynamic magnetic permeability μ _V(H) is calculated for a thin (of thickness d∼λ) high-T _c superconducting wafer in a magnetic field parallel to the surface. The resulting curves are plotted on the basis of an exact numerical analysis of the vortex structures both for the thermodynamic-equilibrium vortex lattice and in the presence of pinning forces and the Bean-Livingston surface barrier. It is shown that the μ _V(H) curves are highly sensitive to the size factor (d/λ) and exhibit abrupt changes corresponding to a change in the number of vortex rows. The equilibrium μ _V(H) curve is found to be similar in its general behavior and absolute value (obtained with allowance for the distribution of grain sizes and with appropriate values of λ and ϰ) to the experimental μ _V(H) curve plotted at nitrogen temperature for fine-grained YBa₂Cu₃O_x with grain diameters 〈D〉∼λ in an increasing magnetic field. It is established that the main cause of the experimentally observed irreversible behavior of the μ _V(H) curves during cyclic variation of the applied magnetic field is the existence of a surface barrier to the exit of vortices from the superconductor. The lower limit H _min(B) of stability of the mixed state in the presence of an ideal surface barrier in a thin, high-T _c superconducting wafer (d∼λ) is determined, along with the range of the vortex state (H _max-H _min) for a fixed number of vortices in micrometer-size grains of the investigated YBaCuO samples. Fiz. Tverd. Tela (St. Petersburg) 39, 1943–1947 (November 1997)     

Diffusion-controlled first-order surface reaction in turbulent flow

The effect of advective diffusion on the rate of reactant consumption by a first-order surface reaction is analyzed in the fast-reaction limit. The decay of reactant concentration is described by the function n(t) ～ exp(?λt). In the limit of well-developed turbulence, the scaling estimates λ ～ L ^?1κ^3/4μ^1/4 and λ ～ fκ^3/4μ^1/4 are obtained, respectively, for a confined flow with characteristic length scale L and in the case when the reactants are contained near the surface by an external field with potential U/T = fx, where κ is molecular diffusivity and μ is the constant parameter in the eddy diffusivity D _adv = μx ⁴ (x is distance to the wall). The coefficients in the scaling laws are evaluated by a variational method and by numerical solution of the governing equations.     

Viscous dark cosmology with account of quantum effects

The analytic properties of the energy density ρ(t) of the cosmic fluid and the Hubble parameter H(t) are investigated close to the future singularity t=t _s assuming different forms for the equation of state. First, it is shown that the inclusion of quantum effects coming from the conformal anomaly modifies the singularity. Thereafter, we consider the effect coming from the bulk viscosity in the fluid. The viscosity tends to reduce the magnitude of t _s, but does not alter the singularity itself (the exponent). The main emphasis is laid on the simple case when the equation of state is p=w ρ, with w a constant.     

Large-scale amplification of a magnetic field by turbulent helicity fluctuations

In this paper the procedure of large-scale averaging of the magnetic-field diffusion equation with the α-term curlα(r,t)B(r,t) is used to show that a nonuniform distribution of the turbulent helicity fluctuations (more precisely, the fluctuations of the coefficient α) with a zero average value gives rise to large-scale amplification of the initial magnetic field. A detailed study is carried out of the dependence of the resulting large-scale α effect on the characteristics of the correlator 〈〈α(r, t)α(r″,t″)〉〉 in a rotating medium with a nonuniform distribution of the angular velocity ω=ω(ρ,z) (ρ is the distance for the rotation axis z). The effect of helicity fluctuations and the diffusion coefficient on the turbulent diffusion process is also investigated. Zh. éksp. Teor. Fiz. 116, 85–104 (July 1999)     

On mean flow universality of turbulent wall flows. II. Asymptotic flow analysis

Understanding of the structure of turbulent flows at extreme Reynolds numbers (Re) is relevant because of several reasons: almost all turbulence theories are only valid in the high Re limit, and most turbulent flows of practical relevance are characterized by very high Re. Specific questions about wall-bounded turbulent flows at extreme Re concern the asymptotic laws of the mean velocity and turbulence statistics, their universality, the convergence of statistics towards their asymptotic profiles, and the overall physical flow organization. In extension of recent studies focusing on the mean flow at moderate and relatively high Re, the latter questions are addressed with respect to three canonical wall-bounded flows (channel flow, pipe flow, and the zero-pressure gradient turbulent boundary layer). Main results reported here are the asymptotic logarithmic law for the mean velocity and corresponding scale-separation laws for bulk flow properties, the Reynolds shear stress, the turbulence production and turbulent viscosity. A scaling analysis indicates that the establishment of a self-similar turbulence state is the condition for the development of a strict logarithmic velocity profile. The resulting overall physical flow structure at extreme Re is discussed.     

On the Large Time Asymptotics¶of Decaying Burgers Turbulence

The decay of Burgers turbulence with compactly supported Gaussian “white noise” initial conditions is studied in the limit of vanishing viscosity and large time. Probability distribution functions and moments for both velocities and velocity differences are computed exactly, together with the “time-like” structure functions T _n (t,τ)≡< (u(t+τ) -u(t)) ⁿ >. The analysis of the answers reveals both well known features of Burgers turbulence, such as the presence of dissipative anomaly, the extreme anomalous scaling of the velocity structure functions and self similarity of the statistics of the velocity field, and new features such as the extreme anomalous scaling of the “time-like” structure functions and the non-existence of a global inertial scale due to multiscaling of the Burgers velocity field. We also observe that all the results can be recovered using the one point probability distribution function of the shock strength and discuss the implications of this fact for Burgers turbulence in general. Received: 4 October 1999 / Accepted: 4 February 2000