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


