|
107 Порядок перехода от одного расчетного координатного слоя к другому показан на рис 5.11 пунктирной линией и стрелками. Заметим, что возможна и другая организация расчета с применением (3.108). Ввиду того, что в выражения V, V и прежде всего IV на (.7+1 )-м временном слое входят неизвестные значения концентрации , на каждом временном слое 3 необходимо проводить итерационное уточнение распределения концентрации в сеточной области 2<,К 2 $ Ы $ М ' . При этом в качестве начального приближения временного слоя, т.е. значения р „ к. Рис. 3.10. К построению разностного уравнения Рис. 3.11. Схема расчета В сеточной области, в узлах которой расход жидкости ц подчиняется неравенствам 0 и д у Ы К >0, разностный аналог уравнения массопереноса (3.82) может быть представлен в виде рис.3.12. _ -У _./+1 _ _./♦! _./+1 _ ./♦! Рн.К Ры.К . . .У+1 Ры.К РН-\.К ж, 3+1 Р.М.К Р,Ч.К+\ _ +1 а 7 уы.к--------7 ”ы,к Ат Ах Ау (3.115) Выразив из (3.115) р„*1 к и дополнив полученное соотношение начальным условием (3.110), граничным условиям (3.91), а так же приведенным к дискретной форме граничным условиям (3.86), приходим к разностной схеме |
201 дено значение р[ . в узле разностной сетки, в котором последний раз вы-2 .К.Ч г полняется неравенство ду,чК >0. После этого подобным же образом вычисляются значения концентрации руАЛ для всех К = 2,3,...,К*ы при фиксированной величине индекса N 3, N = 4 и т. д. Расчет по (6.8.45) на первом временном слое У = 1 заканчивается, если встречается координатный слой (обозначим его номер через +1), во всех узлах которого д; составляющая расхода ^x подчиняется неравенству ^xN. к <0. После этого аналогичные вычисления с использованием полученных значений концентрации р\ к проводят для следующего временного слоя У = 2, затем для У = 3 и так вплоть до временного слоя У = У,, соответствующего конечному моменту времени т( моделирования исследуемого процесса массопереноса. Рис. 6.9.1. К построению разностного уравнения (6.9.12)^ Порядок перехода от одного расчетного координатного слоя к другому показан на рис. 6.9.2 пунктирной линией и стрелками. Заметим, что возможна и другая организация расчета с применением (6.9.7). Ввиду того, что в выражения для и, V и прежде всего для Ж на У + 1-м временном слое входят неизвестные значения концентрации р^к, на каждом временном слое У необходимо проводить итерационное уточнение распреде 202 ления концентрации в сеточной области 2 < К < К9 Ы, 2 < N < N'. При этом в качестве начального приближения к р^] к используются значения концентрации с предыдущего временного слоя, т. е. значения р# к. В сеточной области, в узлах которой расход жидкости # подчиняется неравенствам 0 и >0, разностный аналог уравнения массопереноса (6.8.39) может быть представлен в виде рис. 6.9.3. .7+1 _7 + 74 Л7+1 Л7+1 Л7+1 _ ]/•>*У \',К Рл'Ж Рн.К . ..7+1 Рн.К Рн-\.К т/7+1 Ряж РN,К+1_ _ уу 7 +1 (6 9 14) Дг ' Ах Ау .7+1 Выразив из (6.9.14) ры к и дополнив полученное соотношение начальным условием (6.9.9), граничным условием (6.8.48), а также приведенным к дискретной форме граничным условием (6.8.43), приходим к разностной схеме /У -—(о3 +Е о3 1 -Е о3 )+Мг3~хАгРх.К г \Р\\К ' ^хРя-ЪК ЬуРХ'К+ХГ (6.9.15) Р\,к 0; (6.9.16) 7+1 _ Л Р\,К ~ Ро> (6.9.17) |