A non-negative, non-linear, Petrov-Galerkin method for bilinear discontinuous differencing of the SN equations Conference Paper uri icon

abstract

  • We have developed a new, non-negative, non-linear, Petrov-Galerkin bilinear discontinuous (BLD) finite element differencing of the 2-D Cartesian geometry SN equations for quadrilaterals on an unstructured mesh. This work is an extension of a scheme we previously developed for use with linear discontinuous (LD) differencing of the 2-D Sn equations for rectangular mesh cells. We present the theory and equations that describe the new method. Additionally, we numerically compare the accuracy of our proposed method to the accuracy of BLD without lumping and the subcell corner balance method (equivalent to a "fully" lumped BLD scheme) for a test problem that causes BLD scheme to generate negative angular flux solutions.

published proceedings

  • Mathematics and Computations, Supercomputing in Nuclear Applications and Monte Carlo International Conference, M and C+SNA+MC 2015

author list (cited authors)

  • Maginot, P. G., Ragusa, J. C., & Morel, J. E.

complete list of authors

  • Maginot, PG||Ragusa, JC||Morel, JE

publication date

  • January 2015