Researcher profile

Benjamin Graille

Benjamin Graille contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

6 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.

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.