Meshless local Petrov–Galerkin solution of the neutron transport equation with streamline-upwind Petrov–Galerkin stabilization

Abstract

The meshless local Petrov–Galerkin (MLPG) method is applied to the steady-state and k-eigenvalue neutron transport equations, which are discretized in energy using the multigroup approximation and in angle using the discrete ordinates approximation. To prevent oscillations in the neutron flux, the MLPG transport equation is stabilized by the streamline upwind Petrov–Galerkin (SUPG) method. Global neutron conservation is enforced by using moving least squares basis and weight functions and appropriate SUPG parameters. The cross sections in the transport equation are approximated in accordance with global particle balance and without constraint on their spatial dependence or the location of the basis and weight functions. The equations for the strong-form meshless collocation approach are derived for comparison to the MLPG equations. The method of manufactured solutions is used to verify the resulting MLPG method in one, two and three dimensions. Results for realistic problems, including two-dimensional pincells, a reflected ellipsoid and a three-dimensional problem with voids, are verified by comparison to Monte Carlo simulations.