A Web services accessible database of turbulent channel flow and its use for testing a new integral wall model for LES

The output from a direct numerical simulation (DNS) of turbulent channel flow at Re_τ ≈ 1000 is used to construct a publicly and Web services accessible, spatio-temporal database for this flow. The simulated channel has a size of 8πh × 2h × 3πh, where h is the channel half-height. Data are stored at 2048 × 512 × 1536 spatial grid points for a total of 4000 time samples every 5 time steps of the DNS. These cover an entire channel flow-through time, i.e. the time it takes to traverse the entire channel length 8πh at the mean velocity of the bulk flow. Users can access the database through an interface that is based on the Web services model and perform numerical experiments on the slightly over 100 terabytes (TB) DNS data on their remote platforms, such as laptops or local desktops. Additional technical details about the pressure calculation, database interpolation, and differentiation tools are provided in several appendices. As a sample application of the channel flow database, we use it to conduct an a-priori test of a recently introduced integral wall model for large eddy simulation of wall-bounded turbulent flow. The results are compared with those of the equilibrium wall model, showing the strengths of the integral wall model as compared to the equilibrium model.     

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.     

Large eddy simulation of channel flow with and without periodic constrictions using the explicit algebraic subgrid-scale model

We analyse the performance of the explicit algebraic subgrid-scale (SGS) stress model (EASSM) in large eddy simulation (LES) of plane channel flow and the flow in a channel with streamwise periodic hill-shaped constrictions (periodic hill flow) which induce separation. The LESs are performed with the Code_Saturne which is an unstructured collocated finite volume solver with a second-order spatial discretisation suitable for LES of incompressible flow in complex geometries. At first, performance of the EASSM in LES of plane channel flow at two different resolutions using the Code_Saturne and a pseudo-spectral method is analysed. It is observed that the EASSM predictions of the mean velocity and Reynolds stresses are more accurate than the conventional dynamic Smagorinsky model (DSM). The results with the pseudo-spectral method were, in general, more accurate. In the second step, LES with the EASSM of flow separation in the periodic hill flow is compared to LES with the DSM, no SGS model and a highly resolved LES data using the DSM. Results show that the mean velocity profiles, the friction and pressure coefficients, the length and shape of the recirculation bubble, as well as the Reynolds stresses are considerably better predicted by the EASSM than the DSM and the no SGS model simulations. It was also observed that in some parts of the domain, the resolved strain-rate and SGS shear stress have the same sign. The DSM cannot produce a correct SGS stress in this case, in contrast to the EASSM.     

Periodic Filtering as a subgrid-scale model for LES of laminar separation bubble flows

Laminar separation bubbles develop over many blades and airfoils at moderate angles of attack and Reynolds numbers ranging from 10⁴ to 10⁵. More accurate simulation tools are necessary to enable higher fidelity design optimisation for these airfoils and blades as well as to test flow control schemes. Following previous investigators, an equivalent problem is formulated by imposing suitable boundary conditions for flow over a flat plate which allows to use a high accuracy spectral solver. Large eddy simulation (LES) of such a flow were performed at drastically reduced resolution to assess the accuracy of several LES modelling approaches: the explicit dynamic Smagorinsky model, implicit LES, and the truncated Navier–Stokes approach (TNS). To mimic dissipation that occurs in implicit LES, the solution on a coarse mesh is filtered at every time step and two different filter strengths are used. In the TNS approach, the solution is filtered periodically, every few hundred time steps. The performance of each approach is evaluated against benchmark direct numerical simulation (DNS) data focusing on pressure and skin friction distributions, which are critical to airfoil designers. TNS results confirm that periodic filtering can act as an apt substitute for explicit subgrid-scale models, whereas filtering at every time step demonstrates the dependence of implicit LES on details of numerics.     

Use of DES in mildly separated internal flow: dimples in a turbulent channel

Detached eddy simulation (DES) is investigated as a means to study an array of shallow dimples with depth to diameter ratios of 1.5% and 5% in a turbulent channel. The DES captures large-scale flow features relatively well, but is unable to predict skin friction accurately due to flow modelling near the wall. The current work instead relies on the accuracy of DES to predict large-scale flow features, as well as its well-documented reliability in predicting flow separation regions to support the proposed mechanism that dimples reduce drag by introducing spanwise flow components near the wall through the addition of streamwise vorticity. Profiles of the turbulent energy budget show the stabilising effect of the dimples on the flow. The presence of flow separation however modulates the net drag reduction. Increasing the Reynolds number can reduce the size of the separated region and experiments show that this increases the overall drag reduction.     

A Reynolds stress model for turbulent flows of viscoelastic fluids

A second-order closure is developed for predicting turbulent flows of viscoelastic fluids described by a modified generalised Newtonian fluid model incorporating a nonlinear viscosity that depends on a strain-hardening Trouton ratio as a means to handle some of the effects of viscoelasticity upon turbulent flows. Its performance is assessed by comparing its predictions for fully developed turbulent pipe flow with experimental data for four different dilute polymeric solutions and also with two sets of direct numerical simulation data for fluids theoretically described by the finitely extensible nonlinear elastic – Peterlin model. The model is based on a Newtonian Reynolds stress closure to predict Newtonian fluid flows, which incorporates low Reynolds number damping functions to properly deal with wall effects and to provide the capability to handle fluid viscoelasticity more effectively. This new turbulence model was able to capture well the drag reduction of various viscoelastic fluids over a wide range of Reynolds numbers and performed better than previously developed models for the same type of constitutive equation, even if the streamwise and wall-normal turbulence intensities were underpredicted.     

Drag reduction in turbulent channel flow using bidirectional wavy Lorentz force

Turbulent control and drag reduction in a channel flow via a bidirectional traveling wave induced by spanwise oscillating Lorentz force have been investigated in the paper. The results based on the direct numerical simulation (DNS) indicate that the bidirectional wavy Lorentz force with appropriate control parameters can result in a regular decline of near-wall streaks and vortex structures with respect to the flow direction, leading to the effective suppression of turbulence generation and significant reduction in skin-friction drag. In addition, experiments are carried out in a water tunnel via electro-magnetic (EM) actuators designed to produce the bidirectional traveling wave excitation as described in calculations. As a result, the actual substantial drag reduction is realized successfully in these experiments.     

高聚物减阻溶液对壁湍流输运过程的影响

基于相同雷诺数下清水和高分子聚合物溶液壁湍流的高时间分辨率粒子图像测速技术(time-resolved particle image velocimetry, TRPIV)的对比实验, 从高聚物溶液对湍流边界层动量能量输运影响的角度分析其减阻的机理. 对比两者的雷诺应力发现高聚物的存在抑制了湍流输运过程. 这一影响与高聚物对壁湍流中占主导地位的涡旋运动和低速条带等相干结构的作用密切相关. 运用条件相位平均、相关函数和线性随机估计(linear stochastic estimation, LSE)等方法, 分析提取了高聚物溶液流场中的发卡涡和发卡涡包等典型相干结构的空间拓扑形态. 相比于清水, 高聚物溶液中相干结构的流向尺度增大, 涡旋运动的发展及低速流体喷射的强度受到削弱, 表明了添加的高聚物阻碍了湍流原有的能量传递和自维持的机理. 正是通过影响相干结构, 高聚物抑制了湍流边界层中近壁区与外区之间的动量和能量输运, 使得湍流的无序性降低, 从而减小了湍流流动的阻力.     

A priori evaluation of large eddy simulation subgrid-scale scalar flux models in isotropic passive-scalar and anisotropic buoyancy-driven homogeneous turbulence

Direct numerical simulation (DNS) of passive (non-buoyant) and active (buoyant) scalar homogeneous turbulence is carried out using a standard pseudo-spectral numerical method. The flow settings simulated include stationary forced and decaying passive-scalar turbulence, as well as decaying anisotropic active-scalar turbulence. The Schmidt number is unity in all cases. The results are compared with, and are found to be in very good agreement with, previous similar DNS studies. The well-validated DNS data are divided into 19 sets, and are employed to study different large eddy simulation (LES) subgrid-scale (SGS) models for the SGS scalar flux. The models examined include three eddy-viscosity-type models (Smagorinsky, Vreman and Sigma with a constant SGS Schmidt number), a Dynamic Structure model and two versions of the Gradient (Gradient and Modulated Gradient) model. The models are investigated with respect to their ability to predict the orientation, and the magnitude, of the SGS scalar flux. Eddy-viscosity models are found to predict the magnitude of the SGS scalar flux accurately, but are poor at predicting the orientation of the SGS scalar flux. The Dynamic Structure and Gradient models are better than eddy-viscosity models at predicting both the magnitude and direction. However, neither of them can be realised in an actual LES, without carrying additional transport equations. Based on these observations, four new models are proposed – combining directions from Dynamic Structure and Gradient models, and magnitudes from Smagorinsky and Vreman eddy-viscosity models. These models are expected to be better than eddy-viscosity and Modulated Gradient models, and this is confirmed by preliminary a posteriori tests.     

Particle subgrid scale modelling in large-eddy simulations of particle-laden turbulence

This study is concerned with particle subgrid scale (SGS) modelling in large-eddy simulations (LESs) of particle-laden turbulence. Although many particle-laden LES studies have neglected the effect of the SGS on the particles, several particle SGS models have been proposed in the literature. In this research, the approximate deconvolution method (ADM) and the stochastic models of Fukagata et al. (Dynamics of Brownian particles in a turbulent channel flow, Heat Mass Transf. 40 (2004), 715–726) Shotorban and Mashayek (A stochastic model for particle motion in large-eddy simulation, J. Turbul. 7 (2006), 1–13) and Berrouk et al. (Stochastic modelling of inertial particle dispersion by subgrid motion for LES of high Reynolds number pipe flow, J. Turbul. 8 (2007), pp. 1–20) are analysed. The particle SGS models are assessed using both a priori and a posteriori simulations of inertial particles in a periodic box of decaying, homogeneous and isotropic turbulence with an initial Reynolds number of Re_λ = 74. The model results are compared with particle statistics from a direct numerical simulation (DNS). Particles with a large range of Stokes numbers are tested using various filter sizes and stochastic model constant values. Simulations with and without gravity are performed to evaluate the ability of the models to account for the crossing trajectory and continuity effects. The results show that ADM improves results but is only capable of recovering a portion of the SGS turbulent kinetic energy. Conversely, the stochastic models are able to recover sufficient SGS energy, but show a large range of results dependent on the Stokes number and filter size. The stochastic models generally perform best at small Stokes numbers, but are unable to predict preferential concentration.     

Relevance of the subgrid-scales for large eddy simulations of turbulence-radiation interactions in a turbulent plane jet

An a priori study based on direct numerical simulation (DNS) of a non-isothermal turbulent plane jet has been carried out in order to analyse the role of the small-scales of turbulence on thermal radiation. Filtered DNS and large eddy simulation (LES) without subgrid-scale (SGS) model have been estimated for the radiative heat transfer. The comparison of the results highlights the subgrid-scale influence over the filtered radiation quantities, such as the radiative intensity and the radiative emission. The influence of the optical thickness is also studied. It is shown that the subgrid-scales are not significant near the centerline of the jet, where the radiative heat transfer is more important, and therefore that the SGS can be neglected in this configuration. However, when the optical thickness increases, the SGS become relevant and SGS modeling may be needed.     

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

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

Visualisation and analysis of large-scale vortex structures in three-dimensional turbulent lid-driven cavity flow

In this paper, large eddy simulation (LES) of a three-dimensional turbulent lid-driven cavity (LDC) flow at Re = 10,000 has been performed using the multiple relaxation time lattice Boltzmann method. A Smagorinsky eddy viscosity model was used to represent the sub-grid scale stresses with appropriate wall damping. The prediction for the flow field was first validated by comparing the velocity profiles with previous experimental and LES studies, and then subsequently used to investigate the large-scale three-dimensional vortical structures in the LDC flow. The instantaneous three-dimensional coherent structures inside the cavity were visualised using the second invariant (Q), Δ criterion, λ₂ criterion, swirling strength (λ_ci) and streamwise vorticity. The vortex structures obtained using the different criteria in general agree well with each other. However, a cleaner visualisation of the large vortex structures was achieved with the λ_ci criterion and also when the visualisation is based on the vortex identification criteria expressed in terms of the swirling strength parameters. A major objective of the study was to perform a three-dimensional proper orthogonal decomposition (POD) on the fluctuating velocity fields. The higher energy POD modes efficiently extracted the large-scale vortical structures within the flow which were then visualised with the swirling strength criterion. Reconstruction of the instantaneous fluctuating velocity field using a finite number of POD modes indicated that the large-scale vortex structures did effectively approximate the large-scale motion. However, such a reduced order reconstruction of the flow based on the large-scale vortical structures was clearly not as effective in predicting the small-scale details of the fluctuating velocity field which relate to the turbulent transport.     

Studies of the flow and turbulence fields in a turbulent pulsed jet flame using LES/PDF

A turbulent piloted jet flame subject to a rapid velocity pulse in its fuel jet inflow is proposed as a new benchmark case for the study of turbulent combustion models. In this work, we perform modelling studies of this turbulent pulsed jet flame and focus on the predictions of its flow and turbulence fields. An advanced modelling strategy combining the large eddy simulation (LES) and the probability density function (PDF) methods is employed to model the turbulent pulsed jet flame. Characteristics of the velocity measurements are analysed to produce a time-dependent inflow condition that can be fed into the simulations. The effect of the uncertainty in the inflow turbulence intensity is investigated and is found to be very small. A method of specifying the inflow turbulence boundary condition for the simulations of the pulsed jet flame is assessed. The strategies for validating LES of statistically transient flames are discussed, and a new framework is developed consisting of different averaging strategies and a bootstrap method for constructing confidence intervals. Parametric studies are performed to examine the sensitivity of the predictions of the flow and turbulence fields to model and numerical parameters. A direct comparison of the predicted and measured time series of the axial velocity demonstrates a satisfactory prediction of the flow and turbulence fields of the pulsed jet flame by the employed modelling methods.     

Large eddy simulation of orientation and rotation of ellipsoidal particles in isotropic turbulent flows

The rotational motion and orientational distribution of ellipsoidal particles in turbulent flows are of significance in environmental and engineering applications. Whereas the translational motion of an ellipsoidal particle is controlled by the turbulent motions at large scales, its rotational motion is determined by the fluid velocity gradient tensor at small scales, which raises a challenge when predicting the rotational dispersion of ellipsoidal particles using large eddy simulation (LES) method due to the lack of subgrid scale (SGS) fluid motions. We report the effects of the SGS fluid motions on the orientational and rotational statistics, such as the alignment between the long axis of ellipsoidal particles and the vorticity, the mean rotational energy at various aspect ratios against those obtained with direct numerical simulation (DNS) and filtered DNS. The performances of a stochastic differential equation (SDE) model for the SGS velocity gradient seen by the particles and the approximate deconvolution method (ADM) for LES are investigated. It is found that the missing SGS fluid motions in LES flow fields have significant effects on the rotational statistics of ellipsoidal particles. Alignment between the particles and the vorticity is weakened; and the rotational energy of the particles is reduced in LES. The SGS-SDE model leads to a large error in predicting the alignment between the particles and the vorticity and over-predicts the rotational energy of rod-like particles. The ADM significantly improves the rotational energy prediction of particles in LES.     

A priori analysis of subgrid scale pressure and heat flux in high pressure mixing and reacting shear layers

Direct Numerical Simulation (DNS) data on high pressure H₂/O₂ and H₂/air flames using the compressible flow formulation, detailed kinetics, a real fluid equation of state, and generalised diffusion are analysed. The DNS is filtered over a range of filter widths to provide exact terms in the Large Eddy Simulation (LES) governing equations, including unclosed terms. The filtered pressure and the filtered heat flux vector are extensively compared with the pressure and the heat flux vector calculated as a function of the filtered primitive variables (i.e. the exact LES term is compared with its form available within an actual LES). The difference between these forms defines the subgrid pressure and the subgrid heat flux vector. The analyses are done both globally across the entire flame, as well as by conditionally averaging over specific regions of the flame; including regions of large subgrid kinetic energy, subgrid scalar dissipation, subgrid temperature variance, flame temperature, etc. In this work, although negligible for purely mixing cases, the gradient of the subgrid pressure is shown to be of the same order as, and larger than, the corresponding divergence of the turbulent subgrid stresses for reacting cases. This is despite the fact that all species behave essentially as ideal gases for this flame and holds true even when the ideal gas law is used to calculate the pressure. The ratio of the subgrid pressure gradient to the subgrid stress tensor divergence is shown to increase with increasing Reynolds number. Both the subgrid heat flux vector and its divergence are found to be substantially larger in reacting flows in comparison with mixing due to the associated larger temperature gradients. However, the divergence of the subgrid heat flux vector tends to be significantly smaller than other unclosed terms in the energy equation with decreasing significance with increasing Reynolds number.     

Use of Lagrangian statistics for the analysis of the scale separation hypothesis in turbulent channel flow

Turbulence models often involve Reynolds averaging, with a closure providing the Reynolds stress tensor as function of mean velocity gradients, through a turbulence constitutive equation. The main limitation of this linear closure is that it rests on an analogy with kinetic theory. For this analogy to be valid there has to be a scale separation between the mean velocity variations and the turbulent Lagrangian free path whose mean value is the turbulent mixing length. The aim of this work is to better understand this hypothesis from a microscopic point of view. Therefore, fluid elements are tracked in a turbulent channel flow. The flow is resolved by direct numerical simulation (DNS). Statistics on particle trajectories ending on a certain distance y₀ from the wall are computed, leading to estimations of the turbulent mixing length scale and the Knudsen number. Comparing the computed values to the Knudsen number in the case of scale separation, we may know in which region of the flow and to what extent the turbulence constitutive equation is not verified. Finally, a new non-local formulation for predicting the Reynolds stress is proposed.     

On minimum aspect ratio for duct flow facilities and the role of side walls in generating secondary flows

To the surprise of some of our colleagues, we recently recommended aspect ratios of at least 24 (instead of accepted values over last few decades ranging from 5 to 12) to minimise effects of side walls in turbulent duct flow experiments, in order to approximate the two-dimensional channel flow. Here we compile available results from hydraulics and civil engineering literature, where this was already documented in the 1980s. This is of great importance due to the large amount of computational studies (mainly direct numerical simulations, DNSs) for spanwise-periodic turbulent channel flows, and the extreme complexity of constructing a fully developed duct flow facility with aspect ratio of 24 for high Reynolds numbers with adequate probe resolution. Results from this non-traditional literature for the turbulence community are compared to our recent database of DNS of turbulent duct flows with aspect ratios ranging from 1 to 18 at Re_{τ, c} values of 180 and 330, leading to very good agreement between their experimental and our computational results at these low Reynolds numbers. The DNS results also reveal the complexity of a multitude of streamwise vortical structures in addition to the secondary corner flows (which extend up to z ? 5h). These time-dependent and meandering streamwise structures are located at the core of the duct and scale with its half-height. Comparisons of these structures with the vortical motions found in spanwise-periodic channels reveal similitudes in their time-averages and the same rate of decay of their mean kinetic energy ～ T^{? 1}_A, with T_A being the averaging time. However, differences between the two flows are identified and ideas for their future analysis are proposed.