Multiple-Relaxation-Time Lattice Boltzmann scheme for fractional advection–diffusion equation

Alain Cartalade (1) , Amina Younsi (1, 2) , Marie-Christine Neel (3)

1 DM2S - Département de Modélisation des Systèmes et Structures
2 FRAMATOME-ANP
3 EMMAH - Environnement Méditerranéen et Modélisation des Agro-Hydrosystèmes

Partial differential equations (p.d.e) equipped with spatial derivatives of fractional order capture anomalous transport behaviors observed in diverse fields of Science. A number of numerical methods approximate their solutions in dimension one. Focusing our effort on such p.d.e. in higher dimension with Dirichlet boundary conditions, we present an approximation based on Lattice Boltzmann Method with Bhatnagar-Gross-Krook (BGK) or Multiple-Relaxation-Time (MRT) collision operators. First, an equilibrium distribution function is defined for simulating space-fractional diffusion equations in dimensions 2 and 3. Then, we check the accuracy of the solutions by comparing with i) random walks derived from stable Lévy motion, and ii) exact solutions. Because of its additional freedom degrees, the MRT collision operator provides accurate approximations to space-fractional advection-diffusion equations, even in the cases which the BGK fails to represent because of anisotropic diffusion tensor or of flow rate destabilizing the BGK LBM scheme.

Mots clés: Lattice Boltzmann method; Multiple-Relaxation-Time; Random Walk; Stable Process; Fractional Advection-Diffusion Equation

Domaines: Equations aux dérivées partielles [math.AP]; Physique Numérique [physics.comp-ph]

Alain Cartalade, Amina Younsi, Marie-Christine Neel. Multiple-Relaxation-Time Lattice Boltzmann scheme for fractional advection–diffusion equation. Computer Physics Communications, 2019, 234, pp.40 - 54. ⟨10.1016/j.cpc.2018.08.005⟩. ⟨cea-01900173⟩