A Fast Semi-implicit Method for Anisotropic Diffusion

Abstract

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.

accepted in the Journal of Computational Physics; 13 pages, 7 figures; updated to the accepted version

References in corpus (3)

• Preserving Monotonicity in Anisotropic Diffusion
• Anisotropic Thermal Conduction and the Cooling Flow Problem in Galaxy Clusters
• Numerical simulations of buoyancy instabilities in galaxy cluster plasmas with cosmic rays and anisotropic thermal conduction

Cited by in corpus (15)

• Anisotropic Diffusion in Mesh-Free Numerical Magnetohydrodynamics
• A New Numerical Scheme for Cosmic Ray Transport
• Semi-implicit anisotropic cosmic ray transport on an unstructured moving mesh
• Simulations of cosmic ray propagation
• Thermal Conduction and Multiphase Gas in Cluster Cores
• Accurately simulating anisotropic thermal conduction on a moving mesh
• The impact of cosmic rays on dynamical balance and disk-halo interaction in Lstar disk galaxies
• Scalable explicit implementation of anisotropic diffusion with Runge-Kutta-Legendre super-time-stepping
• Asymptotic Preserving scheme for strongly anisotropic parabolic equations for arbitrary anisotropy direction
• A self-consistent hybrid model of kinetic striations in low-current Argon discharges
• Modeling the thermal conduction in the solar atmosphere with the code MANCHA3D
• Polarization of thermal bremsstrahlung emission due to electron pressure anisotropy
• A Simple Sub-Grid Model For Cosmic Ray Effects on Galactic Scales
• Propagation of weak shocks in cool-core galaxy clusters in two-temperature magnetohydrodynamics with anisotropic thermal conduction
• Comparing different closure relations for cosmic ray hydrodynamics