Численные методы Часть 2
..pdf100 200 с - снижение концентрации примесей на входе в область до
нулевого значения, что соответствует прекращению работы источника примеси (рис. 6.7,а - г);
200 ... 300 с - свободное рассеивание примеси при отсутствии ветра (рис. 6.8,а - г);
300 с ... - поступление с противоположной стороны относительно чистого воздушного потока (рис. 6.9,6.10).
Высота, и
Хорошо прослеживаются тенденции, отмеченные ранее при анализе решения одномерной задачи. Виден процесс заполнения лесного массива примесью при обдувании его воздушным потоком, происходящий с некоторым "отставанием" по сравнению с распределением концентрации примесей в воздушном потоке, обтекающем этот лесной участок (рис. 6.6 и 6.7).
Также отчетливо прослеживается этап понижения концентрации примесей вокруг леса непосредственно при воздействии чистого воздуха, движущегося с направления, противоположного первоначальному. Результаты отчетливо демонстрируют, что при рассмотренном ветровом режиме лес начинает играть
роль вторичного источника примесей (рис. 6.9 и 6.10) с постепенно снижающейся интенсивностью.
При указанных допущениях задача определения распределения скоростей и глубин потока, концентраций примесей для каждого прямолинейного участка сводится к интегрированию системы дифференциальных уравнений в частных производных [19] неразрывности
д а |
dQ |
п |
— + — «0. |
||
a |
ds |
’ |
движения |
|
|
т 1 av |
vav . . |
a s + g a |
+ g d s e | " ,T?' |
сохранения массы (для каждого компонента раствора)
afac,) |
<*Q C ,) |
и - |
а |
& |
’ |
с граничными |
|
|
H(t,0) = H0(t), |
||
V(t,0) = V0(t), |
||
c,(t.0)= c1°(t). |
1 = й |
|
и начальными условиями |
|
|
(6.26)
(6.27)
(6.28)
(629)
(6.30) (6 31)
H(0,s) = R(s), |
<(>33) |
|
V(0,s) = V(s), |
(<• |
V') |
C|(0,s) = C,(s), l = lTk. |
(0 |
31) |
Здесь обозначено: V(t,s), H(t,s) - скорость течения жидкости и |
глубина |
потока; ©(H) - площадь поперечного сечения потока; Q(t,V) = ©V |
расход |
жидкости; i9iyp - уклон трубы (синус угла ее наклона), уклон трения, пщ |
коэффициент шероховатости труб; C,(t,s) - концентрация 1-й примеси, к количество рассматриваемых примесей; t - время; s - продольная координата.
Уклон трения |
вычисляется |
в соответствии |
с |
[19] |
по |
формуле |
|
= vV(sl(H)R(H>). Здесь S(H) - коэффициент Шези, |
S(H) = R>""(H) n„, |
||||||
R(H) - гидравлический радиус трубы. Показатель степени у(Н) вычисляется |
|||||||
согласно [19] по формуле Павловского |
|
|
|
|
|
||
у(Н) = 2 . 5 ^ - 0 ,7 5 # ( Н } ( ^ - 0.1) - 0.13 |
|
|
|
||||
Уравнения (6.26) - (6.28) содержат к+2 дифференциальных уравнений |
с |
||||||
к+2 неизвестными |
функциями |
V(t,s), H(t,s), |
C^t.s), |
i = 1,2— k |
В |
||
граничных условиях |
(6.29) - (6.31) |
предполагается, что жидкость |
вi екает в |
трубу с левого края (s = 0) и вытекает из правого (s = L)