Moment-preserving particle merging via non-negative least squares

A novel particle merging algorithm for rarefied gas dynamics simulations is proposed that can conserve arbitrary velocity and spatial moments of the particle distribution via solving a non-negative least squares problem. An extension that preserves both exact and approximate collision rates is also derived. The algorithm is applied to the simulation of several model rarefied gas dynamics problems, where it exhibits noticeably lower merging-induced error in key macroscopic quantities.