Proper orthogonal decomposition methods for noise reduction in particle-based transport calculations

Proper orthogonal decomposition techniques to reduce noise in the reconstruction of the distribution function in particle-based transport calculations are explored. For two-dimensional steady-state problems, the method is based on low rank truncations of the singular value decomposition of a coarse-grained representation of the particle distribution function. For time-dependent two-dimensional problems or three-dimensional time-independent problems, the use of a generalized low-rank approximation of matrices technique is proposed. The methods are illustrated and tested with Monte Carlo particle simulation data of plasma collisional relaxation and guiding-center transport with collisions in a magnetically confined plasma in toroidal geometry. It is observed that the proposed noise reduction methods achieve high levels of smoothness in the particle distribution function by using significantly fewer particles in the computations.