Earth
Мерзлота і відтавання
Стежте за ґрунтом під будинком у тундрі зі зміною сезонів: щоліта верх відтає на глибину діяльного шару, щозими знову замерзає, а нижче ґрунт лишається мерзлим при температурі плавлення свого льоду. Товстий сніг не пускає зимовий холод і тримає ґрунт на п’ять градусів теплішим; згорніть його — і мерзлота зміцніє. Нагрівайте клімат на градус за десятиліття — діяльний шар поглиблюється рік у рік, доки не відтане льодистий шар і будинок не осяде на його надлишковий лід. Прихована теплота льоду стримує фронт відтавання точно за розв’язком Неймана, а тепло, що пройшло крізь поверхню й основу, — це ентальпія, набрана колоною.

Що відбувається
Колона ґрунту завглибшки 15 м у 75 комірках, 35 % води з льодистим шаром 70 % між 1,5 і 3 м. Кожна комірка зберігає ентальпію відносно мерзлого ґрунту при 0 °C: нижче нуля вона мерзла, між нулем і прихованою теплотою її води — суміш при 0 °C, вище — тала; теплопровідність бере гармонічні середні мерзлої і талої. Повітря (середнє −6 °C, ±18 °C, тренд потепління) діє на верх крізь 0,3 м снігу взимку; 60 мВт/м² тепла Землі входять в основу. Надлишковий лід дає осідання при першому відтаванні. Явні кроки по 3 год; тепло крізь верх і основу дорівнює набраній ентальпії, а розв’язок Неймана перевіряє фронт відтавання.
Що спробувати
Подивіться річне відтавання й промерзання; згорніть сніг або поглибте його; нагрівайте клімат на градус за десятиліття і дочекайтеся осідання будинку; вимкніть приховану теплоту і дивіться, як швидко біжить фронт.
Модель та її обмеження
Одновимірна теплопровідність, різка точка замерзання, фіксовані властивості ґрунту, без міграції води, морозного здимання й рослинності; осідання — заявлений облік надлишкового льоду. Глибини — цієї ідеалізованої колони.
Рівняння модуля
Рівняння H = C_f·T (T < 0), 0…L (T = 0), L + C_u·T (T > 0) dH/dt = ∂/∂z(k·∂T/∂z), q0 = (T_air − T1)/(R_snow + Δz/2k) X(t) = 2λ·√(α·t), λ·e^(λ²)·erf λ = St/√π ΔH_column = ∫q0 dt + ∫q_geo dt, s = Σ(w_ice − w)·Δz


