J.2 Lecture des tables eos

Lors du lancement d'une itération de Mach2 , les tables d'équations d'état utilisées par la namelist sont envoyées dans un fichier appelé eos. Ce fichier se présente de la façon suivante:

1matid =  3333    znbar =  29.0000    anbar =  63.5400    dens0 =  8.9300E+03   ifile = sesame               1

  f(x,y) = total pressure (Gj/m3 / fcnvt)                 x = dens (Mg/m3 / xcnvt)    y = temp (kelvin / ycnvt)
  fcnvt =  1.0000E+09                                     xcnvt =  1.0000E+03         ycnvt =  8.6173E-05

         y =  0.0000E+00  2.5850E-02  3.1600E-02  4.0000E-02  5.0000E-02  6.3000E-02  8.0000E-02  1.0000E-01
 x =
  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
  1.0000E-03  6.0678E-59  2.0535E-37  8.4946E-29  1.1681E-20  6.0582E-17  1.0468E-16  1.5289E-07  7.8443E-04
  1.6000E-03  6.0678E-59  2.0535E-37  8.4946E-29  1.1681E-20  6.0582E-17  1.0468E-16  1.5289E-07  7.8443E-04
  2.5000E-03  6.0678E-59  2.0535E-37  8.4946E-29  1.1681E-20  6.0582E-17  1.0468E-16  1.5289E-07  7.8443E-04
  4.0000E-03  6.0678E-59  2.0535E-37  8.4946E-29  1.1681E-20  6.0582E-17  1.0468E-16  1.5289E-07  7.8443E-04
  6.3000E-03  6.0678E-59  2.0535E-37  8.4946E-29  1.1681E-20  6.0582E-17  1.0468E-16  1.5289E-07  7.8443E-04
  1.0000E-02  6.0678E-59  2.0535E-37  8.4946E-29  1.1681E-20  6.0582E-17  1.0468E-16  1.5289E-07  7.8443E-04
  1.6000E-02  6.0678E-59  2.0535E-37  8.4946E-29  1.1681E-20  6.0582E-17  1.0468E-16  1.5289E-07  7.8443E-04

La première ligne comporte l'identificatif du matériau : 1matid = 3333. Celui-ci doit être comparé avec la liste de matériaux donnée en annexe J.4^J.2 en sachant que les catégories sont de trois types:

propriétés générales

telles que pression, ...c'est la catégorie 1;

opacités ...

c'est la catégorie 2;

conductivité électrique

, c'est la catégorie 3.

L'identificatif est de la forme cmmmv (dans le cas d'une table personnelle, la forme est c9mmmv), avec:

c

= catégorie - 1 ^J.3;

mmm

= numéro du matériau;

v

= numéro de version.

Par exemple, 3333 rerésente une propriété de catégorie 1 du cuivre (matériau 333), version 3. Suivent le numéro atomique et la masse atomique (respectivement znbar = 29.0000 et anbar = 63.5400) ainsi que la densié initiale présente dans la namelist. Le paragraphe suivant comporte les unités et titres de colonnes pour le tableau qui occupe le reste de la section. On y apprend donc que la première colonne (celle des x) représente la densité, en $\, \mathrm{Mg/m^3} \cdot \frac{1}{1 \cdot {10}^3}$, soit des $\, \mathrm{kg/m^3}$; que la première ligne représente la température, en $\frac{\, \mathrm{K}}{8.6173 \cdot {10}^{-5}}$, c'est à dire en $\, \mathrm{eV}$. De plus, la grandeur calculée est la pression (f(x,y) = total pressure). Ce`` bloc logique '' se répète plusieurs fois à la suite: il faut imaginer que c'est la représentation d'un tableau $Y = \mathop{F}(\rho,T)$ très large que l'on a choisi de découper à une largeur de page, donc qui se poursuit pour des températures plus élevées en dessous du bloc que l'on est en train de lire. D'après le critère défini en section J.1, il faut chercher à annuler la pression pour le matériau à la densité initiale fixée dans la namelist Mach2 . Il faut donc:

1. rechercher le matériau pour lequel on veut la température magique;
2. rechercher la table correspondant à sa pression, indiquée par f(x,y) = total pressure;
3. rechercher dans celle-ci la ligne correspondant à la densité initiale du matériau.

Lorsque ceci est fait, on cherche dans cette ligne le changement de signe de la pression:

  8.9300E+03 -2.1128E+09  1.0000E+05  4.0654E+08  1.0823E+09  1.8876E+09  2.9358E+09  4.3087E+09  5.9272E+09

Le premier nombre est la densité, les suivants sont les valeurs de la pression. L'annulation de la pression a donc lieux entre $T_1 = 0.0000E+00$ et $T_2 = 2.5850E-02$ (passage de $P_1 = -2.1128E+09$ à $P_2 = 1.0000E+05$).

Notes

...annexe:mat_Tmag^J.2

Cette liste est issue de [76, section 2.22]

... 1^J.3

rien si c-1=0