2.2.1.1 Mach2

Initialement développé par l'US Air Force Research Laboratory: Phillips Research Site, puis maintenu par Numerex , ce code était initialement destiné à la simulation de plasma focus. Son nom signifie Multiblock Arbitrary Coordinate Hydromagnetic. C'est donc un code de MHD non-idéale^2.4 2.5D , ALE ^2.5, écrit en Fortran, dont une version parallélisée existe, qui possède les caractéristiques suivantes:

monofluide ;
multi-températures ;
structure multiblocs ;
maillage non-orthogonal ;
remaillage possible ;

Les équations résolues par le code sont donc les suivantes (en supprimant le modèle de matériau élastoplastique):

$\begin{displaymath} \index{équation!conservation de la masse} \index{équation!... ... \, \rho}{dt} + \mathop{div}(\rho \overrightarrow{v}) = 0 \\ \end{displaymath}$

$\begin{displaymath} \index{équation!conservation de la quantité de mouvement} ... ...overrightarrow{B} - \mathop{div}(\stackrel{\Rightarrow}{\pi}) \end{displaymath}$

$\begin{displaymath} \index{équation!thermique} \left\{ \begin{array}{rcl} \f... ...ghtarrow{grad}}({E_{rad.}})) + \Phi_{eR} \end{array} \right. \end{displaymath}$

$\begin{displaymath} \index{équation!diffusion magnétique} \frac{\partial \, \o... ...( \frac{\mathop{\overrightarrow{grad}}({P_e})}{e n_e} \Biggr) \end{displaymath}$

$\begin{displaymath} \index{équation!puissance rayonnée} dV \cdot \Phi_{eR} = 4 \sigma \rho \lambda_{planck} ( T_e^4 -T_{rad.}^4 ) \cdot dV \end{displaymath}$

Un terme de viscosité numérique

sous forme d'une pression^2.6 est ajouté aux termes de pression

(le terme $\stackrel{\Rightarrow}{\pi}$ représentant le tenseur de viscosité). Les équations d'états peuvent être définies par des modèles analytiques simples ou bien par des tables SESAME , ainsi que les résistivités ( $\sigma$ ), opacités (libre parcours moyen de ROSSELAND $\lambda_{ross.}$ ), coefficients de conductibilité thermique électronique et ionique (respectivement $\lambda_e$ et $\lambda_i$ ).... Pour l'établissement des équations données au dessus, se référer à la section C.2 page

. Des détails quant à la résolution numérique de certaines équations sont donnés en section H.1 page

, les explications des paramètres définissant une simulation (paramètres regroupés dans une namelist) sont largement décrits dans le manuel [76].

Notes

... non-idéale ^2.4: et non relativiste
... ALE ^2.5: vient de l'anglais Arbitrary Lagragian Eulerian.
... pression ^2.6: voir en section H.1.7 page pour plus de précisions.

Mathias.Bavay_at_ingenieurs-supelec.org - juillet 2002