Approximate Wall Boundary Conditions in the Large-Eddy Simulation of High Reynolds Number Flow

The near-wall regions of high Reynolds numbers turbulent flows must be modelled to treat many practical engineering and aeronautical applications. In this review we examine results from simulations of both attached and separated flows on coarse grids in which the near-wall regions are not resolved and are instead represented by approximate wall boundary conditions. The simulations use the dynamic Smagorinsky subgrid-scale model and a second-order finite-difference method. Typical results are found to be mixed, with acceptable results found in many cases in the core of the flow far from the walls, provided there is adequate numerical resolution, but with poorer results generally found near the wall. Deficiencies in this approach are caused in part by both inaccuracies in subgrid-scale modelling and numerical errors in the low-order finite-difference method on coarse near-wall grids, which should be taken into account when constructing models and performing large-eddy simulation on coarse grids. A promising new method for developing wall models from optimal control theory is also discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

A ghost-cell immersed boundary method for large eddy simulation of flows in complex geometries

An efficient ghost-cell immersed boundary (IB) method is proposed for large eddy simulations of three-dimensional incompressible flow in complex geometries. In the framework of finite volume method, the Navier–Stokes equations are integrated using an explicit time advancement scheme on a collocated mesh. Since the IB method is known to generate an unphysical velocity field inside the IB that violates the mass conservation of the cells near the IB, a new IB treatment is devised to eliminate the unphysical velocity generated near the IB and to improve the pressure distribution on the body surface. To validate the proposed method, both laminar and turbulent flow cases are presented. In particular, large eddy simulations were performed to simulate the turbulent flows over a circular cylinder and a sphere at subcritical Reynolds numbers. The computed results show good agreements with the published numerical and experimental data.     

Application of an Explicit Algebraic Reynolds Stress Model within a Hybrid LES–RANS Method

Large-eddy simulations (LES) still suffer from extremely large resources required for the resolution of the near-wall region, especially for high-Re flows. That is the main motivation for setting up hybrid LES–RANS methods. Meanwhile a variety of different hybrid concepts were proposed mostly relying on linear eddy-viscosity models. In the present study a hybrid approach based on an explicit algebraic Reynolds stress model (EARSM) is suggested. The model is applied in the RANS mode with the aim of accounting for the Reynolds stress anisotropy emerging especially in the near-wall region. For the implementation into a CFD code this anisotropy-resolving closure can be formally expressed in terms of a non-linear eddy-viscosity model (NLEVM). Its extra computational effort is small, still requiring solely the solution of one additional transport equation for the turbulent kinetic energy. In addition to this EARSM approach, a linear eddy-viscosity model (LEVM) is used in order to verify and emphasize the advantages of the non-linear model. In the present formulation the predefinition of RANS and LES regions is avoided and a gradual transition between both methods is assured. A dynamic interface criterion is suggested which relies on the modeled turbulent kinetic energy and the wall distance and thus automatically accounts for the characteristic properties of the flow. Furthermore, an enhanced version guaranteeing a sharp interface is proposed. The interface behavior is thoroughly investigated and it is shown how the method reacts on dynamic variations of the flow field. Both model variants, i.e. LEVM and EARSM, have been tested on the basis of the standard plane channel flow and even more detailed on the flow over a periodic arrangement of hills using fine and coarse grids.     

大涡模拟的壁模型及其应用

大涡模拟是研究湍流的非定常特性的重要方法. 但解析壁面层的大涡模拟所需的计算量与直接数值模拟相当,是大涡模拟在高雷诺数壁湍流数值模拟中所面临的主要困难. 解析壁面层所需的网格尺度与壁面黏性长度同量级,是引起壁湍流大涡模拟计算量增加的主要原因. 壁模型通过模化近壁流动避免了完全解析壁面层,可以显著地降低壁湍流大涡模拟的计算量,是克服上述困难的有效方法. 本文介绍了大涡模拟壁模型的主要类型;详细讨论了常用的壁面应力模型,特别是平衡层模型和双层模型的构建思路和特点;基于近壁流动的特征讨论了应力边界条件的必要性和适用性;指出了壁面应力模型的局限性以及考虑非平衡效应修正的各种方法;讨论了壁面应力模型的研究历史、最新进展和发展趋势,给出了常用的壁面应力模型的分支与发展关系图;并基于Werner-Wengle模型实现了周期山状流的大涡模拟.      

LES-based filter-matrix lattice Boltzmann model for simulating fully developed turbulent channel flow

In this paper, a three-dimensional filter-matrix lattice Boltzmann (FMLB) model based on large eddy simulation (LES) was verified for simulating wall-bounded turbulent flows. The Vreman subgrid-scale model was employed in the present FMLB–LES framework, which had been proved to be capable of predicting turbulent near-wall region accurately. The fully developed turbulent channel flows were performed at a friction Reynolds number Re_τ of 180. The turbulence statistics computed from the present FMLB–LES simulations, including mean stream velocity profile, Reynolds stress profile and root-mean-square velocity fluctuations greed well with the LES results of multiple-relaxation-time (MRT) LB model, and some discrepancies in comparison with those direct numerical simulation (DNS) data of Kim et al. was also observed due to the relatively low grid resolution. Moreover, to investigate the influence of grid resolution on the present LES simulation, a DNS simulation on a finer gird was also implemented by present FMLB–D3Q19 model. Comparisons of detailed computed various turbulence statistics with available benchmark data of DNS showed quite well agreement.     

An advanced switching parameter for a hybrid LES/RANS model considering the characteristics of near-wall turbulent length scales

In this study, we proposed an idea for an advanced switching parameter used in a hybrid approach connecting large eddy simulation (LES) with Reynolds-averaged Navier–Stokes modeling [the hybrid LES/RANS (HLR) model]. Although the HLR model is promising way to predict engineering turbulent flows, an important problem is that RANS is always adopted in the near-wall region, even if the grid resolution is fine enough for LES. To overcome this difficulty, the switching parameter proposed here introduced knowledge of the Kolmogorov microscale that is thought to be reasonable for representing the near-wall turbulence. This parameter enabled the present HLR model to be smoothly replaced by a full LES if a grid resolution was fine enough in the near-wall region. To confirm model performance, the present HLR model was applied to numerical simulations of a periodic hill flow as well as fundamental plane channel flows. The model generally provided reasonable predictions for these test cases that include complex turbulence with massive flow separation.     

LES and RANS for Turbulent Flow over Arrays of Wall-Mounted Obstacles

Large-eddy simulation (LES) has been applied to calculate the turbulent flow over staggered wall-mounted cubes and staggered random arrays of obstacles with area density 25%, at Reynolds numbers between 5 × 10³ and 5 10⁶, based on the free stream velocity and the obstacle height. Re = 5 × 10³ data were intensively validated against direct numerical simulation (DNS) results at the same Re and experimental data obtained in a boundary layer developing over an identical roughness and at a rather higher Re. The results collectively confirm that Reynolds number dependency is very weak, principally because the surface drag is predominantly form drag and the turbulence production process is at scales comparable to the roughness element sizes. LES is thus able to simulate turbulent flow over the urban-like obstacles at high Re with grids that would be far too coarse for adequate computation of corresponding smooth-wall flows. Comparison between LES and steady Reynolds-averaged Navier-Stokes (RANS) results are included, emphasising that the latter are inadequate, especially within the canopy region.     

Compound wall treatment with low‐Re turbulence model

A generalized treatment for the wall boundary conditions relating to turbulent flows is developed that blends the integration to a solid wall with wall functions. The blending function ensures a smooth transition between the viscous and turbulent regions. An improved low Reynolds number k?ε model is coupled with the proposed compound wall treatment to determine the turbulence field. The eddy viscosity formulation maintains the positivity of normal Reynolds stresses and Schwarz' inequality for turbulent shear stresses. The model coefficients/functions preserve the anisotropic characteristics of turbulence. Computations with fine and coarse meshes of a few flow cases yield appreciably good agreement with the direct numerical simulation and experimental data. The method is recommended for computing the complex flows where computational grids cannot satisfy a priori the prerequisites of viscous/turbulence regions. Copyright © 2011 John Wiley & Sons, Ltd.     

Extension of the local domain-free discretization method to large eddy simulation of turbulent flows

In this work, an immersed boundary method, called the local domain-free discretization (DFD) method, is extended to large eddy simulation (LES) of turbulent flows. The discrete form of partial differential equations at an interior node may involve some nodes outside the solution domain. The flow variables at these exterior dependent nodes are evaluated via linear extrapolation along the direction normal to the wall. To alleviate the requirement of mesh resolution in the near-wall region, a wall model based on the turbulence boundary layer equations is introduced. The wall shear stress yielded by the wall model and the no-penetration condition are enforced at the immersed boundary to evaluate the velocity components at an exterior dependent node. For turbulence closure, a dynamic subgrid scale (SGS) model is adopted and the Lagrangian averaging procedure is used to compute the model coefficient. The SGS eddy viscosity at an exterior dependent node is set to be equal to that at the outer layer. To maintain the mass conservation near the immersed boundary, a mass source/sink term is added into the continuity equation. Numerical experiments on relatively coarse meshes with stationary or moving solid boundaries have been conducted to verify the ability of the present LES-DFD method. The predicted results agree well with the published experimental or numerical data.     

Lagrangian particle dispersion in turbulent flow over a wall mounted obstacle

Large-eddy simulations (LES) of particle-laden turbulent flows are presented in order to investigate the effects of particle response time on the dispersion patterns of a space developing flow with an obstruction, where solid particles are injected inside the wake of an obstacle [Vincont, J.Y., Simoens, S., Ayrault M., Wallace, J.M., 2000. Passive scalar dispersion in a turbulent boundary layer from a line source at the wall and downstream of an obstacle. J. Fluid Mech. 424, 127–167]. The numerical method is based on a fully explicit fractional step approach and finite-differences on Cartesian grids, using the immersed boundary method (IBM) to represent the existence of solid obstacles. Two different turbulence models have been tested, the classical Smagorinsky turbulence model and the filtered structure function model. The dispersed phase was modelled either by an Eulerian approach or a Lagrangian particle tracking scheme of solid particles with Stokes numbers in the range St = 0–25, assuming one-way coupling between the two phases. A very good agreement was observed between the Lagrangian and Eulerian approaches. The effect of particle size was found to significantly differentiate the dispersion pattern for the inhomogeneous flow over the obstacle. Although in homogeneous flows like particle-laden turbulent channels near-wall particle clustering increases monotonically with particle size, for the examined flow over an obstacle, preferential concentration effects were stronger only for an intermediate range of Stokes numbers.     

Spectral-Dynamic Model for Large-Eddy Simulations of Turbulent Rotating Channel Flow

A new subgrid-scale model called the spectral-dynamic model is proposed. It consists of a refinement of spectral eddy-viscosity models taking into account nondeveloped turbulence in the subgrid-scales. The proposed correction, which is derived from eddy-damped quasi-normal Markovian statistical theory, is based on an adjustment of the turbulent eddy-viscosity coefficient to the deviation of the spectral slope (at small scales) with respect to the standard Kolmogorov law. The spectral-dynamic model is applied to large eddy simulation (LES) of rotating and nonrotating turbulent plane channel flows. It is shown that the proposed refinement allows for clear improvement of the statistical predictions due to a correct prediction of the near-wall behavior. Cases of rotating and nonrotating low (DNS) and high Reynolds (LES) numbers are then compared. It is shown that the principal structural features of the rotating turbulent channel flow are reproduced by the LES, such as the presence of the near-zero mean absolute vorticity region, the modification of the anisotropic character of the flow (with respect to the nonrotating case), the enhancement of flow organization, and the inhibition of the high- and low-speed streaks near the anticyclonic wall. Only a moderate Reynolds number dependence is exhibited, resulting in a more unstable character of the longitudinal large-scale roll cells at high Reynolds number, and a slight increase of the laminarization tendency on the cyclonic side of the channel. Received 16 October 1997 and accepted 1 October 1998     

An Attempt to Combine Large Eddy Simulation with the k-ε Model in a Channel-Flow Calculation

Large eddy simulation (LES) is combined with the Reynolds-averaged Navier–Stokes (RANS) equation in a turbulent channel-flow calculation. A one-equation subgrid-scale model is solved in a three-dimensional grid in the near-wall region whereas the standard k–ε model is solved in a one-dimensional grid in the outer region away from the wall. The two grid systems are overlapped to connect the two models smoothly. A turbulent channel flow is calculated at Reynolds numbers higher than typical LES and several statistical quantities are examined. The mean velocity profile is in good agreement with the logarithmic law. The profile of the turbulent kinetic energy in the near-wall region is smoothly connected with that of the turbulent energy for the k–ε model in the outer region. Turbulence statistics show that the solution in the near-wall region is as accurate as a usual LES. The present approach is different from wall modeling in LES that uses a RANS model near the wall. The former is not as efficient as the latter for calculating high-Reynolds-number flows. Nevertheless, the present method of combining the two models is expected to pave the way for constructing a unified turbulence model that is useful for many purposes including wall modeling. Received 11 June 1999 and accepted 15 December 2000     

A low Reynolds number variant of partially-averaged Navier-Stokes model for turbulence

A low Reynolds number (LRN) formulation based on the Partially Averaged Navier-Stokes (PANS) modelling method is presented, which incorporates improved asymptotic representation in near-wall turbulence modelling. The effect of near-wall viscous damping can thus be better accounted for in simulations of wall-bounded turbulent flows. The proposed LRN PANS model uses an LRN k-ε model as the base model and introduces directly its model functions into the PANS formulation. As a result, the inappropriate wall-limiting behavior inherent in the original PANS model is corrected. An interesting feature of the PANS model is that the turbulent Prandtl numbers in the k and ε equations are modified compared to the base model. It is found that this modification has a significant effect on the modelled turbulence. The proposed LRN PANS model is scrutinized in computations of decaying grid turbulence, turbulent channel flow and periodic hill flow, of which the latter has been computed at two different Reynolds numbers of Re = 10,600 and 37,000. In comparison with available DNS, LES or experimental data, the LRN PANS model produces improved predictions over the standard PANS model, particularly in the near-wall region and for resolved turbulence statistics. Furthermore, the LRN PANS model gives similar or better results - at a reduced CPU time - as compared to the Dynamic Smagorinsky model.     

Large Eddy Simulations of convergent–divergent channel flows at moderate Reynolds numbers

The paper is focused on the study of fully turbulent channel flows, using Large Eddy Simulations (LES), in order to address the effects of adverse pressure gradient regions. Analyses of the effects of streak instabilities, which have been shown to be relevant in such regions, are extended to moderate Reynolds numbers. The work considers two different channel geometries in order to further separate influences from wall curvature, flow separation and adverse pressure gradients. Turbulent kinetic energy and Reynolds stress budgets are investigated at separation and re-attachment points. The numerical approach used in the present work is based on the incompressible Navier–Stokes equations, which are solved by a pseudo-spectral methodology for structured grids. Wall-resolved LES calculations are performed using the WALE subgrid scale model. The study shows that the streak instability mechanism persists at higher Reynolds numbers with and without wall curvature in the adverse pressure gradient regions. Moreover, the observed effects are also present regardless of the existence of flow separation regions. Finally, the study of turbulent kinetic energy budgets indicates that, independently of the flow condition, there are well-defined patterns for such turbulent properties at separation and re-attachment points.     

Development of a Two-velocities Hybrid RANS-LES Model and its Application to a Trailing Edge Flow

The flow around a trailing edge is computed with a new hybrid method designed to more clearly separate the effects of total and sub-grid turbulent stress-modelling on the time-averaged and instantaneous velocity fields, and in turn, mean momentum and kinetic energy balances. These two velocity fields independently define Reynolds averaged and sub-grid-scale viscosities, and distinct stresses, at the same location. In particular, resolved eddies can emerge, or sweep in and out of the Reynolds averaged near wall layer, without being dampened by higher levels of the viscosity in this RANS dominated layer. The two-field hybrid model, first tested on channel flows, gives accurate predictions of mean velocities and stresses for different Reynolds numbers and coarse meshes. For the trailing edge flow the results of the hybrid model are close to the reference fine LES for mean velocity and turbulent content, whereas the DES-SST on the same coarse mesh gives too early separation.     

基于全局动态Vreman模型的复杂非定常湍流的大涡模拟

在湍流数值模拟方法中,大涡模拟方法可以提供丰富的大涡旋信息,已逐渐成为复杂湍流问题数值研究的重要方法。而大涡模拟中,最重要的一环是尽量准确地构建能反映流场物理本质特征的亚格子应力模型。基于该思想,将一种新型的大涡模拟亚格子应力模型-Vreman亚格子应力模型用于高雷诺数三维后台阶流动的求解,计算结果与实验结果进行对比分析结果较吻合,验证了该模型的可靠性。这是对该模型用于无任何均匀流动方向的高雷诺数复杂湍流非定常流动的首次检验,计算结果优于基于传统的Smagorinsky涡粘性的动态亚格子模型。     

基于滑移速度壁模型的复杂边界湍流大涡模拟

采用滑移速度壁模型实现了浸入边界方法与壁模型相结合的大涡模拟.本文首先分别采用平衡层模型和非平衡壁模型对周期山状流进行数值模拟,以考查在壁模型中考虑切向压力梯度的作用.数值结果表明,流场的压力对本文所采用的壁模型形式并不敏感,但是考虑切向压力梯度可以显著改进壁面摩擦力的计算结果,并且能够准确的预测强压力梯度区以及分离区内的流动平均统计特性.不考虑压力梯度效应的平衡层模型显著低估了壁面摩擦力的分布,同时无法准确预测分离区内的平均速度剖面.非平衡模型的修正项正比于切向压力梯度和壁面法向距离,因此在强压力梯度区或者网格较粗时,计算得到的平均压力和摩擦力分布以及流动的低阶统计量均与参考的实验和计算结果吻合.在此基础上,通过回转体绕流的大涡模拟考查了该方法用于模拟高雷诺数壁湍流的适用性,非平衡壁模型可以准确地捕捉流动的物理结构并较准确地预测其水动力学特性.结果表明,将浸入边界方法与非平衡滑移速度壁模型相结合的大涡模拟,有望成为数值模拟复杂边界高雷诺数壁湍流的工具.      

A multiscale approach to hybrid RANS/LES wall modeling within a high-order discontinuous Galerkin scheme using function enrichment

We present a novel approach to hybrid Reynolds-averaged Navier-Stokes (RANS)/ large eddy simulation (LES) wall modeling based on function enrichment, which overcomes the common problem of the RANS-LES transition and enables coarse meshes near the boundary. While the concept of function enrichment as an efficient discretization technique for turbulent boundary layers has been proposed in an earlier article by Krank & Wall (A new approach to wall modeling in LES of incompressible flow via function enrichment. J Comput Phys. 2016;316:94-116), the contribution of this work is a rigorous derivation of a new multiscale turbulence modeling approach and a corresponding discontinuous Galerkin discretization scheme. In the near-wall area, the Navier-Stokes equations are explicitly solved for an LES and a RANS component in one single equation. This is done by providing the Galerkin method with an independent set of shape functions for each of these two methods; the standard high-order polynomial basis resolves turbulent eddies, where the mesh is sufficiently fine and the enrichment automatically computes the ensemble-averaged flow if the LES mesh is too coarse. As a result of the derivation, the RANS model is applied solely to the RANS degrees of freedom, which effectively prevents the typical issue of a log-layer mismatch in attached boundary layers. As the full Navier-Stokes equations are solved in the boundary layer, spatial refinement gradually yields wall-resolved LES with exact boundary conditions. Numerical tests show the outstanding characteristics of the wall model regarding grid independence, superiority compared to equilibrium wall models in separated flows, and achieve a speed-up by two orders of magnitude compared to wall-resolved LES.     

Mesh Node Distribution in Terms of Wall Distance for Large-eddy Simulation of Wall-bounded Flows

In this note, basic turbulent statistics in a pipe flow are computed accurately by large-eddy simulation using a mesh resolution coarser than the viscous sublayer. These results are obtained when a regular Cartesian mesh is used for the spatial discretization of the circular pipe thanks to an immersed boundary method combined with high-order schemes. In this particular computational configuration, the near-wall features of mean velocity and Reynolds stress profiles are found to be correctly captured at a scale significantly smaller than the mesh size. Comparisons between channel and pipe flow configurations suggest that an irregular mesh distribution in terms of wall distance may be a favourable condition to explicitly compute by large-eddy simulation reliable wall turbulence without any extra-modelling in the near-wall region.     

Validation of a novel very large eddy simulation method for simulation of turbulent separated flow

The paper describes the validation of a newly developed very LES (VLES) method for the simulation of turbulent separated flow. The new VLES method is a unified simulation approach that can change seamlessly from Reynolds‐averaged Navier–Stokes to DNS depending on the numerical resolution. Four complex test cases are selected to validate the performance of the new method, that is, the flow past a square cylinder at Re = 3000 confined in a channel (with a blockage ratio of 20%), the turbulent flow over a circular cylinder at Re = 3900 as well as Re = 140,000, and a turbulent backward‐facing step flow with a thick incoming boundary layer at Re = 40,000. The simulation results are compared with available experimental, LES, and detached eddy simulation‐type results. The new VLES model performs well overall, and the predictions are satisfactory compared with previous experimental and numerical results. It is observed that the new VLES method is quite efficient for the turbulent flow simulations; that is, good predictions can be obtained using a quite coarse mesh compared with the previous LES method. Discussions of the implementation of the present VLES modeling are also conducted on the basis of the simulations of turbulent channel flow up to high Reynolds number of Re_τ = 4000. The efficiency of the present VLES modeling is also observed in the channel flow simulation. From a practical point of view, this new method has considerable potential for more complex turbulent flow simulations at relative high Reynolds numbers. Copyright © 2013 John Wiley & Sons, Ltd.