Effects of resolution inhomogeneity in large-eddy simulation

Abstract

Large Eddy Simulation (LES) of turbulence in complex geometries is often conducted using strongly inhomogeneous resolution. The issues associated with resolution inhomogeneity are related to the noncommutativity of the filtering and differentiation operators, which introduces a commutation term into the governing equations. Neglect of this commutation term gives rise to commutation error. While the commutation error is well recognized, it is often ignored in practice. Moreover, the commutation error arising from the implicit part of the filter (i.e., projection onto the underlying discretization) has not been well investigated. Modeling the commutator between numerical projection and differentiation is crucial for correcting errors induced by resolution inhomogeneity in practical LES settings, which typically rely solely on implicit filtering. Here, we employ a multiscale asymptotic analysis to investigate the characteristics of the commutator. This provides a statistical description of the commutator, which can serve as a target for the statistical characteristics of a commutator model. Further, we investigate how commutation error manifests in simulation and demonstrate its impact on the convection of a packet of homogeneous isotropic turbulence through an inhomogeneous grid. A connection is made between the commutation error and the propagation properties of the underlying numerics. A modeling approach for the commutator is proposed that is applicable to LES with filters that include projections to the discrete solution space and that respects the numerical properties of the LES evolution equation. It may also be useful in addressing other LES modeling issues such as discretization error.