Discrete filter operators for large‐eddy simulation using high‐order spectral difference methods

The combination of a high‐order unstructured spectral difference (SD) spatial discretization scheme with sub‐grid scale (SGS) modeling for large‐eddy simulation is investigated with particular focus on the consistent implementation of a structural mixed model based on the scale similarity hypothesis. The difficult task of deriving a consistent formulation for the discrete filter within the SD element of arbitrary order led to the development of a new class of three‐dimensional constrained discrete filters. The discrete filters satisfy a set of selected criteria and are completely local within the SD element. Their weights can be automatically computed at run time from the number of solution points within each element and the expected filter cutoff length scale. The novel discrete filters can be applied to any SGS model involving explicit filtering and to a broad class of high‐order discontinuous finite element numerical schemes. The code is applied to the computation of turbulent channel flows at three Reynolds numbers, namely Re_τ = 180, 395, and 590 (based on the friction velocity u_τ and channel half‐width δ). Results from computations with and without the SGS model are compared against results from direct numerical simulation. The numerical experiments suggest that the results are sensitive to the use of the SGS model, even when a high‐order numerical scheme is used, especially when the grid resolution is kept relatively low and mostly in terms of resolved Reynolds stresses. Results obtained using existing filters based on the projection of the solution over lower‐order polynomial bases are also shown and demonstrate that these filters are inadequate for SGS modeling purposes, mostly because of their inability to enforce the selected cutoff length scale with sufficient accuracy. The use of the similarity mixed formulation proved to be particularly accurate in reproducing SGS interactions, confirming that its well‐known potential can be realized in conjunction with state‐of‐the‐art high‐order numerical schemes.Copyright © 2012 John Wiley & Sons, Ltd.     

An improved anisotropy-resolving subgrid-scale model for flows in laminar–turbulent transition region

Some types of mixed subgrid-scale (SGS) models combining an isotropic eddy-viscosity model and a scale-similarity model can be used to effectively improve the accuracy of large eddy simulation (LES) in predicting wall turbulence. Abe (2013) has recently proposed a stabilized mixed model that maintains its computational stability through a unique procedure that prevents the energy transfer between the grid-scale (GS) and SGS components induced by the scale-similarity term. At the same time, since this model can successfully predict the anisotropy of the SGS stress, the predictive performance, particularly at coarse grid resolutions, is remarkably improved in comparison with other mixed models. However, since the stabilized anisotropy-resolving SGS model includes a transport equation of the SGS turbulence energy, k_SGS, containing a production term proportional to the square root of k_SGS, its applicability to flows with both laminar and turbulent regions is not so high. This is because such a production term causes k_SGS to self-reproduce. Consequently, the laminar–turbulent transition region predicted by this model depends on the inflow or initial condition of k_SGS. To resolve these issues, in the present study, the mixed-timescale (MTS) SGS model proposed by Inagaki et al. (2005) is introduced into the stabilized mixed model as the isotropic eddy-viscosity part and the production term in the k_SGS transport equation. In the MTS model, the SGS turbulence energy, k_es, estimated by filtering the instantaneous flow field is used. Since the k_es approaches zero by itself in the laminar flow region, the self-reproduction property brought about by using the conventional k_SGS transport equation model is eliminated in this modified model. Therefore, this modification is expected to enhance the applicability of the model to flows with both laminar and turbulent regions. The model performance is tested in plane channel flows with different Reynolds numbers and in a backward-facing step flow. The results demonstrate that the proposed model successfully predicts a parabolic velocity profile under laminar flow conditions and reduces the dependence on the grid resolution to the same degree as the unmodified model by Abe (2013) for turbulent flow conditions. Moreover, it is shown that the present model is effective at transitional Reynolds numbers. Furthermore, the present model successfully provides accurate results for the backward-facing step flow with various grid resolutions. Thus, the proposed model is considered to be a refined anisotropy-resolving SGS model applicable to laminar, transitional, and turbulent flows.     

Multi-scale analysis of subgrid stress and energy dissipation in turbulent channel flow

In present study, the subgrid scale （SGS） stress and dissipation for multiscale formulation of large eddy simulation are analyzed using the data of turbulent channel flow at Ret = 180 obtained by direct numerical simulation. It is found that the small scale SGS stress is much smaller than the large scale SGS stress for all the stress components. The dominant contributor to large scale SGS stress is the cross stress between small scale and subgrid scale motions, while the cross stress between large scale and subgrid scale motions make major contributions to small scale SGS stress. The energy transfer from resolved large scales to subgrid scales is mainly caused by SGS Reynolds stress, while that between resolved small scales and subgrid scales are mainly due to the cross stress. The multiscale formulation of SGS models are evaluated a priori, and it is found that the small- small model is superior to other variants in terms of SGS dissipation.     

Large-Eddy and Direct Numerical Simulations of the Bachalo-Johnson Flow with Shock-Induced Separation

The Bachalo-Johnson experiment on an axisymmetric bump has been a primary validation case for turbulence models in shock-boundary-layer interactions since the 1980’s. In the present work, Wall-Modelled Large-Eddy Simulations (WMLES) of this flow were conducted using Improved Delayed Detached-Eddy Simulation (IDDES) as the sub-grid-scale (SGS) and wall model, with a synthetic turbulence generator, expecting close enough agreement with experiment. However, the WMLES results are disappointing, even in terms of the shock position, even though the results from two grids with 4.7 × 10⁸ and 1.6 × 10⁹ cells respectively agree well with each other. This failure of grid refinement to warn of an inaccurate simulation is of great concern, and the reasons for it are explored. We then conducted a Direct Numerical Simulation (DNS) embedded in the LES over a reduced domain, with 8 × 10⁹ grid cells. The DNS has a far more accurate shock position and overall pressure distribution. The skin friction in the favourable pressure gradient is also much higher than in the LES; thus, wide differences appear upstream of the shock wave, most probably caused by the rapid acceleration which leads to atypical shear-stress profiles. Other SGS models were tried, and performed worse than IDDES. The DNS essentially fulfils the initial expectations although in a reduced domain and provides data for turbulence-modelling studies, for instance by extracting an effective eddy viscosity from it. The most noticeable remaining disagreement with experiment is over the Reynolds shear stress.     

Effects of the Similarity Model in Finite-Difference LES of Isotropic Turbulence Using a Lagrangian Dynamic Mixed Model

A Lagrangian dynamic formulation of the mixed similarity subgrid (SGS) model for large-eddy simulation (LES) of turbulence is proposed. In this model, averaging is performed over fluid trajectories, which makes the model applicable to complex flows without directions of statistical homogeneity. An alternative version based on a Taylor series expansion (nonlinear mixed model) is also examined. The Lagrangian models are implemented in a finite difference code and tested in forced and decaying isotropic turbulence. As comparison, the dynamic Smagorinsky model and volume-averaged formulations of the mixed models are also tested. Good results are obtained, except in the case of low-resolution LES (32³) of decaying turbulence, where the similarity coefficient becomes negative due to the fact that the test-filter scale exceeds the integral scale of turbulence. At a higher resolution (64³), the dynamic similarity coefficient is positive and good agreement is found between predicted and measured kinetic energy evolution. Compared to the eddy viscosity term, the similarity or the nonlinear terms contribute significantly to both SGS dissipation of kinetic energy and SGS force. In order to dynamically test the accuracy of the modeling, the error incurred in satisfying the Germano identity is evaluated. It is found that the dynamic Smagorinksy model generates a very large error, only 3% lower than the worst-case scenario without model. Addition of the similarity or nonlinear term decreases the error by up to about 50%, confirming that it represents a more realistic parameterization than the Smagorinsky model alone.     

On the logarithmic-law constants and the turbulent boundary layer at low Reynolds numbers

It is shown that a family of formally derived similarity solutions describe to leading order the outer region of a turbulent boundary layer for all Reynolds numbers for which the layer satisfies the logarithmic law-of-the-wall. The family includes Coles' [1] hypothesis. For consistency with this hypothesis and the logarithmic law-of-the-wall, it is further shown that the constants in the latter form the product κC=2+O(ε), suggesting the logarithmic law of the wall be written $${U \mathord{\left/ {\vphantom {U {U_\tau = \kappa ^{ - 1} }}} \right. \kern-\nulldelimiterspace} {U_\tau = \kappa ^{ - 1} }}\ln \left( {e^2 U_\tau {y \mathord{\left/ {\vphantom {y \nu }} \right. \kern-\nulldelimiterspace} \nu }} \right) + O\left( \in \right).$$ A range of data are reprocessed to determine the skin friction coefficientC _f using κC = 2 and these collapse well when plotted against momentum thickness Reynolds number, Re _θ . It is also shown that the form parameter, Π, in Coles hypothesis is not unique but is determined by history effects peculiar to the boundary layer. Expressions are derived forC _f (Re _θ ) and the shape factorH (Re _θ ); both agree closely with the data and are valid over all Reynolds numbers for which the logarithmic law of the wall is satisfied.     

Development and Analysis of Wall Models for Internal Combustion Engine Simulations Using High-speed Micro-PIV Measurements

The performance, efficiency and emissions of internal combustion (IC) engines are affected by the thermo-viscous boundary layer region and heat transfer. Computational models for the prediction of engine performance typically rely on equilibrium wall-function models to overcome the need for resolving the viscous boundary layer structure. The wall shear stress and heat flux are obtained as boundary conditions for the outer flow calculation. However, these equilibrium wall-function models are typically derived by considering canonical flow configurations, introducing substantial modeling assumptions that are not necessarily justified for in-cylinder flows. The objective of this work is to assess the validity of several model approximations that are commonly introduced in the development of wall-function models for IC-engine applications. This examination is performed by considering crank-angle resolved high-resolution micro-particle image velocimetry (µ-PIV) measurements in a spark-ignition direct-injection single cylinder engine. Using these measurements, the performance of an algebraic equilibrium wall-function model commonly used in RANS and LES IC-engine simulations is evaluated. By identifying shortcomings of this model, a non-equilibrium differential wall model is developed and two different closures are considered for the determination of the turbulent viscosity. It is shown that both wall models provide adequate predictions if applied inside the viscous sublayer. However, the equilibrium wall-function model consistently underpredicts the shear stress if applied in the log-layer. In contrast, the non-equilibrium wall model provides improved predictions of the near-wall region and shear stress irrespective of the wall distance and the piston location. By utilizing the experimental data, significant adverse pressure gradients due to the large vortical motion inside the cylinder (induced by tumble, swirl and turbulence) are observed and included in the non-equilibrium wall model to further improve the model performance. These investigations are complemented by developing a consistent wall heat transfer model, and simulation results are compared against the equilibrium wall-function model and Woschni’s empirical correlation.     

Development and application of wall-function treatments for turbulent forced and mixed convection flows

CFD calculations of turbulent flow near smooth walls generally employ one of two broad strategies to resolve the very influential, complex, but thin near-wall viscosity-affected sub-layer. One approach uses a fine numerical mesh and a turbulence model incorporating viscous influences; the other employs “wall functions”—formulae that attempt to account for the overall resistance of the sublayer to momentum and heat transport. The latter requires only a fraction of the computational effort of the former and is thus strongly favoured for industrial calculations. However, the wall-function performance is often poor, partly because of inappropriate implementations and partly because the schemes themselves have inherent limitations.The present paper reviews the evolution of wall-function strategies. Attention is then given to two new schemes developed by the authors, one based on an analytical treatment and the other on a numerical resolution of the near-wall sub-layer. Several applications are shown of mixed and forced convection.     

A general optimal formulation for the dynamic Smagorinsky subgrid‐scale stress model

In this paper, a general optimal formulation for the dynamic Smagorinsky subgrid‐scale (SGS) stress model is reported. The Smagorinsky constitutive relation has been revisited from the perspective of functional variation and optimization. The local error density of the dynamic Smagorinsky SGS model has been minimized directly to determine the model coefficient C_S. A sufficient and necessary condition for optimizing the SGS model is obtained and an orthogonal condition (OC), which governs the instantaneous spatial distribution of the optimal dynamic model coefficient, is formulated. The OC is a useful general optimization condition, which unifies several classical dynamic SGS modelling formulations reported in the literature. In addition, the OC also results in a new dynamic model in the form of a Picard's integral equation. The approximation tensorial space for the projected Leonard stress is identified and the physical meaning for several basic grid and test‐grid level tensors is systematically discussed. Numerical simulations of turbulent Couette flow are used to validate the new model formulation as represented by the Picard's integral equation for Reynolds numbers ranging from 1500 to 7050 (based on one half of the velocity difference of the two plates and the channel height). The relative magnitudes of the Smagorinsky constitutive parameters have been investigated, including the model coefficient, SGS viscosity and filtered strain rate tensor. In general, this paper focuses on investigation of fundamental mathematical and physical properties of the popular Smagorinsky constitutive relation and its related dynamic modelling optimization procedure. Copyright © 2005 John Wiley & Sons Ltd.     

Near-wall flow in a three-dimensional boundary layer on the endwall of a 30° bend

 Turbulence measurements are reported on the three-dimensional turbulent boundary layer along the centerline of the flat endwall in a 30° bend. Profiles of mean velocities and Reynolds stresses were obtained down to y ⁺≈2 for the mean flow and y ⁺≈8 for the turbulent stresses. Mean velocity data collapsed well on a simple law-of-the-wall based on the magnitude of the resultant velocity. The turbulence intensity and turbulent shear stress magnitude both increased with increased three-dimensionality. The ratio of these two quantities, the a ₁ structure parameter, decreased in the central regions of the boundary layer and showed profile similarity for y ⁺<50. The shear stress vector angle lagged behind the velocity gradient vector angle in the outer region of the boundary layer, however there was an indication that the shear stress vector tends to lead the velocity gradient vector close to the wall. Received: 16 July 1996/Accepted: 14 July 1997     

LES and DES Investigations of Turbulent Flow over a Sphere at Re = 10,000

Large Eddy Simulation (LES) using a dynamic Smagorinsky type subgridstress (SGS) model and Detached Eddy Simulation (DES) are applied toprediction and investigation of the flow around a sphere at a Reynoldsnumber of 10⁴ in the subcritical regime. In this regime the boundarylayers at separation are laminar, and transition to turbulence occursfarther downstream in the separated shear layers via Kelvin–Helmholtz(K–H) instabilities. The dynamic eddy viscosity model of Germano et al.(Physics of Fluids 3 (1991) 1760–1765) is used in the LES, while the current implementation of the DESemploys a formulation based on the Spalart–Allmaras (S–A) model. DES isa hybrid approach in which the closure is a modification to theproduction/destruction term of the original Reynolds-AveragedNavier–Stokes (RANS) model, reducing to RANS in the attached regions,and to LES away from the wall. In the present work where we simulate theflow over a sphere in the subcritical regime in which the boundarylayers at separation are laminar, DES can be viewed as LES with adifferent SGS model. Effects of the discretization scheme used toapproximate the convective terms are considered, along with sensitivityof predictions to changes in the additional model coefficient, C _DES, in the DES formulation. DES and LES yield similar predictions of the wakestructure, large-scale vortex shedding and the Strouhal numberassociated with the low frequency mode in the wake. Predictions ofquantities such as the drag coefficient, wake frequencies, position oflaminar separation on the sphere, and the mean pressure andskin-friction distributions along the sphere are in good agreement withthe measurements of Achenbach (Journal of Fluid Mechanics 54 (1972) 565–575). Predictions of the primaryReynolds shear stress, turbulent kinetic energy, eddy viscosity, andturbulent dissipation for the two models are also similar. In addition,both models successfully resolve the formation of the vortex tubes inthe detached shear layers along with the value of the Strouhal numberassociated with the high frequency instability mode, provided that thelevel of numerical dissipation introduced by the discretization schemeremains sufficiently low. Flow physics investigations are focused onunderstanding the wake structure in the subcritical regime.     

A new wall-damping function for large eddy simulation employing Kolmogorov velocity scale

A new wall-damping function, based on the Kolmogorov velocity scale, for large eddy simulation (LES) is proposed, which accounts for the near-wall effect. To calculate the Kolmogorov velocity scale, u_ε, the dissipation rate of turbulent energy, ε, is needed. In LES, however, the dissipation rate is generally not solved, unlike in the Reynolds averaged Navier-Stokes (RANS) simulations, e.g., k-ε models. Although, in some previous studies, the dissipation rate of the subgrid-scale (SGS) turbulent energy, ε_SGS, is used instead of ε in calculating the Kolmogorov velocity scale, the scale obtained using such a method overly depends on the grid resolution employed and is generally inappropriate. Accordingly, the wall-damping function using the incorrect velocity scale also depends on the grid resolution and gives an inadequate wall effect. This is because ε_SGS contains only the components in the scale smaller than the grid-filter width, which obviously varies with the grid resolution employed. In this study, to overcome this problem, we propose a method for estimating the Kolmogorov velocity scale with a technique of conversion in LES, and the estimated one is utilized in the wall-damping function. The revised wall-damping function for LES is tested in channel flows and a backward-facing step flow. The results show that it yields a proper near-wall effect in all test cases which cover a wide range of grid resolution and Reynolds numbers. It is also shown that all three kinds of SGS models incorporating the present wall-damping function provide good predictions, and it is effective both in one-equation and 0-equation SGS models. These results suggest that the use of the proposed wall-damping function is a refined and versatile near-wall treatment in LES with various kinds of SGS models.     

DNS of premixed turbulent V-flame: coupling spectral and finite difference methods

To allow for a reliable examination of the interaction between velocity fluctuations, acoustics and combustion, a novel numerical procedure is discussed in which a spectral solution of the Navier–Stokes equations is directly associated to a high-order finite difference fully compressible DNS solver (sixth order PADE). Using this combination of high-order solvers with accurate boundary conditions, simulations have been performed where a turbulent premixed V-shape flame develops in grid turbulence. In the light of the DNS results, a sub-model for premixed turbulent combustion is analyzed. To cite this article: R. Hauguel et al., C. R. Mecanique 333 (2005).     

Constrained large-eddy simulation of turbulent flow over inhomogeneous rough surfaces

In this work we extend the method of the constrained large-eddy simulation(CLES) to simulate the turbulent flow over inhomogeneous rough walls. In the original concept of CLES, the subgrid-scale(SGS) stress is constrained so that the mean part and the fluctuation part of the SGS stress can be modelled separately to improve the accuracy of the simulation result. Here in the simulation of the rough-wall flows, we propose to interpret the extra stress terms in the CLES formulation as the roughness-induced stress so that the roughness inhomogeneity can be incorporated by modifying the formulation of the constrained SGS stress. This is examined with the simulations of the channel flow with the spanwise alternating high/low roughness strips. Then the CLES method is employed to investigate the temporal response of the turbulence to the change of the wall condition from rough to smooth. We demonstrate that the temporal development of the internal boundary layer is just similar to that in a spatial rough-tosmooth transition process, and the spanwise roughness inhomogeneity has little impact on the transition process.     

A dynamic two-equation Sub Grid Scale model

In this paper we demonstrate that the transport equation of the generalised subgrid scale (SGS) turbulent stress tensor is form-invariant but not frame-indifferent under Euclidean transformations of the frame. A new closure equation between the generalized SGS turbulent stress tensor and the resolved kinematic quantities is proposed. The closure equation at the basis of the proposed model (Two-Equation Model, TEM): a) respects the principle of the turbulence frame indifference [1]; b) takes into account both the anisotropy of the turbulence velocity scales and turbulence length scales; c) removes any balance assumption between the production and dissipation of SGS turbulent kinetic energy; d) assumes scale similarity in the definition of the second-order tensor representing the turbulent velocity scales. In the proposed model: a) the closure coefficient C which appears in the constitutive equation is uniquely determined without using Germanos dynamic procedure [2]; b) the generalized SGS turbulent stress tensor is related exclusively to the generalized SGS turbulent kinetic energy (which is calculated by means of its balance equation) and the modified Leonard tensor; c) the viscous dissipation of the generalized SGS turbulent kinetic energy is calculated by solving the balance equation. The proposed model is tested for a turbulent channel flow at Reynolds numbers (based on friction velocity and channel half-width) ranging from 180 to 2340.Received: 11 February 2004, Accepted: 20 August 2004, Published online: 22 February 2005PACS: 02.60.Cb, 47.27.Eq, 47.11. + j Correspondence to: F. Gallerano     

Evaluations of SGS Combustion,Scalar Flux and Stress Models in a Turbulent Jet Premixed Flame

A newly developed fractal dynamic SGS (FDSGS) combustion model and a scale self-recognition mixed (SSRM) SGS stress model are evaluated along with other SGS combustion, scalar flux and stress models in a priori and a posteriori manners using DNS data of a hydrogen-air turbulent plane jet premixed flame. A posteriori tests reveal that the LES using the FDSGS combustion model can predict the combustion field well in terms of mean temperature distributions and peak positions in the transverse distributions of filtered reaction progress variable fluctuations. A priori and a posteriori tests of the scalar flux models show that a model proposed by Clark et al. accurately predicts the counter-gradient transport as well as the gradient diffusion, and introduction of the model of Clark et al. into the LES yields slightly better predictions of the filtered progress variable fluctuations than that of a gradient diffusion model. Evaluations of the stress models reveal that the LES with the SSRM model predicts the velocity fluctuations well compared to that with the Smagorinsky model.     

Similarity solution of mixed convection over a horizontal moving plate

Investigation to the mixed convective heat and mass transfer over a horizontal plate has been carried out. By applying transformation group theory to analysis of the governing equations of continuity, momentum, energy and diffusion, we show the existence of similarity solution for the problem provided that the temperature and concentration at the wall are proportional to x ^4/(7-5n) and that the moving speed of the plate is proportional to x ^(3-n)/(7-5n), and further obtain a similarity representation of the problem. The similarity equations have been solved numerically by a fourth-order Runge–Kutta scheme. The numerical results obtained for Pr=0.72 and various values of the parameters Sc, K ₁, K ₂ and K ₃ reveals the influence of the parameters on the flow, heat and mass transfer behavior.     

A Non-Linear SGS Model Based On The Spatial Velocity Increment

A new subgrid scale (SGS) modelling concept for large-eddy simulation (LES) of incompressible flow is proposed based on the three-dimensional spatial velocity increment δ u _i . The new model is inspired by the structure function formulation developed by Métais and Lesieur [39] and applied in the context of the scale similarity type formulation. First, the similarity between the SGS stress tensor τ _ij and the velocity increment tensor Q _ij = δ u _i δ u _j is analyzed analytically and numerically using a priori tests of fully developed pipe flow at Re _τ = 180. Both forward and backward energy transfers between resolved and unresolved scales of the flow are well predicted with a SGS model based on Q _ij . Secondly, a posteriori tests are performed for two families of turbulent shear flows. LES of fully developed pipe flow up to Re _τ = 520 and LES of round turbulent jet at Re _D = 25000 carried out with a dynamic version of the model provide promising results that confirm the power of this approach for wall-bounded and free shear flows.     

Parallel FEM LES with one-equation subgrid-scale model for incompressible flows

This article develops a parallel large-eddy simulation (LES) with a one-equation subgrid-scale (SGS) model based on the Galerkin finite element method and three-dimensional (3D) brick elements. The governing filtered Navier–Stokes equations were solved by a second-order accurate fractional-step method, which decomposed the implicit velocity–pressure coupling in incompressible flow and segregated the solution to the advection and diffusion terms. The transport equation for the SGS turbulent kinetic energy was solved to calculate the SGS processes. This FEM LES model was applied to study the turbulence of the benchmark open channel flow at a Reynolds number Re_τ = 180 (based on the friction velocity and channel height) using different model constants and grid resolutions. By comparing the turbulence statistics calculated by the current model with those obtained from direct numerical simulation (DNS) and experiments in literature, an optimum set of model constants for the current FEM LES model was established. The budgets of turbulent kinetic energy and vertical Reynolds stress were then analysed for the open channel flow. Finally, the flow structures were visualised to further reveal some important characteristics. It was demonstrated that the current model with the optimum model constants can predict well the organised structure near the wall and free surface, and can be further applied to other fundamental and engineering applications.