Simulations of isotropic turbulent flows using lattice Boltzmann method with different forcing functions

Abstract

Three different forcing functions are used with the lattice Boltzmann method (LBM) to simulate the forced isotropic turbulence in periodic boxes at different resolutions ranging from [Formula: see text] to [Formula: see text] grid points using the D3Q19 model. The aims of this study are to examine the effect of using different forcing functions on the LBM stability; to track the development of the turbulent fields at several resolutions, to investigate the effect of the weak compressibility due to change of fluid density on the flow simulations, and to identify the effective force type. The injection is performed through adding the force randomly to the collision term. The three forcing methods depend on sine and cosine as functions of the wave numbers and space. The forcing amplitude values of [Formula: see text] and the relaxation time [Formula: see text] are fixed in all cases. The single relaxation time model is found stable at such values of the forcing amplitude and the relaxation time. However, the development of the turbulent data at the different resolutions needs about 10000 time-steps to reach the required statistical state including clear visualizations of fine scale vortices. Many simulations have been tested using different values of the relaxation time [Formula: see text] and the development of the turbulent fields is found faster with fewer time-steps but the stability of the LBM is broken at some resolutions (not necessary the higher resolution). The statistical features of all fields, such as the Taylor and the Kolmogorov micro-scales, the Taylor Reynolds number, the flatness and the skewness, are calculated and compared with the previous efforts. The worm-like vortices are visualized at all cases and it is found that more fine vortices can be extracted as the resolution increases. The energy spectrum has a reasonable Kolmogorov power law at the resolutions of [Formula: see text] and [Formula: see text], respectively. Results show that the third forcing method that uses a cosine disturbance function has the best statistical features and the finest visualized vortical structures especially at higher resolutions. Extensive discussions about the density field and its evolution with time at different forcing functions, comparison to Navier–Stokes solutions and the time development of the energy spectra for all cases are also carried out.

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

“…In this study, the Guo-method [1] is used with the SRT and MRT models, respectively. The three components of the forcing function F which is used in this study can be written as [6] F ≡…”

“…Simulations of forced homogeneous and isotropic turbulent flows are not popular because adding the force to the LBM is still under development due to the stability issue and the nature of the method. In a recent study [6], the forcing function and its amplitude, initialization of the flow field, the injection cutoff wavenumber, and the relaxation time are carefully chosen to strength the stability process at several resolutions and with different forcing functions. Also, comparison between the variable forcing techniques of the lattice Boltzmann method for simulations of box turbulent is discussed in [7] with resolutions of 128 3 and 256 3 , respectively.…”

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

“…In this study, the Guo-method [13] is used with the SRT and MRT models, respectively. The three components of the forcing function F which is used in this study can be written as [4]…”

A modified lattice Boltzmann model for microcylindrical Couette gas flows

“…Macroscale flows are usually characterized by the Reynolds number Re, and so the magnitude of t in LB models is dependent on the value of Re [33,34]. By contrast, microscale gas flows are generally featured by Kn.…”