Comparison of the lattice Boltzmann and pseudo-spectral methods for decaying turbulence: Low-order statistics

We conduct a detailed comparison of the lattice Boltzmann equation (LBE) and the pseudo-spectral (PS) methods for direct numerical simulations (DNS) of the decaying homogeneous isotropic turbulence in a three-dimensional periodic cube. We use a mesh size of N-3 = 128(3) and the Taylor micro-scale Reynolds number 24.35 <= Re-lambda <= 72.37, and carry out all simulations to t approximate to 30 tau(0), where tau(0) is the turbulence turnover time. In the PS method, the second-order Adam-Bashforth scheme is used to numerically integrate the nonlinear term while the viscous term is treated exactly. We compare the following quantities computed by the LBE and PS methods: instantaneous velocity omega and vorticity omega fields, and statistical quantities such as, the total energy K(t) and the energy spectrum E(k, t), the dissipation rate epsilon(t), the root-meansquared (rms) pressure fluctuation delta p(t) and the pressure spectrum P(k, t), and the skewness and flatness of the velocity derivative. Our results show that the LBE method performs very well when compared to the PS method in terms of accuracy and efficiency: the instantaneous flow fields, u and omega, and all the statistical quantities - except the rms pressure fluctuation delta p(t) and the pressure spectrum P(k, t) - computed from the LBE and PS methods agree well with each other, provided that the initial flow field is adequately resolved by both methods. We note that delta p(t) and P(k, t) computed from the two methods agree with each other in a period of time much shorter than that for other quantities, indicating that the pressure field p computed by using the LBE method is less accurate than other quantities. The skewness and flatness computed from the LBE method contain high-frequency oscillations due to acoustic waves in the system, which are absent in PS methods. Our results indicate that the resolution requirement for the LBE method is delta x/eta(0) < 1.0, approximately twice of the requirement for PS methods, where delta x and eta(0) are the grid spacing and the initial Kolmogorov length, respectively. Overall, the LBE method is shown to be a reliable and accurate method for the DNS of decaying turbulence. (C) 2010 Elevier Ltd. All rights reserved.