MODELISATION D'UNE SUBDUCTION GRACE AU LOGICIEL THERMOMECANIQUE ADELI

Paramétrisation de la durée du modèle et du nombre d’itération

Avec TFIN le temps de subduction en seconde, NTIME le nombre d’itération effectuée durant le calcul.

Il faut choisir un temps de subduction d’une durée réaliste c'est-à-dire dans le cas de la plaque nazca pas plus de 23 millions d’années et un nombre d’itération pas trop important afin que la durée de calcul soit raisonnable.

 

Paramétrisation de la vitesse des plaques

 Nous désirons conserver la même longueur de plaque subduite. Il faut donc d’adapter la vitesse de poussée sur le slab.

avec L la longueur de plaque subduite en mètre, TFIN durée de la simulation en seconde, V vitesse de poussée en.

Paramétrisation de la viscosité

Pour choisir la viscosité, il est également nécessaire de prendre en compte la durée de subduction simulée. Cette durée doit être supérieure aux temps de relaxation du matériau considéré pour que le système puisse revenir à l’équilibre.

Le temps de relaxation est donné par

d’où

Avec  la viscosité en Pa.s,  le coefficient de poisson , E le module d’Young en Pa.  

La viscosité n’est pas directement présente dans le fichier d’entré iessai, c’est un paramètre qui lui est relié, la fluidité, qui est utilisé pour paramétrer les simulations.

 

 

 

Avec  en , Ea l’énergie d’activation en , T la température en Kelvin, R constante des gaz parfait en .

Afin de pouvoir négliger les effets du gradient de température terrestre  et par conséquent avoir une viscosité qui ne dépend pas de la profondeur dans le modèle que nous utilisons, nous avons pris une énergie d’activation nulle d’où