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.
Quantification of cervical elasticity during pregnancy based on transvaginal ultrasound imaging and stress measurement
Authors:
Peng Hu,
Peinan Zhao,
Yuan Qu,
Konstantin Maslov,
Jessica Chubiz,
Methodius G. Tuuli,
Molly J. Stout,
Lihong V. Wang
Abstract:
Objective: Strain elastography and shear wave elastography are two commonly used methods to quantify cervical elasticity; however, they have limitations. Strain elastography is effective in showing tissue elasticity distribution in a single image, but the absence of stress information causes difficulty in comparing the results acquired from different imaging sessions. Shear wave elastography is ef…
▽ More
Objective: Strain elastography and shear wave elastography are two commonly used methods to quantify cervical elasticity; however, they have limitations. Strain elastography is effective in showing tissue elasticity distribution in a single image, but the absence of stress information causes difficulty in comparing the results acquired from different imaging sessions. Shear wave elastography is effective in measuring shear wave speed (an intrinsic tissue property correlated with elasticity) in relatively homogeneous tissue, such as in the liver. However, for inhomogeneous tissue in the cervix, the shear wave speed measurement is less robust. To overcome these limitations, we develop a quantitative cervical elastography system by adding a stress sensor to an ultrasound imaging system.
Methods: In an imaging session for quantitative cervical elastography, we use the transvaginal ultrasound imaging system to record B-mode images of the cervix showing its deformation and use the stress sensor to record the probe-surface stress simultaneously. We develop a correlation-based automatic feature tracking algorithm to quantify the deformation, from which the strain is quantified. After each imaging session, we calibrate the stress sensor and transform its measurement to true stress. Applying a linear regression to the stress and strain, we obtain an approximation of the cervical Young's modulus.
Results: We validate the accuracy and robustness of this elastography system using phantom experiments. Applying this system to pregnant participants, we observe significant softening of the cervix during pregnancy (p-value < 0.001) with the cervical Young's modulus decreasing 3.95% per week. We estimate that geometric mean values of cervical Young's moduli during the first (11 to 13 weeks), second, and third trimesters are 13.07 kPa, 7.59 kPa, and 4.40 kPa, respectively.
△ Less
Submitted 9 March, 2023;
originally announced March 2023.