近壁区湍流脉动压力相关项的理论模型

提出在湍流边界层近壁区采用三维波的理论模型描述湍流相干结构,根据理论模型计算了Ｒｅｙｎｏｌｄｓ应力输运方程中的脉动速度与脉动压力梯度的相关项,理论计算结果与直接数值模拟（ＤＮＳ）符合很好。表明该理论方法不仅有益于对湍流机制的了解,而且可能为湍流的近壁模型化开辟一条新的途径。     

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

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

Near-wall modelling in computations of turbulent flows on unstructured grids

The questions related to the formulation and numerical realization of boundary conditions on a wall at the computation of turbulent flows on unstructured grids are considered. A technique is proposed for realization of weak boundary conditions assuming a non-zero value of the tangent velocity on the wall at a discretization of the Reynolds-averaged Navier — Stokes equations by the control volume method. The capabilities of the developed approach are demonstrated by the example of computing the flow in the inter-blade channel of a low-velocity compressor. The influence of the near-wall grid step on the accuracy of computations, in particular, the pressure distribution near the profile trailing edge is shown, and the solution grid dependence is investigated when using the method of near-wall functions and weak boundary conditions.     

High-order compact schemes for incompressible flows: A simple and efficient method with quasi-spectral accuracy

In this paper, a finite difference code for Direct and Large Eddy Simulation (DNS/LES) of incompressible flows is presented. This code is an intermediate tool between fully spectral Navier–Stokes solvers (limited to academic geometry through Fourier or Chebyshev representation) and more versatile codes based on standard numerical schemes (typically only second-order accurate). The interest of high-order schemes is discussed in terms of implementation easiness, computational efficiency and accuracy improvement considered through simplified benchmark problems and practical calculations. The equivalence rules between operations in physical and spectral spaces are efficiently used to solve the Poisson equation introduced by the projection method. It is shown that for the pressure treatment, an accurate Fourier representation can be used for more flexible boundary conditions than periodicity or free-slip. Using the concept of the modified wave number, the incompressibility can be enforced up to the machine accuracy. The benefit offered by this alternative method is found to be very satisfactory, even when a formal second-order error is introduced locally by boundary conditions that are neither periodic nor symmetric. The usefulness of high-order schemes combined with an immersed boundary method (IBM) is also demonstrated despite the second-order accuracy introduced by this wall modelling strategy. In particular, the interest of a partially staggered mesh is exhibited in this specific context. Three-dimensional calculations of transitional and turbulent channel flows emphasize the ability of present high-order schemes to reduce the computational cost for a given accuracy. The main conclusion of this paper is that finite difference schemes with quasi-spectral accuracy can be very efficient for DNS/LES of incompressible flows, while allowing flexibility for the boundary conditions and easiness in the code development. Therefore, this compromise fits particularly well for very high-resolution simulations of turbulent flows with relatively complex geometries without requiring heavy numerical developments.     

Generalized wall functions for turbulent flows with strong adverse pressure gradient

The generalized wall functions for turbulent flows with strong adverse pressure gradients are derived on the basis of the asymptotic theory of near-wall turbulence. The generalized wall functions have a correct asymptotic behavior in the limit of zero friction velocity and can be applied to computations of flows under a strong adverse pressure gradient and with separation or reattachment. Calculations of a turbulent boundary layer in a strong adverse pressure gradient with the aid of the developed modified k-ɛ model of turbulence and comparison with the experimental data validate the advantages of the generalized wall functions over traditional wall functions based on the logarithmic law of the wall.     

A novel iterative direct-forcing immersed boundary method and its finite volume applications

We present a novel iterative immersed boundary (IB) method in which the body force updating is incorporated into the pressure iterations. Because the body force and pressure are solved simultaneously, the boundary condition on the immersed boundary can be fully verified. The computational costs of this iterative IB method is comparable to those of conventional IB methods. We also introduce an improved body force distribution function which transfers the body force in the discrete volume of IB points to surrounding Cartesian grids totally. To alleviate the demanding computational requirements of a full-resolved direct numerical simulation, a wall-layer model is presented. The accuracy and capability of the present method is verified by a variety of two- and three-dimensional numerical simulations, ranging from laminar flow past a cylinder and a sphere to turbulent flow around a cylinder. The improvement of the iterative IB method is fully demonstrated and the influences of different body force distribution strategies is discussed.     

Finite-size particles in turbulent channel flow: quadrant analysis and acceleration statistics

The interaction between finite-size particles and turbulent channel flow in the absence of gravity is studied by direct numerical simulations (DNS). The study is motivated by DNS observations of a turbulent channel flow with high-density, pointwise particles, that cluster in regions of high streamwise root-mean-square (RMS) acceleration close to the wall, contrary to what is observed in homogeneous isotropic turbulence. The aim of the present study is to explore if this is still the case when size effects are taken into account in the DNS. Based on the analysis of the velocity and acceleration statistics, the present DNS shows that, close to the wall, particles with ρ_p/ρ_f ranging from 2 to 4 are surrounded by regions with low streamwise RMS velocity but high streamwise RMS acceleration. According to the normalised particle acceleration probability density functions (PDFs), size effects become important in the near-wall region. As particle inertia increases, the normalised PDFs of particle acceleration tend to a Gaussian distribution. The tails of the normalised PDFs of the fluid conditioned by the presence of particles are higher than that of the unconditioned fluid close to the wall.     

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary(IB) method has been developed to simulate laminar flows. However,its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work,we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case,we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.     

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.     

A differentially interpolated direct forcing immersed boundary method for predicting incompressible Navier–Stokes equations in time-varying complex geometries

A dispersion-relation-preserving dual-compact scheme developed in Cartesian grids is applied together with the immersed boundary method to solve the flow equations in irregular and time-varying domains. The artificial momentum forcing term applied at certain points in cells containing fluid and solid allows an imposition of velocity condition to account for the motion of solid body. We develop in this study a differential-based interpolation scheme which can be easily extended to three-dimensional simulation. The results simulated from the proposed immersed boundary method agree well with other numerical and experimental results for the chosen benchmark problems. The accuracy and fidelity of the IB flow solver developed to predict flows with irregular boundaries are therefore demonstrated.     

确定分布的展向Lorentz力调制下的槽道湍流涡结构

采用直接数值模拟方法,对槽道湍流中确定分布的Lorentz力的流动控制与减阻问题进行研究.讨论了Lorentz力作用于槽道湍流后,流场的特性和涡结构的特性,并对此类Lorentz力对槽道湍流的控制与减阻机理进行了讨论.研究发现:1)Lorentz力诱导的层流流场壁面附近存在梯度极大的展向速度剪切层,该剪切层容易形成流向涡结构;2)在给定合适参数的确定分布的Lorentz力作用下,湍流流场仅剩周期分布的准流向涡;3)与未控制流场相比,控制后的流场中,准流向涡的抬升高度大大降低,从而减小猝发强度,使壁面阻力下降.     

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.     

A parallel solution – adaptive method for three-dimensional turbulent non-premixed combusting flows

A parallel adaptive mesh refinement (AMR) algorithm is proposed and applied to the prediction of steady turbulent non-premixed compressible combusting flows in three space dimensions. The parallel solution-adaptive algorithm solves the system of partial-differential equations governing turbulent compressible flows of reactive thermally perfect gaseous mixtures using a fully coupled finite-volume formulation on body-fitted multi-block hexahedral meshes. The compressible formulation adopted herein can readily accommodate large density variations and thermo-acoustic phenomena. A flexible block-based hierarchical data structure is used to maintain the connectivity of the solution blocks in the multi-block mesh and to facilitate automatic solution-directed mesh adaptation according to physics-based refinement criteria. For calculations of near-wall turbulence, an automatic near-wall treatment readily accommodates situations during adaptive mesh refinement where the mesh resolution may not be sufficient for directly calculating near-wall turbulence using the low-Reynolds-number formulation. Numerical results for turbulent diffusion flames, including cold- and hot-flow predictions for a bluff-body burner, are described and compared to available experimental data. The numerical results demonstrate the validity and potential of the parallel AMR approach for predicting fine-scale features of complex turbulent non-premixed flames.     

On mean flow universality of turbulent wall flows. I. High Reynolds number flow analysis

ABSTRACT

The universality and mathematical physical structure of wall-bounded turbulent flows is a topic of discussions over many decades. There is no agreement about questions like what is the physical mean flow structure, how universal is it, and how universal are theoretical concepts for local and global flow variations. These questions are addressed by using latest direct numerical simulation (DNS) data at moderate Reynolds numbers Re and experimental data up to extreme Re. The mean flow structure is explained by analytical models for three canonical wall-bounded turbulent flows (channel flow, pipe flow, and the zero-pressure gradient turbulent boundary layer). Thorough comparisons with DNS and experimental data provide support for the validity of models. Criteria for veritable physics derived from observations are suggested. It is shown that the models presented satisfy these criteria. A probabilistic interpretation of the mean flow structure shows that the physical constraints of equal entropies and equally likely mean velocity values in a region unaffected by boundary effects impose a universal log-law structure. The structure of wall-bounded turbulent flows is much more universal than previously expected. There is no discrepancy between local logarithmic velocity variations and global friction law and bulk velocity variations. Flow effects are limited to the minimum: the difference of having a bounded or unbounded domain, and the variation range of mean velocity values allowed by the geometry.     

Velocity,temperature, and Reynolds-stress scaling in the wall region of turbulent boundary layer on a permeable surface

A finite total number of flow parameters in the wall region of a turbulent boundary layer points to universal behavior of turbulent shear stress as a function of mean-velocity gradient and turbulent heat flux as a function of both mean-velocity and mean-temperature gradients. Combined with dimensional arguments, this fact is used to reduce the momentum and heat equations to first-order ordinary differential equations for temperature and velocity profiles amenable to general analysis. Scaling laws for velocity and temperature in boundary layer flows with transpiration are obtained as generalizations of well-known logarithmic laws. Scaling relations are also established for shear stress and rms transverse velocity fluctuation. The proposed method has substantial advantages as compared to the classical approach (which does not rely on fluid-dynamics equations [1–3]). It can be applied to establish scaling laws for a broader class of near-wall turbulence problems without invoking closure hypotheses.     

Simulation of Incompressible Viscous Flows by Local DFD-Immersed Boundary Method

A local domain-free discretization-immersed boundary method (DFD-IBM) is presented in this paper to solve incompressible Navier-Stokes equations in the primitive variable form. Like the conventional immersed boundary method (IBM), the local DFD-IBM solves the governing equations in the whole domain including exterior and interior of the immersed object. The effect of immersed boundary to the surrounding fluids is through the evaluation of velocity at interior and exterior dependent points. To be specific, the velocity at interior dependent points is computed by approximate forms of solution and the velocity at exterior dependent points is set to the wall velocity. As compared to the conventional IBM, the present approach accurately implements the non-slip boundary condition. As a result, there is no flow penetration, which is often appeared in the conventional IBM results. The present approach is validated by its application to simulate incompressible viscous flows around a circular cylinder. The obtained numerical results agree very well with the data in the literature.     

Direct numerical simulation of shock/turbulent boundary layer interaction in a supersonic compression ramp

A direct numerical simulation of the shock/turbulent boundary layer interaction flow in a supersonic 24-degree compression ramp is conducted with the free stream Mach number 2.9.The blow-and-suction disturbance in the upstream wall boundary is used to trigger the transition.Both the mean wall pressure and the velocity profiles agree with those of the experimental data,which validates the simulation.The turbulent kinetic energy budget in the separation region is analyzed.Results show that the turbulent production term increases fast in the separation region,while the turbulent dissipation term reaches its peak in the near-wall region.The turbulent transport term contributes to the balance of the turbulent conduction and turbulent dissipation.Based on the analysis of instantaneous pressure in the downstream region of the mean shock and that in the separation bubble,the authors suggest that the low frequency oscillation of the shock is not caused by the upstream turbulent disturbance,but rather the instability of separation bubble.     

Sharp interface immersed-boundary/level-set method for wave–body interactions

A sharp interface Cartesian grid method for the large-eddy simulation of two-phase turbulent flows interacting with moving bodies is presented. The overall approach uses a sharp interface immersed boundary formulation and a level-set/ghost–fluid method for solid–fluid and fluid–fluid interface treatments, respectively. A four-step fractional-step method is used for velocity–pressure coupling, and a Lagrangian dynamic Smagorinsky subgrid-scale model is adopted for large-eddy simulations. A simple contact angle boundary condition treatment that conforms to the immersed boundary formulation is developed. A variety of test cases of different scales ranging from bubble dynamics, water entry and exit, landslide-generated waves, to ship hydrodynamics are performed for validation. Extensions for high Reynolds number ship flows using wall-layer models are also considered.     

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.     

Boundary conditions on a compliant wall in the turbulent flow

It is shown that with an increase in elasticity modulus of the coating material from 0.01 MPa to 100 MPa the pulsation velocity of coating surface changes not more than by 30 % and equals 0.17–0.24 of the value of dynamic flow velocity, and this can change significantly Reynolds stresses in the near-wall area. According to performed analysis, it was found out that the deformation value of the compliant coatings within the velocity range, optimal for their interaction with the turbulent flow, is only several units of the viscous scale. Moreover, these deformations are very gentle: the ratio of amplitude of deformation to wavelength is less than 10^?3. It is assumed that while modeling the interaction between the compliant coating and turbulent flow it is not necessary to transfer the boundary conditions to the moving coating surface. Perhaps, it will be sufficient to determine the velocity of wall motion over the undisturbed coating surface.