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
