An unconventional divergence free Finite Volume discretization of Lagrangian MHD - Direction des applications militaires Accéder directement au contenu
Pré-Publication, Document De Travail Année : 2023

An unconventional divergence free Finite Volume discretization of Lagrangian MHD

Résumé

We construct an unconventional divergence free discretization of updated Lagrangian ideal MHD over simplicial grids. The cell-centered FV method employed to discretize the conservation laws of volume, momentum and total energy is rigorously the same than the one developed to simulate hyperelasticity equations. By construction this moving mesh method ensures the compatibility between the mesh displacement and the approximation of the volume flux by means of the nodal velocity and the attached unit corner normal vector which is nothing but the partial derivative of the cell volume with respect to the node coordinate under consideration. This is precisely the definition of the compatibility with the Geometrical Conservation Law which is the cornerstone of any proper multi-dimensional moving mesh FV discretization. The momentum and the total energy fluxes are approximated utilizing the partition of cell faces into sub-faces and the concept of sub-face force which is the traction force attached to each sub-face impinging at a node. We observe that the time evolution of the magnetic field might be simply expressed in terms of the deformation gradient which characterizes the Lagrange-to-Euler mapping. In this framework the divergence of the magnetic field is conserved with respect to time thanks to Piola An unconventional divergence free FV discretization of Lagrangian MHD formula. Therefore, we solve the fully compatible updated Lagrangian discretization of the deformation gradient tensor for updating in a simple manner the cell-centered value of the magnetic field. Finally, the sub-face traction force is expressed in terms of the nodal velocity to ensure a semi-discrete entropy inequality within each cell. The conservation of momentum and total energy is recovered prescribing the balance of all the sub-face forces attached to the sub-faces impinging at a given node. This balance corresponds to a vectorial system satisfied by the nodal velocity. It always admits a unique solution which provides the nodal velocity. The robustness and the accuracy of this unconventional FV scheme has been demonstrated employing various representative test cases. Finally, it is worth emphasizing that once you have an updated Lagrangian code for solving hyperelasticity you also get an almost free updated Lagrangian code for solving ideal MHD ensuring an exact divergence free constraint of the magnetic field at the discrete level.
Fichier principal
Vignette du fichier
Manuscript.pdf (37.98 Mo) Télécharger le fichier
Origine : Fichiers produits par l'(les) auteur(s)

Dates et versions

hal-04010654 , version 1 (01-03-2023)
hal-04010654 , version 2 (13-08-2023)

Identifiants

  • HAL Id : hal-04010654 , version 1

Citer

Walter Boscheri, Raphaël Loubère, Pierre-Henri Maire. An unconventional divergence free Finite Volume discretization of Lagrangian MHD. 2023. ⟨hal-04010654v1⟩
182 Consultations
6 Téléchargements

Partager

Gmail Facebook X LinkedIn More