A new development of the dynamic procedure in large-eddy simulation based on a Finite Volume integral approach. Application to stratified turbulence

A Finite Volume-based large-eddy simulation method is proposed along with a suitable extension of the dynamic modelling procedure that takes into account for the integral formulation of the governing filtered equations. Discussion about the misleading interpretation of FV in some literature is addressed. Then, the classical Germano identity is congruently rewritten in such a way that the determination of the modelling parameters does not require any arbitrary averaging procedure and thus retains a fully local character. The numerical modelling of stratified turbulence is the specific problem considered in this study, as an archetypal of simple geophysical flows. The original scaling formulation of the dynamic sub-grid scale model proposed by Wong and Lilly (Phys. Fluids 6(6), 1994) is suitably extended to the present integral formulation. This approach is preferred with respect to traditional ones since the eddy coefficients can be independently computed by avoiding the addition of unjustified buoyancy production terms in the constitutive equations. Simple scaling arguments allow us not to use the equilibrium hypothesis according to which the dissipation rate should equal the sub-grid scale energy production. A careful a priori analysis of the relevance of the test filter shape as well as the filter-to-grid ratio is reported. Large-eddy simulation results are a posteriori compared with a reference pseudo-spectral direct numerical solution that is suitably post-filtered in order to have a meaningful comparison. In particular, the spectral distribution of kinetic and thermal energy as well as the viscosity and diffusivity sub-grid scale profiles are illustrated. The good performances of the proposed method, in terms of both evolutions of global quantities and statistics, are very promising for the future development and application of the method.