Earth
Gradient géothermique
Chauffez par le bas une colonne crustale de 40 km à 30 mW/m² pendant que sa propre radioactivité en ajoute 45 : le profil se courbe, raide dans les sédiments et plus doux dans la croûte inférieure conductrice, et après quelques constantes de temps de plaque de 23 Ma il tombe à un centième de degré de l’état stationnaire fermé tandis que le flux de surface atteint 75 mW/m². Divisez par deux la conductivité des sédiments et toute la croûte dessous tourne 180 K plus chaud ; doublez le flux basal et le Moho dépasse 1300 degrés. Un demi-espace homogène après un saut de surface suit la solution en erf. Chaleur fournie, produite et évacuée égale la chaleur stockée à l’arrondi près.

Ce qui se passe
Conduction unidimensionnelle à travers une croûte stratifiée : quatre-vingts volumes finis d’un demi-kilomètre avec la conductivité et la production de chaleur des sédiments (2 W/m/K, 1 µW/m³), de la croûte supérieure (2,5, 2) et de la croûte inférieure (2,2, 0,5), conductivités de face harmoniques, la surface tenue à 10 °C et la base alimentée à 30 mW/m². Des pas explicites de 2 ka comptabilisent chaque flux de face et la production, de sorte que la chaleur stockée égale ce qui est entré, produit et sorti à l’arrondi près ; le profil stationnaire est la forme fermée quadratique par morceaux, et un demi-espace homogène après un saut de surface la solution en erf. Une seconde d’expérience vaut un million d’années.
À essayer
Regardez le profil se courber vers sa forme stationnaire et lisez la température du Moho contre la forme fermée ; divisez la conductivité de la couverture par deux, doublez le flux basal, ou démarrez la colonne chaude et suivez erf.
Le modèle et ses limites
Conduction pure à propriétés constantes par couche, pas d’érosion, d’advection ni de production variable dans le temps, flux basal fixé. Les températures sont celles de cette colonne.
Équations du module
Équations ρc·∂T/∂t = ∂/∂z(k·∂T/∂z) + A, q = −k·∂T/∂z k_face = 2k_a·k_b/(k_a + k_b), q(z) = q_b + ∫_z^L A·dz′ T(z) = T_s + ∫_0^z q/k·dz′, q_surface = q_b + Σ A·Δz T(z, t) = T_s + (T_0 − T_s)·erf(z/2√(κt)), τ = 4L²/π²κ


