Criterion of numerical instability of liquid state in LBE simulations

The numerical stability of the lattice Boltzmann equation (LBE) method in simulations of a fluid described by an equation of state with possible vapor-liquid phase transitions is considered. The Courant-Friedrichs-Lewy number defined by the advection term in the Boltzmann equation is exactly equal to unity in classical LBE models. However, this condition does not ensure the numerical stability of LBE simulations with the equation of state. In our numerical LBE simulations, we find out that instability arises initially in the liquid phase, even if the vapor phase and, consequently, the vapor-liquid interface are absent. We demonstrate both in numerical tests and theoretically that the numerical stability of LBE simulations requires the criterion [email protected][email protected][email protected]?"c"r to be fulfilled for the liquid phase, where [email protected]?=c"[email protected]/h is the hydrodynamic Courant number. The hydrodynamic Courant number is proportional to the speed of sound c"s, obtained from an equation of state of a fluid. This criterion is very similar to the well-known criteria of numerical stability of explicit finite difference schemes for a compressible fluid. The critical value of the Courant number [email protected]?"c"r depends neither on the temperature T, nor on the fluid velocity, nor on the form of the equation of state. This critical value is equal to [email protected]?"c"r=1.1547 for the kinetic temperature of LBE pseudo-particles @[email protected]?=1/3.

[1]  I. Karlin,et al.  Lattice Boltzmann method for thermal flow simulation on standard lattices. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[2]  G. Doolen,et al.  Discrete Boltzmann equation model for nonideal gases , 1998 .

[3]  Laura Schaefer,et al.  Equations of state in a lattice Boltzmann model , 2006 .

[4]  Li-Shi Luo,et al.  Unified Theory of Lattice Boltzmann Models for Nonideal Gases , 1998 .

[5]  Raoyang Zhang,et al.  Lattice Boltzmann method for simulations of liquid-vapor thermal flows. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[6]  P. Adler,et al.  Boundary flow condition analysis for the three-dimensional lattice Boltzmann model , 1994 .

[7]  D. A. Medvedev,et al.  On equations of state in a lattice Boltzmann method , 2009, Comput. Math. Appl..

[8]  Baoming Li,et al.  Mean-field free-energy approach to the lattice Boltzmann method for liquid-vapor and solid-fluid interfaces. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[9]  D. P. Agoris,et al.  Stochastic models of partial discharge activity in solid and liquid dielectrics , 2007 .

[10]  Yeomans,et al.  Lattice Boltzmann simulation of nonideal fluids. , 1995, Physical review letters.

[11]  P. Roache.  QUANTIFICATION OF UNCERTAINTY IN COMPUTATIONAL FLUID DYNAMICS , 1997 .

[12]  R. D. Richtmyer,et al.  Difference methods for initial-value problems , 1959 .

[13]  Shan,et al.  Lattice Boltzmann model for simulating flows with multiple phases and components. , 1993, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[14]  Y. Qian,et al.  Lattice BGK Models for Navier-Stokes Equation , 1992 .