известным уравнением К1сВагф$
д0 _ д 
д 
дй
К(@ 
— | К(д) — 
||-О(х, 
! 
1
я дд[ ( )( 
Л+дх[ ( )[дхл О(х,х,Ё) 
()
где @(х х, ©) - объёмная влажность, (1/1/); No(@0)= \у + 7 — гидравлический
напор воды почвы, (1); у — 
капиллярный потенциал почвенной влаги
потенциал, (1.); х - вертикальная координата, с положительным направлением
вниз, (1.); х - горизонтальное координата, (1-); К(@) - функция гидравлической
проводимости при неполном насыщении почвы, (1./Т); & - время (Т); О(х,7,ё)
функция источника, характеризующая расход капельниц(в единицу времени)
(1/1).
Однозначность рассматриваемого решения уравнения (1) определяется
заданием краевых условий (начальные + граничные). Начальные условия
(начальный профиль) приняты однородными в горизонтальном направлении,
по вертикали возможно наличие слоёв почвогрунта, принимая во внимание
их переменную 
мощность 
и 
однородные 
воднофизические 
свойства 
в
пределах 
слоя. 
Заданный 
— объем — воды, 
подаваемой — капельницей,
преобразуется в слой воды в лунке, таким образом, на каждом шаге по
времени, получаем граничное условие первого рода — переменный во времени
напор воды в условиях инфильтрации. По вертикальным границам слева и
справа расчетного профиля задаётся условия симметрии, отражая таким
образом, «системы капельниц». Снизу рассматриваемого профиля задается
«свободное стекание».
Решение
Из-за нелинейности уравнения Ричарда, могут быть получены только
относительно немногие упрощенные аналитические решения. Практические
применения (1) требуют численного решения, которое может быть получено,
используя разнообразие методов, типа конечных разностей или конечных
элементов. Мы не обсуждаем использования — численных методов, при
решении уравнения (1), 
но 
отмечаем важность развития графического
интерфейса пользователя, чтобы гарантировать принятие и использование
пакета программ. Решение уравнения (1) проводится методом конечных
разностей. Расчетная область покрывается нерегулярной сеткой, т.е. ис-
пользовался 
переменный 
шаг 
по 
пространственным 
координатам, 
что
повышает точность расчетов. Чтобы соблюдать механизм физики процесса
(подача 
воды 
в 
лунку, 
заполнение 
лунки 
и 
сработка 
уровня 
при
инфильтрация) вводится набор правил:
1. Разбивка 
области 
расчета 
в 
поверхностной 
зоне 
почвогрунта
проводится с шагом А по горизонтальной и вертикальной координатам <=
0.01м. При этом обеспечивается, в том числе, реальное преставление лунки с
заданным радиусом (Клунки^ 0.4; 
0.2;0.12;0.05м) при — математическом
моделировании.
2. Лунка представляющая (рис. 1) часть круга (в плоскости), покрыта
квадратной сеткой. В лунке выделяются горизонтальные слои, начиная с
нижней части лунки (основание лунки).