Journal
COMPUTERS & FLUIDS
Volume 39, Issue 4, Pages 568-591Publisher
PERGAMON-ELSEVIER SCIENCE LTD
DOI: 10.1016/j.compfluid.2009.10.002
Keywords
Lattice Boltzmann equation; Pseudo-spectral method; Incompressible flow; Decaying turbulence; Low-order statistics; Instantaneous flow fields
Funding
- US Department of Defense
- National Science Foundation [ATM-0527140]
- National Natural Science Foundation of China [10628206]
Ask authors/readers for more resources
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.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available