# American Institute of Mathematical Sciences

doi: 10.3934/krm.2021008

## Projective integration schemes for hyperbolic moment equations

 Department of Computer Science, NUMA, KU Leuven, Celestijnenlaan 200A, 3001 Leuven, Belgium

* Corresponding author: Julian Koellermeier

Received  May 2020 Revised  November 2020 Published  January 2021

In this paper, we apply projective integration methods to hyperbolic moment models of the Boltzmann equation and the BGK equation, and investigate the numerical properties of the resulting scheme. Projective integration is an explicit scheme that is tailored to problems with large spectral gaps between slow and (one or many) fast eigenvalue clusters of the model. The spectral analysis of a linearized moment model clearly shows spectral gaps and reveals the multi-scale nature of the model for which projective integration is a matching choice. The combination of the non-intrusive projective integration method with moment models allows for accurate, but efficient simulations with significant speedup, as demonstrated using several 1D and 2D test cases with different collision terms, collision frequencies and relaxation times.

Citation: Julian Koellermeier, Giovanni Samaey. Projective integration schemes for hyperbolic moment equations. Kinetic & Related Models, doi: 10.3934/krm.2021008
##### References:

show all references

##### References:
Stability domain (37) of PFE method for $\delta t = 10^{-4}$, $\Delta t = 10^{-3}$, $K = 1$
Increasing spectral gap in eigenvalue spectra of HSM4 for constant collision frequency $\nu = 1$ and varying $\tau$ is ideally suited for the application of projective integration
Additional intermediate cluster in eigenvalue spectra of HSM4 for piecewise constant collision frequency $\nu \in \{0.1,1\}$ and varying $\tau$ is ideally suited for the application of a two-level telescopic projective integration
Extended fast eigenvalue spectra of HSM4 for space-dependent collision frequency $\nu = \rho(x) \in [1,7]$ and varying $\tau$ requires a connected stability region of a TPI method
Shock tube for constant collision frequency $\nu = 1$ and varying $\tau$
Shock tube for space-dependent collision frequency $\nu = \rho(x)$ and varying $\tau$
Model comparison for shock tube using space-dependent collision frequency $\nu = \rho(x)$ and varying $\tau$
Two-beam test for QBME, constant collision frequency $\nu = 1$, and varying $\tau$
Two-beam test for QBME, piecewise constant collision frequency $\nu \in \{0.01,1\}$, third order, and varying $\tau$
Computational domain for the forward facing step test case, taken from [41]
Forward facing step for QBME, BGK collision operator, and varying $\tau$
Forward facing step for QBME, Boltzmann collision operator, and varying $\tau$
Maximum extrapolation factors $N+K+1$ for connected stability region depending on $K$ according to [48]
 $K$ 1 2 3 4 5 6 7 $N+K+1$ 4 6 10.66 13.32 18.21 21.24 26.21 $N$ 2 3 6.66 8.32 12.21 14.24 18.21
 $K$ 1 2 3 4 5 6 7 $N+K+1$ 4 6 10.66 13.32 18.21 21.24 26.21 $N$ 2 3 6.66 8.32 12.21 14.24 18.21
Stability of different parameter settings for PFE. QBME model, $\nu = 1, \tau = 10^{-5}$. Base parameters $K = 1$, $\delta = 1 \cdot 10^{-5}$. Parameters predicted by linear stability analysis indicated by gray column. Instable simulation indicated by red numbers
 $\delta t / 10^{-5}$ $1.5$ $1.1$ $1$ $0.9$ $0.5$
 $\delta t / 10^{-5}$ $1.5$ $1.1$ $1$ $0.9$ $0.5$
Stability of different parameter settings for TPFE. QBME model, $\nu = \rho(x), \tau = 10^{-5}$. Each line changes only one parameter from the chosen parameters predicted by linear stability analysis $K = 6$, $\delta_0 = 1.4 \cdot 10^{-6}$, $\delta_1 = 3 \cdot 10^{-5}$ indicated by gray column. Instable simulation indicated by {red} numbers
 $K$ $8$ $7$ $6$ $5$ $4$ $\delta t_0 / 10^{-6}$ $2.5$ $2$ $1.4$ $1.3$ $1.2$ $\delta t_1 / 10^{-5}$ $5$ $4$ $3$ $2$ $1$
 $K$ $8$ $7$ $6$ $5$ $4$ $\delta t_0 / 10^{-6}$ $2.5$ $2$ $1.4$ $1.3$ $1.2$ $\delta t_1 / 10^{-5}$ $5$ $4$ $3$ $2$ $1$
Speedup of (T)PI schemes in comparison to standard FE scheme
 relaxation time $\tau$ $10^{-2}$ $10^{-3}$ $10^{-4}$ $10^{-5}$ $10^{-6}$ shock tube $\nu = 1$ 1 1 1.925 19.25 192.25 shock tube $\nu = \rho$ 2 1.375 3.93 5.61 8.02 two-beam $\nu = 1$ 1 1 1.925 19.25 192.25 two-beam $\nu = \nu_i$ 1 1 1.925 9.625 96.25 forward facing step 2 2.5 6.25 15.63 39.06
 relaxation time $\tau$ $10^{-2}$ $10^{-3}$ $10^{-4}$ $10^{-5}$ $10^{-6}$ shock tube $\nu = 1$ 1 1 1.925 19.25 192.25 shock tube $\nu = \rho$ 2 1.375 3.93 5.61 8.02 two-beam $\nu = 1$ 1 1 1.925 19.25 192.25 two-beam $\nu = \nu_i$ 1 1 1.925 9.625 96.25 forward facing step 2 2.5 6.25 15.63 39.06
 [1] Alina Chertock, Alexander Kurganov, Mária Lukáčová-Medvi${\rm{\check{d}}}$ová, Șeyma Nur Özcan. An asymptotic preserving scheme for kinetic chemotaxis models in two space dimensions. Kinetic & Related Models, 2019, 12 (1) : 195-216. doi: 10.3934/krm.2019009 [2] M. Grasselli, V. Pata. Asymptotic behavior of a parabolic-hyperbolic system. Communications on Pure & Applied Analysis, 2004, 3 (4) : 849-881. doi: 10.3934/cpaa.2004.3.849 [3] Ronald E. Mickens. Positivity preserving discrete model for the coupled ODE's modeling glycolysis. Conference Publications, 2003, 2003 (Special) : 623-629. doi: 10.3934/proc.2003.2003.623 [4] Kin Ming Hui, Soojung Kim. Asymptotic large time behavior of singular solutions of the fast diffusion equation. Discrete & Continuous Dynamical Systems - A, 2017, 37 (11) : 5943-5977. doi: 10.3934/dcds.2017258 [5] Ka Luen Cheung, Man Chun Leung. Asymptotic behavior of positive solutions of the equation $\Delta u + K u^{\frac{n+2}{n-2}} = 0$ in $IR^n$ and positive scalar curvature. Conference Publications, 2001, 2001 (Special) : 109-120. doi: 10.3934/proc.2001.2001.109 [6] Guillaume Bal, Wenjia Jing. Homogenization and corrector theory for linear transport in random media. Discrete & Continuous Dynamical Systems - A, 2010, 28 (4) : 1311-1343. doi: 10.3934/dcds.2010.28.1311 [7] Alexander A. Davydov, Massimo Giulietti, Stefano Marcugini, Fernanda Pambianco. Linear nonbinary covering codes and saturating sets in projective spaces. Advances in Mathematics of Communications, 2011, 5 (1) : 119-147. doi: 10.3934/amc.2011.5.119 [8] W. Cary Huffman. On the theory of $\mathbb{F}_q$-linear $\mathbb{F}_{q^t}$-codes. Advances in Mathematics of Communications, 2013, 7 (3) : 349-378. doi: 10.3934/amc.2013.7.349 [9] Graziano Crasta, Philippe G. LeFloch. Existence result for a class of nonconservative and nonstrictly hyperbolic systems. Communications on Pure & Applied Analysis, 2002, 1 (4) : 513-530. doi: 10.3934/cpaa.2002.1.513 [10] Teddy Pichard. A moment closure based on a projection on the boundary of the realizability domain: 1D case. Kinetic & Related Models, 2020, 13 (6) : 1243-1280. doi: 10.3934/krm.2020045 [11] Petra Csomós, Hermann Mena. Fourier-splitting method for solving hyperbolic LQR problems. Numerical Algebra, Control & Optimization, 2018, 8 (1) : 17-46. doi: 10.3934/naco.2018002 [12] Lakmi Niwanthi Wadippuli, Ivan Gudoshnikov, Oleg Makarenkov. Global asymptotic stability of nonconvex sweeping processes. Discrete & Continuous Dynamical Systems - B, 2020, 25 (3) : 1129-1139. doi: 10.3934/dcdsb.2019212 [13] Vladimir Georgiev, Sandra Lucente. Focusing nlkg equation with singular potential. Communications on Pure & Applied Analysis, 2018, 17 (4) : 1387-1406. doi: 10.3934/cpaa.2018068 [14] Daoyin He, Ingo Witt, Huicheng Yin. On the strauss index of semilinear tricomi equation. Communications on Pure & Applied Analysis, 2020, 19 (10) : 4817-4838. doi: 10.3934/cpaa.2020213 [15] Vassili Gelfreich, Carles Simó. High-precision computations of divergent asymptotic series and homoclinic phenomena. Discrete & Continuous Dynamical Systems - B, 2008, 10 (2&3, September) : 511-536. doi: 10.3934/dcdsb.2008.10.511 [16] Marian Gidea, Rafael de la Llave, Tere M. Seara. A general mechanism of instability in Hamiltonian systems: Skipping along a normally hyperbolic invariant manifold. Discrete & Continuous Dynamical Systems - A, 2020, 40 (12) : 6795-6813. doi: 10.3934/dcds.2020166 [17] Diana Keller. Optimal control of a linear stochastic Schrödinger equation. Conference Publications, 2013, 2013 (special) : 437-446. doi: 10.3934/proc.2013.2013.437 [18] Simone Cacace, Maurizio Falcone. A dynamic domain decomposition for the eikonal-diffusion equation. Discrete & Continuous Dynamical Systems - S, 2016, 9 (1) : 109-123. doi: 10.3934/dcdss.2016.9.109 [19] Naeem M. H. Alkoumi, Pedro J. Torres. Estimates on the number of limit cycles of a generalized Abel equation. Discrete & Continuous Dynamical Systems - A, 2011, 31 (1) : 25-34. doi: 10.3934/dcds.2011.31.25 [20] Wentao Huang, Jianlin Xiang. Soliton solutions for a quasilinear Schrödinger equation with critical exponent. Communications on Pure & Applied Analysis, 2016, 15 (4) : 1309-1333. doi: 10.3934/cpaa.2016.15.1309

2019 Impact Factor: 1.311

## Tools

Article outline

Figures and Tables