Kerstin Küpper

Kerstin Küpper, Martin Frank, Shi Jin

We propose a two-dimensional asymptotic preserving scheme for linear transport equations with diffusive scalings. It is an extension of the time splitting developed by Jin, Pareschi and Toscani [SINUM,2000], but uses spatial discretizations on staggered grids, which preserves the discrete diffusion limit with a more compact stencil. The first novelty of this paper is that we propose a staggering in two dimensions that requires fewer unknowns than one could have naively expected. The second contribution of this paper is that we rigorously analyze the scheme of Jin, Pareschi, and Toscani [SINUM,2000] We show that the scheme is AP and obtain an explicit CFL condition, which couples a hyperbolic and a parabolic condition. This type of condition is common for asymptotic preserving schemes and guarantees uniform stability with respect to the mean free path. In addition, we obtain an upper bound on the relaxation parameter, which is the crucial parameter of the used time discretization. Several numerical examples are provided to verify the accuracy and asymptotic property of the scheme.

Thomas Camminady, Martin Frank, Kerstin Küpper et al.

Solving the radiation transport equation is a challenging task, due to the high dimensionality of the solution's phase space. The commonly used discrete ordinates (S$_N$) method suffers from ray effects which result from a break in rotational symmetry from the finite set of directions chosen by S$_N$. The spherical harmonics (P$_N$) equations, on the other hand, preserve rotational symmetry, but can produce negative particle densities. The discrete ordinates (S$_N$) method, in turn, by construction ensures non-negative particle densities. In this paper we present a modified version of the S$_N$ method, the rotated S$_N$ (rS$_N$) method. Compared to S$_N$, we add a rotation and interpolation step for the angular quadrature points and the respective function values after every time step. Thereby, the number of directions on which the solution evolves is effectively increased and ray effects are mitigated. Solution values on rotated ordinates are computed by an interpolation step. Implementation details are provided and in our experiments the rotation/interpolation step only adds 5% to 10% to the runtime of the S$_N$ method. We apply the rS$_N$ method to the line-source and a lattice test case, both being prone to ray-effects. Ray effects are reduced significantly, even for small numbers of quadrature points. The rS$_N$ method yields qualitatively similar solutions to the S$_N$ method with less than a third of the number of quadrature points, both for the line-source and the lattice problem. The code used to produce our results is freely available and can be downloaded.