Ray effect mitigation for the discrete ordinates method through quadrature rotation
Abstract
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 (SN) method suffers from ray effects which result from a break in rotational symmetry from the finite set of directions chosen by SN. The spherical harmonics (PN) equations, on the other hand, preserve rotational symmetry, but can produce negative particle densities. By construction, the discrete ordinates (SN) method ensures non-negative particle densities.
In this paper we present a modified version of the SN method, the rotated SN (rSN) method. Compared to SN, 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 and interpolation steps only add 5% to 10% to the runtime of the SN method. We apply the rSN 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 rSN method yields qualitatively similar solutions to the SN 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 [4].- Publication:
-
Journal of Computational Physics
- Pub Date:
- April 2019
- DOI:
- 10.1016/j.jcp.2019.01.016
- arXiv:
- arXiv:1808.05846
- Bibcode:
- 2019JCoPh.382..105C
- Keywords:
-
- Discrete ordinates method;
- Ray effects;
- Radiative transfer;
- Quadrature;
- Mathematics - Numerical Analysis;
- 85A25;
- 70F45;
- 82D75;
- 65M70;
- 41A55
- E-Print:
- doi:10.1016/j.jcp.2019.01.016