Single Relaxation Time and Multiple Revised Matrix Lattice Boltzmann Simulations of Forced Isotropic Turbulence

TL;DR

This study compares single relaxation time and revised matrix lattice Boltzmann models for simulating forced isotropic turbulence, confirming the MRT model's superiority and analyzing stability and turbulence characteristics at different resolutions.

Contribution

It introduces a comparison of SRT and MRT lattice Boltzmann models for turbulence simulation, highlighting the stability improvements and turbulence features of MRT.

Findings

MRT model shows superior turbulence simulation performance.

Turbulent velocity fields exhibit universal characteristics.

MRT maintains near-unity density PDF indicating better incompressibility.

Abstract

The single relaxation time (SRT) and the revised matrix (RM) lattice Boltzmann models are compared for simulations of three dimensional forced isotropic turbulence with resolutions of 128^3 and 256^3, respectively. The forcing technique by Guo et al. (2002) is applied with the two models using the same parameters and conditions. Some new aspects and results have been confirmed such as the superiority of the MRT model to simulate forced turbulence and using the Courant-Friedichs-Lewy condition (CFL) (Courant et al., 1967) by multiplying the velocity input with the coefficient CFL < 1.0 to overcome the stability problem and divide the output velocity data by the same CFL. The initial velocity field is chosen as u(x; 0) = 0 and the force is injected at low wave-numbers with a fixed forcing amplitude to 10^-4 for all cases. The single relaxation time is set to 0.503 in all SRT simulations. Results show that the obtained turbulent velocity fields yield universal characteristics as proven in previous theoretical, experimental and numerical studies. The Taylor Reynolds number for the simulations are found as SRT: R = 62 and MRT : R = 65 for 1283 and SRT: R= 107 and MRT: R = 82 for the case of 2563, respectively. To test the weak incompressibility for the SRT model in comparison to the MRT model, the density probability distribution function (PDF) is depicted and it is found that ? is almost unity at all timesteps for the MRT case, while a clear disturbance about unity is observed for the SRT case.