On the numerical simulation of particle dynamics in the radiation belt: 2. Procedure based on the diagonalization of the diffusion tensor

被引：12

作者：

Camporeale, E. ^{[1
]}

Delzanno, G. L. ^{[1
]}

Zaharia, S. ^{[2
]}

Koller, J. ^{[2
]}

机构：

[1] Los Alamos Natl Lab, Los Alamos, NM 87545 USA

[2] Los Alamos Natl Lab, ISR 1, Los Alamos, NM 87545 USA

来源：

JOURNAL OF GEOPHYSICAL RESEARCH-SPACE PHYSICS | 2013年 / 118卷 / 06期

关键词：

radiation belt; diffusion code; numerical simulation; QUASI-LINEAR DIFFUSION; MODEL;

D O I：

10.1002/jgra.50278

中图分类号：

P1 [天文学];

学科分类号：

0704 ;

摘要：

In this paper, we conclude the survey and comparison of different numerical methods used to solve the diffusion equation for particle dynamics in the Earth's radiation belt, initiated in Camporeale et al. (2013). Here we focus on the diagonalization procedure introduced by Albert and Young (2005) that, by performing a change of coordinates, solves the diffusion equation in a space where the mixed diffusion terms are null. We describe the diagonalization procedure and its numerical implementation, which is not as straightforward as the implementation of a traditional solver in a rectangular domain. We compare the computing times with and without the diagonalization procedure, and we conclude that this procedure is generally not advantageous from the point of view of computational efficiency.

引用

页码：3476 / 3484

页数：9

共 15 条

[1] Three-dimensional diffusion simulation of outer radiation belt electrons during the 9 October 1990 magnetic storm [J].