109 разностной сетки с номером К3-\ и переход от узла к узлу ведется в направлении, противоположном принятому в разностной схеме (3.108), (3.109), (3.112), (3.113), т.е. последовательно перебираются значения К-К3-1, К3-2,..., вплоть до номера такого узла разностной сетки, за которым происходит изменение знака д> ы к. В узлах сеточной области, характеризуемым условием дхЫ К, справедливо равенство т.к. в плоскости водоупорах=х3ур(х3%у,т)=0. Рис. 3.13. К построению разностного уравнения (3.201) Для каждого значения коэффициента кинетики / при шаге Ах, большем некоторого в силу того, что <0, значения р^] к, вычисленные по (3.108) или (3.115), оказываются отрицательными. Это приводит к дополнительному ограничению на шаг Ат, которое при больших к может быть достаточно жестким. Если в уравнении кинетики (3.28) То в (3.87) р(р4)=1р(1-<Г), В этом случае, разрешив (3.107) относительно имеем (рис.3.13) где Е' =\ + Ех + Еу + (\-с)ЛАт. (3.201) Аналогично из (3.94) получаем |
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) г(р,е)=лр([-с). |