WikiSort.ru - Не сортированное

Уравнение Ричардса, описывающее влагоперенос в зоне аэрации (в ненасыщенной зоне), было сформулировано Лоренцо А. Ричардсом в 1931 году^[1]. Оно представляет собой нелинейное дифференциальное уравнение в частных производных, основная трудность решения которого заключается в отсутствии точных аналитических решений.

Закон Дарси разработан для описания потока влаги в водонасыщенной пористой среде, адаптируя его, Ричардс использовал предложение Букингема (англ.)русск. (1907) и вывел общее дифференциальное уравнение в частных производных, которое описывает ненасыщенный влагоперенос в зоне аэрации. Наиболее известной формой записи является:

${\displaystyle {\frac {\partial \theta }{\partial t}}={\frac {\partial }{\partial z}}\left[K(\theta )\left({\frac {\partial \psi }{\partial z}}+1\right)\right]\ }$

где

${\displaystyle K}$  — коэффициент фильтрации в ненасыщенной зоне,

${\displaystyle \psi }$  — высота всасывающего давления,

${\displaystyle z}$  — высота над плоскостью сравнения,

${\displaystyle \theta }$  — объемная влажность, и

${\displaystyle t}$  — время.

Уравнение Ричардса аналогично уравнению влагопереноса в насыщенной зоне, переход от одного к другому обусловлен представлением напора в виде h = ψ + z, и заменой насыщенного потока ненасыщенным. Использование формы записи приведенной выше обусловлено удобством описания граничных условий (часто описываемых в терминах напора, например для использования атмосферного условия ψ = 0).

Вывод уравнения

В этом разделе будет показан вывод уравнения Ричардса для вертикального влагопереноса в очень упрощенном виде. Закон сохранения массы гласит, что величина изменения водонасыщения в закрытом объеме равна величине изменения суммы всех (отрицательных и положительных) потоков влаги в данном объеме. Опишем это математическим языком:

${\displaystyle {\frac {\partial \theta }{\partial t}}={\vec {\nabla }}\cdot \left(\sum _{i=1}^{n}{{\vec {q}}_{i,\,{\text{in}}}}-\sum _{j=1}^{m}{{\vec {q}}_{j,\,{\text{out}}}}\right)}$

Введем одномерную форму записи для направления ${\displaystyle {\hat {k}}}$ :

${\displaystyle {\frac {\partial \theta }{\partial t}}=-{\frac {\partial }{\partial z}}q}$

Горизонтальный поток описан эмпирическим Законом Дарси:

${\displaystyle q=-K{\frac {\partial h}{\partial z}}}$

Подставив q в уравнение выше, получим:

${\displaystyle {\frac {\partial \theta }{\partial t}}={\frac {\partial }{\partial z}}\left[K{\frac {\partial h}{\partial z}}\right]}$

И используем выражение для напора h = ψ + z:

${\displaystyle {\frac {\partial \theta }{\partial t}}={\frac {\partial }{\partial z}}\left[K\left({\frac {\partial \psi }{\partial z}}+{\frac {\partial z}{\partial z}}\right)\right]={\frac {\partial }{\partial z}}\left[K\left({\frac {\partial \psi }{\partial z}}+1\right)\right]}$

Таким образом, получено уравнение Ричардса^[2] .

Формулировка

Уравнение Ричардса используется во многих статьях, посвященных вопросам экологии, потому что оно описывает поток на границе подземных и поверхностных вод (атмосферных осадков, водоемов, рек и т. д.), что является чрезвычайно важным при моделировании миграции различных компонентов. Его так же рассматривают и в чисто математических журналах, потому что его решение не тривиально. Обычно, оно выражается в трех разных формах. Смешанная форма (mixed form), рассмотренная выше, включает в себя описание в терминах и напора и влажности. А также две другие формы записи, в терминах напора (head-based) и водонасыщения (saturation-based).

Ограничения применимости

Численное решение уравнения влагопереноса Ричардса критиковалось за большие вычислительные затраты и непредсказуемость^[4][5], обусловленную тем, что нет гарантии, что вычислительный комплекс (solver) найдет решение (сойдется) для конкретных почвенных характеристик. Также отмечалось, что в данном методе преувеличена  роль капиллярных сил^[6] и в некоторой степени его критиковали за излишнюю простоту^[7]. При моделировании одномерного влагопереноса метод требует разбивки высокой дискретности для области граничащей с поверхностью (величина дискретизации должна быть не более одного сантиметра).  При трехмерном моделировании численное решение уравнения Ричардса определяется соотношением сторон, где отношение горизонтального и вертикального размеров ячейки модели в расчетной области должно быть не больше 7.

Ссылки

См. также

• Гидрогеология
• Основная гидрофизическая характеристика
• Finite water-content vadose zone flow method