A rational length scale for large-eddy simulation of turbulence on anisotropic grids

Abstract

Due to the prohibitive cost of resolving all relevant scales, direct numerical simulations of turbulence remain unfeasible for most real-world applications. Consequently, dynamically simplified formulations are needed for coarse-grained simulations. In this regard, eddy-viscosity models for large-eddy simulation (LES) are widely used in both academia and industry. These models require a subgrid characteristic length, typically linked to the local grid size. While this length scale corresponds to the mesh step for isotropic grids, its definition for unstructured or anisotropic Cartesian meshes, such as the pancake-like meshes commonly used to capture near-wall turbulence or shear layers, remains an open question. Despite its significant influence on LES model performance, no consensus has been reached on its proper formulation. In this work, we introduce a novel subgrid characteristic length. This length scale is derived from the analysis of the entanglement between the numerical discretization and the filtering in LES. Its mathematical properties and simplicity make it a robust choice for reducing the impact of mesh anisotropies on simulation accuracy. The effectiveness of the proposed subgrid length is demonstrated through simulations of decaying isotropic turbulence and a turbulent channel flow using different codes.

F.X.T. and J.R. are supported by SIMEX project (PID2022142174OB-I00) of the Ministerio de Ciencia e Innovacion MCIN/AEI/ 10.13039/501100011033 and the European Union Next GenerationEU. The work of A.G. and A.D. has been financially supported by the Russian Science Foundation, Project No. 24-11-00287 (implementation and simulations with the NOISEtte code). Calculations were carried out on the MareNostrum 5-GPP supercomputer at BSC. We thankfully acknowledge these institutions.

Peer Reviewed

Postprint (author's final draft)