108 (3.116) р°ы.к= о ; „./♦I _ „ Р\.к ~ Ро-> (3.117) (3.118) ЛЙ.-0, (З.П9) где Е\ = \ + Ех Еу. В частности, для узла N=2, К=К3-1 при .1=0 разностное выражение (3.95) принимает вид л*,.,= . . + м * , . *Лс,) ■ (3.200) Рис. 3.72 Последовательность вычисления В (3.200) известно из начального условия (3.116), р\к^ из граничного условия (3.117), а р\к^ из граничного условия (3.118). Алгоритм расчета поля концентрации, использующий разностные уравнения (3.115)-(3.118), подобен описанному выше. Отличие заключается в том, что для каждого координатного слоя N вычисления начинают с узла |
204 Алгоритм расчета поля концентрации, использующий разностные уравнения (6.9.14)-(6.9.17), подобен описанному выше. Отличие заключается в том, что для каждого координатного слоя N вычисления начинают с узла разностной сетки с номером -1 и переход от узла к узлу ведется в направлении, противоположном принятому в разностной схеме (6.9.7), (6.9.8), (6.9.11), (6.9.12), т. е. последовательно перебираются значения К = К3 -1, К3 -2,..., вплоть до номера такого узла разностной сетки, за которым происходит изменение знака дуЛ>к. В узлах сеточной области, характеризуемых условием Л. к, справедливо равенство 7 + 1г\ РЛ’,К ~ т. к. в плоскости водоупора х = х3, р(*3,.у>г) = 0. Рис. 6.9.4. К построению разностного уравнения (6.9.20) Для каждого значения коэффициента кинетики Я при шаге Дг, большем некоторого в силу того, что IV^ <0, значения р^к, вычисленные по (6.9.7) или (6.9.15), оказываются отрицательными. Это приводит к дополнительному ограничению на шаг Дг, которое при больших Я может быть достаточно жестким. Если в уравнении кинетики (6.7.6) г(р,е)=лр([-с). |