Direct SN Treatment of the Scattering Source with Phase Functions
In deterministic transport, the scattering operator is commonly expressed using a spherical harmonics (pn) expansion of the scattering source. Not only does this require minimal modifications if the problem dimensionality is changed, but it also contains little to no approximation whenever the scattering anisotropy degree is small. An alternative method relying on a discrete ordinates (sn) approach -- essentially computing the angular integral of the phase function with a finite number of weighted directions -- is theoretically possible and has shown some advantages. In particular, this method was introduced for neutron elastic scattering and was shown to be more accurate for highly anisotropic group-to-group transfer cross-sections. Benefits were further shown in the case of a collimated beam in highly forward peaked scattering media, particularly reducing the memory usage by storing an angle-to-angle collision matrix rather than a large number of angular moments. In addition, the scattering source was expressed using that same method, the phase function having been previously expanded using Legendre polynomials, in the context of achieving extreme precision while solving the 1D radiative transfer equation. Based on these previous works, we aim at exploring two questions. First, we wish to derive an expression of the phase function in reduced geometries that does not require any more angular quadrature points than for a standard sn calculation. This is particularly relevant in 2-D problems where the reduced phase function does not solely depend on the product of the incoming and outgoing directions. Second, we wish to better understand the relationship between the sn and pn treatments. We in particular show equivalence if the phase function in the former case is expanded Legendre polynomials. This equivalence is verified numerically and the computational costs of the methods are compared on a test problem involving a Dirac delta phase function.
- Research Organization:
- Idaho National Lab. (INL), Idaho Falls, ID (United States)
- Sponsoring Organization:
- USDOE
- DOE Contract Number:
- AC07-05ID14517
- OSTI ID:
- 1499689
- Report Number(s):
- INL/CON-16-40750-Rev000
- Resource Relation:
- Conference: 2017 ANS Annual Meeting, San Francisco, CA, 06/11/2017 - 06/15/2017
- Country of Publication:
- United States
- Language:
- English
Similar Records
Discontinouous finite element equivalent to nonstandard S{sub N} methods for slab geometry transport problems
Why are the P{sub n} and S{sub n} Methods Equivalent