Direct numerical simulation of compressible turbulent flows

This paper reviews the authors＇ recent studies on compressible turbulence by using direct numerical simulation （DNS）,including DNS of isotropic（decaying） turbulence, turbulent mixing-layer,turbulent boundary-layer and shock/boundary-layer interaction.Turbulence statistics, compressibility effects,turbulent kinetic energy budget and coherent structures are studied based on the DNS data.The mechanism of sound source in turbulent flows is also analyzed. It shows that DNS is a powerful tool for the mechanistic study of compressible turbulence.     

Direct numerical simulation of wall turbulent flows with microbubbles

The marker‐density‐function (MDF) method has been developed to conduct direct numerical simulation (DNS) for bubbly flows. The method is applied to turbulent bubbly channel flows to elucidate the interaction between bubbles and wall turbulence. The simulation is designed to clarify the structure of the turbulent boundary layer containing microbubbles and the mechanism of frictional drag reduction. It is deduced from the numerical tests that the interaction between bubbles and wall turbulence depends on the Weber and Froude numbers. The reduction of the frictional resistance on the wall is attained and its mechanism is explained from the modulation of the three‐dimensional structure of the turbulent flow. Copyright © 2001 John Wiley & Sons, Ltd.     

Direct numerical simulation of particle interaction with ejections in turbulent channel flows

Direct numerical simulations (DNS) of incompressible turbulent channel flows coupled with Lagrangian particle tracking are performed to study the characteristics of ejections that surround solid particles. The behavior of particles in dilute turbulent channel flows, without particle collisions and without feedback of particles on the carrier fluid, is studied using high Reynolds number DNS (Re = 12,500). The results show that particles moving away from the wall are surrounded by ejections, confirming previous studies on this issue. A threshold value separating ejections with only upward moving particles is established. When normalized by the square root of the Stokes number and the square of the friction velocity, the threshold profiles follow the same qualitative trends, for all the parameters tested in this study, in the range of the experiments. When compared to suspension thresholds proposed by other studies in the Shields diagram, our simulations predict a much larger value because of the measure used to characterize the fluid and the criterion chosen to decide whether particles are influenced by the surrounding fluid. However, for intermediate particle Reynolds numbers, the threshold proposed here is in fair agreement with the theoretical criterion proposed by Bagnold (1966) [Bagnold, R., 1966. Geological Survey Professional Paper, vol. 422-1]. Nevertheless, further studies will be conducted to understand the normalization of the threshold.     

气固两相流中颗粒弥散的拉格朗日模拟

本文提出了一种对于均匀,稳定及各向同性气固两相紊流场中圆形固体颗粒弥散的拉格朗日模拟计算方法,应用该方法对带有网栅的垂直与水平管道中均匀,稳定的气固两相流模拟计算结果与Snyder及Wells等人所做的相同情况下的试验结果进行了比较,以证明该模拟计算方法的有效性,。     

Reynolds stress and vorticity in turbulent wall flows

Experimental measurements in a boundary layer and a large-eddy simulation of plane channel flow have been used to study the dynamics of vorticity and mass transport in the nearwall region. It was found that Reynolds stress generation occurs in the vicinity of quasi-streamwise vortices, and that smoke particles tend to be ejected from the wall near these vortical structures.     

Direct numerical simulation of particle dispersion in a turbulent jet considering inter-particle collisions

To investigate the behaviour of inter-particle collision and its effects on particle dispersion, direct numerical simulation of a three-dimensional two-phase turbulent jet was conducted. The finite volume method and the fractional-step projection algorithm were used to solve the governing equations of the gas phase fluid and the Lagrangian method was applied to trace the particles. The deterministic hard-sphere model was used to describe the inter-particle collision. In order to allow an analysis of inter-particle collisions independent of the effect of particles on the flow, two-way coupling was neglected. The inter-particle collision occurs frequently in the local regions with higher particle concentration of the flow field. Under the influence of the local accumulation and the turbulent transport effects, the variation of the average inter-particle collision number with the Stokes number takes on a complex non-linear relationship. The particle distribution is more uniform as a result of inter-particle collisions, and the lateral and the spanwise dispersion of the particles considering inter-particle collision also increase. Furthermore, for the case of particles with the Rosin–Rammler distribution (the medial particle size is set d₅₀ = 36.7 μm), the collision number is significantly larger than that of the particles at the Stokes number of 10, and their effects on calculated results are also more significant.     

Differential and algebraic models for the second moments of particle velocity and temperature fluctuations in turbulent flows

Differential and algebraic models are constructed for the dispersed-phase turbulent stresses and heat fluxes and for the mixed moments of the velocity and temperature fluctuations in the continuous and dispersed phases. The models are based on a kinetic equation for the joint probability density of the particle velocity and temperature in an anisotropic turbulent flow. The results are compared with the available direct numerical simulation (DNS) data.     

湍流加速火焰的三维数值模拟

火焰在设有障碍物的管内传播时会自身加速，并可能导致爆炸。本文基于湍流κ-ε模型和改进的EBU—Arrhenius反应模型，对该现象进行了三维空间的数值模拟。计算结果反映了障碍物、湍流和火焰之间相互作用的正反馈机理，描绘了火焰在管内加速传播的三维图像。     

Direct numerical simulation of modulation of isotropic turbulence by poly‐dispersed particles

A direct numerical simulation technique based on two‐way coupling is presented to study a particle‐laden, decaying isotropic turbulent flow. Physical characteristics of turbulence modulation because of the mono‐dispersed (i.e., particles with single Stokes number) and poly‐dispersed particles (i.e., particles with more than one Stokes number) were investigated. A scale dependent effective viscosity that summarizes the aspects of the interaction between the velocity field and particles is defined in the study. Particles of Stokes number (St) 3.2,6.4 and 12.8 were used in performing the simulations. Poly‐dispersed particles were acquired by mixing particles of two different Stokes numbers at a time. As a whole, decay of turbulence because of the poly‐dispersed particles is observed to be larger than that of the decay of turbulence because of the mono‐dispersed particles. Simulations of poly‐dispersed particle indicate nonlinear characteristics in the modification of the temporal evolution of turbulence energy and dissipation. The scale dependent effective viscosity, which correlates with the energy spectrum plot, indicates that the decay of turbulence is mostly observed at the intermediate scales of turbulence. The effective viscosity for the simulations of the poly‐dispersed particles was calculated to be higher than that of the simulations of the mono‐dispersed particles. Copyright © 2011 John Wiley & Sons, Ltd.     

Direct numerical simulation of turbulent pipe flow

Fully developed turbulent pipe flow at low Re-number is studied by means of direct numerical simulation (DNS). In contrast to many previous DNS's of turbulent flows in rectangular geometries, the present DNS code, developed for a cylindrical geometry, is based on the finite volume technique rather than being based on a spectral method. The statistical results are compared with experimental data obtained with two different experimental techniques. The agreement between numerical and experimental results is found to be good which indicates that the present DNS code is suited for this kind of numerical simulations.     

Large eddy simulation of turbulent statistical and transport properties in stably stratified flows

Three dimensional large eddy simulation （LES） is performed in the investigation of stably stratified turbulence with a sharp thermal interface. Main results are focused on the turbulent characteristic scale, statistical properties, transport properties, and temporal and spatial evolution of the scalar field. Results show that the buoyancy scale increases first, and then goes to a certain constant value. The stronger the mean shear, the larger the buoyancy scale. The overturning scale increases with the flow, and the mean shear improves the overturning scale. The flatness factor of temperature departs from the Gaussian distribution in a fairly large region, and its statistical properties are clearly different from those of the velocity fluctuations in strong stratified cases. Turbulent mixing starts from small scale motions, and then extends to large scale motions.     

Direct numerical simulation of driven cavity flows

Direct numerical simulations of 2D driven cavity flows have been performed. The simulations exhibit that the flow converges to a periodically oscillating state at Re=11,000, and reveal that the dynamics is chaotic at Re=22,000. The dimension of the attractor and the Kolmogorov entropy have been computed. Explicit time-integration techniques are discussed.     

Direct numerical simulation of the turbulent energy balance and the shear stresses in power-law fluid flows in pipes

The results of direct numerical simulation of turbulent flows of non-Newtonian pseudoplastic fluids in a straight pipe are presented. The data on the distributions of the turbulent stress tensor components and the shear stress and turbulent kinetic energy balances are obtained for steady turbulent flows at the Reynolds numbers of 10⁴ and 2×10⁴. As distinct from Newtonian fluid flows, the viscous shear stresses turn out to be significant even far from the wall. In power-law fluid flows the mechanism of the energy transport from axial to transverse component fluctuations is suppressed. It is shown that with decrease in the fluid index the turbulent transfer of the momentum and the velocity fluctuations between the wall layer and the flow core reduces, while the turbulent energy flux toward the wall increases. The earlier-proposed models for the average viscosity and the non-Newtonian one-point correlations are in good agreement with the data of direct numerical simulation.     

Direct numerical simulation of homogeneous stratified rotating turbulence

The effects of the Prandtl number on stratified rotating turbulence have been studied in homogeneous turbulence by using direct numerical simulations and a rapid distortion theory. Fluctuations under strong stable-density stratification can be theoretically divided into the WAVE and the potential vorticity (PV) modes. In low-Prandtl-number fluids, the WAVE mode deteriorates, while the PV mode remains. Imposing rotation on a low-Prandtl-number fluid makes turbulence two-dimensional as well as geostrophic; it is found from the instantaneous turbulent structure that the vortices merge to form a few vertically-elongated vortex columns. During the period toward two-dimensionalization, the vertical vortices become asymmetric in the sense of rotation. Communicated by S. Obi PACS 47.55.Hd     

Direct numerical simulation of preferential particle concentration in decaying turbulence under the influence of magnetic field

This paper presents a direct numerical simulation of an isotropic turbulence at Reynolds number R_λ=71.4. The statistics of heavy particles in such turbulence has been investigated. The numerical results indicate the merit of method in solving the magnetohydrodynamic turbulence. The revision of the paper was made based on the reviewers' comments. The paper was improved significantly. Hence, review suggests accepting the submission in current form. Copyright © 2009 John Wiley & Sons, Ltd.     

Direct numerical simulation of confined swirling jets

This study investigates the Lagrangian acceleration and velocity of fluid particles in swirling flows via direct numerical simulation. The intermittency characteristics of acceleration and velocity of fluid particles are investigated at different swirl numbers and Reynolds numbers. The flatness factor and trajectory curvature are used to analyse the effect of Lagrangian intermittency. The joint probability density function of Lagrangian acceleration and turbulence intensity is shown to explain the augmentation effect of Lagrangian intermittency by the strongly swirling levels under the relatively low intensity of turbulence. In addition, the correlation between the Lagrangian acceleration and the turbulence intensity is enhanced as the swirl level increases. It shows the important effect of swirl on the motion behaviour of fluid particles in the strongly swirling flows.     

Numerical simulation of turbulent flows in trapezoidal meandering compound open channels

Three‐dimensional (3D) numerical study is presented to investigate the turbulent flow in meandering compound open channels with trapezoidal cross‐sections. The flow simulation is carried out by solving the 3D Reynolds‐averaged continuity and Navier–Stokes equations with Reynolds stress equation model (RSM) for steady‐state flow. Finite volume method (FVM) is applied to numerically solve the governing equations of fluid flow. The velocity magnitude, tangential velocity, transverse velocity and Reynolds stresses are calculated for various flow conditions. Good agreement between the simulated and available laboratory measurements was obtained, indicating that the RSM can accurately predict the complicated flow phenomenon. Comparison of the calculated secondary currents of four cases (one being inbank flow and other three being overbank flow) with different water depths reveals that (i) the inbank flow exhibits different flow behaviors from that of the overbank flow does and (ii) the water depth has significant effects on the magnitude and direction of secondary currents. Copyright © 2010 John Wiley & Sons, Ltd.     

Applications of a higher‐order bounded numerical scheme to turbulent flows

This paper applies the higher‐order bounded numerical scheme Weighted Average Coefficients Ensuring Boundedness (WACEB) to simulate two‐ and three‐dimensional turbulent flows. In the scheme, a weighted average formulation is used for interpolating the variables at cell faces and the weighted average coefficients are determined from a normalized variable formulation and total variation diminishing (TVD) constraints to ensure the boundedness of the solution. The scheme is applied to two turbulent flow problems: (1) two‐dimensional turbulent flow around a blunt plate; and (2) three‐dimensional turbulent flow inside a mildly curved U‐bend. In the present study, turbulence is evaluated by using a low‐Reynolds number version of the k–ω model. For the flow simulation, the QUICK scheme is applied to the momentum equations while either the WACEB scheme (Method 1) or the UPWIND scheme (Method 2) is used for the turbulence equations. The present study shows that the WACEB scheme has at least second‐order accuracy while ensuring boundedness of the solutions. The present numerical study for a pure convection problem shows that the ‘TVD’ slope ranges from 2 to 4. For the turbulent recirculating flow, two different mixed procedures (Method 1 and Method 2) produce a substantial difference for the mean velocities as well as for the turbulence kinetic energy. Method 1 predicts better results than Method 2 does, comparing the analytical solution and the experimental data. For the turbulent flow inside the mildly curved U‐bend, although the predictions of velocity distributions with two procedures are very close, a noticeable difference of turbulence kinetic energy is exhibited. It is noticed that the discrepancy exists between numerical results and the experimental data. The reason is the limit of the two‐equation turbulence model to such complex turbulent flows with extra strain‐rates. Copyright © 2001 John Wiley & Sons, Ltd.