VLASIX
Codes and methods
Theory and analysis
Past events
Publications People
edit SideBar

Waterbags and related models
We are interested in the mathematical analysis of physical models which consists in showing that they are wellposed, i.e. they admit a unique solution, in a wellsuited functional framework. The notion of solution can change according to the nature of the nonlinearities in the model. We have to study the particular properties of the solution such as its regularity or its longtime asymptotic behavior. Moreover spectral analysis of the models is very useful since it supplies a lot of information about the behavior of the system such as instability growth rates and spatial structure of the unstable global eigenmodes. For example, we have obtained a set of results about the derivation and the mathematical analysis of waterbagtype models in different physical configurations (Besse 2011a, 2011b, 2012; Bardos & Besse 2013, 2015).
Waterbag models (Figure 1) come from the exact geometric reduction of the Vlasov equation thanks to geometric Liouville invariants. Their mathematical structures are closely related to hyperbolic systems of conservation laws, in finite or infinite dimension with nonlocal fluxes.
Figure 1: Schematic representation of phasespace with waterbags. In each of the region with a given color, the phasespace distribution function is approximated by a constant and is conserved during motion as a consequence of Liouville theorem. From Besse & Coulette (2016).
Figure 2: Toroidal geometry (left panel) and the structure of the magnetic field lines (right panels). From Besse & Coulette (2016).
We have also made the mathematical analysis of the VlasovDiracBenney model, and show its intrinsic links with plasma physics, fluid mechanics and quantum mechanics (Bardos & Besse 2013, 2015). More precisely, in Bardos & Besse (2016), we have studied the semiclassical limit of
an infinite dimensional system of coupled nonlinear Schrödinger equations
towards weak solutions of the VlasovDiracBenney equation, for initial
data with analytical regularity in space. Regarding to spectral analysis, we have performed a numerical study of the eigenvalue problem for the gyrowaterbag model in cylindrical geometry (Coulette & Besse 2013a, 2013b; Morel et al., 2014). An asymptotic and spectral analysis of the gyrowaterbag integrodifferential operator in toroidal geometry (Figure 2; Besse 2016, 2017; Besse & Coulette 2016) has led to the development of a parallel numerical code whose results are very promising and in good agreement with those of a quasilinear code and standard gyrokineticVlasov code (Coulette & Besse 2017). In addition we often have to build and justify new reduced models, starting from original base ones, in order to take into account their own properties. This modeling effort is necessary to reduce the computational cost of original base models while it allows to retain the essential physics of the studied phenomenon. For example, using an Eulerian variational principle of least action (EulerPoincaré), asymptotic perturbation methods and Lagrangian averaging techniques, we have obtained a system of PDEs called « Lagrangian averaged gyrowaterbag continuum » (LAGWBC) designed to accurately capture the dynamics of the turbulent plasma flow at length scales larger than a typical length scale while averaging the motion at scales smaller than this typical length scale. In the isotropic setting, we have shown that this system has unique strong solutions locally in time. Regarding to the original gyrowaterbag continuum, the LAGWBC equations show some additional properties and several advantages from the mathematical (regularizing terms, wellposedness) and physical (modeling of small scales, anisotropic tensor of turbulent fluctuations, dispersion nonlinear terms) viewpoints, which make this model a good candidate for describing accurately gyrokinetic turbulence in magnetically confined plasma (Besse 2016).
We also aim at performing a numerical analysis of the physical models. In other words, we design efficient numerical schemes for solving nonlinear PDEs (kinetic and multifluid transport equations, elliptic and waves equations). We also make the mathematical analysis of these numerical schemes, i.e. we prove in a wellsuited mathematical framework, the convergence of the approximated numerical solution of the discrete problem towards the exact solution of the continuous problem and obtain a priori error estimates. For instance, we constructed and analyzed several discontinuousGalerkin numerical schemes which do not increase a priori the total energy for the gyrokineticwaterbag equations (Besse 2017).
Finally, we are currently developing very highorder schemes for both the 3D incompressible Euler equation in bounded domains and the VlasovPoisson
equation in the whole space or with periodic boundary conditions. These schemes are based on a new Lagrangian formulation of these equations and make a crucial use of the invariants preserved by the flow. Moreover these numerical schemes give constructive proofs for showing high regularity of the flow launched by initial conditions with moderate smoothness.
References
 C. Bardos, N. Besse, The Cauchy problem for the VlasovDiracBenney equation and related issues in fluid mechanics and semiclassical limits, Kinetic and Related Models, 6 (2013), 893917.
 C. Bardos, N. Besse, Hamiltonian structure, fluid representation and stability for the VlasovDiracBenney equation Fields Inst. Commun. 75 (2015), 130.
 C. Bardos, N. Besse, Semiclassical limit of an infinite dimensional system of nonlinear Schrödinger equations, Bardos C., Besse N., Bull. Inst. Math. Acad. Sin. 11 (2016), 4361
 N. Besse, On the Cauchy problem for the gyrowaterbag model, Math. Mod. Meth. Appl. Sci. 21 (2011a), 18391869.
 N. Besse, On the waterbag continuum, Arch. Rational Mech. Anal., 199 (2011b), 453491.
 N. Besse, Global weak solution for the relativistic waterbag continuum, Math. Mod. Meth. Appl. Sci. 22 (2012), 1150001.
 N. Besse, Lagrangian averaged gyrokineticwaterbag continuum, Commun. Math. Sci. 14 (2016), 593626
 N. Besse, Discontinuous Galerkin finite element methods for the gyrokineticwaterbag equations, IMA Journal of Numerical Analysis 37 (2017), 9851040
 N. Besse, D. Coulette, Asymptotic and spectral analysis of the gyrokineticwaterbag integrodifferential operator in toroidal geometry, J. Math. Phys. 57 (2016), 081518
 D. Coulette, N. Besse, Numerical comparisons of gyrokinetic multiwaterbag models, J. Comput. Phys., 248 (2013a), 132.
 D. Coulette, N. Besse, Multiwaterbag models of ion temperature gradient instability in cylindrical geometry, Phys. Plasmas, 20 (2013b), 052107.
 P. Morel, F. Dreydemy Ghiro, V. Berionni, D. Coulette, N. Besse, O. Gurcan, A Multiwaterbag model of driftkinetic electron, Eur. Phys. J. D, 68 (2014), 220.
 D. Coulette, N. Besse, Numerical resolution of the globaleigenvalue problem for the gyrokineticwaterbag model in toroidal geometry, J. Plasma Phys. 83 (2017), 905830207
