Large Eddy Simulation Of Incompressible Flows Research Articles

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

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

A conservative immersed interface method for Large-Eddy Simulation of incompressible flows

We propose a conservative, second-order accurate immersed interface method for representing incompressible fluid flows over complex three dimensional solid obstacles on a staggered Cartesian grid. The method is based on a finite-volume discretization of the incompressible Navier–Stokes equations which is modified locally in cells that are cut by the interface in such a way that accuracy and conservativity are maintained. A level-set technique is used for description and tracking of the interface geometry, so that an extension of the method to moving boundaries and flexible walls is straightforward. Numerical stability is ensured for small cells by a conservative mixing procedure. Discrete conservation and sharp representation of the fluid–solid interface render the method particularly suitable for Large-Eddy Simulations of high-Reynolds number flows. Accuracy, second-order grid convergence and robustness of the method is demonstrated for several test cases: inclined channel flow at Re = 20, flow over a square cylinder at Re = 100, flow over a circular cylinder at Re = 40, Re = 100 and Re = 3900, as well as turbulent channel flow with periodic constrictions at Re = 10,595.

Meshless approach for wall treatment in Large-Eddy Simulation

An extended meshless approach for wall-modelling in Large-Eddy Simulation of incompressible flows that may take into account several terms of Navier–Stokes equations (the pressure gradient, the unsteady term, buoyancy effects, …) is presented. The method is based on an analytical integration of the momentum equations along the normal direction to the wall. After a detailed description of the algorithm, the present approach is applied to a two-dimensional periodic hill flow at a bulk Reynolds number of 7100. It is shown that taking into account the instantaneous pressure gradient has no effect on the mean values of the velocity. However, a substantial improvement is obtained concerning the rms values.

Spectral Vanishing Viscosity Method for Large-Eddy Simulation of Turbulent Flows

An efficient spectral vanishing viscosity method for the large-eddy simulation of incompressible flows is proposed, both for standard spectral and spectral element approximations. The approach is integrated in a collocation spectral Chebyshev-Fourier solver and then used to compute the turbulent wake of a cylinder in a crossflow confined geometry (Reynolds number Re=3900)

High-order LES modeling of turbulent incompressible flow

This article presents the high-order algorithms that we have developed for large-eddy simulation of incompressible flows, and the results that have been obtained for the 3D turbulent wake of a cylinder at a Reynolds number of Re = 3900 . To cite this article: R. Pasquetti, C. R. Mecanique 333 (2005).

An approximate deconvolution model for large-eddy simulation with application to incompressible wall-bounded flows

The approximate deconvolution model (ADM) for the large-eddy simulation of incompressible flows is detailed and applied to turbulent channel flow. With this approach an approximation of the unfiltered solution is obtained by repeated filtering. Given a good approximation of the unfiltered solution, the nonlinear terms of the filtered Navier–Stokes equations can be computed directly. The effect of nonrepresented scales is modeled by a relaxation regularization involving a secondary filter operation. Large-eddy simulations are performed for incompressible channel flow at Reynolds numbers based on the friction velocity and the channel half-width of Reτ=180 and Reτ=590. Both simulations compare well with direct numerical simulation (DNS) data and show a significant improvement over results obtained with classical subgrid scale models such as the standard or the dynamic Smagorinsky model. The computational cost of ADM is lower than that of dynamic models or the velocity estimation model.

A physical-space version of the stretched-vortex subgrid-stress model for large-eddy simulation

A physical-space version of the stretched-vortex subgrid-stress model is presented and applied to large-eddy simulations of incompressible flows. This version estimates the subgrid-kinetic energy required for evaluation of the subgrid-stress tensor using local second-order structure-function information of the resolved velocity field at separations of order the local cell size. A relation between the structure function and the energy spectrum is derived using the kinematic assumptions of the stretched-vortex model for locally homogeneous anisotropic turbulence. Results of large-eddy simulations using this model are compared to experimental and direct numerical simulation data. Comparisons are shown for the decay of kinetic energy and energy spectra of decaying isotropic turbulence and for mean velocities, root-mean-square velocity fluctuations and turbulence-kinetic energy balances of channel flow at three different Reynolds numbers.