Two methods to analyse radial diffusion ensembles: the peril of space- and time- dependent diffusion
Authors:
Sarah N. Bentley,
J. Stout,
Rhys Thompson,
Daniel J. Ratliff,
Clare E. J. Watt
Abstract:
Particle dynamics in Earth's outer radiation belt can be modelled using a diffusion framework, where large-scale electron movements are captured by a diffusion equation across a single adiabatic invariant, $L^{*}$ $``(L)"$. While ensemble models are promoted to represent physical uncertainty, as yet there is no validated method to analyse radiation belt ensembles. Comparisons are complicated by th…
▽ More
Particle dynamics in Earth's outer radiation belt can be modelled using a diffusion framework, where large-scale electron movements are captured by a diffusion equation across a single adiabatic invariant, $L^{*}$ $``(L)"$. While ensemble models are promoted to represent physical uncertainty, as yet there is no validated method to analyse radiation belt ensembles. Comparisons are complicated by the domain dependent diffusion, since diffusion coefficient $D_{LL}$ is dependent on $L$. We derive two tools to analyse ensemble members: time to monotonicity $t_m$ and mass/energy moment quantities $\mathcal{N}, \mathcal{E}$. We find that the Jacobian ($1/L^2$) is necessary for radiation belt error metrics. Components of $\partial\mathcal{E}/\partial t$ are explicitly calculated to compare the effects of outer and inner boundary conditions, and loss, on the ongoing diffusion. Using $t_m$, $\mathcal{N}$ and $\mathcal{E}$, we find that: (a) different physically motivated choices of outer boundary condition and location result in different final states and different rates of evolution; (b) the gradients of the particle distribution affect evolution more significantly than $D_{LL}$; (c) the enhancement location, and the amount of initial background particles, are both significant factors determining system evolution; (d) loss from pitch-angle scattering is generally dominant; it mitigates but does not remove the influence of both initial conditions and outer boundary settings, which are due to the $L$-dependence of $D_{LL}$. We anticipate this study will promote renewed focus on the distribution gradients, on the location and nature of the outer boundary in radiation belt modelling, and provide a foundation for systematic ensemble modelling.
△ Less
Submitted 5 July, 2024;
originally announced July 2024.
Reappraisal of Whitham's 1967 theory for wave-mean flow interaction in shallow water
Authors:
Thomas J. Bridges,
Daniel J. Ratliff
Abstract:
The modulation equations for Stokes waves in shallow water coupled to wave generated mean flow, derived in Whitham (1967), based on an averaged Lagrangian are revisited. Firstly, it is shown that they can be recast into two coupled classical shallow water equations, with modified gravity having the sign of the Whitham index: sign($ω_0'' ω_2$). Secondly, it is shown that the amplitude of the mean f…
▽ More
The modulation equations for Stokes waves in shallow water coupled to wave generated mean flow, derived in Whitham (1967), based on an averaged Lagrangian are revisited. Firstly, it is shown that they can be recast into two coupled classical shallow water equations, with modified gravity having the sign of the Whitham index: sign($ω_0'' ω_2$). Secondly, it is shown that the amplitude of the mean flow and amplitude of the wave are, in general, independent. Thirdly, the implications of the coalescing characteristics, whose unfolding is associated with the Benjamin-Feir instability, are studied.
△ Less
Submitted 8 March, 2022;
originally announced March 2022.
