Source author record

Benjamin Graille

Benjamin Graille appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

Does the multiresolution lattice Boltzmann method allow to deal with waves passing through mesh jumps?

We consider an adaptive multiresolution-based lattice Boltzmann scheme, which we have recently introduced and studied from the perspective of the error control and the theory of the equivalent equations. This numerical strategy leads to high compression rates, error control and its high accuracy has been explained on uniform and dynamically adaptive grids. However, one key issue with non-uniform meshes within the framework of lattice Boltzmann schemes is to properly handle acoustic waves passing through a level jump of the grid. It usually yields spurious effects, in particular reflected waves. In this paper, we propose a simple mono-dimensional test-case for the linear wave equation with a fixed adapted mesh characterized by a potentially large level jump. We investigate this configuration with our original strategy and prove that we can handle and control the amplitude of the reflected wave, which is of fourth order in the space step of the finest mesh. Numerical illustrations show that the proposed strategy outperforms the existing methods in the literature and allow to assess the ability of the method to handle the mesh jump properly.

preprint2022arXiv

Multiresolution-based mesh adaptation and error control for lattice Boltzmann methods with applications to hyperbolic conservation laws

Lattice Boltzmann Methods (LBM) stand out for their simplicity and computational efficiency while offering the possibility of simulating complex phenomena. While they are optimal for Cartesian meshes, adapted meshes have traditionally been a stumbling block since it is difficult to predict the right physics through various levels of meshes. In this work, we design a class of fully adaptive LBM methods with dynamic mesh adaptation and error control relying on multiresolution analysis. This wavelet-based approach allows to adapt the mesh based on the regularity of the solution and leads to a very efficient compression of the solution without loosing its quality and with the preservation of the properties of the original LBM method on the finest grid. This yields a general approach for a large spectrum of schemes and allows precise error bounds, without the need for deep modifications on the reference scheme. An error analysis is proposed. For the purpose of assessing the approach, we conduct a series of test-cases for various schemes and scalar and systems of conservation laws, where solutions with shocks are to be found and local mesh adaptation is especially relevant. Theoretical estimates are retrieved while a reduced memory footprint is observed. It paves the way to an implementation in a multi-dimensional framework and high computational efficiency of the method for both parabolic and hyperbolic equations, which is the subject of a companion paper.

preprint2020arXiv

Unexpected convergence of lattice Boltzmann schemes

In this work, we study numerically the convergence of the scalar D2Q9 lattice Boltzmann scheme with multiple relaxation times when the time step is proportional to the space step and tends to zero. We do this by a combination of theory and numerical experiment. The classical formal analysis when all the relaxation parameters are fixed and the time step tends to zero shows that the numerical solution converges to solutions of the heat equation, with a constraint connecting the diffusivity, the space step and the coefficient of relaxation of the momentum. If the diffusivity is fixed and the space step tends to zero, the relaxation parameter for the momentum is very small, causing a discrepency between the previous analysis and the numerical results. We propose a new analysis of the method for this specific situation of evanescent relaxation, based on the dispersion equation of the lattice Boltzmann scheme. A new asymptotic partial differential equation, the damped acoustic system, is emergent as a result of this formal analysis. Complementary numerical experiments establish the convergence of the scalar D2Q9 lattice Boltzmann scheme with multiple relaxation times and acoustic scaling in this specific case of evanescent relaxation towards the numerical solution of the damped acoustic system.

preprint2015arXiv

Lattice Boltzmann schemes with relative velocities

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.

preprint2015arXiv

On the stability of a relative velocity lattice Boltzmann scheme for compressible Navier-Stokes equations

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.

preprint2015arXiv

Stability of a bidimensional relative velocity lattice Boltzmann scheme

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.

preprint2015arXiv

Third order equivalent equation for the relative velocity lattice Boltzmann schemes with one conservation law

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.

preprint2014arXiv

On rotational invariance of lattice Boltzmann schemes

We propose the derivation of acoustic-type isotropic partial differential equations that are equivalent to linear lattice Boltzmann schemes with a density scalar field and a momentum vector field as conserved moments. The corresponding linear equivalent partial differential equations are generated with a new "Berliner version" of the Taylor expansion method. The details of the implementation are presented. These ideas are applied for the D2Q9, D2Q13, D3Q19 and D3Q27 lattice Boltzmann schemes. Some limitations associated with necessary stability conditions are also presented.

preprint2011arXiv

Isotropy conditions for lattice Boltzmann schemes. Application to D2Q9

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.

preprint2011arXiv

Linear Lattice Boltzmann Schemes for Acoustic: parameters choices and isotropy properties

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.

preprint2007arXiv

Kinetic Theory of Plasmas: Translational Energy

In the present contribution, we derive from kinetic theory a unified fluid model for multicomponent plasmas by accounting for the electromagnetic field influence. We deal with a possible thermal nonequilibrium of the translational energy of the particles, neglecting their internal energy and the reactive collisions. Given the strong disparity of mass between the electrons and heavy particles, such as molecules, atoms, and ions, we conduct a dimensional analysis of the Boltzmann equation. We then generalize the Chapman-Enskog method, emphasizing the role of a multiscale perturbation parameter on the collisional operator, the streaming operator, and the collisional invariants of the Boltzmann equation. The system is examined at successive orders of approximation, each of which corresponding to a physical time scale. The multicomponent Navier-Stokes regime is reached for the heavy particles, which follow a hyperbolic scaling, and is coupled to first order drift-diffusion equations for the electrons, which follow a parabolic scaling. The transport coefficients exhibit an anisotropic behavior when the magnetic field is strong enough. We also give a complete description of the Kolesnikov effect, i.e., the crossed contributions to the mass and energy transport fluxes coupling the electrons and heavy particles. Finally, the first and second principles of thermodynamics are proved to be satisfied by deriving a total energy equation and an entropy equation. Moreover, the system of equations is shown to be conservative and the purely convective system hyperbolic, thus leading to a well-defined structure.