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


  • 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.

author list (cited authors)

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

publication date

  • January 2015