Numerical approximations of the 10-moment Gaussian closure

We propose a numerical scheme to approximate the weak solutions of the 10-moment Gaussian closure. The moment Gaussian closure for gas dynamics is governed by a conservative hyperbolic system supplemented by entropy inequalities whose solutions satisfy positiveness of density and tensorial pressure. We consider a Suliciu-type relaxation numerical scheme to approximate the solutions. These methods are proved to satisfy all the expected positiveness properties and all the discrete entropy inequalities. The scheme is illustrated by several numerical experiments.