Pierre Lallemand

Pierre Lallemand, François Dubois

We show that a single particle distribution for the D2Q13 lattice Boltzmann scheme can simulate coupled effects involving advection and diffusion of velocity and temperature. We consider various test cases: non-linear waves with periodic boundary conditions, a test case with buoyancy, propagation of transverse waves, Couette and Poiseuille flows. We test various boundary conditions and propose to mix bounce-back and anti-bounce-back numerical boundary conditions to take into account velocity and temperature Dirichlet conditions. We present also first results for the de Vahl Davis heated cavity. Our results are compared with the coupled D2Q9-D2Q5 lattice Boltzmann approach for the Boussinesq system and with an elementary finite differences solver for the compressible Navier-Stokes equations.

François Dubois, Pierre Lallemand, Christian Obrecht et al.

We show that the asymptotic properties of the link-wise artificial compressibility method are not compatible with a correct approximation of fluid properties. We propose to adapt the previous method through a framework suggested by the Taylor expansion method and to replace first order terms in the expansion by appropriate three or five points finite differences and to add non linear terms. The "FD-LBM" scheme obtained by this method is tested in two dimensions for shear wave, Stokes modes and Poiseuille flow. The results are compared with the usual lattice Boltzmann method in the framework of multiple relaxation times.

Bruce Boghosian, François Dubois, Benjamin Graille et al.

We consider the D1Q3 lattice Boltzmann scheme with an acoustic scale for the simulation of diffusive processes. When the mesh is refined while holding the diffusivity constant, we first obtain asymptotic convergence. When the mesh size tends to zero, however, this convergence breaks down in a curious fashion, and we observe qualitative discrepancies from analytical solutions of the heat equation. In this work, a new asymptotic analysis is derived to explain this phenomenon using the Taylor expansion method, and a partial differential equation of acoustic type is obtained in the asymptotic limit. We show that the error between the D1Q3 numerical solution and a finite-difference approximation of this acoustic-type partial differential equation tends to zero in the asymptotic limit. In addition, a wave vector analysis of this asymptotic regime demonstrates that the dispersion equation has nontrivial complex eigenvalues, a sign of underlying propagation phenomena, and a portent of the unusual convergence properties mentioned above.

François Dubois, Benjamin Graille, Pierre Lallemand

We consider multi relaxation times lattice Boltzmann scheme with two particle distributions for the thermal Navier Stokes equations formulated with conservation of mass and momentum and dissipation of volumic entropy.Linear stability is taken into consideration to determine a coupling between two coefficients of dissipation.We present interesting numerical results for one-dimensional strong nonlinear acoustic waves with shocks.

François Dubois, Pierre Lallemand

In this contribution we extend the Taylor expansion method proposed previously by one of us and establish equivalent partial differential equations of DDH lattice Boltzmann scheme at an arbitrary order of accuracy. We derive formally the associated dynamical equations for classical thermal and linear fluid models in one to three space dimensions. We use this approach to adjust relaxation parameters in order to enforce fourth order accuracy for thermal model and diffusive relaxation modes of the Stokes problem. We apply the resulting scheme for numerical computation of associated eigenmodes and compare our results with analytical references.

François Dubois, Pierre Lallemand

Using the Taylor expansion method, we show that it is possible to improve the lattice Boltzmann method for acoustic applications. We derive a formal expansion of the eigenvalues of the discrete approximation and fit the parameters of the scheme to enforce fourth order accuracy. The corresponding discrete equations are solved with the help of symbolic manipulation. The solutions are explicited in the case of D3Q27 lattice Boltzmann scheme. Various numerical tests support the coherence of this approach.