Earth
Осадочный бассейн
Смотрите, как за десять миллионов лет заполняется бассейн длиной 100 км: фундамент опускается тем быстрее, чем дальше от шарнира, осадок приходит у шарнира и расползается вниз по склону диффузией, быстрее на суше, чем под водой, и каждые четверть миллиона лет ложится маркирующий пласт. Песок остаётся у источника, ил уходит дальше; пористость выжимается с глубиной, так что уплотнённая толща выше, чем твёрдое, что она содержит. Поднимите море — берег отступит; утройте снос — бассейн переполнится и перельётся за дальний край. Поданный, вынесенный и удержанный осадок сходятся до округления, а диффундирующий холм растит свою дисперсию ровно на 2 K t.
Что происходит
Разрез бассейна из пятидесяти колонок по 2 км: фундамент опускается вокруг шарнира у левого края, 0,2 км/млн лет на дальнем конце; осадок входит у шарнира со скоростью 10 км² за млн лет и расползается диффузией высоты поверхности с 1000 км²/млн лет на суше и 300 под морем; дальний край сбрасывает всё, что стоит выше уровня моря. Твёрдая толща сохраняется; уплотнённая толща каждой колонки следует из пористости 0,5 у поверхности, убывающей с масштабом 2 км; маркирующий горизонт ложится каждые 0,25 млн лет, а доля песка в отложенном убывает с расстоянием от источника. Явный шаг 1 тыс. лет, исходящие потоки колонки масштабируются вместе, чтобы она не отдала больше, чем держит; одна секунда опыта — один миллион лет.
Что попробовать
Смотрите, как ложатся пласты и движется берег; поднимите море — берег отступит; поднимайте снос, пока бассейн не перельётся; прочитайте песок у шарнира против ила вдали и уплотнённую толщу против твёрдой.
Модель и её ограничения
Диффузионный перенос с двумя постоянными коэффициентами, кинематический шарнир вместо изгиба и термического погружения, уплотнение по фиксированному закону пористости, без волн, русел и карбонатов. Учебный разрез, не бассейновая модель.
Уравнения модуля
Уравнения ∂h/∂t = −∂q/∂x − s(x), q = −K·∂h/∂x, s = s_max·x/L T − φ_0·λ·(1 − e^(−T/λ)) = S_solid, φ(z) = φ_0·e^(−z/λ) f_sand(x) = e^(−x/λ_s) Σ S·Δx = Q_in·t − Q_out, σ² = σ_0² + 2K·t
