1,151 research outputs found
A Fast Semi-implicit Method for Anisotropic Diffusion
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.Comment: accepted in the Journal of Computational Physics; 13 pages, 7
figures; updated to the accepted versio
The cutoff method for the numerical computation of nonnegative solutions of parabolic PDEs with application to anisotropic diffusion and lubrication-type equations
The cutoff method, which cuts off the values of a function less than a given
number, is studied for the numerical computation of nonnegative solutions of
parabolic partial differential equations. A convergence analysis is given for a
broad class of finite difference methods combined with cutoff for linear
parabolic equations. Two applications are investigated, linear anisotropic
diffusion problems satisfying the setting of the convergence analysis and
nonlinear lubrication-type equations for which it is unclear if the convergence
analysis applies. The numerical results are shown to be consistent with the
theory and in good agreement with existing results in the literature. The
convergence analysis and applications demonstrate that the cutoff method is an
effective tool for use in the computation of nonnegative solutions. Cutoff can
also be used with other discretization methods such as collocation, finite
volume, finite element, and spectral methods and for the computation of
positive solutions.Comment: 19 pages, 41 figure
An implicit scheme for solving the anisotropic diffusion of heat and cosmic rays in the RAMSES code
Astrophysical plasmas are subject to a tight connection between magnetic
fields and the diffusion of particles, which leads to an anisotropic transport
of energy. Under the fluid assumption, this effect can be reduced to an
advection-diffusion equation augmenting the equations of magnetohydrodynamics.
We introduce a new method for solving the anisotropic diffusion equation using
an implicit finite-volume method with adaptive mesh refinement and adaptive
time-stepping in the RAMSES code. We apply this numerical solver to the
diffusion of cosmic ray energy, and diffusion of heat carried by electrons,
which couple to the ion temperature. We test this new implementation against
several numerical experiments and apply it to a simple supernova explosion with
a uniform magnetic field.Comment: 11 pages, 10 figures, A&
Rank-preserving geometric means of positive semi-definite matrices
The generalization of the geometric mean of positive scalars to positive
definite matrices has attracted considerable attention since the seminal work
of Ando. The paper generalizes this framework of matrix means by proposing the
definition of a rank-preserving mean for two or an arbitrary number of positive
semi-definite matrices of fixed rank. The proposed mean is shown to be
geometric in that it satisfies all the expected properties of a rank-preserving
geometric mean. The work is motivated by operations on low-rank approximations
of positive definite matrices in high-dimensional spaces.Comment: To appear in Linear Algebra and its Application
- …