Benjamin Graille

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, Tony Fevrier, Benjamin Graille

This paper studies the stability properties of a two dimensional relative velocity scheme for the Navier-Stokes equations. This scheme inspired by the cascaded scheme has the particularity to relax in a frame moving with a velocity field function of space and time. Its stability is studied first in a linear context then on the non linear test case of the Kelvin-Helmholtz instability. The link with the choice of the moments is put in evidence. The set of moments of the cascaded scheme improves the stability of the d'Humières scheme for small viscosities. On the contrary, a relative velocity scheme with the usual set of moments deteriorates the stability.

François Dubois, Tony Fevrier, Benjamin Graille

In this contribution, a new class of lattice Boltzmann schemes is introduced and studied. These schemes are presented in a framework that generalizes the multiple relaxation times method of d'Humières. They extend also the Geier's cascaded method. The relaxation phase takes place in a moving frame involving a set of moments depending on a given relative velocity field. We establish with the Taylor expansion method that the equivalent partial differential equations are identical to the ones obtained with the multiple relaxation times method up to the second order accuracy. The method is then performed to derive the equivalent equations up to third order accuracy.

Adeline Augier, François Dubois, Benjamin Graille

In this paper, we investigate the numerous parameters choices for linear lattice Boltzmann schemes according to the definition of the isotropic order given in \cite{ADG11}. This property---written in a general framework including all of the \ddqq schemes---can be read through a group operation. It implies some relations on the parameters of the scheme (equilibrium states and relaxation times) that give rigorous methodology to select them according to the desired order of isotropy. For acoustic applications in two spaces dimensions (namely \ddqn and \ddqt schemes) this methodology is used to propose a full description of the sets of parameters that involve isotropy of order $m$ ($m\in\{1,2,3,5\}$ for \ddqn and $m\in\{1,2\}$ for \ddqt). We then propose numerical illustrations for the \ddqn scheme.

Adeline Augier, François Dubois, Benjamin Graille

In this paper, we recall the linear version of the lattice Boltzmann schemes in the framework proposed by d'Humiéres. According to the equivalent equations we introduce a definition for a scheme to be isotropic at some order. This definition is chosen such that the equivalent equations are preserved by orthogonal transformations of the frame. The property of isotropy can be read through a group operation and then implies a sequence of relations on relaxation times and equilibrium states that characterizes a lattice Boltzmann scheme. We propose a method to select the parameters of the scheme according to the desired order of isotropy. Applying it to the D2Q9 scheme yields the classical constraints for the first and second orders and some non classical for the third and fourth orders.

Benjamin Graille, François Dubois, Tony Fevrier

We study the formal precision of the relative velocity lattice Boltzmann schemes. They differ from the d'Humières schemes by their relaxation phase: it occurs for a set of moments parametrized by a velocity field function of space and time. We deal with the asymptotics of the relative velocity schemes for one conservation law: the third order equivalent equation is exposed for an arbitrary number of dimensions and velocities.

François Dubois, Tony Février, Benjamin Graille

In this contribution, we study the theoretical and numerical stability of a bidimensional relative velocity lattice Boltzmann scheme. These relative velocity schemes introduce a velocity field parameter called "relative velocity" function of space and time. They generalize the d'Humières multiple relaxation times scheme and the cascaded automaton. This contribution studies the stability of a four velocities scheme applied to a single linear advection equation according to the value of this relative velocity. We especially compare when it is equal to 0 (multiple relaxation times scheme) or to the advection velocity ("cascaded like" scheme). The comparison is made in terms of L1 and L2 stability. The L1 stability area is fully described in terms of relaxation parameters and advection velocity for the two choices of relative velocity. These results establish that no hierarchy of these two choices exists for the L1 notion. Instead, choosing the parameter equal to the advection velocity improves the numerical L2 stability of the scheme. This choice cancels some dispersive terms and improve the numerical stability on a representative test case. We theoretically strengthen these results with a weighted L2 notion of stability.