A new pressure-based lattice-Boltzmann method (HRR-p) is proposed for the simulation of flows for Mach numbers ranging from 0 to 1.5. Compatible with nearest-neighbor lattices (e.g., D3Q19), the model consists of a predictor step comparable to classical athermal lattice-Boltzmann methods, appended with a fully local and explicit correction step for the pressure. Energy conservation—for which the Hermitian quadrature is not accurate enough on such a lattice—is solved via a classical finite volume MUSCL-Hancock scheme based on the entropy equation. The Euler part of the model is then validated for the transport of three canonical modes (vortex, entropy, and acoustic propagation), while its diffusive/viscous properties are assessed via thermal Couette flow simulations. All results match the analytical solutions with very limited dissipation. Last, the robustness of the method is tested in a one-dimensional shock tube and a two-dimensional shock–vortex interaction.

International audienceAn original penalization method is applied to model the interaction of magnetically confined plasma with limiter in the frame of minimal transport model for ionic density and parallel momentum. The limiter is considered as a pure particle sink for the plasma and consequently the density and the momentum are enforced to be zero inside. Comparisons of the numerical results with one dimensional analytical solutions show a very good agreement. In particular, presented method provides a plasma velocity which is almost sonic at the boundaries obstacles as expected from the sheath conditions through the Bohm criterion. The new system being solved in an obstacle free domain, an efficient pseudo-spectral algorithm based on a Fast Fourier transform is also proposed, and associated with an exponential filtering of the unphysical oscillations due to Gibbs phenomenon. Finally, the efficiency of the method is illustrated by investigating the flow spreading from the plasma core to the Scrape Off Layer at the wall in a two-dimensional system with one then two limiters neighboring