Sur la structure mathématique et l'approximation numérique de l'hydrodynamique lagrangienne bidimensionnelle

Abstract

Ce travail étudie une nouvelle formulation des équations d’Euler compressibles écrites en coordonnées lagrangiennes multidimensionnelles, sous la structure d’un système de lois de conservation associé à une contrainte de divergence nulle. Cette structure s’applique également à une physique plus large, incluant par exemple la magnétohydrodynamique. Elle permet l’étude mathématique du problème global couplant les inconnues physiques avec les inconnues géométriques associées au déplacement de la matière. Nous montrons que la partie physique du système, dont l’ensemble est faiblement hyperbolique, est symétrisable sous cette contrainte, alors que la perte de régularité des inconnues géométriques est caractéristique des cisaillements. Nous construisons ensuite une méthode d’approximation originale, de type Volumes-Finis sur maillage mobile, avec degrés de liberté aux noeuds. Le schéma repose uniquement sur des considérations physiques (principe de conservation, croissance de l’entropie), qui assurent sa stabilité ainsi qu’un résultat de positivité et de non-croisement sur maillage triangulaire. On montre théoriquement qu’il converge avec un ordre 1/2 sur les équations linéaires de l’acoustique. Enfin, une extension à la géométrie axisymétrique généralise la structure symétrisée et la contrainte différentielle, ainsi que la méthode d’approximation associée.This work studies a nex formulation of compressible Euler equations written in multidimensional Lagrangian coordinates, as a system of conservation laws linked to a free divergence constraint : It applies also to an extended physics, including for instance magnetohydrodynamics. This structure allows the mathematical study of the whole Lagrangian problematics, coupling physical unknowns with geometrical one’s, associated to the displacement of matter. We prove that the physical part of the system, whose entire formulation is known to be only weakly hyperbolic, is symmetrizable under the differential constraint, although the loss of regularity of geometrical unknowns is characteristic of shear discontinuities. We then derive an original approximate method, of Finite Volumes kind on moving mesh, whose degree of freedom are placed on nodes. The scheme relies only on physical considerations (conservation principle, entropy production), which ensure its stability as well as a result of positivity and non-crossing on triangular mesh. We prove theoretically its convergence with a rate of 1/2 on the linearized equations of acoustics. Eventually, we extent the symmeterized structure and numerical method to the case of axisymmetric geometry

    Similar works