Multiscale Capability in Rattlesnake using Contiguous Discontinuous Discretization of SelfAdjoint Angular Flux Equation
Abstract
In this report, we present a new upwinding scheme for the multiscale capability in Rattlesnake, the MOOSE based radiation transport application. Comparing with the initial implementation of multiscale utilizing Lagrange multipliers to impose strong continuity of angular flux on interface of inbetween subdomains, this scheme does not require the particular domain partitioning. This upwinding scheme introduces discontinuity of angular flux and resembles the classic upwinding technique developed for solving first order transport equation using discontinuous finite element method (DFEM) on the subdomain interfaces. Because this scheme restores the causality of radiation streaming on the interfaces, significant accuracy improvement can be observed with moderate increase of the degrees of freedom comparing with the continuous method over the entire solution domain. Hybrid SNPN is implemented and tested with this upwinding scheme. Numerical results show that the angular smoothing required by Lagrange multiplier method is not necessary for the upwinding scheme.
 Authors:
 Texas A & M Univ., College Station, TX (United States)
 Idaho National Lab. (INL), Idaho Falls, ID (United States)
 Publication Date:
 Research Org.:
 Idaho National Lab. (INL), Idaho Falls, ID (United States)
 Sponsoring Org.:
 USDOE Office of Nuclear Energy (NE)
 OSTI Identifier:
 1364491
 Report Number(s):
 INL/EXT1639793
TRN: US1703367
 DOE Contract Number:
 AC0705ID14517
 Resource Type:
 Technical Report
 Country of Publication:
 United States
 Language:
 English
 Subject:
 22 GENERAL STUDIES OF NUCLEAR REACTORS; FINITE ELEMENT METHOD; COMPARATIVE EVALUATIONS; DEGREES OF FREEDOM; RADIATION STREAMING; RADIATION TRANSPORT; multiscale; Rattlesnake; TREAT
Citation Formats
Zheng, Weixiong, Wang, Yaqi, and DeHart, Mark D. Multiscale Capability in Rattlesnake using Contiguous Discontinuous Discretization of SelfAdjoint Angular Flux Equation. United States: N. p., 2016.
Web. doi:10.2172/1364491.
Zheng, Weixiong, Wang, Yaqi, & DeHart, Mark D. Multiscale Capability in Rattlesnake using Contiguous Discontinuous Discretization of SelfAdjoint Angular Flux Equation. United States. doi:10.2172/1364491.
Zheng, Weixiong, Wang, Yaqi, and DeHart, Mark D. Thu .
"Multiscale Capability in Rattlesnake using Contiguous Discontinuous Discretization of SelfAdjoint Angular Flux Equation". United States.
doi:10.2172/1364491. https://www.osti.gov/servlets/purl/1364491.
@article{osti_1364491,
title = {Multiscale Capability in Rattlesnake using Contiguous Discontinuous Discretization of SelfAdjoint Angular Flux Equation},
author = {Zheng, Weixiong and Wang, Yaqi and DeHart, Mark D.},
abstractNote = {In this report, we present a new upwinding scheme for the multiscale capability in Rattlesnake, the MOOSE based radiation transport application. Comparing with the initial implementation of multiscale utilizing Lagrange multipliers to impose strong continuity of angular flux on interface of inbetween subdomains, this scheme does not require the particular domain partitioning. This upwinding scheme introduces discontinuity of angular flux and resembles the classic upwinding technique developed for solving first order transport equation using discontinuous finite element method (DFEM) on the subdomain interfaces. Because this scheme restores the causality of radiation streaming on the interfaces, significant accuracy improvement can be observed with moderate increase of the degrees of freedom comparing with the continuous method over the entire solution domain. Hybrid SNPN is implemented and tested with this upwinding scheme. Numerical results show that the angular smoothing required by Lagrange multiplier method is not necessary for the upwinding scheme.},
doi = {10.2172/1364491},
journal = {},
number = ,
volume = ,
place = {United States},
year = {Thu Sep 01 00:00:00 EDT 2016},
month = {Thu Sep 01 00:00:00 EDT 2016}
}

Nonlinear acceleration of a continuous finite element (CFE) discretization of the transport equation requires a modification of the transport solution in order to achieve local conservation, a condition used in nonlinear acceleration to define the stopping criterion. In this work we implement a coarsemesh finite difference acceleration for a CFE discretization of the secondorder selfadjoint angular flux (SAAF) form of the transport equation and use a postprocessing to enforce local conservation. Numerical results are given for onegroup source calculations of onedimensional slabs. We also give a novel formal derivation of the boundary conditions for the SAAF.

Nonlinear acceleration of a continuous finite element discretization of the selfadjoint angular flux form of the transport equation
Nonlinear acceleration of a continuous finite element (CFE) discretization of the transport equation requires a modification of the transport solution in order to achieve local conservation, a condition used in nonlinear acceleration to define the stopping criterion. In this work we implement a coarsemesh finite difference acceleration for a CFE discretization of the secondorder self adjoint angular flux (SAAF) form of the transport equation and use a post processing to enforce local conservation. Numerical results are given for onegroup source calculations of onedimensional slabs. We also give a formal derivation of the boundary conditions for the SAAF. (authors) 
Spatial and angular variation and discretization of the selfadjoint transport operator
This mathematical treatise begins with a variational derivation of a secondorder, selfadjoint form of the transport equation. Next, a space variational functional whose minimization solves the transport equation is derived. A onedimensional example is given. Then, {ital S{sub N}} and {ital P{sub N}} discretized functionals are expressed. Next, the surface contributions to the functionals are discretized. Finally, the explicit forms of the {rvec D} and {rvec H} matrices are given for four different geometries: hexahedron, wedge, tetrahedron, and pyramid. 
LeastSquares PN Formulation of the Transport Equation Using SelfAdjointAngularFlux Consistent Boundary Conditions.
In this paper, we study the LeastSquares (LS) PN form of the transport equation compatible with voids in the context of Continuous Finite Element Methods (CFEM).We first deriveweakly imposed boundary conditions which make the LS weak formulation equivalent to the SelfAdjoint Angular Flux (SAAF) variational formulation with a void treatment, in the particular case of constant crosssections and a uniform mesh. We then implement this method in Rattlesnake with the Multiphysics Object Oriented Simulation Environment (MOOSE) framework using a spherical harmonics (PN) expansion to discretize in angle. We test our implementation using the Method of Manufactured Solutions (MMS) and findmore » 
LeastSquares PN Formulation of the Transport Equation Using SelfAdjointAngularFlux Consistent Boundary Conditions
In this paper, we study the LeastSquares (LS) PN form of the transport equation compatible with voids [1] in the context of Continuous Finite Element Methods (CFEM).We first deriveweakly imposed boundary conditions which make the LS weak formulation equivalent to the SelfAdjoint Angular Flux (SAAF) variational formulation with a void treatment [2], in the particular case of constant crosssections and a uniform mesh. We then implement this method in Rattlesnake with the Multiphysics Object Oriented Simulation Environment (MOOSE) framework [3] using a spherical harmonics (PN) expansion to discretize in angle. We test our implementation using the Method of Manufactured Solutionsmore »