scholarly journals Linear and brute force stability of orthogonal moment multiple-relaxation-time lattice Boltzmann methods applied to homogeneous isotropic turbulence

Author(s):  
Stephan Simonis ◽  
Marc Haussmann ◽  
Louis Kronberg ◽  
Willy Dörfler ◽  
Mathias J. Krause

Multiple-relaxation-time (MRT) lattice Boltzmann methods (LBM) based on orthogonal moments exhibit lattice Mach number dependent instabilities in diffusive scaling. The present work renders an explicit formulation of stability sets for orthogonal moment MRT LBM. The stability sets are defined via the spectral radius of linearized amplification matrices of the MRT collision operator with variable relaxation frequencies. Numerical investigations are carried out for the three-dimensional Taylor–Green vortex benchmark at Reynolds number 1600. Extensive brute force computations of specific relaxation frequency ranges for the full test case are opposed to the von Neumann stability set prediction. Based on that, we prove numerically that a scan over the full wave space, including scaled mean flow variations, is required to draw conclusions on the overall stability of LBM in turbulent flow simulations. Furthermore, the von Neumann results show that a grid dependence is hardly possible to include in the notion of linear stability for LBM. Lastly, via brute force stability investigations based on empirical data from a total number of 22 696 simulations, the existence of a deterministic influence of the grid resolution is deduced. This article is part of the theme issue ‘Progress in mesoscale methods for fluid dynamics simulation’.

2016 ◽  
Vol 20 (2) ◽  
pp. 374-404 ◽  
Author(s):  
S. Karimi ◽  
K. B. Nakshatrala

AbstractThe Lattice Boltzmann Method (LBM) has established itself as a popular numerical method in computational fluid dynamics. Several advancements have been recently made in LBM, which include multiple-relaxation-time LBM to simulate anisotropic advection-diffusion processes. Because of the importance of LBM simulations for transport problems in subsurface and reactive flows, one needs to study the accuracy and structure preserving properties of numerical solutions under the LBM. The solutions to advective-diffusive systems are known to satisfy maximum principles, comparison principles, the non-negative constraint, and the decay property. In this paper, using several numerical experiments, it will be shown that current single- and multiple-relaxation-time lattice Boltzmann methods fail to preserve these mathematical properties for transient diffusion-type equations. We will also show that these violations may not be removed by simply refining the discretization parameters. More importantly, it will be shown that meeting stability conditions alone does not guarantee the preservation of the aforementioned mathematical principles and physical constraints in the discrete setting. A discussion on the source of these violations and possible approaches to avoid them is included. A condition to guarantee the non-negativity of concentration under LBM in the case of isotropic diffusion is also derived. The impact of this research is twofold. First, the study poses several outstanding research problems, which should guide researchers to develop LBM-based formulations for transport problems that respect important mathematical properties and physical constraints in the discrete setting. This paper can also serve as a good source of benchmark problems for such future research endeavors. Second, this study cautions the practitioners of the LBM for transport problems with the associated numerical deficiencies of the LBM, and provides guidelines for performing predictive simulations of advective-diffusive processes using the LBM.


2007 ◽  
Vol 18 (04) ◽  
pp. 635-643 ◽  
Author(s):  
XIAOWEN SHAN ◽  
HUDONG CHEN

We formulate a simple extension to the Bhatnagar-Gross-Krook collision model by expanding the distribution function in Hermite polynomials and assigning a relaxation time to each hydrodynamic moment. By discretizing the velocity space, multiple-relaxation-time lattice Boltzmann models can be constructed. The transport coefficients are analytically calculated and numerically verified. At the lowest order, allowing different relaxation rates for the second and third Hermite components results in a variable Prandtl number. Comparing with the previously proposed multiple-relaxation-time lattice Boltzmann models, the present formulation is general in the sense that it is independent of the underlying lattice structure and does not require a procedure for transformation of base vectors.


Sign in / Sign up

Export Citation Format

Share Document