Influence of polymer additives on turbulent energy cascading in forced homogeneous isotropic turbulence studied by direct numerical simulations

Direct numerical simulations(DNS) were performed for the forced homogeneous isotropic turbulence(FHIT) with/without polymer additives in order to elaborate the characteristics of the turbulent energy cascading influenced by drag-reducing effects.The finite elastic non-linear extensibility-Peterlin model(FENE-P) was used as the conformation tensor equation for the viscoelastic polymer solution.Detailed analyses of DNS data were carried out in this paper for the turbulence scaling law and the topological dynamics of FHIT as well as the important turbulent parameters,including turbulent kinetic energy spectra,enstrophy and strain,velocity structure function,small-scale intermittency,etc.A natural and straightforward definition for the drag reduction rate was also proposed for the drag-reducing FHIT based on the decrease degree of the turbulent kinetic energy.It was found that the turbulent energy cascading in the FHIT was greatly modified by the drag-reducing polymer additives.The enstrophy and the strain fields in the FHIT of the polymer solution were remarkably weakened as compared with their Newtonian counterparts.The small-scale vortices and the small-scale intermittency were all inhibited by the viscoelastic effects in the FHIT of the polymer solution.However,the scaling law in a fashion of extended self-similarity for the FHIT of the polymer solution,within the presently simulated range of Weissenberg numbers,had no distinct differences compared with that of the Newtonian fluid case.     

Assessment of subgrid-scale stress statistics in non-premixed turbulent wall-jet flames

We investigate the heat-release effects on the characteristics of the subgrid-scale (SGS) stress tensor and SGS dissipation of kinetic energy and enstrophy. Direct numerical simulation data of a non-premixed reacting turbulent wall-jet flow with and without substantial heat release is employed for the analysis. This study comprises, among others, an analysis of the eigenvalues of the resolved strain rate and SGS stress tensors, to identify the heat-release effects on their topology. An assessment of the alignment between the eigenvectors corresponding to the largest eigenvalues of these two tensors is also given to provide further information for modelling of the SGS stress tensor. To find out the heat-release effects on the dynamics of the turbulent kinetic energy and enstrophy dissipation, probability density functions (PDFs) and mean values are analysed. The mean SGS shear stress and turbulent kinetic energy both slightly increase in the buffer layer and substantially decrease further away from the wall, due to the heat-release effects. Contrary to the kinetic energy, heat release decreases the mean SGS dissipation of enstrophy in the near-wall region. Moreover, differences in the shapes of the PDFs between the isothermal and exothermic cases indicate changes in the intermittency level of both SGS dissipations. Heat release also increases the SGS stress anisotropy in the near-wall region. Although, the structure of the mean resolved strain-rate tensor only marginally differs between the isothermal and exothermic cases in the near-wall region, substantial differences are observed in the jet area, where compressibility effects are important and heat-release effects are found to promote compression states. The differences in the relative alignment between the SGS stress and resolved strain-rate tensors in the isothermal and exothermic cases are discussed in connection with the differences in the SGS dissipation of kinetic energy.     

Curvature Effect on Compressible Turbulent Flow over a Wavy Wall

A fully developed compressible turbulent flow in a channel with a lower wavy wall and a upper plane wall is studied using large eddy simulation. We mainly attempt to deal with the curvature effect on compressible turbulent flow over the wavy wall. Some typical quantities including the mean turbulence statistics, dilatation and baroclinic terms in the enstrophy equation, turbulent kinetic energy budgets and the near-wall turbulent structures are analysed. The results obtained in this study provide physical insight into the understanding of the effects of curvature and compressibility on wall-bounded compressible turbulent flow.     

Convection and correlation of coherent structure in turbulent boundary layer using tomographic particle image velocimetry

The present experimental work focuses on a new model for space–time correlation and the scale-dependencies of convection velocity and sweep velocity in turbulent boundary layer over a flat wall. A turbulent boundary layer flow at Reθ= 2460 is measured by tomographic particle image velocimetry(tomographic PIV). It is demonstrated that arch, cane,and hairpin vortices are dominant in the logarithmic layer. Hairpins and hairpin packets are responsible for the elongated low-momentum zones observed in the instantaneous flow field. The conditionally-averaged coherent structures systemically illustrate the key roles of hairpin vortice in the turbulence dynamic events, such as ejection and sweep events and energy transport. The space–time correlations of instantaneous streamwise fluctuation velocity are calculated and confirm the new elliptic model for the space–time correlation instead of Taylor hypothesis. The convection velocities derived from the space–time correlation and conditionally-averaged method both suggest the scaling with the local mean velocity in the logarithmic layer. Convection velocity result based on Fourier decomposition(FD) shows stronger scale- dependency in the spanwise direction than in streamwise direction. Compared with FD, the proper orthogonal decomposition(POD) has a distinct distribution of convection velocity for the large- and small-scales which are separated in light of their contributions of turbulent kinetic energy.     

Interaction between strain and vorticity in compressible turbulent boundary layer

The interaction of strain and vorticity in compressible turbulent boundary layers at Mach number 2.0 and 4.9 is studied by direct numerical simulation(DNS)of the compressible Navier-Stokes equations.Some fundamental characteristics have been studied for both the enstrophy producing and destroying regions.It is found that large enstrophy production is associated with high dissipation and high enstrophy,while large enstrophy destruction with moderate ones.The enstrophy production and destruction are also correlated with the dissipation production and destruction.Moreover,the enstrophy producing region has a distinct tendency to be‘sheet-like’structures and the enstrophy destroying region tends to be‘tube-like’in the inner layer.Correspondingly,the tendency to be‘sheet-like’or‘tube-like’structures is no longer obvious in the outer layer.Further,the alignment between the vorticity vector and the strain-rate eigenvector is analyzed in the flow topologies.It is noticed that the enstrophy production rate depends mainly on the alignment between the vorticity vector and the intermediate eigenvector in the inner layer,and the enstrophy production(destruction)mainly on the alignment between the vorticity vector and the extensive(compressive)eigenvector in the outer layer.     

Improved formulation of electron kinetic theory approach for laser shortpulse heating: Thermal stress consideration

Nonequilibrium energy transport between excited electrons and lattice site is re-formulated after considering the ballistic contribution of the electron energy to the energy transport process. The improved formulation of the electron kinetic theory predictions are compared with the previously obtained electron kinetic and two-equation models. Thermal stress developed in the region irradiated by a laser beam is formulated during the heating pulse. Copper with variable properties is used in the simulations. It is found that improved electron kinetic theory model predicts less temperature rise than that corresponding to previously formulated electron kinetic theory and two equation models in the surface region; in this case, electron temperature attains high values. Thermal stress developed is compressive and attains the maximum at some depth below the surface. The thermal stress level is well below the yielding limit of the substrate material.     

Physical mechanism of the two-dimensional enstrophy cascade

In two-dimensional turbulence, irreversible forward transfer of enstrophy requires anticorrelation of the turbulent vorticity transport vector and the inertial-range vorticity gradient. We investigate the basic mechanism by numerical simulation of the forced Navier-Stokes equation. In particular, we obtain the probability distributions of the local enstrophy flux and of the alignment angle between vorticity gradient and transport vector. These are surprisingly symmetric and cannot be explained by a local eddy-viscosity approximation. The vorticity transport tends to be directed along streamlines of the flow and only weakly aligned down the fluctuating vorticity gradient. All these features are well explained by a local nonlinear model. The physical origin of the cascade lies in steepening of inertial-range vorticity gradients due to compression of vorticity level sets by the large-scale strain field.     

近壁区理论耗散率k-ε模型的研究

提出在湍流边界层近壁区采用共振三波的理论模型描述湍流相干结构,根据理论模型计算了ε的分布。并且在传统k-ε模式基础上依照理论ε值计算了平均速度分布。在粘性作用层理论值与直接数值模拟符合很好。表明该理论方法不仅有益于对湍流机制的了解,而且可能为湍流的近壁模型化开辟一条新的途径。     

A theoretical study of premixed turbulent flame development

Flame development in a statistically stationary and uniform, planar, one-dimensional turbulent flow is theoretically studied. A generalized balance equation for the mean combustion progress variable, which includes turbulent diffusion and pressure-driven transport terms, as well as the mean rate of product creation, is introduced and analyzed by invoking the sole assumption of a self-similar flame structure, well-supported by numerous experiments. The assumption offers the opportunity to simplify the problem by splitting the aforementioned partial differential equation into two ordinary differential equations, which separately model spatial variations of the progress variable and time variations of flame speed and thickness. The self-similar profile of the progress variable, obtained in numerous experiments, is theoretically predicted. Closures of the normalized pressure-driven transport term and mean rate of product creation are obtained. The closed balance equation shows that turbulent diffusion dominates during the initial stage of flame development, followed by the transition to counter-gradient transport in a sufficiently developed flame. A criterion of the transition is derived. The transition is promoted by the heat release and pressure-driven transport. Fully developed mean flame brush thickness and speed are shown to decrease when either density ratio or pressure-driven transport increases. Solutions for the development of the thickness are obtained. The development is accelerated by the pressure-driven transport and heat release.     

湍流两相流动有燃烧颗粒相概率密度函数输运方程理论

由有燃烧的湍流气粒两相流动的瞬态方程和统计力学概率密度函数概念出发,推导了有燃烧颗粒相的质量－动量－能量联合概率密度函数（ＰＤＦ）输运方程,并对方程中条件期望项用梯度模拟概念进行了模拟封闭。封闭后的ＰＤＦ方程可作为建立颗粒拟流体模型方程和封闭二阶矩模型的基础,也可以通过Ｍｏｎｔｅ－Ｃａｒｌｏ 法求解用以直接计算颗粒雷诺应力和湍流动能,以便和二阶 矩模型的结果相对照,改善二阶矩模型。     

A one-equation turbulence model for recirculating flows

A one-equation turbulence model which relies on the turbulent kinetic energy transport equation has been developed to predict the flow properties of the recirculating flows. The turbulent eddy-viscosity coefficient is computed from a recalibrated Bradshaw’s assumption that the constant a₁ = 0.31 is recalibrated to a function based on a set of direct numerical simulation (DNS) data. The values of dissipation of turbulent kinetic energy consist of the near-wall part and isotropic part, and the isotropic part involves the von Karman length scale as the turbulent length scale. The performance of the new model is evaluated by the results from DNS for fully developed turbulence channel flow with a wide range of Reynolds numbers. However, the computed result of the recirculating flow at the separated bubble of NACA4412 demonstrates that an increase is needed on the turbulent dissipation, and this leads to an advanced tuning on the self-adjusted function. The improved model predicts better results in both the non-equilibrium and equilibrium flows, e.g. channel flows, backward-facing step flow and hump in a channel.     

基于层流动能湍流模型的数值模拟方法

从层流动能出发提出数值模拟原则;综合考虑自然、旁路和分离流转捩的因素构建实用的层流动能湍流模型,结合预处理和基本求解技术发展出适于转捩流动的数值模拟方法和程序.针对预处理技术,以Weiss-Smith矩阵为基础,考虑湍流粘性的影响;针对基本离散格式和边界条件,结合模型方程进行对角占优强化等特殊处理.最后通过平板边界层和典型翼型,特别是低雷诺数翼型的数值模拟,验证数值方法的有效性和鲁棒性.算例表明本文的方法能够为求解更复杂的流动提供参考.     

湍流扩散燃烧的数值研究-PDF方法和火焰面模型的性能比较

采用标量概率密度函数(PDF)方法、稳态和非稳态火焰面模型三种方法对一个值班湍流CH_4/O_2/N_2射流扩散火焰(Sandia Flame D)进行数值计算,以比较不同燃烧模型的性能。PDF方法通过计算反应标量的PDF输运方程来得到标量分布,而火焰面模型只求解单标量混合物分数的PDF方程,组分和温度分布通过火焰面方程的求解或者火焰面数据库的插值得到。计算结果和实验数据对比表明PDF方法计算结果最好但计算量相当大,稳态火焰面模型则反之。综合而言,非稳态火焰面模型的预测结果相对稳态模型有了非常大的改进,而计算量仍然容易接受,非常适合工程应用。     

Phonon and electron transport in aluminum thin film: Influence of film thickness on electron and lattice temperatures

Influence of film thickness on non-equilibrium energy transport in the aluminum thin film is examined. The solutions of Boltzmann equation and the modified two-equation model are presented to predict electron and phonon temperatures in the film for various film thicknesses. It is found that electron and phonon temperatures predicted from the Boltzmann equation differ from the solution of two-equation model in the film for small film thickness. As the film thickness increases, this difference becomes negligibly small. Two-equation model predicts higher electron and phonon temperatures than those obtained from the solutions of the Boltzmann equation in the vicinity of the high temperature edge. This becomes opposite in the region of the low temperature edge.     

Self-diffusion in fluids with weak long-range forces

Diffusion of a test particle in a homogeneous classical fluid with weak long-range forces is studied. The dominant mean-field effect (Vlasov's theory) vanishes for symmetry reasons. Dynamical phenomena follow then from fluctuations of the effective potential energy felt by the propagating particle. The kinetic equation corresponding to this mechanism is derived with the use of the multiple-time-scale method. Its structure resembles very much that of the (linearized) Balescu-Lenard equation of hot plasma theory. It is shown that the kinetic equation holds only if no phase transition occurs in the system. The thermalization of the diffusing particle and the high-temperature and Lorentz gas limits are discussed.     

Turbulent scalar flux transport in head-on quenching of turbulent premixed flames: a direct numerical simulations approach to assess models for Reynolds averaged Navier Stokes simulations

The statistical behaviour and the modelling of turbulent scalar flux transport have been analysed using a direct numerical simulation (DNS) database of head-on quenching of statistically planar turbulent premixed flames by an isothermal wall. A range of different values of Damköhler, Karlovitz numbers and Lewis numbers has been considered for this analysis. The magnitudes of the turbulent transport and mean velocity gradient terms in the turbulent scalar flux transport equation remain small in comparison to the pressure gradient, molecular dissipation and reaction-velocity fluctuation correlation terms in the turbulent scalar flux transport equation when the flame is away from the wall but the magnitudes of all these terms diminish and assume comparable values during flame quenching before vanishing altogether. It has been found that the existing models for the turbulent transport, pressure gradient, molecular dissipation and reaction-velocity fluctuation correlation terms in the turbulent scalar flux transport equation do not adequately address the respective behaviours extracted from DNS data in the near-wall region during flame quenching. Existing models for transport equation-based closures of turbulent scalar flux have been modified in such a manner that these models provide satisfactory prediction both near to and away from the wall.     

Phenomenology of two-dimensional stably stratified turbulence under large-scale forcing

In this paper, we characterise the scaling of energy spectra, and the interscale transfer of energy and enstrophy, for strongly, moderately and weakly stably stratified two-dimensional (2D) turbulence, restricted in a vertical plane, under large-scale random forcing. In the strongly stratified case, a large-scale vertically sheared horizontal flow (VSHF) coexists with small scale turbulence. The VSHF consists of internal gravity waves and the turbulent flow has a kinetic energy (KE) spectrum that follows an approximate k^?3 scaling with zero KE flux and a robust positive enstrophy flux. The spectrum of the turbulent potential energy (PE) also approximately follows a k^?3 power-law and its flux is directed to small scales. For moderate stratification, there is no VSHF and the KE of the turbulent flow exhibits Bolgiano–Obukhov scaling that transitions from a shallow k^?11/5 form at large scales, to a steeper approximate k^?3 scaling at small scales. The entire range of scales shows a strong forward enstrophy flux, and interestingly, large (small) scales show an inverse (forward) KE flux. The PE flux in this regime is directed to small scales, and the PE spectrum is characterised by an approximate k^?1.64 scaling. Finally, for weak stratification, KE is transferred upscale and its spectrum closely follows a k^?2.5 scaling, while PE exhibits a forward transfer and its spectrum shows an approximate k^?1.6 power-law. For all stratification strengths, the total energy always flows from large to small scales and almost all the spectral indicies are well explained by accounting for the scale-dependent nature of the corresponding flux.     

A theoretical model of turbulent fiber suspension and its application to the channel flow

A theoretical model of turbulent fiber suspension is developed by deriving the equations of Reynolds averaged Navier-Stokes,turbulence kinetic energy and turbulence dissipation rate with the additional term of fibers.In order to close the above equations,the equation of probability distribution function for mean fiber orientation is also derived.The theoretical model is applied to the turbulent channel flow and the corresponding equations are solved numerically.The numerical results are verified by comparisons with the experimental ones.The effects of Reynolds number,fiber concentration and fiber aspect-ratio on the velocity profile,turbulent kinetic energy and turbulent dissipation rate are analyzed.Based on the numerical data,the expression for the velocity profile in the turbulent fiber suspension channel flow,which includes the effect of Reynolds number,fiber concentration and aspect-ratio,is proposed.     

用求解概率密度函数输运方程的方法模拟湍流自由射流

用求解概率密度函数输运方程的方法模拟湍流自由射流周向阳，郑楚光，马毓义（华中理工大学煤燃烧国家重点实验室武汉４３００７４）关键词概率密度函数；湍流；自由射流１前言和通常的湍流统计矩模型（如ｋ－ε模型或Ｒｅｙｎｏｌｄｓ应力模型等）相比，用求解速度和标量...