Для концентраций мнгранта в водоносных и разделяющих слоях Л. С. Рыбниковой составлена система конечно-разностных уравнений при следующих долущениях: в проводящих слоях учи- тывается только конвективный перенос, а в разделяющих слоях. принимается схема сосредоточенной емкости, т. е представляется возможным сосредоточить всю емкость этого слоя в его середине. Для некоторого упрощения выкладок рассматривается профиль- ный (линейный в плане) поток. В блоке номера { / (см. рис. 31) ‘удельный расход мигранта, поступающего конвективным путем из блока #, /—-1: туди-1[быу-лвны--Сы (1— в) ], где 6у-1==!, если оуна 2>0 и ву-1==0, если @;у-1<<0, а мигранта, уходящего в блок, , 1 т [ееубун+-сылн (1—вун) ], где 6у== 1, если о:2>0 н 6у== 0, если ©:/<-0. Расходы мигранта, поступающего в расчетный блок из верхнего и нижнего разделяющих ‘слоев, согласно приня- той схемы сосредоточенной емкости составят: ®“елуАж [е ла, в 61 ) (1 — &) + 20% 4 А (° л, -б / т* ль —ча [ -) (! —) + 20 А (с*, у — 6/ тде 6-==1, если 0/-1*2>0, и бу-==0, если: о/-1у* <; 6и==1, если ©7^ 2>0 н бу==0, если 01/*<<0; Э* — коэффициент дисперсии в со- ответствующем разделяющем слое. Учитывая также изменение количества мигранта в расчетном блоке, запишем конечно-разностное уравнение массопереноса: вч уа вна -6) (1 — &) — т ) |е, у - вуы (1—%, ) + Н0а [с* , Вна-Нд, у (1—%-1)]-20" А (с“ уа -б т* — — ') аке с ) (1 — &)] + 20% Ах, (с*,, у — , )т == == луц х ас ЛМ, (7.23) тде Ас‚ — изменение концентрации мигранта в расчетном блоке за время А/. Записывая левую часть этого уравнения по явной схеме, т. е для начала расчетного интервала времени, получаем возмож- мность_ определить величину Ас,у непосредственно из уравне- ния (7.23). Для устойчивости разностной схемы необходнмо, чтобы коэф- фищиент при ©./ был положительным. Отсюда следует условие ‘устойчивости: УУери В. А В.0 АГА с (7.24) пХ лт лт тде |9| и |9*| — максимальные (по модулую) значения 0 н ©*. Аналогично записывается конечно-разностное уравнение для расчетного блока #, / в разделяющем слое: 9% )А) в; -- е у (1 — 9% — в у (1 — %)] - 20)* Ах, (сны, у - в) — Эе“, |ту* = п т Ах) Ас* , ЛАЬ (7.25) тде Асуу* — изменение величины с.,* за время А. 138°