Positivity-preserving scheme for two-dimensional advection–diffusion equations including mixed derivatives

Abstract

In this work, we propose a positivity-preserving scheme for solving two-dimensional advection-diffusion equations including mixed derivative terms, in order to improve the accuracy of lower-order methods. The solution to these equations, in the absence of mixed derivatives, has been studied in detail, while positivity-preserving solutions to mixed derivative terms have received much less attention. A two-dimensional diffusion equation, for which the analytical solution is known, is solved numerically to show the applicability of the scheme. It is further applied to the Fokker-Planck collision operator in two-dimensional cylindrical coordinates under the assumption of local thermal equilibrium. For a thermal equilibration problem, it is shown that the scheme conserves particle number and energy, while the preservation of positivity is ensured and the steady-state solution is the Maxwellian distribution.

A fully implicit, scalable, conservative nonlinear relativistic Fokker-Planck 0D-2P solver for runaway electron

“…The adaptive spline technique results in a small additional exponent in the algorithmic scaling of O(N 0.1 ). Positivity of the distribution function is ensured using a continuum-based reformulation approach [18] combined with robust positivity-preserving discretizations schemes [17]. Using an Anderson Acceleration fixed-point iteration scheme for our nonlinear solves, also preconditioned with multigrid techniques, we obtain an algorithm that overall scales as O(N 1.1 log N ).…”

“…To address this issue, we reformulate the off-diagonal components as effective friction forces as proposed in Ref. [18]:…”

MALTA: A Zonally Averaged Global Atmospheric Transport Model for Long‐Lived Trace Gases

“…Instead, we adapt the approach proposed by du Toit et al. (2018) to ensure positive definite tracer concentrations following the diffusion terms containing a mixed derivative. The diagonal diffusion terms (using D zz and D yy ) are always computed using the Eulerian scheme first described in this section.…”

DREAM: a fluid-kinetic framework for tokamak disruption runaway electron simulations

“…(58) An alternative to the above manner for discretizing off-diagonal diffusion terms was given in [61] where the terms were rewritten as advection terms and combined with a flux limiter scheme to ensure the preservation of monotonicity.…”