Browsing by Author "Wang, C."

Now showing 1 - 2 of 2
  • Thumbnail Image
    Item
    An efficient and positivity-preserving layer method for modeling radiation belt diffusion processes
    (Wiley, 2016) Tao, X.; Zhang, L.; Wang, C.; Li, X.; Albert, J.M.; Chan, A.A.
    An efficient and positivity-preserving layer method is introduced to solve the radiation belt diffusion equation and is applied to study the bounce resonance interaction between relativistic electrons and magnetosonic waves. The layer method with linear interpolation, denoted by LM-L (layer method-linear), requires the use of a large number of grid points to ensure accurate solutions. We introduce a monotonicity- and positivity-preserving cubic interpolation method to be used with the Milstein-Tretyakov layer method. The resulting method, called LM-MC (layer method-monotone cubic), can be used to solve the radiation belt diffusion equation with a much smaller number of grid points than LM-L while still being able to preserve the positivity of the solution. We suggest that LM-MC can be used to study long-term dynamics of radiation belts. We then develop a 2-D LM-MC code and use it to investigate the bounce resonance diffusion of radiation belt electrons by magnetosonic waves. Using a previously published magnetosonic wave model, we demonstrate that bounce resonance with magnetosonic waves is as important as gyroresonance; both can cause several orders of magnitude increase of MeV electron fluxes within 1ﾠday. We conclude that bounce resonance with magnetosonic waves should be taken into consideration together with gyroresonance.
  • Thumbnail Image
    Item
    Modeling radiation belt dynamics using a 3-D layer method code
    (Wiley, 2017) Wang, C.; Ma, Q.; Tao, X.; Zhang, Y.; Teng, S.; Albert, J.M.; Chan, A.A.; Li, W.; Ni, B.; Lu, Q.; Wang, S.
    A new 3-D diffusion code using a recently published layer method has been developed to analyze radiation belt electron dynamics. The code guarantees the positivity of the solution even when mixed diffusion terms are included. Unlike most of the previous codes, our 3-D code is developed directly in equatorial pitch angle (α0), momentum (p), and L shell coordinates; this eliminates the need to transform back and forth between (α0,p) coordinates and adiabatic invariant coordinates. Using (α0,p,L) is also convenient for direct comparison with satellite data. The new code has been validated by various numerical tests, and we apply the 3-D code to model the rapid electron flux enhancement following the geomagnetic storm on 17 March 2013, which is one of the Geospace Environment Modeling Focus Group challenge events. An event-specific global chorus wave model, an AL-dependent statistical plasmaspheric hiss wave model, and a recently published radial diffusion coefficient formula from Time History of Events and Macroscale Interactions during Substorms (THEMIS) statistics are used. The simulation results show good agreement with satellite observations, in general, supporting the scenario that the rapid enhancement of radiation belt electron flux for this event results from an increased level of the seed population by radial diffusion, with subsequent acceleration by chorus waves. Our results prove that the layer method can be readily used to model global radiation belt dynamics in three dimensions.