Validating the Boltzmann approach to the Large-Eddy simulations of forced homogeneous incompressible turbulence

TL;DR

This study validates that lattice Boltzmann method-based large-eddy simulations of forced homogeneous turbulence operate within the hydrodynamic regime, supporting the use of Smagorinsky closure and confirming the kinetic approach's effectiveness.

Contribution

It demonstrates that LBM-Smagorinsky LES remains in the hydrodynamic regime for forced homogeneous turbulence, validating its use and the Smagorinsky model in this context.

Findings

Kinetic turbulent Knudsen number is around 10^{-3}

Spectra and flatness follow canonical LES behavior

Lattice Boltzmann LES operates in the hydrodynamic regime

Abstract

The simulation of turbulent flows remains a central challenge, as even our most powerful computers cannot resolve the finest scales of motion in many flows of practical interest. As a result, the effects of unresolved scales on large eddies must be modelled via closures and coarse-graining procedures. Large-eddy simulation (LES) traditionally coarse-grains Navier-Stokes equations using Smagorinsky's effective viscosity model. This has the merit of simplicity but fails to account for strong non-equilibrium effects, as they typically arise in most flows in the vicinity of solid walls, the reason being that the notion of eddy viscosity assumes scale separation between small and large eddies, an assumption that fails for high-Reynolds flows far from equilibrium. The lattice Boltzmann method (LBM) offers an alternative by coarse-graining at the kinetic level, potentially capturing non-equilibrium effects beyond reach of hydrodynamic closures. This paper addresses whether LBM-Smagorinsky LES of forced homogeneous isotropic turbulence (FHIT) exhibits kinetic behavior. We test whether the turbulent Knudsen number $K_t$, measuring scale separation, reaches order one (kinetic regime) or remains asymptotically small (hydrodynamic regime). Using reference DNS ($800^3$) and iso-Reynolds LES ($100^3$) at $Re = 2 \times 10^4$, we quantify $K_t$ via spatial maps, temporal statistics, energy spectra, and higher-order moments. Results show $K_t \sim O(10^{-3})$, strictly positive without negative excursions, with spectra and flatness following canonical LES behavior. We conclude that despite its kinetic formulation, LBM-Smagorinsky LES operates in the hydrodynamic regime, with small FHIT eddies remaining in local equilibrium with larger ones, validating Smagorinsky viscosity and confirming that LBM-LES functions as conventional hydrodynamic LES while preserving LBM efficiency and locality.