Numerical model of Earth ionosphere F region based on three-dimensional transport and ambipolar diffusion equations
-
Dmitry V. Kulyamin
,
Sergey V. Kostrykin
Abstract
The paper provides a detailed description of the numerical implementation of the transport scheme in the Earth’s ionosphere three-dimensional dynamical model of the Institute of Numerical Mathematics (INMIM). The presented version of INM-IM model takes into account the global dynamical processes of ion transport and ambipolar diffusion in the altitude range between 100 and 500 km (corresponding to F region). Based upon the splitting method the model solves the equations of ambipolar diffusion on the first split step and incorporates specifically designed advective transport scheme on the second step. The accuracy of transport scheme implementation in the model has been investigated through analytical solutions. The stability of the numerical algorithm has been demonstrated even for cases when transport velocities approaching extreme values.
Funding statement: This work was carried out in Fedorov Institute of Applied Geophysics and supported by Russian Federation research and technical development program in ecological strategy and climate change through project ‘Development of an extended version of the Earth system INM RAS model within a new computational framework’, registration No. 1023082900022-7-1.5.1 (Sections 1, 2.2), the testing of transport scheme (Section 2) was carried out in INM RAS with support by Russian Science Foundation project No. 20-17-00190.
References
[1] V. P. Dymnikov, D. V. Kulyamin, and P. A. Ostanin, Coupled model of Earth’s thermosphere and ionosphere global dynamics. Izvestiya, Atmospheric and Oceanic Physics 56 (2020), No. 3, 241–252.10.1134/S0001433820030068Suche in Google Scholar
[2] G. S. Golitsyn, Principle of the fastest response in hydrodynamics, geophysics, and astrophysics. Physics. Doklady 42 (1997), No. 9, 479–482.Suche in Google Scholar
[3] V. M. Goloviznin, S. A. Karabasov, T. K. Kozubskaya, and N. V. Maksimov, CABARET scheme for the numerical solution of aeroacoustics problems: Generalization to linearized one-dimensional Euler equations. Comput. Math. Math. Phys. 49 (2009), No. 12, 2168–2182.10.1134/S096554250912015XSuche in Google Scholar
[4] O. A. Kovyrkina and V. V. Ostapenko, Monotonicity of the CABARET scheme approximating a hyperbolic system of conservation laws. Comput. Math. Math. Phys. 58 (2018), No. 9, 1435–1450.10.1134/S0965542518090129Suche in Google Scholar
[5] D. V. Kulyamin and V. P. Dymnikov, A three-dimensional model of general thermospheric circulation. Russ. J. Numer. Anal. Math. Modelling 28 (2013), No. 4, 353-380.10.1515/rnam-2013-0021Suche in Google Scholar
[6] D. V. Kulyamin, V. P. Dymnikov, and P. A. Ostanin, Modeling the F Layer of the Earth’s ionosphere: Solution of the ambipolar diffusion equations. Math. Models and Comp. Simul. 11 (2019), No. 6, 940–950.10.1134/S2070048219060115Suche in Google Scholar
[7] D. V. Kulyamin, V. P. Dymnikov, and P. A. Ostanin, INM-IM: INM RAS Earth ionosphere F region dynamical model. Russ. J. Numer. Anal. Math. Modelling 37 (2022), No. 6, 349–362.10.1515/rnam-2022-0028Suche in Google Scholar
[8] P. H. Lauritzen, W. C. Skamarock, M. J. Prather, and M. A. Taylor, A standard test case suite for two-dimensional linear transport on the sphere. Geosci. Model Dev. 5 (2012), No. 3, 887–901.10.5194/gmd-5-887-2012Suche in Google Scholar
[9] P. H. Lauritzen et al., A standard test case suite for two-dimensional linear transport on the sphere: results from a collection of state-of-the-art schemes. Geosci. Model Dev. 7 (2014), No. 1, 105–145.Suche in Google Scholar
[10] P. A. Ostanin, D. V. Kulyamin, and V. P. Dymnikov, Numerical modelling of the Earth’s ionosphere F region. IOP Conference Series: Earth and Environmental Science 96 (2017), No. 1, 012011.10.1088/1755-1315/96/1/012011Suche in Google Scholar
[11] R. W. Schunk and A. Nagy, Ionospheres: Physics, Plasma Physics, and Chemistry, 5th ed. Cambridge University Press, 2009.10.1017/CBO9780511635342Suche in Google Scholar
© 2023 Walter de Gruyter GmbH, Berlin/Boston
Artikel in diesem Heft
- Frontmatter
- Explaining breakthrough behaviour in shale rock: influence of capillary effects and geomechanics
- Non-local discretization of the isoneutral diffusion operator in a terrain-following climate ocean model
- Numerical model of Earth ionosphere F region based on three-dimensional transport and ambipolar diffusion equations
- Extracting connectivity paths in digital core images using solution of partial minimum eigenvalue problem
- The study of the local sensitivity of functionals of the optimal solution to observational data and the heat flux input data in a variational assimilation problem for the sea thermodynamics model
- Multiresolution approximation for shallow water equations using summation-by-parts finite differences
Artikel in diesem Heft
- Frontmatter
- Explaining breakthrough behaviour in shale rock: influence of capillary effects and geomechanics
- Non-local discretization of the isoneutral diffusion operator in a terrain-following climate ocean model
- Numerical model of Earth ionosphere F region based on three-dimensional transport and ambipolar diffusion equations
- Extracting connectivity paths in digital core images using solution of partial minimum eigenvalue problem
- The study of the local sensitivity of functionals of the optimal solution to observational data and the heat flux input data in a variational assimilation problem for the sea thermodynamics model
- Multiresolution approximation for shallow water equations using summation-by-parts finite differences
