Quadratic Finite Element Method for 1D Deterministic Transport
In the discrete ordinates, or SN, numerical solution of the transport equation, both the spatial ({und r}) and angular ({und {Omega}}) dependences on the angular flux {psi}{und r},{und {Omega}}are modeled discretely. While significant effort has been devoted toward improving the spatial discretization of the angular flux, we focus on improving the angular discretization of {psi}{und r},{und {Omega}}. Specifically, we employ a Petrov-Galerkin quadratic finite element approximation for the differencing of the angular variable ({mu}) in developing the one-dimensional (1D) spherical geometry S{sub N} equations. We develop an algorithm that shows faster convergence with angular resolution than conventional S{sub N} algorithms.
- Research Organization:
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
- Sponsoring Organization:
- USDOE
- DOE Contract Number:
- W-7405-ENG-48
- OSTI ID:
- 15013672
- Report Number(s):
- UCRL-CONF-201715; TRN: US0801236
- Resource Relation:
- Conference: Presented at: 2004 American Nuclear Society Annual Meeting, Pittsburgh, PA, United States, Jun 13 - Jun 17, 2004
- Country of Publication:
- United States
- Language:
- English
Similar Records
A comparison of angular difference schemes for one-dimensional spherical geometry S{sub N} equations
Nonnegative methods for bilinear discontinuous differencing of the SN equations on quadrilaterals