# American Institute of Mathematical Sciences

April  2018, 11(2): 409-439. doi: 10.3934/krm.2018019

## Numerical schemes for kinetic equation with diffusion limit and anomalous time scale

 ENS Lyon, UMPA UMR 5669 CNRS, and Inria Rhône-Alpes, projet NUMED, 46, allée d'Italie, 69364 Lyon Cedex 07, France

* Corresponding author

Received  October 2016 Revised  May 2017 Published  January 2018

Fund Project: The author is supported by ERC starting grant MESOPROBIO.

In this work, we propose numerical schemes for linear kinetic equation, which are able to deal with a diffusion limit and an anomalous time scale of the form ${\varepsilon ^2}\left( {1 + \left| {\ln \left( \varepsilon \right)} \right|} \right)$. When the equilibrium distribution function is a heavy-tailed function, it is known that for an appropriate time scale, the mean-free-path limit leads either to diffusion or fractional diffusion equation, depending on the tail of the equilibrium. This work deals with a critical exponent between these two cases, for which an anomalous time scale must be used to find a standard diffusion limit. Our aim is to develop numerical schemes which work for the different regimes, with no restriction on the numerical parameters. Indeed, the degeneracy $\varepsilon\to0$ makes the kinetic equation stiff. From a numerical point of view, it is necessary to construct schemes able to undertake this stiffness to avoid the increase of computational cost. In this case, it is crucial to capture numerically the effects of the large velocities of the heavy-tailed equilibrium. Moreover, we prove that the convergence towards the diffusion limit happens with two scales, the second being very slow. The schemes we propose are designed to respect this asymptotic behavior. Various numerical tests are performed to illustrate the efficiency of our methods in this context.

Citation: Hélène Hivert. Numerical schemes for kinetic equation with diffusion limit and anomalous time scale. Kinetic & Related Models, 2018, 11 (2) : 409-439. doi: 10.3934/krm.2018019
##### References:

show all references

##### References:
For $\Delta t = 10^{-2}$, the solutions of (34) at time $T = 0.1$ for different values of $\varepsilon$, when $a_\varepsilon$ is computed with (35), and the solution of the limit scheme (37)
For $\Delta t = 10^{-2}$, the solutions of (34) at time $T = 0.1$ for different values of $\varepsilon$, when $a_\varepsilon$ is computed with (36), and the solution of the limit scheme (37)
For $\Delta t = 10^{-2}$, the relative error between the solution of the scheme (34) and the limit scheme (37) at time $T = 0.1$, in function of $\varepsilon$ (log scale)
For $\Delta t = 10^{-2}$, the solutions of IS, quasi-diff and diff schemes for different values of $\varepsilon$
For $\Delta t = 10^{-2}$, the error (71) between the solution of IS and quasi-diff scheme in function of $\varepsilon$ (log scale)
For $\Delta t = 10^{-4}$ and $\varepsilon = 1$, the solutions of the MMS scheme and of the explicit scheme (69)
The relative consistency error (70) for the MMS scheme (log scale)
For $\Delta t = 10^{-4}$, the solutions of MMS, quasi-diff and diff schemes for different values of $\varepsilon$
For $\Delta t = 10^{-2}$ and $\varepsilon = 1$, the solutions of the IS scheme and of the explicit scheme (69)
The relative consistency error (70) for the IS scheme (log scale)
For $\Delta t = 10^{-4}$, the error (71) between the solution of MMS and quasi-diff scheme in function of $\varepsilon$ (log scale)
For $\Delta t = 10^{-2}$ and $\varepsilon = 1$, the solutions of the DS scheme and of the explicit scheme (69)
The relative consistency error (70) for the DS scheme (log scale)
For $\Delta t = 10^{-2}$, the solutions of DS, quasi-diff and diff schemes for different values of $\varepsilon$
For $\Delta t = 10^{-2}$, the error (71) between the solution of DS and quasi-diff scheme in function of $\varepsilon$ (log scale)
The error (70) as a function of $\varepsilon$. The density $\rho_{reference}$ is the density given by the DS scheme for $\Delta t_{ref} = 5\cdot 10^{-5}$, and $\rho_{\Delta t}$ are the densities given by the DS scheme for different values of $\Delta t$ (log scale)
 [1] Alina Chertock, Changhui Tan, Bokai Yan. An asymptotic preserving scheme for kinetic models with singular limit. Kinetic & Related Models, 2018, 11 (4) : 735-756. doi: 10.3934/krm.2018030 [2] Jingwei Hu, Shi Jin, Li Wang. An asymptotic-preserving scheme for the semiconductor Boltzmann equation with two-scale collisions: A splitting approach. Kinetic & Related Models, 2015, 8 (4) : 707-723. doi: 10.3934/krm.2015.8.707 [3] 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 [4] Nicolas Crouseilles, Mohammed Lemou, SV Raghurama Rao, Ankit Ruhi, Muddu Sekhar. Asymptotic preserving scheme for a kinetic model describing incompressible fluids. Kinetic & Related Models, 2016, 9 (1) : 51-74. doi: 10.3934/krm.2016.9.51 [5] Nicolas Crouseilles, Mohammed Lemou. An asymptotic preserving scheme based on a micro-macro decomposition for Collisional Vlasov equations: diffusion and high-field scaling limits. Kinetic & Related Models, 2011, 4 (2) : 441-477. doi: 10.3934/krm.2011.4.441 [6] Stéphane Brull, Pierre Degond, Fabrice Deluzet, Alexandre Mouton. Asymptotic-preserving scheme for a bi-fluid Euler-Lorentz model. Kinetic & Related Models, 2011, 4 (4) : 991-1023. doi: 10.3934/krm.2011.4.991 [7] Makoto Okumura, Daisuke Furihata. A structure-preserving scheme for the Allen–Cahn equation with a dynamic boundary condition. Discrete & Continuous Dynamical Systems, 2020, 40 (8) : 4927-4960. doi: 10.3934/dcds.2020206 [8] Adrian Viorel, Cristian D. Alecsa, Titus O. Pinţa. Asymptotic analysis of a structure-preserving integrator for damped Hamiltonian systems. Discrete & Continuous Dynamical Systems, 2021, 41 (7) : 3319-3341. doi: 10.3934/dcds.2020407 [9] Jisheng Kou, Huangxin Chen, Xiuhua Wang, Shuyu Sun. A linear, decoupled and positivity-preserving numerical scheme for an epidemic model with advection and diffusion. Communications on Pure & Applied Analysis, , () : -. doi: 10.3934/cpaa.2021094 [10] Tetsuya Ishiwata, Kota Kumazaki. Structure preserving finite difference scheme for the Landau-Lifshitz equation with applied magnetic field. Conference Publications, 2015, 2015 (special) : 644-651. doi: 10.3934/proc.2015.0644 [11] Makoto Okumura, Takeshi Fukao, Daisuke Furihata, Shuji Yoshikawa. A second-order accurate structure-preserving scheme for the Cahn-Hilliard equation with a dynamic boundary condition. Communications on Pure & Applied Analysis, 2022, 21 (2) : 355-392. doi: 10.3934/cpaa.2021181 [12] Mikaela Iacobelli. Asymptotic analysis for a very fast diffusion equation arising from the 1D quantization problem. Discrete & Continuous Dynamical Systems, 2019, 39 (9) : 4929-4943. doi: 10.3934/dcds.2019201 [13] Andrea Signori. Optimal treatment for a phase field system of Cahn-Hilliard type modeling tumor growth by asymptotic scheme. Mathematical Control & Related Fields, 2020, 10 (2) : 305-331. doi: 10.3934/mcrf.2019040 [14] Chuntian Wang, Yuan Zhang. A multiscale stochastic criminal behavior model under a hybrid scheme. Electronic Research Archive, 2021, 29 (4) : 2741-2753. doi: 10.3934/era.2021011 [15] Carmen Cortázar, Manuel Elgueta, Fernando Quirós, Noemí Wolanski. Asymptotic behavior for a nonlocal diffusion equation on the half line. Discrete & Continuous Dynamical Systems, 2015, 35 (4) : 1391-1407. doi: 10.3934/dcds.2015.35.1391 [16] Luis Caffarelli, Juan-Luis Vázquez. Asymptotic behaviour of a porous medium equation with fractional diffusion. Discrete & Continuous Dynamical Systems, 2011, 29 (4) : 1393-1404. doi: 10.3934/dcds.2011.29.1393 [17] Yuto Miyatake, Tai Nakagawa, Tomohiro Sogabe, Shao-Liang Zhang. A structure-preserving Fourier pseudo-spectral linearly implicit scheme for the space-fractional nonlinear Schrödinger equation. Journal of Computational Dynamics, 2019, 6 (2) : 361-383. doi: 10.3934/jcd.2019018 [18] Casimir Emako, Farah Kanbar, Christian Klingenberg, Min Tang. A criterion for asymptotic preserving schemes of kinetic equations to be uniformly stationary preserving. Kinetic & Related Models, 2021, 14 (5) : 847-866. doi: 10.3934/krm.2021026 [19] Roland Zweimüller. Asymptotic orbit complexity of infinite measure preserving transformations. Discrete & Continuous Dynamical Systems, 2006, 15 (1) : 353-366. doi: 10.3934/dcds.2006.15.353 [20] Jyoti Mishra. Analysis of the Fitzhugh Nagumo model with a new numerical scheme. Discrete & Continuous Dynamical Systems - S, 2020, 13 (3) : 781-795. doi: 10.3934/dcdss.2020044

2020 Impact Factor: 1.432