H.3 Modélisation du pseudo vide

Les codes MHD actuels ont des difficultés à traiter le cas du vide idéal. Ce vide est alors remplacé par un matériau à une certaine densité, en général très supérieure à la densité du vide expérimental du fait de la faible tolérance des codes face à une trop grande amplitude de densités. Il est alors nécessaire d'introduire des `` modèles physiques '' particuliers permettant à ce pseudo vide de se comporter le plus possible comme le vide expérimental. Dans Mach2 , ces paramètres sont principalement accessibles dans la section matmdl (ou bien dans la section ezphys pour les paramètres généraux). On modélise le rôle du vide par sa diffusivité, par sa densité minimale, ainsi que par divers seuils : certains modèles physiques sont ainsi activés ou désactivés en fonction de la densité de la cellule courante par rapport à un seuil de densité fixé dans la namelist. Les seuils accessibles sont donnés en table H.8.

Tableau H.8: seuils du pseudo-vide

paramètre	valeur par défaut	description
rof	$10^{-99}$	densité minimale possible.
rofvac	$0$	densité en dessous de laquelle la résistivité devient celle du vide, eta0 si celle ci est constante.
rofanom	$-1$	en dessous de ce seuil, la résistivité anormale devient aresvac.
rofjoule	$-1$	pas de chauffage joule en dessous de ce seuil.
rofhyd	$-1$	pas de résolution de l'hydrodynamique en dessous de ce seuil.
rofrad	$-1$	seuil en dessous duquel il n'y a plus de refroidissement radiatif.
rofcrad	$-1$	seuil au dessus duquel il n'y a plus de refroidissement radiatif.
rofsiecp	$-1$	densité utilisée pour le calcul du seuil de l'énergie interne.

L'impact de ces paramètres peut être énorme, il convient donc de fixer leurs valeurs sur des considérations physiques ou numériques.