The State of the Art of Hybrid RANS/LES Modeling for the Simulation of Turbulent Flows

This review presents the state of the art of hybrid RANS/LES modeling for the simulation of turbulent flows. After recalling the modeling used in RANS and LES methodologies, we propose in a first step a theoretical formalism developed in the spectral space that allows to unify the RANS and LES methods from a physical standpoint. In a second step, we discuss the principle of the hybrid RANS/LES methods capable of representing a RANS-type behavior in the vicinity of a solid boundary and an LES-type behavior far away from the wall boundary. Then, we analyze the principal hybrid RANS/LES methods usually used to perform numerical simulation of turbulent flows encountered in engineering applications. In particular, we investigate the very large eddy simulation (VLES), the detached eddy simulation (DES), the partially integrated transport modeling (PITM) method, the partially averaged Navier-Stokes (PANS) method, and the scale adaptive simulation (SAS) from a physical point of view. Finally, we establish the connection between these methods and more precisely, the link between PITM and PANS as well as DES and PITM showing that these methods that have been built by different ways, practical or theoretical manners have common points of comparison. It is the opinion of the author to consider that the most appropriate method for a particular application will depend on the expectations of the engineer and the computational resources the user is prepared to expend on the problem.     

Wall modeling via function enrichment within a high‐order DG method for RANS simulations of incompressible flow

We present a novel approach to wall modeling for the Reynolds‐averaged Navier‐Stokes equations within the discontinuous Galerkin method. Wall functions are not used to prescribe boundary conditions as usual, but they are built into the function space of the numerical method as a local enrichment, in addition to the standard polynomial component. The Galerkin method then automatically finds the optimal solution among all shape functions available. This idea is fully consistent and gives the wall model vast flexibility in separated boundary layers or high adverse pressure gradients. The wall model is implemented in a high‐order discontinuous Galerkin solver for incompressible flow complemented by the Spalart‐Allmaras closure model. As benchmark examples, we present turbulent channel flow starting from Re_τ=180 and up to Re_τ=100000 as well as flow past periodic hills at Reynolds numbers based on the hill height of Re_H=10595 and Re_H=19000.     

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     

Generation of inflow data for inhomogeneous turbulence

Inflow boundary conditions for turbulent plane channel flow are generated by solving evolution equations only for the most energetic eddies. The dynamical systems are derived by Galerkin projecting the Navier-Stokes equations onto the subspaces spanned by various sets of the most energetic modes from a proper orthogonal decomposition (POD) of the same flow. Low-energy small-scale POD-modes are added randomly in order to impose some energy in the high wave number range. This is found to be crucial in order to more rapidly establish the correct level of dissipation and achieve a more realistic distribution of energy between the velocity components. The method is tested on a DNS of R_*=180 and a LES of R_*=400. Statistics such as mean velocity, rms-profiles, turbulent shear-stress and energy spectra become close to the fully developed state within 1500 wall units downstream the inlet. PACS 47.27.Eq     

Wall modeling for implicit large-eddy simulation and immersed-interface methods

We propose and analyze a wall model based on the turbulent boundary layer equations (TBLE) for implicit large-eddy simulation (LES) of high Reynolds number wall-bounded flows in conjunction with a conservative immersed-interface method for mapping complex boundaries onto Cartesian meshes. Both implicit subgrid-scale model and immersed-interface treatment of boundaries offer high computational efficiency for complex flow configurations. The wall model operates directly on the Cartesian computational mesh without the need for a dual boundary-conforming mesh. The combination of wall model and implicit LES is investigated in detail for turbulent channel flow at friction Reynolds numbers from Re _τ  = 395 up to Re _τ =100,000 on very coarse meshes. The TBLE wall model with implicit LES gives results of better quality than current explicit LES based on eddy viscosity subgrid-scale models with similar wall models. A straightforward formulation of the wall model performs well at moderately large Reynolds numbers. A logarithmic-layer mismatch, observed only at very large Reynolds numbers, is removed by introducing a new structure-based damping function. The performance of the overall approach is assessed for two generic configurations with flow separation: the backward-facing step at Re _h = 5,000 and the periodic hill at Re _H = 10,595 and Re _H = 37,000 on very coarse meshes. The results confirm the observations made for the channel flow with respect to the good prediction quality and indicate that the combination of implicit LES, immersed-interface method, and TBLE-based wall modeling is a viable approach for simulating complex aerodynamic flows at high Reynolds numbers. They also reflect the limitations of TBLE-based wall models.     

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.     

LES/RANS方法混合函数特性研究

将两方程k-ω SST湍流模型和Sagaut的混合尺度亚格子模型通过一个混合函数相结合, 构造出一种混合大涡/雷诺平均N-S方程模拟方法(hybird large eddy simulation/reynolds-averaged navier-stokes, Hybrid LES/RANS), 采用这种混合模拟方法结合5阶WENO格式对Ma=2.8平板湍流边界层进行了数值模拟, 并在计算区域上游入口处采用“回收/调节”方法生成湍流脉动边界条件, 通过考查RANS区域向LES区域的过渡参数及网格分辨率对这种混合模拟方法进行了评价. 计算结果表明: 该文采用的混合模拟方法可以捕捉到湍流边界层中的大尺度结构且入口边界层平均参数不会发生漂移, 混合函数应当将RANS区域和LES区域的过渡点设置在对数律层和尾迹律层的交界处, 而过渡应当迅速以获得正确的雷诺剪切应力分布, 在该文采用的模型及数值方法的条件下, 流向及展向的网格小至与Escudier混合长相当时, 能够获得可以接受的脉动速度的单点-二阶统计值.     

Hybrid LES-RANS: An approach to make LES applicable at high Reynolds number

The main bottle neck for using large eddy simulations (LES) at high Reynolds number is the requirement of very fine meshes near walls. Hybrid LES-Reynolds-averaged Navier-Stokes (RANS) was invented to get rid of this limitation. In this method, unsteady RANS (URANS) is used near walls and away from walls LES is used. The matching between URANS and LES takes place in the inner log-region. In the present paper, a method to improve standard LES-RANS is evaluated. The improvement consists of adding instantaneous turbulent fluctuations (forcing conditions) at the matching plane in order to provide the equations in the LES region with relevant turbulent structures. The fluctuations are taken from a DNS of a generic boundary layer. Simulations of fully developed channel flow and plane asymmetric diffuser flow are presented. Hybrid LES-RANS is used both with and without forcing conditions.     

DNS,LES and RANS of turbulent heat transfer in boundary layer with suddenly changing wall thermal conditions

The objectives of this study are to investigate a thermal field in a turbulent boundary layer with suddenly changing wall thermal conditions by means of direct numerical simulation (DNS), and to evaluate predictions of a turbulence model in such a thermal field, in which DNS of spatially developing boundary layers with heat transfer can be conducted using the generation of turbulent inflow data as a method. In this study, two types of wall thermal condition are investigated using DNS and predicted by large eddy simulation (LES) and Reynolds-averaged Navier–Stokes equation simulation (RANS). In the first case, the velocity boundary layer only develops in the entrance of simulation, and the flat plate is heated from the halfway point, i.e., the adiabatic wall condition is adopted in the entrance, and the entrance region of thermal field in turbulence is simulated. Then, the thermal boundary layer develops along a constant temperature wall followed by adiabatic wall. In the second case, velocity and thermal boundary layers simultaneously develop, and the wall thermal condition is changed from a constant temperature to an adiabatic wall in the downstream region. DNS results clearly show the statistics and structure of turbulent heat transfer in a constant temperature wall followed by an adiabatic wall. In the first case, the entrance region of thermal field in turbulence can be also observed. Thus, both the development and the entrance regions in thermal fields can be explored, and the effects upstream of the thermal field on the adiabatic region are investigated. On the other hand, evaluations of predictions by LES and RANS are conducted using DNS results. The predictions of both LES and RANS almost agree with the DNS results in both cases, but the predicted temperature variances near the wall by RANS give different results as compared with DNS. This is because the dissipation rate of temperature variance is difficult to predict by the present RANS, which is found by the evaluation using DNS results.     

Numerical prediction of turbulent boundary layer noise from a sharp-edged flat plate

An efficient hybrid uncorrelated wall plane waves–boundary element method (UWPW-BEM) technique is proposed to predict the flow-induced noise from a structure in low Mach number turbulent flow. Reynolds-averaged Navier-Stokes equations are used to estimate the turbulent boundary layer parameters such as convective velocity, boundary layer thickness, and wall shear stress over the surface of the structure. The spectrum of the wall pressure fluctuations is evaluated from the turbulent boundary layer parameters and by using semi-empirical models from literature. The wall pressure field underneath the turbulent boundary layer is synthesized by realizations of uncorrelated wall plane waves (UWPW). An acoustic BEM solver is then employed to compute the acoustic pressure scattered by the structure from the synthesized wall pressure field. Finally, the acoustic response of the structure in turbulent flow is obtained as an ensemble average of the acoustic pressures due to all realizations of uncorrelated plane waves. To demonstrate the hybrid UWPW-BEM approach, the self-noise generated by a flat plate in turbulent flow with Reynolds number based on chord Re_c = 4.9 × 10⁵ is predicted. The results are compared with those obtained from a large eddy simulation (LES)-BEM technique as well as with experimental data from literature.     

A Rapid and Accurate Switch from RANS to LES in Boundary Layers Using an Overlap Region

An efficient recycling algorithm is developed for injecting resolved turbulent content in a boundary layer as it switches from a Reynolds Averaged Navier-Stokes (RANS) type treatment to a Large Eddy Simulation (LES) type treatment inside a generalized Detached-Eddy Simulation (DES). The motivation is to use RANS in the thinnest boundary-layer area, following the original argument in favour of DES, and LES in the thicker boundary-layer areas especially approaching separation, to improve accuracy and possibly obtain unsteady outputs. The algorithm relies on an overlap of the RANS and LES domains and, therefore, the availability of both RANS and LES solutions in the recycling region, which is about 5 boundary-layer thicknesses long. This permits a smooth transfer of the turbulent stresses from this section to the LES inflow. The continuity of the skin-friction distribution is very good, reflecting the excellent viability of the resolved turbulence. The approach is validated in a flat-plate boundary layer and an airfoil near stall, with mild pressure gradient near the interface, and then applied to the compressible flow over an idealized airliner windshield wiper. The pressure fluctuations at reattachment are 12dB more intense than under a simple boundary layer at the same speed, and the output contains all the quantities needed to calculate the transmission of sound through the glass.     

A stress transport equation model for simulating turbulence at any mesh resolution

Prior work has demonstrated the effectiveness of using two-equation closures as the basis for universal, self-adapting turbulence models that are effective at any mesh resolution (Perot and Gadebusch in Phys. Fluids 19:115105, 2007). In order to demonstrate the broad applicability of the fundamental approach, the same behavior is now demonstrated for a second-moment closure (SMC). The SMC has the advantage over the earlier two-equation universal closure of being more accurate in the coarse mesh limit and of having a natural mechanism for backscattering energy from the modeled to the resolved turbulent fluctuations. The mathematical explanation for why Reynolds averaged (RANS) transport equation closures are applicable at any mesh resolution, including the large eddy simulation (LES) regime, is reviewed. It is demonstrated that for the problem of isotropic decaying turbulence, the SMC model produces good predictions at any mesh resolution and with arbitrary initial conditions. In addition, it is shown that the proposed model automatically adapts to the mesh resolution provided. The self-adaptive nature of the method is clearly observed when different initial conditions are used. It is shown that classic RANS models (often thought to produce steady and smooth solutions) can produce three-dimensional, unsteady, and chaotic solutions when generalized correctly and when provided with sufficient mesh resolution. The implications of these observations on the fundamental theories of RANS and LES turbulence modeling are discussed.      

An immersed boundary/wall modeling method for RANS simulation of compressible turbulent flows

In this paper, an immersed boundary (IB) method is developed to simulate compressible turbulent flows governed by the Reynolds‐averaged Navier‐Stokes equations. The flow variables at the IB nodes (interior nodes in the immediate vicinity of the solid wall) are evaluated via linear interpolation in the normal direction to close the discrete form of the governing equations. An adaptive wall function and a 2‐layer wall model are introduced to reduce the near‐wall mesh density required by the high resolution of the turbulent boundary layers. The wall shear stress modified by the wall modeling technique and the no‐penetration condition are enforced to evaluate the velocity at an IB node. The pressure and temperature at an IB node are obtained via the local simplified momentum equation and the Crocco‐Busemann relation, respectively. The SST k ? ω and S‐A turbulence models are adopted in the framework of the present IB approach. For the Shear‐Stress Transport (SST) k ? ω model, analytical solutions in near‐wall region are utilized to enforce the boundary conditions of the turbulence equations and evaluate the turbulence variables at an IB node. For the S‐A model, the turbulence variable at an IB node is calculated by using the near‐wall profile of the eddy viscosity. In order to validate the present IB approach, numerical experiments for compressible turbulent flows over stationary and moving bodies have been performed. The predictions show good agreements with the referenced experimental data and numerical results.     

Towards a new partially integrated transport model for coarse grid and unsteady turbulent flow simulations

A new two-equation model is proposed for large eddy simulations (LESs) using coarse grids. The modeled transport equations are obtained from a direct transposition of well-known statistical models by using multiscale spectrum splitting given by the filtering operation applied to the Navier–Stokes equations. The model formulation is compatible with the two extreme limits that are on one hand a direct numerical simulation and on the other hand a full statistical modeling. The characteristic length scale of subgrid turbulence is no longer given by the spatial discretization step size, but by the use of a dissipation equation. The proposed method is applied to a transposition of the well-known k- statistical model, but the same method can be developed for more advanced closures. This approach is intended to contribute to non-zonal hybrid models that bridge Reynolds-averaged Navier–Stokes (RANS) and LES, by using a continuous change rather than matching zones. The main novelty in the model is the derivation of a new equation for LES that is formally consistent with RANS when the filter width is very large. This approach is dedicated to applications to non-equilibrium turbulence and coarse grid simulations. An illustration is made of large eddy simulations of turbulence submitted to periodic forcing. The model is also an alternative approach to hybrid models. PACS 47.27.Eq     

Joint RANS/LES Approach to Premixed Flame Modelling in the Context of the TFC Combustion Model

We present an original timesaving joint RANS/LES approach to simulate turbulent premixed combustion. It is intended mainly for industrial applications where LES may not be practical. It is based on successive RANS/LES numerical modelling, where turbulent characteristics determined from RANS simulations are used in LES equations for estimation of the subgrid chemical source and viscosity. This approach has been developed using our TFC premixed combustion model, which is based on a generalization of the Kolmogorov’s ideas. We assume existence of small-scale statistically equilibrium structures not only of turbulence but also of the reaction zones. At the same time, non-equilibrium large-scale structures of reaction sheets and turbulent eddies are described statistically by model combustion and turbulence equations in RANS simulations or follow directly without modelling in LES. Assumption of small-scale equilibrium gives an opportunity to express the mean combustion rate (controlled by small-scale coupling of turbulence and chemistry) in the RANS and LES sub-problems in terms of integral or subgrid parameters of turbulence and the chemical time, i.e. the definition of the reaction rate is similar to that of the mean dissipation rate in turbulence models where it is expressed in terms of integral or subgrid turbulent parameters. Our approach therefore renders compatible the combustion and turbulent parts of the RANS and LES sub-problems and yields reasonable agreement between the RANS and averaged LES results. Combining RANS simulations of averaged fields with LES method (and especially coupled and acoustic codes) for simulation of corresponding nonstationary process (and unsteady combustion regimes) is a promising strategy for industrial applications. In this work we present results of simulations carried out employing the joint RANS/LES approach for three examples: High velocity premixed combustion in a channel, combustion in the shear flow behind an obstacle and the impinging flame (a premixed flame attached to an obstacle).     

Numerical study of an oscillatory turbulent flow over a flat plate

Oscillatory turbulent flow over a flat plate is studied using large eddy simulation (LES) and Reynolds-average Navier-Stokes (RANS) methods. A dynamic subgrid-scale model is employed in LES and Saffman's turbulence model is used in RANS. The flow behaviors are discussed for the accelerating and decelerating phases during the oscillating cycle. The friction force on the wall and its phase shift from laminar to turbulent regime are also investigated for different Reynolds numbers. The project supported by the Youngster Funding of Academia Sinica and by the National Natural Science Foundation of China     

Towards a Unified Turbulence Simulation Approach for Wall-Bounded Flows

A hybrid Reynolds-averaged Navier–Stokes/Large-Eddy Simulation (RANS/LES) methodology has received considerable attention in recent years, especially in its application to wall-bounded flows at high-Reynolds numbers. In the conventional zonal hybrid approach, eddy-viscosity-type RANS and subgrid scale models are applied in the RANS and LES zones, respectively. In contrast, the non-zonal hybrid approach uses only a generalized turbulence model, which provides a unified simulation approach that spans the continuous spectrum of modeling/simulation schemes from RANS to LES. A particular realization of the non-zonal approach, known as partially resolved numerical simulation (PRNS), uses a generalized turbulence model obtained from a rescaling of a conventional RANS model through the introduction of a resolution control function F _R , where F _R is used to characterize the degree of modeling required to represent the unresolved scales of turbulent motion. A new generalized functional form for F _R in PRNS is proposed in this study, and its performance is compared with unsteady RANS (URANS) and LES computations for attached and separated wall-bounded turbulent flows. It is demonstrated that PRNS behaves similarly to LES, but outperforms URANS in general.     

Large eddy simulation of turbulence in engineering applications

Flows of engineering interest which have been investigated with the large eddy simulation (LES) technique are discussed. Using this technique the filtered Navier-Stokes equations are solved to simulate the three-dimensional unsteady motion. Not only the statistically-averaged fields, but in addition, the spatial and temporal variations in the flow fields can be examined. Grids are equidistant in the three coordinate directions and typically consist of 150,000 to 200,000 mesh cells. Examples of flows include the turbulent boundary layer with a sudden change in the freestream pressure gradient and the separating and reattaching flow behind a rearward-facing step. Results which illustrate the large-scale structure of these two high Reynolds number turbulent flows are discussed in the present paper.     

Passive Scalar Transport Modeling for Hybrid RANS/LES Simulation

A transport model for hybrid RANS/LES simulation of passive scalars is proposed. It invokes a dynamically computed subgrid Prandtl number. The method is based on computing test-filter fluxes. The formulation proves to be especially effective on coarse grids, as occur in DES. After testing it in a wall resolved LES, the present formulation is applied to the Adaptive DDES model of Yin et al. (Phys. Fluids 27, 025105 2015). It is validated by turbulent channel flow and turbulent boundary layer computations.