GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems
SIAM Journal on Scientific and Statistical Computing
Scale-Space and Edge Detection Using Anisotropic Diffusion
IEEE Transactions on Pattern Analysis and Machine Intelligence
ScaLAPACK user's guide
A Technique for Accelerating the Convergence of Restarted GMRES
SIAM Journal on Matrix Analysis and Applications
Monotonicity of control volume methods
Numerische Mathematik
Preserving monotonicity in anisotropic diffusion
Journal of Computational Physics
Interpolation-free monotone finite volume method for diffusion equations on polygonal meshes
Journal of Computational Physics
Monotone Finite Volume Schemes of Nonequilibrium Radiation Diffusion Equations on Distorted Meshes
SIAM Journal on Scientific Computing
SIAM Journal on Scientific Computing
SIAM Journal on Scientific Computing
Asymptotic-preserving scheme for highly anisotropic non-linear diffusion equations
Journal of Computational Physics
Hi-index | 31.45 |
Simple finite differencing of the anisotropic diffusion equation, where diffusion is only along a given direction, does not ensure that the numerically calculated heat fluxes are in the correct direction. This can lead to negative temperatures for the anisotropic thermal diffusion equation. In a previous paper we proposed a monotonicity-preserving explicit method which uses limiters (analogous to those used in the solution of hyperbolic equations) to interpolate the temperature gradients at cell faces. However, being explicit, this method was limited by a restrictive Courant-Friedrichs-Lewy (CFL) stability timestep. Here we propose a fast, conservative, directionally-split, semi-implicit method which is second order accurate in space, is stable for large timesteps, and is easy to implement in parallel. Although not strictly monotonicity-preserving, our method gives only small amplitude temperature oscillations at large temperature gradients, and the oscillations are damped in time. With numerical experiments we show that our semi-implicit method can achieve large speed-ups compared to the explicit method, without seriously violating the monotonicity constraint. This method can also be applied to isotropic diffusion, both on regular and distorted meshes.