A cost-effective nonlinear extremum-preserving finite volume scheme for highly anisotropic diffusion on Cartesian grids, with application to radiation belt dynamics

by   Nour Dahmen, et al.

We construct a new nonlinear finite volume (FV) scheme for highly anisotropic diffusion equations, that satisfies the discrete minimum-maximum principle. The construction relies on the linearized scheme satisfying less restrictive monotonicity conditions than those of an M-matrix, based on a weakly regular matrix splitting and using the Cartesian structure of the mesh (extension to quadrilateral meshes is also possible). The resulting scheme, obtained by expressing fluxes as nonlinear combinations of linear fluxes, has a larger stencil than other nonlinear positivity preserving or minimum-maximum principle preserving schemes. Its larger "linearized" stencil, closer to the actual complete stencil (that includes unknowns appearing in the convex combination coefficients), enables a faster convergence of the Picard iterations used to compute the solution of the scheme. Steady state dimensionless numerical tests as well as simulations of the highly anisotropic diffusion in electron radiation belts show a second order of convergence of the new scheme and confirm its computational efficiency compared to usual nonlinear FV schemes.



page 17


Long-time behaviour of hybrid finite volume schemes for advection-diffusion equations: linear and nonlinear approaches

We are interested in the long-time behaviour of approximate solutions to...

On the square-root approximation finite volume scheme for nonlinear drift-diffusion equations

We study a finite volume scheme for the approximation of the solution to...

Maximum Principle Preserving Finite Difference Scheme for 1-D Nonlocal-to-Local Diffusion Problems

In a recent paper (see [7]), a quasi-nonlocal coupling method was introd...

Large time behavior of nonlinear finite volume schemes for convection-diffusion equations

In this contribution we analyze the large time behavior of a family of n...

On differentiable local bounds preserving stabilization for Euler equations

This work is focused on the design of nonlinear stabilization techniques...

A Generalised Complete Flux scheme for anisotropic advection-diffusion equations

In this paper, we consider separating the discretisation of the diffusiv...

Monotonicity-preserving finite element schemes with adaptive mesh refinement for hyperbolic problems

This work is focused on the extension and assessment of the monotonicity...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.