A LES/PDF simulator on block-structured meshes

A block-structured mesh large-eddy simulation (LES)/probability density function (PDF) simulator is developed within the OpenFOAM framework for computational modelling of complex turbulent reacting flows. The LES/PDF solver is a hybrid solution methodology consisting of (i) a finite-volume (FV) method for solving the filtered mass and momentum equations (LES solver), and (ii) a Lagrangian particle-based Monte Carlo algorithm (PDF solver) for solving the modelled transport equation of the filtered joint PDF of compositions. Both the LES and the PDF methods are developed and combined to form a hybrid LES/PDF simulator entirely within the OpenFOAM framework. The in situ adaptive tabulation method [S.B. Pope, Computationally efficient implementation of combustion chemistry using in situ adaptive tabulation, Combust. Theory Model. 1 (1997), pp. 41–63; L. Lu, S.R. Lantz, Z. Ren, and B.S. Pope, Computationally efficient implementation of combustion chemistry in parallel PDF calculations, J. Comput. Phys. 228 (2009), pp. 5490–5525] is incorporated into the new LES/PDF solver for efficient computations of combustion chemistry with detailed reaction kinetics. The method is designed to utilise a block-structured mesh and can readily be extended to unstructured grids. The three-stage velocity interpolation method of Zhang and Haworth [A general mass consistency algorithm for hybrid particle/finite-volume PDF methods, J. Comput. Phys. 194 (2004), pp. 156–193] is adapted to interpolate the LES velocity field onto particle locations accurately and to enforce the consistency between LES and PDF fields at the numerical solution level. The hybrid algorithm is fully parallelised using the conventional domain decomposition approach. A detailed examination of the effects of each stage and the overall performance of the velocity interpolation algorithm is performed. Accurate coupling of the LES and PDF solvers is demonstrated using the one-way coupling methodology. Then the fully two-way coupled LES/PDF solver is successfully applied to simulate the Sandia Flame-D, and a turbulent non-swirling premixed flame and a turbulent swirling stratified flame from the Cambridge turbulent stratified flame series [M.S. Sweeney, S. Hochgreb, M.J. Dunn, and R.S. Barlow, The structure of turbulent stratified and premixed methane/air flames I: Non-swirling flows, Combust. Flame 159 (2012), pp. 2896–2911; M.S. Sweeney, S. Hochgreb, M.J. Dunn, and R.S. Barlow, The structure of turbulent stratified and premixed methane/air flames II: Swirling flows, Combust. Flame 159 (2012), pp. 2912–2929]. It is found that the LES/PDF method is very robust and the results are in good agreement with the experimental data for both flames.     

Large-eddy simulation of a bluff-body stabilised turbulent premixed flame using the transported flame surface density approach

A premixed propane–air flame stabilised on a triangular bluff body in a model jet-engine afterburner configuration is investigated using large-eddy simulation (LES). The reaction rate source term for turbulent premixed combustion is closed using the transported flame surface density (TFSD) model. In this approach, there is no need to assume local equilibrium between the generation and destruction of subgrid FSD, as commonly done in simple algebraic closure models. Instead, the key processes that create and destroy FSD are accounted for explicitly. This allows the model to capture large-scale unsteady flame propagation in the presence of combustion instabilities, or in situations where the flame encounters progressive wrinkling with time. In this study, comprehensive validation of the numerical method is carried out. For the non-reacting flow, good agreement for both the time-averaged and root-mean-square velocity fields are obtained, and the Karman type vortex shedding behaviour seen in the experiment is well represented. For the reacting flow, two mesh configurations are used to investigate the sensitivity of the LES results to the numerical resolution. Profiles for the velocity and temperature fields exhibit good agreement with the experimental data for both the coarse and dense mesh. This demonstrates the capability of LES coupled with the TFSD approach in representing the highly unsteady premixed combustion observed in this configuration. The instantaneous flow pattern and turbulent flame behaviour are discussed, and the differences between the non-reacting and reacting flow are described through visualisation of vortical structures and their interaction with the flame. Lastly, the generation and destruction of FSD are evaluated by examining the individual terms in the FSD transport equation. Localised regions where straining, curvature and propagation are each dominant are observed, highlighting the importance of non-equilibrium effects of FSD generation and destruction in the model afterburner.     

LES modelling of turbulent non-premixed jet flames with correlated dynamic adaptive chemistry

Large eddy simulations (LES) for turbulent flames with detailed kinetic mechanisms have received growing interest. However, a direct implementation of detailed kinetic mechanisms in LES modelling of turbulent combustion remains a challenge due to the requirement of huge computational resources. An on-the-fly mechanism reduction method named correlated dynamic adaptive chemistry (CoDAC) is proposed to overcome this issue. A LES was conducted for Sandia Flame-D, with the reaction mechanism of GRI-Mech 3.0 consisting of 53 species and 325 reactions. The reduction threshold used in LES was obtained a priori by using auto-ignition model and partially stirred reactor (PaSR) with pairwise mixing model. LES results with CoDAC are in good agreement with experimental data and those without reduction. The conditional mean of the number of selected species indicates that a large size of locally reduced mechanism is required in the reaction zone where CH₄ is destructed. A computational time analysis shows that the PaSR model predicts better than the auto-ignition model on the wall time reduction with CoDAC in LES.     

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.     

A direct approach to generalised multiple mapping conditioning for selected turbulent diffusion flame cases

This work presents a direct and transparent interpretation of two concepts for modelling turbulent combustion: generalised Multiple Mapping Conditioning (MMC) and sparse-Lagrangian Large Eddy Simulation (LES). The MMC approach is presented as a hybrid between the Probability Density Function (PDF) method and approaches based on conditioning (e.g. Conditional Moment Closure, flamelet, etc.). The sparse-Lagrangian approach, which allows for a dramatic reduction of computational cost, is viewed as an alternative interpretation of the Filtered Density Function (FDF) methods. This work presents simulations of several turbulent diffusion flame cases and discusses the universality of the localness parameter between these cases and the universality of sparse-Lagrangian FDF methods with MMC.     

Large eddy simulation of a turbulent spray jet generated by high-pressure injection: impact of the in-nozzle flow

The turbulence properties of a gas spray jet generated through the injection of a high-pressure atomising spray were studied in the context of the Euler–Lagrange formulation using a large eddy simulation (LES) approach. The study's main aim was to investigate nozzle flow effects on flow and turbulence statistics for a two-phase jet flow in the near and far fields of the spray. The study investigated the injection of an existing in-nozzle flow at the spray inlet. Simulation results were compared with measurements obtained for non-evaporating sprays under quasi-steady conditions, which, in practice, correspond to a long-duration injection. Both simulated high-injection pressure cases showed a good agreement with experimental data. A two-stage evolution of the spray was observed. The flow of the gas jet first developed along the spray cone, followed by the appearance of highly vortical structures around and within the spray. Cavitation in the nozzle flow resulted in an observed asymmetry of the generated gas spray jet: on the gas spray jet side of the nozzle (relative to the cavitation side), stronger jet spreading occurred, and droplets were smaller and slower. The spray jet typically exhibited an ellipsoidal cross-sectional shape. The paper provides a characterisation of the asymmetry that persisted for the duration of the spray evolution (at least for the distances considered here).     

Snapshot POD analysis of transient flow in the pilot stage of a jet pipe servo valve

The spatial evolution of a turbulent flow in the pilot stage of a jet pipe servo valve at the inlet pressure and deflection angle of the jet pipe is investigated using a large eddy simulation (LES). The pressure of the same flow field is measured by a high frequency dynamic pressure sensor in the experiments and is compared with the results of the LES, as well as their root-mean-square (RMS) and fast Fourier transform (FFT) results. The results of experiments and LES are in good agreement, indicating that LES is able to predict the flow dynamics. Velocity datasets based on LES are utilised to conduct the snapshot proper orthogonal decomposition (snapshot POD) technique. The snapshot POD analysis results of the first 4 modes show a full ability to directly visualise details of the coherent structures. The influences of the inlet pressure and deflection angle of the jet pipe are also discussed. Under different inlet pressures, the velocity eigenfunctions of the first mode are similar, while the locations and strengths of the vortices in high modes are different. The Lamb-Oseen vortices that affect the trajectory of jet streams are observed in the vicinity of the entrances of receiver channels only in the first mode, and several spindly vortices appear in the region of ?5?y/n?相似文献   

Comparisons of Lagrangian and Eulerian PDF methods in simulations of non-premixed turbulent jet flames with moderate-to-strong turbulence-chemistry interactions

Transported probability density function (PDF) methods have been applied widely and effectively for modelling turbulent reacting flows. In most applications of PDF methods to date, Lagrangian particle Monte Carlo algorithms have been used to solve a modelled PDF transport equation. However, Lagrangian particle PDF methods are computationally intensive and are not readily integrated into conventional Eulerian computational fluid dynamics (CFD) codes. Eulerian field PDF methods have been proposed as an alternative. Here a systematic comparison is performed among three methods for solving the same underlying modelled composition PDF transport equation: a consistent hybrid Lagrangian particle/Eulerian mesh (LPEM) method, a stochastic Eulerian field (SEF) method and a deterministic Eulerian field method with a direct-quadrature-method-of-moments closure (a multi-environment PDF-MEPDF method). The comparisons have been made in simulations of a series of three non-premixed, piloted methane–air turbulent jet flames that exhibit progressively increasing levels of local extinction and turbulence-chemistry interactions: Sandia/TUD flames D, E and F. The three PDF methods have been implemented using the same underlying CFD solver, and results obtained using the three methods have been compared using (to the extent possible) equivalent physical models and numerical parameters. Reasonably converged mean and rms scalar profiles are obtained using 40 particles per cell for the LPEM method or 40 Eulerian fields for the SEF method. Results from these stochastic methods are compared with results obtained using two- and three-environment MEPDF methods. The relative advantages and disadvantages of each method in terms of accuracy and computational requirements are explored and identified. In general, the results obtained from the two stochastic methods (LPEM and SEF) are very similar, and are in closer agreement with experimental measurements than those obtained using the MEPDF method, while MEPDF is the most computationally efficient of the three methods. These and other findings are discussed in detail.     

Proposed stochastic parameterisations of subgrid turbulence in large eddy simulations of turbulent channel flow

Stochastic and deterministic subgrid parameterisations are developed for the large eddy simulation (LES) of a turbulent channel flow with friction-velocity-based Reynolds number of Re_τ = 950 and centreline-based Reynolds number of Re₀ = 20,580. The subgrid model coefficients (eddy viscosities) are determined from the statistics of truncated reference direct numerical simulations (DNSs). The stochastic subgrid model consists of a mean-field shift, a drain eddy viscosity acting on the resolved field and a stochastic backscatter force of variance proportional to the backscatter eddy viscosity. The deterministic variant consists of a net eddy viscosity acting on the resolved field, which represents the net effect of the drain and backscatter. LES adopting the stochastic and deterministic models is shown to reproduce the time-averaged kinetic energy spectra of the DNS within the resolved scales.     

Transported scalar PDF calculations of autoignition of a hydrogen jet in a heated turbulent co-flow

The joint scalar PDF method, as implemented in FLUENT, was used to simulate the autoignition of a jet of hydrogen in a turbulent co-flow of heated air. While the autoignition phenomenon is intermittent in the experiment, ensemble-averaged data on the effect of the flow on ignition length are available, which enables us to compare them with the steady state calculations.Results of sensitivity tests showed that the choice of chemical mechanism affects the calculation more than the mixing model and model constants. Further calculations for different initial conditions (i.e. temperature and velocity of the jet T _jet and U _jet and the co-flow T _air and U _air) have been done using a set of parameters selected after the sensitivity study. Scatter plots and conditional scalar profiles confirmed that the ignition is always initiated in lean mixture fractions. The ignition length was predicted with good accuracy for the case of U _jet>U _air but not so well for the case of U _jet≈U _air. For the equal velocity case, increasing the velocity resulted in delayed autoignition time (defined as the ignition length divided by the mean velocity), in agreement with the experimental trend. The results give credence to the use of the joint scalar PDF method for autoignition in non-premixed flows.     

基于FGM的湍流射流扩散火焰超大涡模拟

为了研究超大涡模拟(VLES)预测湍流燃烧问题的能力,本文采用VLES结合基于假定概率密度函数的火焰面生成流型(FGM)建表湍流燃烧模型对值班甲烷/空气湍流射流扩散火焰(Sandia Flame D)开展了高精度数值研究,并与实验结果进行了详细比较.结果表明本文发展的VLES-FGM方法可以较准确地预测出湍流射流扩散火...     

Large eddy simulation of turbulent premixed combustion using tabulated detailed chemistry and presumed probability density function

A method of chemistry tabulation combined with presumed probability density function (PDF) is applied to simulate piloted premixed jet burner flames with high Karlovitz number using large eddy simulation. Thermo-chemistry states are tabulated by the combination of auto-ignition and extended auto-ignition model. To evaluate the predictive capability of the proposed tabulation method to represent the thermo-chemistry states under the condition of different fresh gases temperature, a-priori study is conducted by performing idealised transient one-dimensional premixed flame simulations. Presumed PDF is used to involve the interaction of turbulence and flame with beta PDF to model the reaction progress variable distribution. Two presumed PDF models, Dirichlet distribution and independent beta distribution, respectively, are applied for representing the interaction between two mixture fractions that are associated with three inlet streams. Comparisons of statistical results show that two presumed PDF models for the two mixture fractions are both capable of predicting temperature and major species profiles, however, they are shown to have a significant effect on the predictions for intermediate species. An analysis of the thermo-chemical state-space representation of the sub-grid scale (SGS) combustion model is performed by comparing correlations between the carbon monoxide mass fraction and temperature. The SGS combustion model based on the proposed chemistry tabulation can reasonably capture the peak value and change trend of intermediate species. Aspects regarding model extensions to adequately predict the peak location of intermediate species are discussed.     

Large eddy simulation of a turbulent nonpremixed piloted methane jet flame (Sandia Flame D)

Large eddy simulation (LES) is conducted of the Sandia Flame D [Proc. Combust. Inst. 27 (1998) 1087, Sandia National Laboratories (2004)], which is a turbulent piloted nonpremixed methane jet flame. The subgrid scale (SGS) closure is based on the scalar filtered mass density function (SFMDF) methodology [J. Fluid Mech. 401 (1999) 85]. The SFMDF is basically the mass weighted probability density function (PDF) of the SGS scalar quantities [Turbulent Flows (2000)]. For this flame (which exhibits little local extinction), a simple flamelet model is used to relate the instantaneous composition to the mixture fraction. The modelled SFMDF transport equation is solved by a hybrid finite-difference/Monte Carlo scheme. This is the first LES of a realistic turbulent flame using the transported PDF method as the SGS closure. The results via this method capture important features of the flame as observed experimentally.     

Numerical simulations of the flow field ahead of an accelerating flame in an obstructed channel

The development of the unburned gas flow field ahead of a flame front in an obstructed channel was investigated using large eddy simulation (LES). The standard Smagorinsky–Lilly and dynamic Smagorinsky–Lilly subgrid models were used in these simulations. The geometry is essentially two-dimensional. The fence-type obstacles were placed on the top and bottom surfaces of a square cross-section channel, equally spaced along the channel length at the channel height. The laminar rollup of a vortex downstream of each obstacle, transition to turbulence, and growth of a recirculation zone between consecutive obstacles were observed in the simulations. By restricting the simulations to the early stages of the flame acceleration and by varying the domain width and domain length, the three-dimensionality of the vortex rollup process was investigated. It was found that initially the rollup process was two-dimensional and unaffected by the domain length and width. As the recirculation zone grew to fill the streamwise gap between obstacles, the length and width of the computational domain started to affect the simulation results. Three-dimensional flow structures formed within the shear layer, which was generated near the obstacle tips, and the core flow was affected by large-scale turbulence. The simulation predictions were compared to experimental schlieren images of the convection of helium tracer. The development of recirculation zones resulted in the formation of contraction and expansion regions near the obstacles, which significantly affected the centerline gas velocity. Oscillations in the centerline unburned gas velocity were found to be the dominate cause for the experimentally observed early flame-tip velocity oscillations. At later simulation times, regular oscillations in the unburned streamwise gas velocity were not observed, which is contrary to the experimental evidence. This suggests that fluctuations in the burning rate might be the source of the late flame-tip velocity oscillations. The effect of the obstacle blockage ratio (BR) on the development of the unburned gas flow field was also investigated by varying the obstacle height. Simulation predictions show favorable agreement with the experimental results and indicate that turbulence production increases with increasing obstacle BR.     

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.     

Modelling of turbulent lifted jet flames using flamelets: a priori assessment and a posteriori validation

This study focuses on the modelling of turbulent lifted jet flames using flamelets and a presumed Probability Density Function (PDF) approach with interest in both flame lift-off height and flame brush structure. First, flamelet models used to capture contributions from premixed and non-premixed modes of the partially premixed combustion in the lifted jet flame are assessed using a Direct Numerical Simulation (DNS) data for a turbulent lifted hydrogen jet flame. The joint PDFs of mixture fraction Z and progress variable c, including their statistical correlation, are obtained using a copula method, which is also validated using the DNS data. The statistically independent PDFs are found to be generally inadequate to represent the joint PDFs from the DNS data. The effects of Z–c correlation and the contribution from the non-premixed combustion mode on the flame lift-off height are studied systematically by including one effect at a time in the simulations used for a posteriori validation. A simple model including the effects of chemical kinetics and scalar dissipation rate is suggested and used for non-premixed combustion contributions. The results clearly show that both Z–c correlation and non-premixed combustion effects are required in the premixed flamelets approach to get good agreement with the measured flame lift-off heights as a function of jet velocity. The flame brush structure reported in earlier experimental studies is also captured reasonably well for various axial positions. It seems that flame stabilisation is influenced by both premixed and non-premixed combustion modes, and their mutual influences.     

LES studies of the flow in a swirl gas combustor

Environmental and other practical concerns have led to the development of compact gas turbine combustors burning lean mixtures leading to potentially low CO and NO_x emissions. The compact design requires efficient atomization and mixing together with a compact premixed flame. Associated with these requirements are higher temperatures, increased heat transfer, and thermal load, thus increasing the danger of combustion instabilities (causing performance deterioration and excessive mechanical loads), and possible off-design operation. Numerical simulations of reacting flows are well suited to address these issues. To this end, large eddy simulation (LES) is particularly promising. The philosophy behind LES is to explicitly simulate the large scales of the flow and the thermochemistry, affected by boundary conditions whilst modeling only the small scales, including the interaction between the flow and the combustion processes. Here, we examine the flow and the flame in a model gas turbine combustor (General Electric’s lean premixed dry low NO_x LM6000) to evaluate the potential of LES for design studies of engineering applications and to study the effects of the combustor confinement geometry on the flow and on the flame dynamics. Two LES models, a Monotone Integrated LES model with 1 and 2 step Ahrrenius chemistry, and a fractal flame-wrinkling LES model coupled to a conventional one-equation eddy-viscosity subgrid model, are used. Reasonable agreement is found when comparing predictions with experimental data and with other LES computations of the same case. Furthermore, the combustor confinement geometry is found to strongly affect the vortical flow, and hence also the flame and its dynamics.     

Large Eddy Simulation of flame stabilisation dynamics and vortex control in a lifted H2/N2 jet flame

Flame stabilisation in (highly) preheated mixture is common in several industrial applications. When the reactants are injected separately in the device (usually at high-speed), the flame is lifted so that the fuel and oxidant first mix to give an ignitable mixture. If the temperature of the mixture is adequate, it auto-ignites stabilizing the flame. Here we focus on an academic lifted jet flame and Large Eddy Simulation (LES) is used to capture the flame and auto-ignition dynamics. Comparisons with experimental data show that LES simulates accurately high OH fluctuation levels at the stabilisation location. The vortex dynamics linked to these fluctuations is analyzed and it is found that small scale coherent structures play a vital role in the auto-ignition process. These structures are axial vorticity tubes (braids) and are located relatively far (in the radial direction) from the shear-layer. As a consequence, the lift-off height varies dramatically in time leading to OH fluctuations of the order of the mean OH concentration. This scenario is monitored in the compositional space highlighting the simultaneous evolution of OH, HO ₂ and temperature. Further, different strategies for open-loop control of the flame lift-off height are tested. In order to anchor the flame at different positions downstream of the nozzle, the vortex dynamics in the shear-layer was modified. Promoting successively vortex ring and braids, the auto-ignition region was moved significantly. In particular, modified nozzle geometries impacted the formation of braids and ensured a good premixing very close to the nozzle. As a consequence, it was possible to reduce significantly the lift-off height and stabilise the flame few diameters downstream of the nozzle.     

Effects of combined dimension reduction and tabulation on the simulations of a turbulent premixed flame using a large-eddy simulation/probability density function method

A turbulent lean-premixed propane–air flame stabilised by a triangular cylinder as a flame-holder is simulated to assess the accuracy and computational efficiency of combined dimension reduction and tabulation of chemistry. The computational condition matches the Volvo rig experiments. For the reactive simulation, the Lagrangian Large-Eddy Simulation/Probability Density Function (LES/PDF) formulation is used. A novel two-way coupling approach between LES and PDF is applied to obtain resolved density to reduce its statistical fluctuations. Composition mixing is evaluated by the modified Interaction-by-Exchange with the Mean (IEM) model. A baseline case uses In Situ Adaptive Tabulation (ISAT) to calculate chemical reactions efficiently. Its results demonstrate good agreement with the experimental measurements in turbulence statistics, temperature, and minor species mass fractions. For dimension reduction, 11 and 16 represented species are chosen and a variant of Rate Controlled Constrained Equilibrium (RCCE) is applied in conjunction with ISAT to each case. All the quantities in the comparison are indistinguishable from the baseline results using ISAT only. The combined use of RCCE/ISAT reduces the computational time for chemical reaction by more than 50%. However, for the current turbulent premixed flame, chemical reaction takes only a minor portion of the overall computational cost, in contrast to non-premixed flame simulations using LES/PDF, presumably due to the restricted manifold of purely premixed flame in the composition space. Instead, composition mixing is the major contributor to cost reduction since the mean-drift term, which is computationally expensive, is computed for the reduced representation. Overall, a reduction of more than 15% in the computational cost is obtained.     

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.