Direct S{sub N} Treatment of the Scattering Source with Phase Functions
- Reactor Physics Analysis and Design, Idaho National Laboratory, Idaho Falls, ID (United States)
- Nuclear Methods Development, Idaho National Laboratory, Idaho Falls, ID (United States)
In deterministic transport, the angular dependency of the scattering source is described using a so-called phase function. The scattering operator is commonly expressed using a spherical harmonics (P{sub N}) expansion of this phase function. 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 (S{sub N}) 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 S N calculation. This is particularly relevant because 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 S{sub N} and P{sub N} 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.
- OSTI ID:
- 23047504
- Journal Information:
- Transactions of the American Nuclear Society, Vol. 116; Conference: 2017 Annual Meeting of the American Nuclear Society, San Francisco, CA (United States), 11-15 Jun 2017; Other Information: Country of input: France; 5 refs.; available from American Nuclear Society - ANS, 555 North Kensington Avenue, La Grange Park, IL 60526 (US); ISSN 0003-018X
- 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