Simulations of Three-Dimensional Turbulent Mixing for Schmidt Numbers of the Order 1000

We report basic results from new numerical simulations of passive scalar mixing at Schmidt numbers (Sc) of the order of 1000 in isotropic turbulence. The required high grid-resolution is made possible by simulating turbulence at very low Reynolds numbers, which nevertheless possesses universality in dissipative scales of motion. The results obtained are qualitatively consistent with those based on another study (Yeung et al., Phys. Fluids 14 (2002) 4178-4191) with a less extended Schmidt number range and a higher Reynolds number. In the stationary state maintained by a uniform mean scalar gradient, the scalar variance increases slightly with Sc but scalar dissipation is nearly constant. As the Schmidt number increases, there is an increasing trend towards k ^?1 scaling predicted by Batchelor (Batchelor, J. Fluid Mech. 5 (1959) 113-133) for the viscous-convective range of the scalar spectrum; the scalar gradient skewness approaches zero; and the intermittency measured by the scalar gradient flatness approaches its asymptotic state. However, the value of Sc needed for the asymptotic behavior to emerge appears to increase with decreasing Reynolds number of the turbulence. In the viscous-diffusive range, the scalar spectrum is in better agreement with Kraichnan's (Kraichnan., Phys. Fluids 11 (1968) 945-953) result than with Batchelor's.     

How to estimate the low wavenumber turbulence frequency spectrum

 An estimate of the low wavenumber component of surface turbulence shear stress as a function of frequency has been obtained through measurements of the correlations of the longitudinal component of turbulence velocity made close to the surface at y ⁺=7. The data were acquired in a fully-developed turbulent pipe flow at a Reynolds number (based on centreline velocity and pipe diameter) of 268000, using two single hot-wire anemometer probes. A novel data analysis procedure has been introduced to establish the accuracy limits of the low wavenumber turbulence energy estimate for frequencies in the similarity regime of wall turbulence and the results are compared with other measurement techniques. Received: 18 November 1993/Accepted: 21 April 1997     

Large eddy simulation using a dynamic mixing length subgrid‐scale model

A novel dynamic mixing length (DML) subgrid‐scale model for large eddy simulations is proposed in this work to improve the cutoff length of the Smagorinsky model. The characteristic mixing length (or the characteristic wave number) is dynamically estimated for the subgrid‐scale fluctuation of turbulence by the cutoff wave‐number, k_c, and the dissipation wave‐number, k_d. The dissipation wave number is derived from the kinetic energy spectrum equation and the dissipation spectrum equation. To prove the promise of the DML model, this model is used to simulate the lid‐driven cubical cavity with max‐velocity‐based Reynolds numbers 8850 and 12,000, the channel flows with friction‐velocity‐based Reynolds numbers 180, 395, 590, and 950, and the turbulent flow past a square cylinder at the higher Reynolds number 21,400, respectively, compared with the Smagorinsky model and Germano et al.'s dynamic Smagorinsky model. Different numerical experiments with different Reynolds numbers show that the DML model can be used in simulations of flows with a wide range of Reynolds numbers without the occurrence of singular values. The DML model can alleviate the dissipation of the Smagorinsky model without the loss of its robustness. The DML model shows some advantages over Germano et al.'s dynamic Smagorinsky model in its high stability and simplicity of calculation because the coefficient of the DML model always stays positive. The characteristic mixing length in the DML model reflects the subgrid‐scale fluctuation of turbulence in nature and thus the characteristic mixing length has a spatial and temporal distribution in turbulent flow. Copyright © 2011 John Wiley & Sons, Ltd.     

Sharp vorticity gradients in two-dimensional turbulence and the energy spectrum

Formation of sharp vorticity gradients in two-dimensional (2D) hydrodynamic turbulence and their influence on the turbulent spectra are considered. The analog of the vortex line representation as a transformation to the curvilinear system of coordinates moving together with the di-vorticity lines is developed and compressibility of this mapping appears as the main reason for the formation of the sharp vorticity gradients at high Reynolds numbers. In the case of strong anisotropy the sharp vorticity gradients can generate spectra which fall off as k ^?3 at large k, which appear to take the same form as the Kraichnan spectrum for the enstrophy cascade. For turbulence with weak anisotropy the k dependence of the spectrum due to the sharp gradients coincides with the Saffman spectrum: E(k) ~ k ^?4. Numerical investigations of decaying turbulence reveal exponential growth of di-vorticity with a spatial distributed along straight lines. Thus, indicating strong anisotropy and accordingly the spectrum is close to the k ^?3-spectrum.     

The influence of normal stress anisotropy in predicting scalar dispersion with the v–f model

A numerical study of scalar dispersion is presented to investigate the effectiveness of pairing the v²–f turbulence model with algebraic models for the scalar flux. This approach is contrasted with utilizing a full Second Moment Closure (SMC) as the flow field input to the scalar model. Predictions of scalar transport in a turbulent channel and over a wavy wall are compared to available DNS databases. The latter case includes a scalar release from a point source and therefore detailed comparisons of the three-component turbulent scalar flux are reported. It is found that the transported variable v², representing the near wall turbulent velocity fluctuation scale, can be used to increase the level of normal stress anisotropy provided to algebraic scalar models and thereby improve mean scalar prediction over that of the Standard Gradient Diffusion Hypothesis (SGDH). Improvement is most significant in the near wall region. Three specifications of the normal stresses, derived from v², are considered to provide the link from the v²–f model to the algebraic flux models used to close the scalar transport equation. Barycentric maps are used to examine the state of turbulence anisotropy in each case. As the anisotropy in the normal stress specification becomes more accurate, improvements are realized in the prediction of the spanwise flux as well as the mean concentration.     

A modified Lin equation for the energy balance in isotropic turbulence

At sufficiently large Reynolds numbers, turbulence is expected to exhibit scale-invariance in an intermediate (“inertial”) range of wavenumbers, as shown by power law behavior of the energy spectrum and also by a constant rate of energy transfer through wavenumber. However, there is an apparent contradiction between the definition of the energy flux (i.e., the integral of the transfer spectrum) and the observed behavior of the transfer spectrum itself. This is because the transfer spectrum T(k) is invariably found to have a zero-crossing at a single point (at k = k_*), implying that the corresponding energy flux cannot have an extended plateau but must instead have a maximum value at k = k_*. This behavior was formulated as a paradox and resolved by the introduction of filtered/partitioned transfer spectra, which exploited the symmetries of the triadic interactions (J. Phys. A: Math. Theor., 2008). In this paper we consider the more general implications of that procedure for the spectral energy balance equation, also known as the Lin equation. It is argued that the resulting modified Lin equations (and their corresponding Navier–Stokes equations) offer a new starting point for both numerical and theoretical methods, which may lead to a better understanding of the underlying energy transfer processes in turbulence. In particular the filtered partitioned transfer spectra could provide a basis for a hybrid approach to the statistical closure problem, with the different spectra being tackled using different methods.     

Vorticity dynamics in turbulence growth

The mechanisms of vorticity amplification in the formation of turbulence are investigated by means of direct numerical simulations of the Navier–Stokes equations with different initial conditions and Reynolds numbers. The simulations show good universality of the enstrophy evolution, that occurs in two stages. The first stage is dominated by the effect of vortex stretching, and it finishes with a k ^?3 power-law energy spectrum. The second stage is dominated by the action of viscosity on the small scales, and it finishes with a Kolmogorov k ^?5/3 energy spectrum.     

Optimal Gradient Estimates and Asymptotic Behaviour for the Vlasov–Poisson System with Small Initial Data

The Vlasov–Poisson system describes interacting systems of collisionless particles. For solutions with small initial data in three dimensions it is known that the spatial density of particles decays as t ⁻³ at late times. In this paper this statement is refined to show that each derivative of the density which is taken leads to an extra power of decay, so that in N dimensions for N \geqq 3{N \geqq 3} the derivative of the density of order k decays as t ^−N-k . An asymptotic formula for the solution at late times is also obtained.     

Response of mean turbulent energy dissipation rate and spectra to concentrated wall suction

The response of mean turbulent energy dissipation rate and spectra to concentrated suction applied through a porous wall strip has been quantified. Both suction and no suction data of the spectra collapsed reasonably well for Kolmogorov normalised wavenumber k ₁^* > 0.2. Similar results were also observed for second-order structure functions (not shown) for Kolmogorov normalised radius r* < 10. Although, the quality of collapsed is poorer for transverse component, the result highlights that Kolmogorov similarity hypothesis is reasonably well satisfied. However, the suction results shows a significant departure from the no suction case of the Kolmogorov normalised spectra and second-order structure functions for k ₁^* < 0.2 and r* > 20, respectively. The departure at the larger scales with collapse at the small scales suggests that suction induce a change in the small-scale motion. This is also reflected in the alteration of mean turbulent energy dissipation rate and Taylor microscale Reynolds number. This change is a result of the weakening of the large-scale structures. The effect is increased as the suction rate is increased.     

MHD Turbulence at the Laboratory Scale: Established Ideas and New Challenges

The properties of MHD turbulence in the electrically conducting fluids available in the laboratory (where the magnetic Reynolds number is significantly smaller than unity) may be summarised as follows:(1) The Alfven waves, even under their degenerated form at this scale, are responsible for a tendency to two-dimensionality. Eddies tend to become aligned with the applied magnetic field and inertia tends to restore isotropy. The competition between these mechanisms results in a spectral law t^-2k^-3.(2) When insulating walls, perpendicular to the magnetic field, are present and close enough to each other, two-dimensionality can be established with a good approximation within the large scales, and the predominant mechanism is the inverse energy cascade.(3) These columnar eddies are nevertheless submitted to a dissipation within the Hartmann boundary layers present at their ends, whose time scale is independent of the wave number. When this damping effect is negligible, ordinary 2D turbulence is observed with k^-5/3 spectra. On the contrary when this (ohmic and viscous) damping is significant this 2D turbulence exhibits k^-3 spectra.Besides these homogeneous (except within the Hartmann layers) conditions, for instance in shear flows such as mixing layers, almost nothing is known except that two-dimensionality may be well established. The first results of a recent experimental investigation (still in development) are presented. Some challenging questions are raised, such as the interpretation of a surprising difference between the transport of momentum and the transport of a scalar quantity (heat) across that layer. A video was shown during the oral presentation of this paper, illustrating the energy transfer toward the large scales and the weakness of the dissipation suffered by this 2D velocity field.     

Numerical solution of the turbulent flow of an incompressible fluid in curved planar and axially symmetrical channels

The governing equations for axially symmetric flow, where the Reynolds stresses are expressed by scalar turbulent viscosity, are the Reynolds equations. The turbulence model k, ? is used in the well-known form for fully developed turbulent flow.The numerical method, a continuation of the MAC system¹, is adapted so that even for high Reynolds cell numbers precision (δx²) can be achieved for the steady flow. Irregular cells join the rectangular network on the curved surface. Von Neumann's stability condition of the linearised numerical system is investigated. Special problems concerning the numerical solution of the turbulence model equations are stated and a special procedure is worked out to ensure that the fields k, ? do not converge to physically meaningless values. The program for the computer is universal in that the boundary problems can be assigned by input data.As an example, an axially symmetrical diffuser with an area ratio of widening 1.40 is computed. Fields of velocity and pressure at the wall as well as fields v_T and k are assessed. The results are compared with an experiment. The conclusion is that this method is suitable for the problems mentioned in this study as well as for unsteady flow.     

Numerical prediction of turbulent flow in a conical diffuser usingk-ε model

Fully developed incompressible turbulent flow in a conical diffuser having a total divergence angle of 8° and an area ratio of 4∶1 has been simulated by ak-ε turbulence model with high Reynolds number and adverse pressure gradient. The research has been done for pipe entry Reynolds numbers of 1.16×10⁵ and 2.93×10⁵. The mean flow velocity and turbulence energy are predicted successfully and the advantage of Boundary Fit Coordinates approach is discussed. Furthermore, thek-ε turbulence model is applied to a flow in a conical diffuser having a total divergence angle of 30° with a perforated screen. A simplified mathematical model, where only the pressure drop is considered, has been used for describing the effect of the perforated screen. The optimum combination of the resistance coefficient and the location of the perforated screen is predicted for high diffuser efficiency or the uniform velocity distribution.     

Spray Penetration in a Turbulent Flow

Analytical expressions for mass concentration of liquid fuel in a spray are derived taking into account the effects of gas turbulence, and assuming that the influence of droplets on gas is small (intitial stage of spray development). Beyond a certain distance the spray is expected to be fully dispersed. This distance is identified with the maximum spray penetration. Then the influence of turbulence on the spray stopping distance is discussed and the rms spray penetration is computed from a trajectory (Lagrangian) approach. Finally, the problem of spray penetration is investigated in a homogeneous two-phase flow regime taking into account the dispersion of spray away from its axis. It is predicted that for realistic values of spray parameters the spray penetration at large distances from the nozzle is expected to be proportional to t ^2/3 (in the case when this dispersion is not taken into account this distance is proportional to t ^1/2). The t ^2/3 law is supported by experimental observations for a high pressure injector. This revised version was published online in July 2006 with corrections to the Cover Date.     

Measurement of the local wavenumber and frequency spectrum in a plane wake

A digital technique is presented for experimentally measuring the local wavenumber-frequency spectrum S _L(k, f) of a fluctuating velocity field using two probes. From S _L(k, f), the local wavenumber spectrum S _L(k), the averaged dispersion relation, and the broadening of the wavenumber spectrum for a given frequency can be determined. The technique is demonstrated by applying it to the velocity field of a plane wake which is undergoing transition from laminar to turbulent flow. A specially designed two sensor hotwire probe is used to obtain simultaneous records of streamwise velocity fluctuations at two locations having a fixed streamwise separation. New information is obtained concerning: the spatial characteristics of instability waves in the wake; the importance of local wavenumber matching in nonlinear coupling among waves; and broadening of the dispersion relation associated with the transition to turbulence. The usefulness of the S _L(k, f) approach when Taylor's hypothesis is not valid is also discussed.     

Outer layer similarity in fully rough turbulent boundary layers

Turbulent boundary layer measurements were made on a flat plate covered with uniform spheres and also on the same surface with the addition of a finer-scale grit roughness. The measurements were carried out in a closed return water tunnel, over a momentum thickness Reynolds number (Re) range of 3,000–15,000, using a two-component, laser Doppler velocimeter (LDV). The results show that the mean profiles for all the surfaces collapse well in velocity defect form. Using the maximum peak to trough height (R_t) as the roughness length scale (k), the roughness functions (U⁺) for both surfaces collapse, indicating that roughness texture has no effect on U⁺. The Reynolds stresses for the two rough surfaces also show good agreement throughout the entire boundary layer and collapse with smooth wall results outside of the roughness sublayer. Quadrant analysis and the velocity triple products show changes in the rough wall boundary layers that are confined to y<8k_s, where k_s is the equivalent sand roughness height. The present results provide support for Townsends wall similarity hypothesis for uniform three-dimensional roughness. However, departures from wall similarity may be observed for rough surfaces where 5k_s is large compared to the thickness of the inner layer.     

An application of the k-ɛ model with variable Prandtl number to heat transfer computations in air flows

The two-equation `low Reynolds number' k-? model of turbulence with a set of universal constants suggested by Launder and Sharma is modified in the present paper. The variability of the turbulent Prandtl number Pr_t in the energy equation is assumed along with a change of a constant in the dissipation term of the turbulent kinetic energy equation. The turbulent heat transfer is computed for an air flow in a circular pipe for the Reynolds number within the range of 10⁴?4. The modification considerably improves the agreement between the numerical results and the experiment data published in the available literature.     

Two‐dimensional prediction of time dependent,turbulent flow around a square cylinder confined in a channel

This paper presents two‐dimensional and unsteady RANS computations of time dependent, periodic, turbulent flow around a square block. Two turbulence models are used: the Launder–Sharma low‐Reynolds number k–ε model and a non‐linear extension sensitive to the anisotropy of turbulence. The Reynolds number based on the free stream velocity and obstacle side is Re=2.2×10⁴. The present numerical results have been obtained using a finite volume code that solves the governing equations in a vertical plane, located at the lateral mid‐point of the channel. The pressure field is obtained with the SIMPLE algorithm. A bounded version of the third‐order QUICK scheme is used for the convective terms. Comparisons of the numerical results with the experimental data indicate that a preliminary steady solution of the governing equations using the linear k–ε does not lead to correct flow field predictions in the wake region downstream of the square cylinder. Consequently, the time derivatives of dependent variables are included in the transport equations and are discretized using the second‐order Crank–Nicolson scheme. The unsteady computations using the linear and non‐linear k–ε models significantly improve the velocity field predictions. However, the linear k–ε shows a number of predictive deficiencies, even in unsteady flow computations, especially in the prediction of the turbulence field. The introduction of a non‐linear k–ε model brings the two‐dimensional unsteady predictions of the time‐averaged velocity and turbulence fields and also the predicted values of the global parameters such as the Strouhal number and the drag coefficient to close agreement with the data. Copyright © 2009 John Wiley & Sons, Ltd.     

The Effect of Reynolds Number on the Passive Scalar Field in the Turbulent Wake of a Circular Cylinder

In this paper we investigate by experiments the effect of Reynolds number on a passive scalar (temperature) field in the turbulent wake of a slightly heated circular cylinder. The Reynolds number defined by Re≡U _∞ d/ν (see Nomenclature) is varied from Re= 1200 to Re= 8600. Temperature differential above ambient is chosen to be the passive scalar quantity. Present measurements are conducted using a cold wire (0.63 μm) probe. Results obtained suggest that Reynolds number in general has significant influence on the scalar mixing characteristics in the entire wake flow. Specifically, as Re increases, the mean scalar spreads out more rapidly, the scalar fluctuation intensity increases; however, its variance decays at a lower rate with downstream distance. It is also found that an increase of Re accelerates the streamwise evolution of the scalar probability density function from highly non-Gaussianity to near Gaussianity along the wake centreline. This reflects the reduction in length of the Karman-vortex street caused by an increase of Re.     

Reynolds Number Dependence of Energy Spectra in the Overlap Region of Isotropic Turbulence

A Near-Asymptotics analysis of the turbulence energy spectrum is presented that accounts for the effects of finite Reynolds number recently reported by Mydlarski and Warhaft [21]. From dimensional and physical considerations (following Kolmogorov and von Karman), proper scalings are defined for both low and high wavenumbers, but with functions describing the entire range of the spectrum. The scaling for low wavenumbers uses the kinetic energy and the integral scale, L, based on the integral of the correlation function. The fact that the two scaled profiles describe the entire spectrum for finite values of Reynolds number, but reduce to different profiles in the limit, is used to determine their functional forms in the “overlap” region that both retain in the limit. The spectra in the overlap follow a power law, E(k) =Ck _{−5/3 + μ}, where μ and C are Reynolds number dependent. In the limit of infinite Reynolds number, μ → 0 and C → constant, so the Kolmogorov/Obukhov theory is recovered in the limit. Explicit expressions for μ and the other parameters are obtained, and these are compared to the Mydlarski/Warhaft data. To get a better estimate of the exponent from the experimental data, existing models for low and high wavenumbers are modified to account for the Reynolds number dependence. They are then used to build a spectral model covering all the range of wavenumbers at every Reynolds number. Experimental data from grid-generated turbulence are examined and found to be in good agreement with the theory and the model. Finally, from the theory and data, an explicit form for the Reynolds number dependence of φ = ɛL/u ₃ is obtained. This revised version was published online in July 2006 with corrections to the Cover Date.     

Velocity measurements in a plane turbulent air jet at moderate Reynolds numbers

An experimental investigation of the moderate Reynolds number plane air jets was undertaken and the effect of the jet Reynolds number on the turbulent flow structure was determined. The Reynolds number, which was defined by the jet exit conditions, was varied between 1000 and 7000. Other initial conditions, such as the initial turbulence intensity, were kept constant throughout the experiments. Both hot-wire and laser Doppler anemometry were used for the velocity measurements. In the moderate Reynolds number regime, the turbulent flow structure is in transition. The average size and the number of the large scale of turbulence (per unit length of jet) was unaffected by the Reynolds number. A broadening of the turbulent spectra with increasing Reynolds number was observed. This indicated that there is a decrease in the strength of the large eddies resulting from a reduction of the relative energy available to them. This diminished the jet mixing with the ambient as the Reynolds number increased. Higher Reynolds numbers led to lower jet dilution and spread rates. On the other hand, at higher Reynolds numbers the dependence of jet mixing on Reynolds number became less significant as the turbulent flow structure developed into a self-preserving state.List of symbols b _u velocity half-width of the jet - C _u, C _u,0 constants defining the velocity decay rate - D nozzle width - E _u one dimensional power spectrum of velocity fluctuations - f frequency - K _u, K _u,0 constants defining the jet spread rate - k wavenumber (2f/U) - L longitudinal integral scale - R ₁₁ correlation function - r separation distance - Re jet Reynolds number (U ₀ D/v) - St Strouhal number (fD/U ₀) - t time - U axial component of the mean velocity - U _m mean velocity on the jet axis - U ₀ mean velocity at the jet exit - u the rms of u - u fluctuating component of the axial velocity - V lateral component of the mean velocity - fluctuating component of the lateral velocity - x axial distance from the nozzle exit - y lateral distance from the jet axis - z spanwise distance from the jet axis - v kinematic viscosity - time lag A version of this paper was presented as paper no. 86-0038 at the AIAA 24th Aerospace Sciences Meeting, Reno NV, USA, January 1986