Ray effect mitigation for the discrete ordinates method through quadrature rotation

DOI10.1016/J.JCP.2019.01.016zbMATH Open1451.65116arXiv1808.05846OpenAlexW2885039257WikidataQ128504615 ScholiaQ128504615MaRDI QIDQ2214577FDOQ2214577

Authors: Thomas Camminady, Martin Frank, Kerstin Küpper, Jonas Kusch 

Publication date: 9 December 2020

Published in: Journal of Computational Physics (Search for Journal in Brave)

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 (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.

Full work available at URL: https://arxiv.org/abs/1808.05846

Recommendations

• Filtered discrete ordinates equations for radiative transport
• Quasi-singular angular finite element methods in neutron transport problems
• AnF_Nmethod for the radiative transport equation in three dimensions
• The discrete cones method for two-dimensional neutron transport calculations
• The discrete-ordinate method in diffusive regimes

zbMATH Keywords

quadratureradiative transferdiscrete ordinates methodray effects

Mathematics Subject Classification ID

Numerical interpolation (65D05) Integro-partial differential equations (45K05) Numerical methods for integral equations (65R20) Integro-partial differential equations (35R09) Finite volume methods for initial value and initial-boundary value problems involving PDEs (65M08)

Cites Work

• Title not available (Why is that?)
• An implicit Monte Carlo scheme for calculating time and frequency dependent nonlinear radiation transport
• Title not available (Why is that?)
• Title not available (Why is that?)
• A comparison of moment closures for linear kinetic transport equations: the line source benchmark
• Robust and accurate filtered spherical harmonics expansions for radiative transfer
• Two-dimensional time dependent Riemann solvers for neutron transport
• ANGULAR QUADRATURES FOR IMPROVED TRANSPORT COMPUTATIONS

Cited In (9)

• Goal-based angular adaptivity for Boltzmann transport in the presence of ray-effects
• Ray effect in rarefied flow simulation
• A modified \(S_n\) method for calculating the uncollided flux from a point source without ray effects
• KiT-RT: an extendable framework for radiative transfer and therapy
• Moment methods for the radiative transfer equation based on \(\varphi\)-divergences
• Unified gas-kinetic wave-particle methods. III: Multiscale photon transport
• MIXKIP/RAPCAL: A Computational Package for Integrated Simulations of Large-Scale Atomic Kinetics and Radiation Transport in Non-Local Thermodynamic Equilibrium Plasmas
• High order asymptotic preserving discontinuous Galerkin methods for gray radiative transfer equations
• Filtered discrete ordinates equations for radiative transport

Uses Software

• SNSPH

This page was built for publication: Ray effect mitigation for the discrete ordinates method through quadrature rotation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2214577)