DNS of the Turbulent Channel Flow of a Dilute Polymer Solution

A direct numerical simulation of the turbulent channel flow of a dilute polymer solution has been performed in order to compare its turbulence statistics with those obtained in a Newtonian channel flow. The viscoelastic flow has been simulated by solving the whole set of continuity, momentum and constitutive equations for the six independent components of the extra-stress tensor induced by polymer addition. The Finitely Extensible Nonlinear Elastic dumbbell model was adopted in order to simulate a non-linear modulus of elasticity and a finite extendibility of the polymer macromolecules. Simulations were carried out under the narrow channel assumption at a Reynolds number of 169 based on the channel half height and on the friction velocity; they showed a significant reduction in drag, dependent on the influence of the elastic properties of the chains. A qualitative comparison with experiments at a higher Reynolds number has shown that the model here adopted is capable of reproducing all the main features of the polymer solution flow. Analysis of the turbulence statistics suggests that a dilute polymer solution can affect the intensity of the streamwise vortices, leading to an increase in the spacing between low speed streaks and eventually to a turbulent shear stress reduction.     

DNS and Modelling of Passive Scalar Transport in Turbulent Channel Flow with a Focus on Scalar Dissipation Rate Modelling

A direct numerical simulation of turbulent channel flow with an imposed mean scalar gradient is analyzed with a focus on passive scalar flux modelling and in particular the treatment of the passive scalar dissipation equation. The Prandtl number is 0.71 and the Reynolds number based on the wall friction velocity and the channel half width is 265. Budgets are presented for the passive scalar variance and its dissipation rate, as well as for the individual scalar flux components. These form a basis for a discussion of modelling issues related to explicit algebraic scalar flux modelling. This revised version was published online in July 2006 with corrections to the Cover Date.     

Simulation and Modelling of a Skewed Turbulent Channel Flow

A time-dependent three-dimensionally skewed flow is investigated using direct numerical simulations of the incompressible Navier-Stokes equations. The effect on the instantaneous and mean turbulent field is investigated. Instantaneous flowfields reveal that the skewing has the effect of initially reducing the strength and height of quasi-streamwise vortices of both signs of rotation with respect to the skewing. A mechanism for this process is put forward. The mean flowfields show drops in turbulence quantities such as turbulence kinetic energy. In addition to this, two-equation turbulence modelling of the flow is carried out. This highlights a deficiency, in that the standard turbulence models are unable to capture the drop in turbulence intensity due to the skewing. A modification based on the exact dissipation equation is found to significantly improve the model behaviour for this flow. This revised version was published online in July 2006 with corrections to the Cover Date.     

The Use of DNS to Define Stress Producing Events for Turbulent Flow over a Smooth Wall

Over the past 15 years direct numerical simulations (DNS) of turbulent flow and particle image velocimetry (PIV) have provided the opportunity to obtain information about a turbulent velocity field simultaneously at a large number of locations. This paper gives a personal viewpoint of how these techniques are providing new insights about the Reynolds stress producing structures in turbulence generated by flow over a smooth boundary. This revised version was published online in August 2006 with corrections to the Cover Date.     

Numerical Simulation of a Turbulent Flow in a Channel with Surface Mounted Cubes

In this paper we report on a fourth-order, spectro-consistent simulation of a complex turbulent flow. A spatial discretization of a convection-diffusion equation is termed spectro-consistent if the spectral properties of the convective and diffusive operators are preserved, i.e. convection skew-symmetric; diffusion symmetric positive definite. We consider a fully developed flow in a channel, where a matrix of cubes is placed at a wall of the channel. The Reynolds number (based on the channel width and the mean bulk velocity) is equal to Re = 13,000. The three-dimensional flow around the surface mounted cubes has served at a test case at the 6th ERCOFTAC/IAHR/COST workshop on refined flow modeling (Delft, June 1997). Here, mean velocity profiles as well as Reynolds stresses at various locations in the channel have been computed without using any turbulence models. The results agree well with the available experimental data.     

The Relaminarization Mechanisms of Turbulent Channel Flow at Low Reynolds Numbers

The mechanisms of laminarization in wall-bounded flows have been investigated by performing direct numerical simulations (DNS) of turbulent channel flows. By decreasing Reynolds numbers systematically, the effects of the low Reynolds number are studied in connection with the near-wall turbulent structure and turbulent statistics. At approximately the critical Reynolds number, the turbulent skin friction is reduced, and the turbulent structure changes qualitatively in the very near-wall region. Instantaneous turbulent structures reveal that streamwise vortices, the cores of which are at y+ 10, disappear, although low speed streaks and Reynolds shear stress are still produced by larger streamwise vortices located in the buffer region y+ > 10. Sweep motions induced by these vortical structures are shifted toward the center of a channel and also significantly deterred, which may heighten the effects of the viscous sublayer over most of the channel section and suppress the regeneration mechanisms of new streamwise vortices in the very near-wall region. To investigate the details of how large-scale coherent vortices affect the viscous sublayer and the relevant small-scale streamwise vortices, a body force is virtually imposed in the wall-normal direction to enhance the large streamwise vortices. As a result, it is found that when they are sufficiently enhanced, the small-scale vortices reappear, and the sweep events are again dominant in the viscous sublayer.     

各向同性湍流内颗粒碰撞率的直接模拟研究

对 Re_{\lambda } 约为51均匀各向同性湍流内 St_{k}(=\tau_{p}/\tau_{k}) 为 0 ~10.0 的 有限惯性颗粒的碰撞行为进行了直接数值模拟，以研究湍流对有限惯性 颗粒碰撞的影响. 结果表明，具有一定惯性颗粒的湍流碰撞率完全不同于零惯性的轻颗粒 (St_{k}=0) 和可忽略湍流作用的重颗粒 (St{k} \to \infty) , 其变化趋势极其复杂： 在Stk为 0~1.0 之间，颗粒的碰撞率随 St 的增加而近乎线性地剧烈增长，在 Stk≈1.0 3.0(对应的StE=τp/Te≈0.5)附近，颗粒碰撞率出现两个峰值，在Stk>3.0以后，颗粒的碰撞率随惯性增 大而逐渐趋向于重颗粒极限；在峰值处，有限惯性颗粒的平均碰撞率的峰值较轻颗粒增强了 30倍左右. 为进一步分析湍流作用下颗粒碰撞率的影响因素，分别使用可能发生碰撞 的颗粒对的径向分布函数和径向相对速度来量化颗粒的局部富集效应和湍流掺混效应，表明 St_{k} \approx 1.0 时局部富集效应最为强烈，使得颗粒的碰撞率出现第1个峰值； 湍流掺混效应则随着颗粒Stk的增大而渐近增大；局部富集和湍流掺混联合作用的结果， 使得颗粒碰撞率在 St_{k} \approx 3.0 附近出现另一个峰值.     

Wall Versus Centerline Polymer Injection in Turbulent Channel Flows

This experimental study compares the mean and turbulence characteristics of turbulent channel flows with polymer injection at the wall and at the centerline to assess the impact of the injection location on drag reduction. It also contrasts the drag reduction performance of a hydrolyzed polymer versus a non-ionic polymer under the same conditions. Wall injection of non-ionic and hydrolized polymers resulted in 23% and 9% larger drag reduction than corresponding centerline injection, respectively. In all cases, the polymer was structured and the presence of macromolecular polymer structures, even when concentrated mostly away from the wall, seemed to be able to affect the turbulence structure in the flow.     

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.     

A Conditional Moment Closure Study of Chemical Reaction Source Terms in SCCI Combustion

The objective of this study is to evaluate conditional moment closure (CMC) approaches to model chemical reaction rates in compositionally stratified, autoigniting mixtures, in thermochemical conditions relevant to stratified charge compression ignition (SCCI) engines. First-order closure, second-order closure and double conditioning are evaluated and contrasted as options in comparison to a series of direct numerical simulations (DNSs). The two-dimensional (2D) DNS cases simulate ignitions in SCCI-like thermochemical conditions with compositionally stratified n-heptane/air mixtures in a constant volume. The cases feature two different levels of stratification with three mean temperatures in the negative-temperature coefficient (NTC) regime of ignition delay times. The first-order closure approach for reaction rates is first assessed using hybrid DNS-CMC a posteriori tests when implemented in an open source computational fluid dynamics (CFD) package known as OpenFOAM\(^{{\circledR }}\). The hybrid DNS-CMC a posteriori tests are not a full CMC but a DNS-CMC hybrid in that they compute the scalar and velocity fields at the DNS resolution, thus isolating the first-order reaction rate closure model as the main source of modelling error (as opposed to turbulence model, scalar probability density function model, and scalar dissipation rate model). The hybrid DNS-CMC a posteriori test reveals an excellent agreement between the model and DNS for the cases with low levels of stratification, whereas deviations from the DNS are observed in cases which exhibit high level of stratifications. The a priori analysis reveals that the reason for disagreement is failure of the first-order closure hypothesis in the model due to the high level of conditional fluctuations. Second-order and double conditioning approaches are then evaluated in a priori tests to determine the most promising path forwards in addressing higher levels of stratification. The a priori tests use the DNS data to compute the model terms, thus directly evaluating the model assumptions. It is shown that in the cases with a high level of stratification, even the second-order estimation of the reaction rate source term cannot provide a reasonably accurate closure. Double conditioning using mixture-fraction and sensible enthalpy, however, provides an accurate first-order closure to the reaction rate source term.     

A DNS Study of Closure Relations for Convection Flux Term in Transport Equation for Mean Reaction Rate in Turbulent Flow

The present work aims at modeling the entire convection flux \(\overline {\rho \mathbf {u}W}\) in the transport equation for a mean reaction rate \(\overline {\rho W}\) in a turbulent flow, which (equation) was recently put forward by the present authors. In order to model the flux, several simple closure relations are developed by introducing flow velocity conditioned to reaction zone and interpolating this velocity between two limit expressions suggested for the leading and trailing edges of the mean flame brush. Subsequently, the proposed simple closure relations for \(\overline {\rho \mathbf {u}W}\) are assessed by processing two sets of data obtained in earlier 3D Direct Numerical Simulation (DNS) studies of adiabatic, statistically planar, turbulent, premixed, single-step-chemistry flames characterized by unity Lewis number. One dataset consists of three cases characterized by different density ratios and is associated with the flamelet regime of premixed turbulent combustion. Another dataset consists of four cases characterized by different low Damköhler and large Karlovitz numbers. Accordingly, this dataset is associated with the thin reaction zone regime of premixed turbulent combustion. Under conditions of the former DNS, difference in the entire, \(\overline {\rho {u}W}\), and mean, \(\tilde {u}\overline {\rho W}\), convection fluxes is well pronounced, with the turbulent flux, \(\overline {\rho u^{\prime \prime }W^{\prime \prime }}\), showing countergradient behavior in a large part of the mean flame brush. Accordingly, the gradient diffusion closure of the turbulent flux is not valid under such conditions, but some proposed simple closure relations allow us to predict the entire flux \(\overline {\rho \mathbf {u}W}\) reasonably well. Under conditions of the latter DNS, the difference in the entire and mean convection fluxes is less pronounced, with the aforementioned simple closure relations still resulting in sufficiently good agreement with the DNS data.     

蜗壳内紊流流场的有限元法计算

用有限元法计算径流式叶轮机械蜗壳的紊流时均流场.有关紊流模型采用K-ε两方程模型,用关于压力p的罚函数方法求解.所得结果可供分析蜗壳流场用.文中方法也可用于计算其他形状的二维通道流动.     

Large Eddy Simulation of a Turbulent Lifted Flame using Conditional Moment Closure and Rate-Controlled Constrained Equilibrium

The objective of this paper is the simulation of a turbulent flame by employing the Rate-Controlled Constrained Equilibrium (RCCE) approach for the chemistry reduction, and Large-Eddy Simulation (LES) coupled with Conditional Moment Closure (CMC) for the turbulence-chemistry interaction modelling. RCCE is a systematic method for mechanism reduction, based on the concept that certain species characterized by faster time scales are in a constrained equilibrium state, determined by the concentration of the species controlled by the chemical kinetics. A general system of differential equations can be derived, independent on the selection of the fast and slow species (which appears as a parameter). The RCCE system is used to compute the conditional source term in the CMC equation. The flame simulated here is a methane flame issuing into a vitiated co-flow formed by hot combustion products, the ??Cabra?? flame, which is controlled by auto-ignition and is therefore sensitive to the chemical mechanism. The results show an influence of the chosen chemistry in the ignition length.     

Exchange of Chemical Species between the Atmospheric Boundary Layer and the Reservoir Layer: An Analytical Interpretation

The exchange of chemical species between the atmospheric boundary layer and the reservoir layer is investigated by means of an analytical solution of the conservation equation of a decaying chemical species. The exchange mechanism is governed by two parameters: the Damköhler number (the ratio of the turbulence time scale to the chemical time scale) and the ratio of a concentration scale in the atmospheric boundary layer to the concentration in the reservoir layer. Depending on the value of these two parameters, the exchange flux between the two layers can vary in sign and by several orders of magnitude. The study demonstrates to what extent chemical transformation determines the transfer of chemical species between the atmospheric boundary layer and the reservoir layer.     

湍流边界层流场与噪声实验研究

在重力式水洞中进行了水翼及半翼湍流边界层流场与噪声的实验研究。测量了水翼及半翼边界层附近的湍流脉动速度场;测量了半翼翼型表面三点处的压力脉动及其辐射噪声,测量了水翼内部测点的噪声及外部辐射噪声,在不同流速、不同攻角、光滑和粗糙翼面的情况下都进行了测量分析。试验结果发现,上述因素对模型的边界层湍流速度场有显著影响,15度攻角时,翼面附近湍流强度要比0度时大得多,粗糙翼面附近的湍流强度比光滑的大,而湍流强度随来流速度的变化不大,u(来流）方向和v方向的湍流强度量级相当;一般地,翼面压力脉动、翼内部噪声及外部噪声都是随来流速度的增大而增大,随攻角的增大而增大,粗糙翼面时的结果要比光滑翼面的大。从压力脉动与噪声测量结果与相应的流场测量结果比较可知,可以从湍流区域的湍流强度来判断出声源强度的定性变化。     

基于小型电磁动力驱动器的微流道混合控制研究

对于许多微流体应用系统来说,流体的混合是至关重要的.本文研究了基于电磁流体动力混合器的微流道主动混合控制方法,建立了该微流道混合系统的理论模型并进行了数值模拟.在交变Lorentz驱动下,流经混合器流道的流体及其分界面在混合器流道内往复运动,流体界面的反复折叠与流体局部流动使流体混合,流体的横向往复运动使流体分界面反复折叠从而使流体间的接触面积大大增加.对流体混合过程进行了讨论并给出了漉体混合程度评价方法.     

Conditional Source-Term Estimation as a Method for Chemical Closure in Premixed Turbulent Reacting Flow

A Conditional Source-term Estimation (CSE) model is used to close the mean reaction rates for a turbulent premixed flame. A product-based reaction progress variable is introduced as the conditioning variable for the CSE method. Different presumed probability density function (PDF) models are studied and a modified version of a laminar flame-based PDF model is proposed. Improved predictions of the variable distribution are obtained. The conditional means of reactive scalars are evaluated with CSE and compared to the direct numerical simulation (DNS). The mean reaction rates in a turbulent premixed flame are evaluated with the CSE model and the presumed PDFs. Comparison of the CSE closure method to DNS shows promising results. This paper was presented at the 2nd ECCOMAS Thematic Conference on Computational Combustion.     

AGE method simulations of a turbulent far‐wake compared to spectral DNS

Turbulent flow simulation methods based on finite differences are attractive for their simplicity, flexibility and efficiency, but not always for accuracy or stability. This paper demonstrates that a good compromise is possible with the advected grid explicit (AGE) method. Starting from the same initial field as a previous spectral DNS, AGE method simulations of a planar turbulent wake were carried out as DNS, and then at three levels of reduced resolution. The latter cases were in a sense large‐eddy simulations (LES), although no specific sub‐grid‐scale model was used. Results for the two DNS methods, including variances and power spectra, were very similar, but the AGE simulation required much less computational effort. Small‐scale information was lost in the reduced resolution runs, but large‐scale mean and instantaneous properties were reproduced quite well, with further large reductions in computational effort. Quality of results becomes more sensitive to the value chosen for one of the AGE method parameters as resolution is reduced, from which it is inferred that the numerical stability procedure controlled by the parameter is acting in part as a sub‐grid‐scale model. Copyright © 2002 John Wiley & Sons, Ltd.     

Second-Order Conditional Moment Closure Simulations of Autoignition of an n-heptane Plume in a Turbulent Coflow of Heated Air

Autoignition of an n-heptane plume in a turbulent coflow of heated air has been studied using the conditional moment closure (CMC) method with a second-order closure for the conditional chemical source term. Two different methodologies have been considered: (i) the Taylor expansion method, in which the second order correction was based on the solution of the full covariance matrix for the 31 reactive species in the chemical mechanism and hence was not limited to a few selected reactions, and (ii) the conditional PDF method, in which only the temperature conditional variance equation has been solved and its PDF assumed to be a β-function. The results compare favorably with experiment in terms of autoignition location. The structure of the reaction zone in mixture fraction space has been explored. The relative performance of the two methodologies is discussed.     

Use of a Modified Langevin Equation to Describe Turbulent Dispersion of Fluid Particles in a Channel Flow

A stochastic method to represent the positions and velocities of fluid particles in a nonhomogeneous turbulence was pursued. Spatially varying Lagrangian time scales obtained from direct numerical simulations of turbulent flow in a channel and spatially varying joint Gaussian forcing functions were incorporated into a Langevin equation. The model was tested by comparing calculations of the dispersions and velocities of particles originating from point sources with experiments carried out in a DNS of fully-developed turbulent flows in a channel at Re_τ = 150 and 300. The model captured the dispersions, mean velocities and moments of the velocity fluctuations very well. The use of jointly Gaussian (rather than uncorrelated Gaussian) forcing functions greatly improved the calculation of dispersions and mean velocities in the stream wise direction, as well as turbulence quantities which include stream wise velocity fluctuations. The condition of well-mixedness for the model was verified by considering the behavior of a uniform distribution of sources. This revised version was published online in July 2006 with corrections to the Cover Date.