Jean-Philippe Braeunig

Raphaël Loubère, Jean-Philippe Braeunig, Jean-Michel Ghidaglia

This work concerns the simulation of compressible multi-material fluid flows and follows the method FVCF-NIP described in the former paper Braeunig et al (Eur. J. Mech. B/Fluids, 2009). This Cell-centered Finite Volume method is totally Eulerian since the mesh is not moving and a sharp interface, separating two materials, evolves through the grid. A sliding boundary condition is enforced at the interface and mass, momentum and total energy are conserved. Although this former method performs well on 1D test cases, the interface reconstruction suffers of poor accuracy in conserving shapes for instance in linear advection. This situation leads to spurious instabilities of the interface. The method Enhanced-NIP presented in the present paper cures an inconsistency in the former NIP method that improves strikingly the results. It takes advantage of a more consistent description of the interface in the numerical scheme. Results for linear advection and compressible Euler equations for inviscid fluids are presented to assess the benefits of this new method.

Jerome Guterl, Jean-Philippe Braeunig, Nicolas Crouseilles et al.

The purpose of this work is simulation of magnetised plasmas in the ITER project framework. In this context, Vlasov-Poisson like models are used to simulate core turbulence in the tokamak in a toroidal geometry. This leads to heavy simulation because a 6D dimensional problem has to be solved, 3D in space and 3D in velocity. The model is reduced to a 5D gyrokinetic model, taking advantage of the particular motion of particles due to the presence of a strong magnetic field. However, accurate schemes, parallel algorithms need to be designed to bear these simulations. This paper describes a Hermite formulation of the conservative PSM scheme which is very generic and allows to implement different semi-Lagrangian schemes. We also test and propose numerical limiters which should improve the robustness of the simulations by diminishing spurious oscillations. We only consider here the 4D drift-kinetic model which is the backbone of the 5D gyrokinetic models and relevant to build a robust and accurate numerical method.

Fabien Vivodtzev, Florent Nauleau, Jean-Philippe Braeunig et al.

This work presents an approach for the automatic detection of locally turbulent vortices within turbulent 2D flows such as instabilites. First, given a time step of the flow, methods from Topological Data Analysis (TDA) are leveraged to extract the geometry of the vortices. Specifically, the enstrophy of the flow is simplified by topological persistence, and the vortices are extracted by collecting the basins of the simplified enstrophy's Morse complex. Next, the local kinetic energy power spectrum is computed for each vortex. We introduce a set of indicators based on the kinetic energy power spectrum to estimate the correlation between the vortex's behavior and that of an idealized turbulent vortex. Our preliminary experiments show the relevance of these indicators for distinguishing vortices which are turbulent from those which have not yet reached a turbulent state and thus known as laminar.