Development of Compressible Large-Eddy Simulations Combining High-Order Schemes and Wall Modeling

Abstract

Compressible large-eddy simulations combining high-order methods with a wall model have been developed in order to compute wall-bounded flows at high Reynolds numbers. The high-order methods consist of low-dissipation and low-dispersion implicit finite volume schemes for spatial discretization on structured grids. In the first part, the procedure used to apply these schemes in near-wall regions is presented. It is based on a ghost cell reconstruction. Its validity is assessed by performing the large-eddy simulation of a turbulent channel flow at a friction Reynolds number of Re τ 395. In the second part, to consider flows at higher Reynolds numbers, a large-eddy simulation approach using a wall model is proposed. The coupling between the wall model and the high-order schemes is described. The performance of the approach is evaluated by simulating a turbulent channel flow at Re τ 2000 using meshes with grid spacings of Δx 100, 200, and 300 in the streamwise direction and Δ 50, 100, and 150 in the wall-normal and spanwise directions (in wall units). The effects of the choice of the point used for data exchange between the wall model and the large-eddy simulation algorithm, as well as of the positions of the ghost cells used for the coupling, are examined by performing additional computations in which these parameters vary. The results are in agreement with direct numerical simulation data. In particular, the turbulent intensities obtained in the logarithmic region of the boundary layers of the channel flow are successfully predicted.

Exponential distribution functions for positivity-preserving lattice Boltzmann schemes: Application to 2D compressible flow simulations

“…Such an approximation is pretty common when simulating high-Reynolds number flows with NSF solvers. 67,72 In that case, it is sufficient to enforce up to the trace of third-order equilibrium moments, which allows for the use of 39velocity and 21-velocity lattices in three and two dimensions, i.e., D3Q39 and D2Q21, respectively. In two dimensions, the corresponding set of eight equilibrium moments reads…”

A technique of flux reconstruction at the interfaces of nonconforming grids for aeroacoustic simulations

“…The filter is employed on a uniformly spaced grid due to a coordinate transform. For LES computations, the filter also plays the role of a subgrid‐scale model, relaxing turbulent energy at high frequencies . Time integration is performed by applying a low‐storage six‐stage Runge‐Kutta algorithm .…”

Extended integral wall-model for large-eddy simulations of compressible wall-bounded turbulent flows

“…They provide an explicit relation between the mean velocity and temperature and the wall fluxes τ w and φ w . From these two laws, an analytical wall-model 23,56 can be built,…”

“…For high-Reynolds number flows like those considered in this study, the computational domain size is big enough to capture turbulent eddies and has been used in previous studies. 23,56,79 Two uniform grids G1 and G2 are considered: 25 × 21 × 21 and 49 × 41 × 41. On the former mesh, the first off-wall cell is located at y 1 = 0.05δ, while on the latter, y 1 = 0.025δ.…”

“…For finer grids as encountered in resolved LES, vortex rings are added to the initial flow to mimic turbulent structures. 56,82 However, WMLES grids are not fine enough to discretize these kinds of structures.…”