C. Klingenberg

C. Klingenberg, G. Puppo, M. Semplice

This work presents arbitrary high order well balanced finite volume schemes for the Euler equations with a prescribed gravitational field. It is assumed that the desired equilibrium solution is known, and we construct a scheme which is exactly well balanced for that particular equilibrium. The scheme is based on high order reconstructions of the fluctuations from equilibrium of density, momentum and pressure, and on a well balanced integration of the source terms, while no assumptions are needed on the numerical flux, beside consistency. This technique allows to construct well balanced methods also for a class of moving equilibria. Several numerical tests demonstrate the performance of the scheme on different scenarios, from equilibrium solutions to non steady problems involving shocks. The numerical tests are carried out with methods up to fifth order in one dimension, and third order accuracy in 2D.

A. Crestetto, C. Klingenberg, M. Pirner

This work is devoted to the numerical simulation of the \BGK equation for two species in the fluid limit using a particle method. Thus, we are interested in a gas mixture consisting of two species without chemical reactions assuming that the number of particles of each species remains constant. We consider the kinetic two species model proposed by Klingenberg, Pirner and Puppo in 2017, which separates the intra and interspecies collisions. We want to study numerically the influence of the two relaxation term, one corresponding to intra, the other to interspecies collisions. For this, we use the method of micro-macro decomposition. First, we derive an equivalent model based on the micro-macro decomposition (see Bennoune, Lemou and Mieussens, 2007 and Crestetto, Crouseilles and Lemou, 2013). The kinetic micro part is solved by a particle method, whereas the fluid macro part is discretized by a standard finite volume scheme. Main advantages of this approach are: (i) the noise inherent to the particle method is reduced compared to a standard (without micro-macro decomposition) particle method, (ii) the computational cost of the method is reduced in the fluid limit since a small number of particles is then sufficient.