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).     

Computational analysis of impinging shock-wave boundary layer interaction under conditions of incipient separation

The interaction of an oblique shock wave with a turbulent boundary layer under conditions of incipient separation is analyzed by means of large-eddy simulation (LES) and Reynolds-averaged Navier–Stokes (RANS) turbulence models, with the objective to explore their predictive capabilities, in particular with respect to the unsteady features of the interaction. Consistent with earlier direct numerical simulations, we have found that the flow dynamics in the interaction zone is characterized by strong intermittency associated with the formation of scattered spots of flow reversal near the nominal position of the reflected shock. Comparison with experimental results (at much larger Reynolds number) show that the qualitative features of the interaction are predicted reasonably well by both LES and RANS models. RANS models supplemented with a semi-empirical closure are also found to provide reasonable estimate of the fluctuating pressure loads at the wall.     

Non-zonal detached eddy simulation coupled with a steady RANS solver in the wall region

Xiao and Jenny (2012) proposed an interesting hybrid LES/RANS method in which they use two solvers and solve the RANS and LES equations in the entire computational domain. In the present work this method is simplified and used as a hybrid RANS-LES method, a wall-modeled LES. The two solvers are employed in the entire domain. Near the walls, the flow is governed by the steady RANS solver; drift terms are added to the DES equations to ensure that the time-averaged DES fields agree with the steady RANS field. Away from the walls, the flow is governed by the DES solver; in this region, the RANS field is set to the time-averaged LES field. The disadvantage of traditional DES models is that the RANS models in the near-wall region – which originally were developed and tuned for steady RANS – are used as URANS models where a large part of the turbulence is resolved. In the present method – where steady RANS is used in the near-wall region – the RANS turbulence models are used in a context for which they were developed. In standard DES methods, the near-wall accuracy can be degraded by the unsteady agitation coming from the LES region. It may in the present method be worth while to use an accurate, advanced RANS model. The EARSM model is used in the steady RANS solver. The new method is called NZ S-DES . It is found to substantially improve the predicting capability of the standard DES. A great advantage of the new model is that it is insensitive to the location of the RANS-LES interface.     

Hybrid RANS/LES of flow and heat transfer in round impinging jets

Fluid flow and convective heat transfer predictions are presented of round impinging jets for several combinations of nozzle-plate distances H/D = 2, 6 and 13.5 (where D is the nozzle diameter) and Reynolds numbers Re = 5000, 23,000 and 70,000 with the newest version of the k-ω model of Wilcox (2008) and three hybrid RANS/LES models. In the RANS mode of the hybrid RANS/LES models, the k-ω model is recovered. Three formulations are considered to activate the LES mode. The first model is similar to the hybrid models of Davidson and Peng (2003) and Kok et al. (2004). The turbulent length scale is replaced by the grid size in the destruction term of the k-equation and in the definition of the RANS eddy viscosity. As grid size, a maximum measure of the hexahedral grid cell is used. The second model has the same k-equation, but the eddy viscosity is the minimum of the k-ω eddy viscosity and the Smagorinsky eddy viscosity, following a proposal by Batten et al. (2004). The Smagorinsky eddy viscosity is formed with the cube root of the cell volume. The third model has, again, the same k-equation, but has an eddy viscosity which is an intermediate between the eddy viscosities of the first and second models. This is reached by using the cube root of the cell volume in the eddy viscosity formula of the first model.The simulation results are compared with experimental data for the high Reynolds number cases Re = 23,000 and Re = 70,000 and LES data for the low-Reynolds number case Re = 5000. The Reynolds numbers are defined with the nozzle diameter and the bulk velocity at nozzle outlet. At low nozzle-plate distance (the impingement plate is in the core of the jet), turbulent kinetic energy is overpredicted by RANS in the stagnation flow region. This leads to overprediction of the heat transfer rate along the impingement plate in the impact zone. At high nozzle-plate distance (the impingement plate is in the mixed-out region of the jet), the turbulence mixing is underpredicted by RANS in the shear layer of the jet which gives a too high length of the jet core. This also results in overprediction of the heat transfer rate in the impingement zone caused by too big temperature gradients at impingement.All hybrid RANS/LES models are able to correct the heat transfer overprediction of the RANS model. For good predictions at low nozzle-plate distance, it is necessary to sufficiently resolve the formation and development of the near-wall vortices in the jet impingement region. At high nozzle-plate distance, the essence is to capture the evolution and breakup of the flow unsteadiness in the shear layer of the jet, so that accurate mean and fluctuating velocity profiles are obtained in the impingement region. Although the models have a quite different theoretical justification and generate a quite different eddy viscosity in some flow regions, their overall results are very comparable. The reason is that in zones that are crucial for the results, the models behave similarly.     

Large Eddy Simulation of Scalar Mixing in a Coaxial Confined Jet

The highly turbulent flow occurring inside gas-turbine combustors requires accurate simulation of scalar mixing if CFD methods are to be used with confidence in design. This has motivated the present paper, which describes the implementation of a passive scalar transport equation into an LES code, including assessment/testing of alternative discretisation schemes to avoid over/undershoots and excessive smoothing. Both second order accurate TVD and higher order accurate DRP schemes are assessed. The best performance is displayed by a DRP method, but this is only true on fine meshes; it produces similar (or larger) errors to a TVD scheme on coarser meshes, and the TVD approach has been retained for LES applications. The unsteady scalar mixing performance of the LES code is validated against published DNS data for a slightly heated channel flow. Excellent agreement between the current LES predictions and DNS data is obtained, for both velocity and scalar statistics. Finally, the developed methodology is applied to scalar transport in a confined co-axial jet mixing flow, for which experimental data are available. Agreement with statistically averaged fields for both velocity and scalar, is demonstrated to be very good, and a considerable improvement over the standard eddy viscosity RANS approach. Illustrations are presented of predicted time-resolved information e.g. time histories, and scalar pdf predictions. The LES results are shown, even using a simple Smagorinsky SGS model, to predict (correctly) lower values of the turbulent Prandtl number in the free shear regions of the flow, compared to higher values in the wall-affected regions. The ability to predict turbulent Prandtl number variations (rather than input these as in combustor RANS CFD models) is an important and promising feature of the LES approach for combustor flow simulation since it is known to be important in determining combustor exit temperature traverse.     

Highly Resolved LES and URANS of Turbulent Buoyancy-Driven Flow Within Inclined Differentially-Heated Enclosures

In the present study, both LES and unsteady RANS computations are presented, for turbulent natural convection of air inside differentially-heated rectangular tilted cavities using a finite volume code (Code_Saturne), for an aspect ratio of H/L?=?28.6 and Rayleigh number of 0.86×10⁶. Attention is focused on two angles of inclination: 15° to the horizontal with hot lower and cold upper wall, the 15° unstable case, and the mirror image of this case where the angle is the same but with a hot upper and cold lower wall, the 15° stable case. In accordance with recent experimental data, the LES computations for both the stable and unstable tilted cavities returned three-dimensional time-averaged flow fields. In the case of the unstably stratified enclosure, the flow is highly unsteady with coherent turbulent structures in the core of the enclosure. Time-averaged temperature, velocity and resolved turbulence intensities resulting from LES computations show close agreement to measured data. Subsequent comparisons of different URANS schemes with LES are used in order to explore to what extent these models are able to reproduce the large-scale unsteady flow structures. All URANS schemes have been found to be able to reproduce the 3-D unsteady flow features present in the 15° unstable cavity. However, the low-Reynolds-number model tested, as well as requiring a high resolution near-wall grid, also needed a finer grid in the core region than the high-Reynolds-number models, thus making it computationally very expensive. Flow within the 15° stable cavity also shows some 3-D features, although it is significantly less unsteady, and the URANS models tested here have been less successful in reproducing this flow pattern. The overall heat transfer is presented here for both differentially heated enclosures.     

Investigation of the sensitivity of turbulent closures and coupling of hybrid RANS‐LES models for predicting flow fields with separation and reattachment

Hybrid models have found widespread applications for simulation of wall‐bounded flows at high Reynolds numbers. Typically, these models employ Reynolds‐averaged Navier–Stokes (RANS) and large eddy simulation (LES) in the near‐body and off‐body regions, respectively. A number of coupling strategies between the RANS and LES regions have been proposed, tested, and applied in the literature with varying degree of success. Linear eddy‐viscosity models (LEVM) are often used for the closure of turbulent stress tensor in RANS and LES regions. LEVM incorrectly predicts the anisotropy of Reynolds normal stress at the RANS‐LES interface region. To overcome this issue, use of non‐linear eddy‐viscosity models (NLEVM) have started receiving attention. In this study, a generic non‐linear blended modeling framework for performing hybrid simulations is proposed. Flow over the periodic hills is used as the test case for model evaluation. This case is chosen due to complex flow physics with simplified geometry. Analysis of the simulations suggests that the non‐linear hybrid models show a better performance than linear hybrid models. It is also observed that the non‐linear closures are less sensitive to the RANS‐LES coupling and grid resolution. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical study on supersonic mixing and combustion with hydrogen injection upstream of a cavity flameholder

Characteristics of supersonic mixing and combustion with hydrogen injection upstream of a cavity flameholder are investigated numerically using hybrid RANS/LES (Reynolds-Averaged Navier–Stokes/Large-Eddy Simulation) method. Two types of inflow boundary layer are considered. One is a laminar-like boundary layer with inflow thickness of $\delta_{\inf } = 0.0$ and the other is a turbulent boundary layer with inflow thickness of $\delta_{\inf } = 2.5\,{\text{mm}}$ . The hybrid RANS/LES method acts as a DES (Detached Eddy Simulation) model for the laminar-like inflow condition and a wall-modeled LES for the turbulent inflow condition where the recycling/rescaling method is adopted. Although the turbulent inflow seems to have just minor influences on the supersonic cavity flow without fuel injection, its effects on the mixing and combustion processes are great. It is found that the unsteady turbulent structures in upstream incoming boundary layer interact with the injection jet, resulting in fluctuations of the upstream recirculation region and bow shock, and induce quick dispersion of the hydrogen fuel jet, which enhances the mixing as well as subsequent combustion.     

Vortical structures and pressure pulsations in draft tube of a Francis-99 turbine at part load: RANS and hybrid RANS/LES analysis

Recognizing the limitations of the conventional linear-eddy-viscosity (LEVM) Reynolds-averaged Navier–Stokes (RANS) models to reproduce complex three-dimensional unsteady flows in hydraulic machinery, we performed a comparative assessment of a second-moment (Re-stress model, RSM) RANS closure and a hybrid RANS/LES method in capturing the flow and vortical structures in the draft tube of a Francis hydroturbine at off-design conditions. Considered is a case of part load (PL) at a flow rate of only 35% of the best efficiency point (BEP) characterised by multiple unsteady vortex systems. Despite some remaining uncertainties in generating the inflow conditions, both approaches reproduced reasonably well the measured mean velocity and the rms of its fluctuations, as well as the pressure spectrum with peaks detecting the precessing vortex core. In contrast to the common LEVMs, the Re-stress closure showed sufficient receptivity to intrinsic unsteadiness and reproduced well the overall flow and vortical patterns as well as the associated pressure pulsations in accord with the experiments. The hybrid RANS/LES method gave similar predictions as the RSM, but resolving a wider range of scales, which however, showed no significant effect on the dynamics of the dominant processing vortex core and the pressure pulsations.     

Unsteady CFD simulations of a pump in part load conditions using scale-adaptive simulation

The scope of this work is to demonstrate the applicability of an eddy resolving turbulence model in a turbomachinery configuration. The model combines the Large Eddy Simulation (LES) and the Reynolds Averaged Navier Stokes (RANS) approach. The point of interest of the present investigation is the unsteady rotating stall phenomenon occurring at low part load conditions. Since RANS turbulence models often fail to predict separation correctly, a LES like model is expected to give superior results. In this investigation the scale-adaptive simulation (SAS) model is used. This model avoids the grid dependence appearing in the Detached Eddy Simulation (DES) modelling strategy. The simulations are validated with transient measurement data. The present results demonstrate, that both models are able to predict the major stall frequency at part load. Results are similar for URANS and SAS, with advantages in predicting minor stall frequencies for the turbulence resolving model.     

Unsteady Simulation of a Jet-in-crossflow

The present work represents the first attempt at studying the unsteadiness of jet-in-crossflows using the Reynolds-Averaged Navier-Stokes (RANS) equations. In the past, numerical studies of jet-in-crossflows, with the exception of the few attempts of DNS and LES. are always carried out under the assumption that the flowfield is steady and symmetric about the jet center plane. Growing experimental evidence suggests that this flow is unsteady and periodic. The present work verified the existence of periodicity in jet-in-crossflows numerically. The numerically predicted Slrouhal number is in good agreement with experimental observations. In addition, the present simulation of the jet-in-crossflow were able to resolve the details of flow structures such as the oscillating wake structures, the jet shear-layer vortices, and the counter-rotating vortex pairs.     

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.     

Reconstruction of turbulent fluctuations for hybrid RANS/LES simulations using a Synthetic-Eddy Method

A coupling methodology between an upstream Reynolds Averaged Navier–Stokes (RANS) simulation and a Large Eddy Simulation (LES) further downstream is presented. The focus of this work is on the RANS-to-LES interface inside an attached turbulent boundary layer, where an unsteady LES content has to be explicitly generated from a steady RANS solution. The performance of the Synthetic-Eddy Method (SEM), which generates realistic synthetic eddies at the inflow of the LES, is investigated on a wide variety of turbulent flows, from simple channel and square duct flows to the flow over an airfoil trailing edge. The SEM is compared to other existing methods of generation of synthetic turbulence for LES, and is shown to reduce substantially the distance required to develop realistic turbulence downstream of the inlet.     

Detached Eddy Simulation of Flow around A-Airfoil

In this study the flow around the Aérospatiale A-airfoil at maximum lift (α = 13.3°) and Re = 2 × 10⁶ is investigated by RANS and DES to quantify the influence of transient flow patterns on the quality of the flow prediction. This flow features a highly unsteady pressure-induced trailing-edge separation. The fairly high Reynolds number makes LES rather unattractive from a practical point of view as the numerical costs scale with ≈ Re^1.95 [9] and lead to an unreasonable effort to carry out these simulations. Like LES, DES is designed to capture the unsteady flow features and consequently is supposed to give results superior to RANS. The RANS solution captures the separation, but is notable to predict transient vortex shedding. The application of DES to a 2D domain leads to results comparable to those obtained by 2D-RANS. Only the 3D-DES is capable of predicting both the unsteady flow features and the development of the boundary layer on the airfoil. As expected, the spanwise resolution has a strong impact on the vortex shedding and on the size of the separation bubble. Attention must be paid to both accommodating the full extent of the flow structures and providing the adequate spatial resolution in order to accomplish satisfactory results – that compare favourably with LES – at only a fraction of computational cost.     

Stochastic Simulation of Scalar Mixing Capturing Unsteadiness and Small-scale Structure Based on Mean-flow Properties

A key limitation of Reynolds-Average Navier-Stokes (RANS) simulation of mixing and reaction in turbulent flows is the lack of resolution of small-scale structure and associated unsteadiness. Various subgrid models formulated in state space have been developed to complement the RANS solution in this regard. We here introduce a physical-space formulation that captures unsteady advective and diffusive processes at all scales of the turbulent flow. The approach is a 3D construction based on the Linear Eddy Model (LEM), involving three orthogonally intersecting arrays of 1D LEM domains, and coupled so as to capture the 3D character of fluid trajectories. To illustrate the model performance of the 3D LEM-based formulation, here termed LEM3D, multi-stream mixing in a turbulent round jet is simulated using measured mean-flow properties as input. Comparison to scalar cross-correlation coefficients and other measured properties of this mixing configuration indicate that the LEM3D approach, in conjunction with flow properties that are provided by steady-state models, is a useful representation of complex turbulent mixing processes that would otherwise be difficult to capture within a steady-state CFD framework.     

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.     

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.     

超声速燃烧火焰稳定凹腔质量交换特性的数值研究

对超声速冷流条件下用于超燃冲压发动机的凹腔火焰稳定器的质量交换特性进行研 究. 采用混合RANS/LES方法对非定常流场进行数值模拟,考虑了凹腔的长深比和后缘角度两 个关键参数. 计算得到了凹腔剪切层拟序结构的演化过程. 对凹腔压力振荡历程进行幅频分 析,所得到的频率和理论分析结果及一些计算结果十分吻合. 结果表明,凹腔的长深比和后 缘倾角对凹腔质量交换过程都有较大的影响. 凹腔驻留时间随着长深比的增加而增加. 随着 后缘倾角减小,长深比较小的凹腔的驻留时间增加,长深比较大的凹腔则完全相反. 该计算 结果为凹腔设计提供了新认识.     

A Zonal RANS/LES Approach for Noise Sources Prediction

Hybrid CFD/CAA methods have generally to be used for the numerical simulation of trailing-edge noise (see [9, 20] for instance). This study focuses on the first step of such hybrid methods, which is to predict the unsteady aerodynamic sources by the mean of a 3D unsteady simulation of the flow. Such a simulation is however generally still away from the numerical capabilities of ‘usual’ supercomputers. This paper investigates the use of a zonal LES method (based on the NLDE – Non-Linear Disturbance Equations – technique) for the numerical prediction of the aerodynamic noise sources. This method makes it possible to perform only zonal LES close to the main elements responsible of sound generation, while the overall configuration is only treated by a RANS approach. Attention will be paid to the specific boundary treatment at the interface between the RANS and LES regions. More precisely, the problem of the generation of turbulent inflow conditions for the LES region will be carefully addressed. The method is first assessed in the simulation of a flat plate ended by a blunted trailing-edge, and then applied to the simulation of the flow over a NACA0012 airfoil with blunted trailing-edge.