Numerical modelling of the Earth's ionosphere F region

International Conference on Computational Information Technologies for Environmental Sciences 2017

生态环境科学;计算机科学

Ostanin, P.A.^1 ; Kulyamin, D.V.^2 ; Dymnikov, V.P.^3

Moscow Institute of Physics and Technology (MIPT SU), Institutskiy Per., 9, Dolgoprudny, Moscow Region

141701, Russia^1

Lomonosov Moscow State University, Research Computing Center (RCC), Leninskye gory, 1, Moscow

118991, Russia^2

Institute of Numerical Mathematics, Russian Academy of Sciences (INM RAS), Gubkin str., 8, Moscow

118333, Russia^3

关键词: Ambipolar diffusion;    Different structure;    Electron density profiles;    Neutral temperature;    Numerical experiments;    Numerical solution;    Plasma chemistries;    Two-dimensional approximation;   

Others  :  https://iopscience.iop.org/article/10.1088/1755-1315/96/1/012011/pdf
DOI  :  10.1088/1755-1315/96/1/012011

来源: IOP

【 摘 要 】

This paper presents the first version of a new INM RAS Earth's ionosphere F region dynamical model. A complete set of model equations is formulated taking into account all the key physical processes that form the global state of the ionospheric F region (plasma chemistry, ambipolar diffusion, wind transport, drift across magnetic lines). For the numerical solution, a splitting method based on the physical processes and geometric directions is proposed. The first stage of splitting in a quasi-two-dimensional approximation setting with a projection of ambipolar diffusion on the vertical direction is considered. It is numerically implemented stepwise using various difference schemes for three separate model formulations (taking into account diffusion only along the vertical direction, considering a realistic direction of diffusion along the magnetic field excluding and including a mixed derivative term). The applicability, efficiency, conservation, and monotonicity of these numerical methods are analyzed. The first numerical experiments show convergence of the numerical solution to a stationary vertical profile specific to the F region. The greatest consistency with the observed profiles is obtained in the mid-latitudes. Using the thus constructed model it is shown that the electron density profile is most sensitive to the neutral temperature and ionization level with qualitatively different structures of the corresponding modes of variability.