Earth
Gradiente geotermico
Scalda dal basso una colonna crostale di 40 km con 30 mW/m² mentre la sua radioattività ne aggiunge 45: il profilo si incurva, ripido nei sedimenti e più dolce nella crosta inferiore conduttiva, e dopo alcune costanti di tempo di piastra di 23 Ma cade a un centesimo di grado dallo stato stazionario chiuso mentre il flusso superficiale raggiunge 75 mW/m². Dimezza la conducibilità dei sedimenti e tutta la crosta sotto va 180 K più calda; raddoppia il flusso basale e la Moho supera i 1300 gradi. Un semispazio omogeneo dopo un salto superficiale segue la soluzione erf. Calore fornito, generato ed evacuato eguaglia quello immagazzinato a meno dell’arrotondamento.

Cosa succede
Conduzione monodimensionale attraverso una crosta stratificata: ottanta volumi finiti di mezzo chilometro con la conducibilità e la produzione di calore di sedimenti (2 W/m/K, 1 µW/m³), crosta superiore (2,5, 2) e crosta inferiore (2,2, 0,5), conducibilità di faccia armoniche, la superficie tenuta a 10 °C e la base alimentata con 30 mW/m². Passi espliciti di 2 ka contabilizzano ogni flusso di faccia e la produzione, così che il calore immagazzinato eguaglia ciò che è entrato, generato e uscito a meno dell’arrotondamento; il profilo stazionario è la forma chiusa quadratica a tratti, e un semispazio omogeneo dopo un salto superficiale la soluzione erf. Un secondo di esperimento è un milione di anni.
Da provare
Guarda il profilo incurvarsi verso la sua forma stazionaria e leggi la temperatura della Moho contro la forma chiusa; dimezza la conducibilità della copertura, raddoppia il flusso basale, o avvia la colonna calda e segui erf.
Il modello e i suoi limiti
Conduzione pura con proprietà costanti per strato, niente erosione, avvezione né produzione variabile nel tempo, flusso basale fisso. Le temperature sono quelle di questa colonna.
Equazioni del modulo
Equazioni ρ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²/π²κ


