б 
Начальные и граничные усло-
вия следующие:
е(х, #, 0) = с
р(& 0)=ры 
-- рв — Ра!-т; ›
аее 
е(0, 
й, ) =с, 
при 0ж у<,
7[1 
; 
е(0, й, ) 
—р при 
< у<,
ы 
ее 1.0
//Ё 
е 
ей» в д =
дс
ор ( 
0 @
Рис. 5.1.1. Моделируемая область.
1 — ось сныметрии, # — неточник загряз- 
Р (0, ) = ру»
э( д— ро
тде с, и с, — постоянные, представляющие собой концентрацию
загрязнения в источнике и начальную концентрацию загрязнения
в моделируемой области соответственно; ру и рр — постоянные.
Эта система представляется в конечных разностях согласно
неявной схеме метода переменных направлений (НСМПН). Кон-
вективный член представляется разностями назад по отношению
к направлению потока. Плотность р и скорость @ ЯВЛЯЮТСЯ По-
стоянными в течение временного шага. Возможная схема дискре-
тизации имеет вид:
С3Я 
[1'(!—|—1 /)|И(1+1/)|+М‘ ) и(' Л {0'1! 
—( 
!)]
Р(ч 
о@, ))
Р( !)\д(і-і)Н-Р(‘—Ъ/)ШН—!_... е( :*(ійі.п} Ч
2 
Р@) 
р@ -1)
е(& 7 
) -НоМ, ) [ е( 
7+1) _ е( ))
ЕЕ Пк НЕНЕ 
НЕВ
—а( 
р #©.0-н06 
10 [ ор 
— 610 _
«& )| 
э 
[р(‹‚і) 
Ф@ 
|›Н
—н рен( рн ее ) 
_ 
еей, р -ей ).
ах 
: 
а() 
*
ВИЗЕТЬ ЛЬ ОЗЕТЬЛЕЕЕР(Ь 
Э1ЕФ ОЦЕРЕНЕЕ 
Й 
ег /›}
‹іх* 
РЕНЬ) о))
_Е(!-!Н”(!›і)Нр(і—1‚і)|И(і71‘і)| е@р 
_ 
@-) ]
2 
В 
лв 
}] +
@ 
гр 
оС1+НО ННр 
7 
131) — ей, )
*@ЧИ(»‹ > 
2 
е
& 
р| 
2 
Л +о 670 [
— и( ...———‚—'——'{
— #(& ре* (
р 
`@ (п)/2 
ы