# 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] Mohammad Parto-haghighi, Mustafa Inc, Mustafa Bayram, Dumitru Baleanu. Solving the time-fractional Schrodinger equation with the group preserving scheme. Discrete & Continuous Dynamical Systems - S, 2018, 0 (0) : 0-0. doi: 10.3934/dcdss.2020146 [8] 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 [9] Mikaela Iacobelli. Asymptotic analysis for a very fast diffusion equation arising from the 1D quantization problem. Discrete & Continuous Dynamical Systems - A, 2019, 39 (9) : 4929-4943. doi: 10.3934/dcds.2019201 [10] Andrea Signori. Optimal treatment for a phase field system of Cahn-Hilliard type modeling tumor growth by asymptotic scheme. Mathematical Control & Related Fields, 2019, 0 (0) : 0-0. doi: 10.3934/mcrf.2019040 [11] 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 - A, 2015, 35 (4) : 1391-1407. doi: 10.3934/dcds.2015.35.1391 [12] Luis Caffarelli, Juan-Luis Vázquez. Asymptotic behaviour of a porous medium equation with fractional diffusion. Discrete & Continuous Dynamical Systems - A, 2011, 29 (4) : 1393-1404. doi: 10.3934/dcds.2011.29.1393 [13] 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 [14] 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 [15] Maria do Carmo Pacheco de Toledo, Sergio Muniz Oliva. A discretization scheme for an one-dimensional reaction-diffusion equation with delay and its dynamics. Discrete & Continuous Dynamical Systems - A, 2009, 23 (3) : 1041-1060. doi: 10.3934/dcds.2009.23.1041 [16] Roland Zweimüller. Asymptotic orbit complexity of infinite measure preserving transformations. Discrete & Continuous Dynamical Systems - A, 2006, 15 (1) : 353-366. doi: 10.3934/dcds.2006.15.353 [17] 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 [18] Genni Fragnelli, A. Idrissi, L. Maniar. The asymptotic behavior of a population equation with diffusion and delayed birth process. Discrete & Continuous Dynamical Systems - B, 2007, 7 (4) : 735-754. doi: 10.3934/dcdsb.2007.7.735 [19] Bhargav Kumar Kakumani, Suman Kumar Tumuluri. Asymptotic behavior of the solution of a diffusion equation with nonlocal boundary conditions. Discrete & Continuous Dynamical Systems - B, 2017, 22 (2) : 407-419. doi: 10.3934/dcdsb.2017019 [20] Abdelaziz Rhandi, Roland Schnaubelt. Asymptotic behaviour of a non-autonomous population equation with diffusion in $L^1$. Discrete & Continuous Dynamical Systems - A, 1999, 5 (3) : 663-683. doi: 10.3934/dcds.1999.5.663

2018 Impact Factor: 1.38