An efficient numerical method for solving the Boltzmann equation in multidimensions

JOURNAL OF COMPUTATIONAL PHYSICS

卷:353

Article

Dimarco, Giacomo¹  Loubere, Raphael^2,3  Narski, Jacek^2,3  Rey, Thomas⁴ 

[1] Univ Ferrara, Dept Math & Comp Sci, Via Machiavelli 30, I-44122 Ferrara, Italy

[2] Univ Paul Sabatier, CNRS, Toulouse, France

[3] Univ Paul Sabatier, IMT, Toulouse, France

[4] Univ Lille 1, Lab Paul Painleve, Lille, France

关键词: Boltzmann equation;    Kinetic equations;    Semi-Lagrangian schemes;    Spectral schemes;    3D/3D;    GPU;    CUDA;    OpenMP;    MPI;   

DOI  :  10.1016/j.jcp.2017.10.010

来源: Elsevier

【 摘 要 】

In this paper we deal with the extension of the Fast Kinetic Scheme (FKS) (Dimarco and Loubere, 2013 [26]) originally constructed for solving the BGK equation, to the more challenging case of the Boltzmann equation. The scheme combines a robust and fast method for treating the transport part based on an innovative Lagrangian technique supplemented with conservative fast spectral schemes to treat the collisional operator by means of an operator splitting approach. This approach along with several implementation features related to the parallelization of the algorithm permits to construct an efficient simulation tool which is numerically tested against exact and reference solutions on classical problems arising in rarefied gas dynamic. We present results up to the 3D x 3D case for unsteady flows for the Variable Hard Sphere model which may serve as benchmark for future comparisons between different numerical methods for solving the multidimensional Boltzmann equation. For this reason, we also provide for each problem studied details on the computational cost and memory consumption as well as comparisons with the BGK model or the limit model of compressible Euler equations. (C) 2017 Elsevier Inc. All rights reserved.

【 授权许可】

Free