Abstract
Electron radiation belts are regions surrounding Earth, filled with highly energetic electrons and overlapping the majority of satellite orbits. Their multi-scale and rapidly evolving dynamics are modelled by the mean of a diffusion equation involving a highly anisotropic and inhomogeneous diffusion tensor. Finite difference based methods have been the preferred method of discretization in physical codes, starting from ONERA’s Salammbô-Electron model, the pioneering 3-dimensional code in the radiation belts community. This choice however does not prevent several numerically induced constraints impacting the reliability of the code as well as its computational cost. Thus in this paper we present the outcome of our investigation to improve Salammbô’s numerical core. In particular, we present our special diffusion frame and its numerically induced limitations on our finite difference based scheme. Then we test potential alternative finite volume schemes with physically relevant properties and we finally highlight with a real case simulation the contribution of the positivity preserving scheme to evaluate the impact of cross diffusion.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.