Numerical analysis
Multiscale numerical schemes for kinetic equations in the anomalous diffusion limit
[Schémas numériques multi-échelles pour les équations cinétiques dans la limite de diffusion anormale]
Comptes Rendus. Mathématique, Tome 353 (2015) no. 8, pp. 755-760.

Nous construisons des schémas numériques pour résoudre les équations cinétiques dans le régime de diffusion anormale. Lorsque l'équilibre présente une queue lourde ou lorsque la fréquence de collision dégénère pour les petites vitesses, un scaling approprié permet d'obtenir un modèle asymptotique appelé modèle de diffusion anormale ou fractionnaire. Le premier schéma que nous construisons est basé sur une décomposition micro–macro de la fonction de distribution, tandis que le second s'appuie sur une formulation de Duhamel de l'équation de départ. Ces deux schémas sont Asymptotic Preserving (AP) : ils sont consistants avec l'équation cinétique lorsque le paramètre d'échelle ε>0 est fixé et dégénèrent en un schéma consistant avec le modèle limite quand ε tend vers 0. Le deuxième schéma est même uniformément précis (UA) par rapport à ε. Les schémas AP qui sont connus dans le cas de la limite de diffusion classique ne peuvent pas directement s'appliquer au cas de la diffusion anormale, car ils ne permettent pas de capturer les effets importants des petites et des grandes vitesses. Nous présentons des tests numériques pour mettre en évidence l'efficacité des schémas que nous présentons.

We construct numerical schemes to solve kinetic equations with anomalous diffusion scaling. When the equilibrium is heavy-tailed or when the collision frequency degenerates for small velocities, an appropriate scaling should be made and the limit model is the so-called anomalous or fractional diffusion model. Our first scheme is based on a suitable micro–macro decomposition of the distribution function, whereas our second scheme relies on a Duhamel formulation of the kinetic equation. Both are Asymptotic Preserving (AP): they are consistent with the kinetic equation for all fixed value of the scaling parameter ε>0 and degenerate into a consistent scheme solving the asymptotic model when ε tends to 0. The second scheme enjoys the stronger property of being uniformly accurate (UA) with respect to ε. The usual AP schemes known for the classical diffusion limit cannot be directly applied to the context of anomalous diffusion scaling, since they are not able to capture the important effects of large and small velocities. We present numerical tests to highlight the efficiency of our schemes.

Reçu le :
Accepté le :
Publié le :
DOI : 10.1016/j.crma.2015.05.003
Crouseilles, Nicolas 1, 2 ; Hivert, Hélène 2 ; Lemou, Mohammed 2, 3

1 IPSO, INRIA, 263, avenue du Général-Leclerc, 35000 Rennes, France
2 IRMAR, Université de Rennes-1, campus de Beaulieu, 35000 Rennes, France
3 IRMAR, CNRS, 35000 Rennes, France
@article{CRMATH_2015__353_8_755_0,
     author = {Crouseilles, Nicolas and Hivert, H\'el\`ene and Lemou, Mohammed},
     title = {Multiscale numerical schemes for kinetic equations in the anomalous diffusion limit},
     journal = {Comptes Rendus. Math\'ematique},
     pages = {755--760},
     publisher = {Elsevier},
     volume = {353},
     number = {8},
     year = {2015},
     doi = {10.1016/j.crma.2015.05.003},
     language = {en},
     url = {http://www.numdam.org/articles/10.1016/j.crma.2015.05.003/}
}
TY  - JOUR
AU  - Crouseilles, Nicolas
AU  - Hivert, Hélène
AU  - Lemou, Mohammed
TI  - Multiscale numerical schemes for kinetic equations in the anomalous diffusion limit
JO  - Comptes Rendus. Mathématique
PY  - 2015
SP  - 755
EP  - 760
VL  - 353
IS  - 8
PB  - Elsevier
UR  - http://www.numdam.org/articles/10.1016/j.crma.2015.05.003/
DO  - 10.1016/j.crma.2015.05.003
LA  - en
ID  - CRMATH_2015__353_8_755_0
ER  - 
%0 Journal Article
%A Crouseilles, Nicolas
%A Hivert, Hélène
%A Lemou, Mohammed
%T Multiscale numerical schemes for kinetic equations in the anomalous diffusion limit
%J Comptes Rendus. Mathématique
%D 2015
%P 755-760
%V 353
%N 8
%I Elsevier
%U http://www.numdam.org/articles/10.1016/j.crma.2015.05.003/
%R 10.1016/j.crma.2015.05.003
%G en
%F CRMATH_2015__353_8_755_0
Crouseilles, Nicolas; Hivert, Hélène; Lemou, Mohammed. Multiscale numerical schemes for kinetic equations in the anomalous diffusion limit. Comptes Rendus. Mathématique, Tome 353 (2015) no. 8, pp. 755-760. doi : 10.1016/j.crma.2015.05.003. http://www.numdam.org/articles/10.1016/j.crma.2015.05.003/

[1] Ben Abdallah, N.; Mellet, A.; Puel, M. Anomalous diffusion limit for kinetic equations with degenerate collision frequency, Math. Models Methods Appl. Sci., Volume 21 (2011), pp. 2249-2262

[2] Ben Abdallah, N.; Mellet, A.; Puel, M. Fractional diffusion limit for collisional kinetic equations: a Hilbert expansion approach, Kinet. Relat. Models, Volume 4 (2011) no. 4, pp. 873-900

[3] Crouseilles, N.; Hivert, H.; Lemou, M. Numerical schemes for kinetic equations in the anomalous diffusion limit. Part I: the case of heavy-tailed equilibrium, 2015 | arXiv

[4] N. Crouseilles, H. Hivert, M. Lemou, Numerical schemes for kinetic equations in the anomalous diffusion limit. Part II: the case of degenerate collision frequency, in preparation.

[5] Lemou, M.; Mieussens, L. A new asymptotic preserving scheme based on micro–macro formulation for linear kinetic equations in the diffusion limit, SIAM J. Sci. Comput., Volume 31 (2008) no. 1, pp. 334-368

[6] Mellet, A.; Mischler, S.; Mouhot, C. Fractional diffusion limit for collisional kinetic equations, Arch. Ration. Mech. Anal., Volume 199 (2011) no. 2, pp. 493-525

Cité par Sources :