Direct numerical simulation of a fully developed compressible wall turbulence over a wavy wall

Compressible turbulent channel flow over a wavy surface is investigated by direct numerical simulations using high-resolution finite difference schemes. The Reynolds number considered in the present paper is 3380 based on the bulk velocity, the channel half-width and the kinetic viscosity at the wall. Four test cases are simulated and analysed at Ma_m = 0.33, 0.8, 1.2, 1.5 based on the bulk velocity and the speed of sound at the wall. We mainly focus on the curvature and the Mach number effects on the compressible turbulent flows. Numerical results show that although the wavy wall has effects on the mean and fluctuation quantities, log law still exists in the distribution of the wave-averaged streamwise velocity if the roughness effects are taken into consideration in the scaling of it. Near-wall streaks are broken by the wavy surface and near-wall quasi-streamwise vortices mostly begin at the upslope of the wave and pass over the crest of it. The wavy wall makes the turbulence more active and the flow easier to be blended. From the viewpoint of turbulent kinetic budgets, curvature effects strengthen both the diffusion terms and the dissipation terms. At the same time, they change the properties of the compressibility-related terms and promote more inner energy transferring into turbulent kinetic energy. As the Mach number increases, the reattachment of the mean flow is delayed, which indicates the mean separation bubble becomes larger. Concerning the near-wall coherent structures, the vortices are more sparsely distributed with the increasing of the Mach number. For the supersonic cases, shock waves appear. Though they have little effects on the mean turbulent quantities, they change the structures of the flow fields and induce local separations at the upper wall of the channel.     

Large Eddy Simulation of Turbulent Channel Flow with 3D Roughness Using a Roughness Element Model

Large eddy simulation of turbulent channel flow with dense and small 3D roughness elements is carried out using a roughness element model Profiles of mean Reynolds stress, mean velocity and rms velocity as well as turbulent structures near the wall are obtained. The shear stress in the rough wall is larger than that in the smooth wall side and the rough wall has a larger influence on the channel flow. Profiles of mean streamwise velocity near the wall have logarithmic velocity distributions for both smooth and roughness walls, while there is a velocity decrease for the rough wall due to larger fractional drag. All the three components of rms velocities in the rough wall region are larger than that in the smooth wall region, and the roughness elements on the wall increase turbulent intensity in all directions. The s~reak spacing and average diameter of near wall quasi-s~reamwise vortices increase with the presence of roughness elements on the wall and it is shown that the rough wall induces complex and strong streamwise vortices. Results of dense and small 3D roughness elements in both turbulent statistics and structure, obtained with a relatively simple method, are found to be comparable to related experiments.     

Turbulent Boundary Layer Control via a Streamwise Travelling Wave Induced by an External Force

Turbulent boundary layer control via a streamwise travelling wave is investigated based on direct numerical simulation of an incompressible turbulent channel flow. The streamwise travelling wave is induced on one side wall of the channel by a spanwise external force, e.g., Lorenz force, which is con~ned in the viscous sublayer. As the control strategy used in this study has never been examined, we pay our attention to its efficiency of drag control. It is revealed that the propagating direction of the travelling wave, i.e., the downstream or upstream propagating direction with respect to the streamwise flow, has an important role on the drag control, leading to a significant drag reduction or enhancement for the parameters considered. The coherent structures of turbulent boundary layer are altered and the underlying mechanisms are analysed. The results obtained provide physical insight into the understanding of turbulent boundary layer control.     

Density-PDFs and Lagrangian statistics of highly compressible turbulence

In isothermal, highly compressible turbulent flows, density fluctuations follow a log-normal distribution. We establish a connection between these density fluctuations and the probability-density-functions (PDF) of Lagrangian tracer particles advected with the flow. Our predicted particle statistics is tested against large scale numerical simulations, which were performed with 512³ collocation points and 2 million tracer particles integrated over several dynamical times.     

Quasi-stationary states in a circular geometry

The quasi-stationary final states of decaying two-dimensional turbulence on a circular domain, with an initial flow field containing either no or a substantial amount of angular momentum, have been investigated numerically. The production of angular momentum is almost absent for these flows on a circular domain with a no-slip wall. Its presence or absence essentially determines the character of the quasi-stationary final state. Based on a minimum-enstrophy principle a diagram is constructed that provides insight into the development of the typical late-time flow patterns on a circular domain with a no-slip wall. The quasi-stationary final states found in the present numerical study can be understood based on the predictions from the minimum-enstrophy principle.     

Towards lowering dissipation bounds for turbulent flows

We examine the background flow variational principle for calculating bounds on the energy dissipation rate in turbulent shear flow, and suggest to select this principle's test functions such that they comply with the small-scale smoothness of real turbulent velocity fields. A self-consistent algorithm implementing this requirement then yields an upper bound on the dimensionless dissipation coefficient which shows a weak power-law decrease at high Reynolds numbers, instead of approaching a nonzero constant, as it did in previous estimates. Received 26 October 1998     

Evaporating droplets in turbulent reacting flows

Three-dimensional direct numerical simulations are carried out to determine the effects of turbulence on the preferential segregation of an evaporating spray and then to study the evolution of the resulting mixture fraction topology and propagating flame. First, the mixing between an initially randomly dispersed phase and the turbulent gaseous carrier phase is studied with non-evaporating particles. According to their inertia and the turbulence properties, the formation of clusters of particles is analyzed (formation delay, cluster characteristic size and density). Once the particles are in dynamical equilibrium with the surrounding turbulent flow, evaporation is considered through the analysis of the mixture fraction evolution. Finally, to mimic ignition, a kernel of burnt gases is generated at the center of the domain and the turbulent flame evolution is described.     

Lattice Boltzmann model for three-dimensional decaying homogeneous isotropic turbulence

We implement a lattice Boltzmann method (LBM) for decaying homogeneous isotropic turbulence based on an analogous Galerkin filter and focus on the fundamental statistical isotropic property. This regularized method is constructed based on orthogonal Hermite polynomial space. For decaying homogeneous isotropic turbulence, this regularized method can simulate the isotropic property very well. Numerical studies demonstrate that the novel regularized LBM is a promising approximation of turbulent fluid flows, which paves the way for coupling various turbulent models with LBM.     

Frequency Response of Near-Wall Coherent Structures to Localized Periodic Blowing and Suction in Turbulent Boundary Layer

We experimentally investigate the frequency response of near-wall coherent structures to localized periodic blow- ing and suction through a spanwise slot in a turbulent boundary layer by changing the frequency of periodic disturbance at similar velocities of free stream. The effects of blowing and suction disturbance on energy redistri- bution, turbulent intensity U^＋ rums over y^＋ and waveforms of phase-averaged velocity during sweeping process are respectively discussed under three frequencies of periodic blowing and suction in near-wall region of turbulent boundary layer, compared with those in a standard turbulent boundary layer. The most effective disturbance frequency is figured out in this system.     

Large scale correlations for energy injection mechanisms in swirling turbulent flows

We report an experimental study of large scale correlations in the power injected in turbulent swirling flows generated in the gap between two coaxial rotating disks. We measure the pressure fluctuations on the blades of one disk, as well as the pressure drop between the leading and the trailing edges of the rotating blades, i.e. the local drag force. Measurements at different positions on one blade and on two successive blades display a correlation length much larger than the ones usually expected in turbulent flows. The time lag for which the correlation between two points is maximum, strongly depends on the global flow configuration. These results help us to understand the statistical properties of the injected power fluctuations in turbulent swirling flows. Received 2 September 1999     

Nonperturbative renormalization group approach to turbulence

We suggest a new, renormalization group (RG) based, nonperturbative method for treating the intermittency problem of fully developed turbulence which also includes the effects of a finite boundary of the turbulent flow. The key idea is not to try to construct an elimination procedure based on some assumed statistical distribution, but to make an ansatz for possible RG transformations and to pose constraints upon those, which guarantee the invariance of the nonlinear term in the Navier-Stokes equation, the invariance of the energy dissipation, and other basic properties of the velocity field. The role of length scales is taken to be inverse to that in the theory of critical phenomena; thus possible intermittency corrections are connected with the outer length scale. Depending on the specific type of flow, we find different sets of admissible transformations with distinct scaling behaviour: for the often considered infinite, isotropic, and homogeneous system K41 scaling is enforced, but for the more realistic plane Couette geometry no restrictions on intermittency exponents were obtained so far. Received: 28 December 1997 / Accepted: 6 August 1998     

Direct Numerical Simulation of a Spatially Evolving Supersonic Turbulent Boundary Layer at Ma = 6

Direct numerical simulation is carried out for a spatially evolving supersonic turbulent boundary layer at freestream Mach number 6. To overcome numerical instability, the seventh-order WENO scheme is used for the convection terms of Navier-Stokes equations, and fine mesh is adopted to minimize numerical dissipation. Compressibility effects on the near-wall turbulent kinetic energy budget are studied. The cross-stream extended self-similarity and scaling exponents including the near-wall region are studied. In high Mach number flows, the coherence vortex structures are arranged to be smoother and streamwised, and the hair-pin vortices are less likely tO OCCUr.     

Large eddy simulation of compressible turbulent channel flow with spanwise wall oscillation

The influences of the modification of turbulent coherent structures on temperature field and heat transfer in turbulent channel flow are studied using large eddy simulation (LES) of compressible turbulent channel flows with spanwise wall oscillation (SWO). The reliability of the LES on such problems is proved by the comparisons of the drag reduction data with those of other researches. The high consistency of coherent velocity structures and temperature structures is found based on the analyses of the turbulent flow field. When the coherent velocity structures are suppressed, the transportations of momentum and heat are reduced simultaneously, demonstrating the same trend. This shows that the turbulent coherent structures have the same effects on the transportations of momentum and heat. The averaged wall heat flux can be reduced with appropriate oscillating parameters. Supported by the Key Subjects of National Natural Science Foundation of China (Grant No. 10732090), the National Natural Science Foundation of China (Grant No. 50476004), and the 111 Project (Grant No. B08009)     

Multi-fractal formulation of compressible fully developed turbulence: Parabolic-profile approximation for the singularity spectrum

A parabolic-profile approximation (PPA) for the singularity spectrum D (h) in the multi-fractal model for compressible fully developed turbulence (FDT) is considered and is then extrapolated to the Kolmogorov microscale regime. The generalization of Kolmogorov’s “4/5th law” relating the third-order velocity structure function to the mean energy dissipation rate ε to compressible FDT is considered. The PPA is also shown to afford, unlike the generic multi-fractal model, an analytical calculation of probability distribution function (PDF) of velocity gradients and to describe intermittency corrections for this PDF that complement those provided by homogeneous-fractal model.     

一种数值模拟S形进气道内可压缩湍流的隐式迎风增量迭代法

在这篇文章里,我们提出了一个二阶隐式迎风的增量迭代法,用以离散BFC(Body-Fitted Coordinates)坐标变换下的Navier-Stokes方程。使用这种算法结合Lomax两层代数湍流模型数值求解了S形进气道内的可压缩湍流。数值实验表明:计算过程具有较好的稳定性和收敛性。二维Poissuil流动和Couette流动的计算结果和精确解符合良好,两种壁面曲率不同的S形进气道内流场的计算结果也是令人满意的。     

Drag reduction of compressible wall turbulence with active dimples

Direct numerical simulations are carried out to assess the potential drag reduction of compressible turbulent flow between isothermal walls.For the sake of achieving drag reduction,the flow is actively controlled by deformable dimples lying on the bottom wall of the channel.The first stage of the procedure consists in assessing the optimum geometry of the dimples.In this regard,the lower wall is allowed to freely deform itself according to the loop of control.This method is called the smart wall approach in...     

Two Phases of Coherent Structure Motions in Turbulent Boundary Layer

Two phases of coherent structure motion are acquired after obtaining conditional phase-averaged waveforms for longitudinal velocity of coherent structures in turbulent boundary layer based on Harr wavelet transfer. The correspondences of the two phases to the two processes （i.e. ejection and sweep） during a burst are determined.     

Acceleration statistics of inertial particles in turbulent flow

Turbulent transport of material inclusions plays an important role in many natural and industrial situations. Being able to accurately model and predict the dynamics of dispersed particles transported by a turbulent carrier flow, remains a challenge. One critical and difficult point is to develop models which correctly describe the dynamics of particles over a wide range of sizes and densities. Our measurements show that acceleration statistics of particles dispersed in a turbulent flow do exhibit specific, and so far unpredicted, size and density effects and that they preserve an extremely robust turbulent signature with lognormal fluctuations, regardless of particles size and density. This has important consequences in terms of modeling for the turbulent transport of dispersed inclusions.     

An experimental study of impulsively started turbulent axisymmetric jets

An impulsively started turbulent jet injected into quiescent surroundings with a constant inlet velocity has been studied experimentally. Results show that the jet length increases linearly with the square-root of time, over a wide range of Reynolds number calculated with respect to the jet diameter. The celerity factor, x_f/t U, has been found to be nearly constant at 2.47 throughout with a 5% variance. Here, x_f is the jet length, t is the time and U is the jet exit velocity. These results compare favourably with earlier results reported at lower Reynolds numbers. Finally, we present a simple model based on the integral energy balance of the turbulent boundary layer equation for an impulsively started turbulent axisymmetric jet. The model predicts a jet length that scales as, where d is the nozzle diameter and B(≈6.0) is the velocity-decay constant. This gives a celerity factor, in close agreement with the experiments.     

Derivation of a Nonlinear Reynolds Stress Model Using Renormalization Group Analysis and Two-Scale Expansion Technique

Adopting Yoshizawa＇s two-scale expansion technique, the fluctuating field is expanded around the isotropic field. The renormalization group method is applied for calculating the covariance of the fluctuating field at the lower order expansion. A nonlinear Reynolds stress model is derived and the turbulent constants inside are evaluated analytically. Compared with the two-scale direct interaction approximation analysis for turbulent shear flows proposed by Yoshizawa, the calculation is much more simple. The analytical model presented here is close to the Speziale model, which is widely applied in the numerical simulations for the complex turbulent flows.