на каждом шаге по времени из решения уравнения влагопереноса. При этом
автоматически учитывается инфильтрационное питание. Поскольку процедура
счета дискретная, то в формупе 7.65 можно полагать йо = ©. В частности,
используя преобразованную (1949), получим
В = — 880 ® Ть
В (В + #Фд) 
{7.65})
где Ро — напор на середине междренья, вычисляемый на хаждом шаге вре-
мени из уравнения влагопереноса; Ко — коэффициент фильтрации; В — меж-
дренное расстояние; Тор = Гд + Рю/2; Т, — расстояние от дрены до водо-
упора; Фд = 0,73То» 2 Тер/ 
я @ (0 — дмаметр дрены).
6. Краевые условия для уравнения теплопереноса.
Аналогичным образом можно записать в общем виде краевые уславия для
уравнения теплопереноса:
на верхней границе при х = ©
вн = (—л‚%Ё+ 6 ут} = а [& - );
7 
{7,66)
на нижней границе при х = (
Вга о о = в [0 — Т, )- 
(т.6т)
При Вт: = 1 получаем граничные условия 3-го рода. В общем случае
1(0 представляет собой некоторую эквивалентную температуру внешней сре-
ды, ад1 — некоторый обобщенный коаффициент теплообмена почвы ©0 сре-
дой. Т; -— температура на нижней границе области. При яя = #2 = О пюл-
‘учаем условия 2-го рода на верхней и нижней границах, когда поток тепла.
за пределы области отсутствует.
При йт = © получаем условия 1-го рода, наиболее часто используемые на
практике.
в,. Краевые условия для уравнения солепереноса.
Как и в предыдущих случаях, можно скомбинировать краевые условия раз-
личных родзв и записать их в виде
при х =
ва (о;‘‚—{{›гогс%} = ва У( - С 
+ @ - вн) (© - ©ох (,
при х = &
—Ве2 (0 с
3;«»&‹.—%}мто/‹:цт—вса)с—цс‹гу
(7,69)
Варьируя параметрами йы (О или 1)}, ссг {0 или 1) и @: {интенсивность
потока солей или минерализации грунтовых вод), получаем разные краевые
условия для уравнения солепереноса на границах области.
Численная процедура решения. Решение системы уравнений (7.35)...
(7.37} в настоящее время аозможно только численными методами на ЭВМ.
При аппроксимации уравнений и краевых условий на равномерной сетке
были использованы неявные, абоолютно устойчивые конечно-разностные схе-
мы, полученные с помощью интегроинтерлоляционного метода {Самарский,
1977). Последукищее решение проводили методом немонстонной прогонки в
итерационном цикле. Для расчета Была составлена программа (совместно с
А.М.Якиревичем, 1979} на языке ФОРТРАН для ЭВМ «Минск-З2» и серии ЕС
(рис. 7.3). Программа позволяет решать как одновременно всю систему урав-
нений тепло и массопереноса, так и отдельно каждое их них. Для этого вво-
дят следующие признаки, задающие режим счета:
93