Патанкар С.В., Сполдинг Д.Б. - Тепло- и массообмен в пограничных слоях (1062125), страница 27
Текст из файла (страница 27)
В примере, приведенном в П!-5, предусмотрено печатание величин абсциссы х н соответствующих профилей скорости и температуры через каждые пять интегрирований. П(-4. КОНТРОЛЬ ПРАВИЛЬНОСТИ СОСТАВЛЕНИЯ ПРОГРАММЫ Приведенное нами выше описание деталей вычислительной программы н сводка отдельных подпрограмм, представленная в й П!-5, составляет информацшо, достаточную лля самостоятельного выполнения машинного счета. Однако после приооретеиия общего понимания программы расчетчик может не пожелать использовать всю эту информацию целиком всякий раз при составлении заново вычислительной программы для решения конкретной задачи.
В этом случае нужна уверенность в том, что соблюдены следующие десять операций (они выражают ту же самую инфориацюо в сжатой и удобной лля использования форме). !. Полг~рограммы СОЕГГ, 5).1Р, 50(лгЕ и ИЕЛОУ сохраняются без изменений. 98 2. В аодпрограмче МЛ1Х определение длины шага (ха —.сс) н условия для окоачаняя вычнслсчии выбираются в соответствии с поставленными целями. 3. Выбирается форт:ула для эффективной вязкости. В случае приехглемости гила.
телы о пути смешения подпрограммы ЛгЕГГХ 1.ЕХСгТН, ЕХТКХ, %ЛЕС, )РЕ(, Л4г)т2 П!-4, ВЫх!ИСЛИТЕЛЬНАЯ ПРОГРАММА В этом парагрзфе мы приводим но.шый текст вьшлслительяоп программы. Подпрограммы, за исключенном МА!Х, которая приводится первой, расположены в алфавит поч порядке. МЛ!Х СОЛ!МОХ(0! Х,'РЕ!,ЛМ1,ЛМЕ,ОРОХ,РКЕЕ(2),РК(2),Р (2),О!ЗХ,ЛЗ(0, !Н,Х!З,ХР,Х!..ОХ,!ХТО,СВАЕЙ ! ~!)Х,ХР!ЗОР2 ХРЗ ХГэО ХРН КЕХ КЕ;,КЛЗЕ !хКЛ(Э 1(В)ВГТЛ САМА(2) ТЛО).ТЛОЕ,Л)! (2),АЛэ (2) .! ХО! (2),! ХНЕ ГО) 1,''эГ(П(43),Г(2 43) К( (3),КНО(43),ОМ(43) А (43) !(С;ЗС(43),ЛО((3),ВО(43),С(1(43),Л(2 43),В (243),С !243) СОЛ!й!ОХг!1.УЛК,Л~.З!61 !6 Продолжение !ХТО =.0 ХЕ =.
20. СЛЕЕ СОХЗТ СЛЕЕ ВГ01Х ЛМ1--0. ЛМЕ= О. Переход к 25 !б СА( Е КГА1)У 25 Продолжение 1ХТО=.!ХТО+ ! СА! Е ).ЕХОТН САГЕ ЕХТКХ С Выбор предыдущего шага РКЛ =-..05 ОХ=ГКЛаРЕ1/(К(!) аАМ! — К(ХРЗ) аАХ!Е) 1Г(ПХ.ОТ..ЗФУ((ХРЗ) ) ОХ=.5ФТ(ХРЗ) ХО=ХО+ОХ 77 Продолжение СЛ1.1. РКЕ(ХУ,ХО,ОРОХ) 1Р(КЛЗЕ.Е(3 2) Перехол к 26 (Е(К)Х.ЕС). ! ) САП.
эТАЗВ(ХП,ХО,ЛЗ!1) ' Обычно трактовки граничных условий в подпрограмме 5!ЛР оказы статочно, но при рассмотрении сложных граничных условий необходима модификация. Ть вается до- некоторак сохраняются неизменных~и. Для любой другой гипотезы течении этн подпрограммы модифицируются соответствующим образом.
4. Выбираются величины ККАО, ХЕО, К1К, КЕХ, Х и вволгмся в подпрограмму ВЕО)Х совместно с началыгымн значениями х н начальными профилями. Выясняется возможность согласования типа границы (т. с. К!Х и КЕХ) с характером движения в рассматриваемой задаче. Если ихжны изменения, то велячипы К!Х али КГХ вводятся в подходящем месте програмиы. Следует изменлть величины б и у ВЕ0)Х, есля течсмнс нс турб)лентно. Нужно, насколько это возможно, ограничить все команлм КЕАО, относящиеся к этой подпрограмме. 5. Подпрогрзмма ЗОГ ВСЕ приводится в соответствие с используемыми От-урзвне. пнями.
Если гранина совпадает со стенкой, то подпрограммы ЕВС н МЛ65 рассматриваются таким образом, чтобы обеспечать задание желаемых Ф-граничных условий' и шпепсиапостн массаперепоса через стенку. 6. Провалится подготонка для вычислений соответствуюш,ш свойств жидкостк а подпрограммах ОЕХЗТУ л Лг!БСО. ). Подготавливаются подпрограмма РКЕ для задаяия величины продольного градиента давления в соответствии с условиями задачи н потпрогрзх:х|а КАП, вводящая информагппо о геометрии !чраннцы. 8. Нужш1я гля расчета информация о постоянных вводится в СОХЗТ. О. Прел!сх1атрнвается вывод и воспроизведение (печзтаннс) нужных результатоа в ОСТРОТ, !О. Потготавлнвастся перфокарта с ланныхги с указанием последовательных ксг.
хганд КЕАО па всей программе. 100 1Г(КЕХ.Е[4.1) САП. МА55 (Х[],ХО,АХ[В) СА1Л. '«т'А1Л. 26 САЕЕ ОСТРОТ САЬЕ СОЕГГ Подстановка скоростей иа свободных границах 1Г(КЕХ ЕЯ2) [!(МРЗ) =ЗЕЙТ(Б(МРЗ) а[)(МР[3) — 2е(ХΠ— ХН) е ОРОХ/ЙНО(МРЗ]) 1Г(К1МЩ 2) 1 «(1) =-ВЩТ(С 61) еС ((1] — 2 е (ХΠ— ХС) еОРОХ[ЙНО(1) ) САЕЕ ЗОРЧЕ(АС,В[),СО,[),МРЗ) Подстановка скоростей на линии симметрии 1Г(К[ММЕ.З) Переход к 7! [Л (1) = [) (2) [Г(КЙАО ГЯ О) Н(1) = 75е[)(2) «- 25еУ[3) 71 1Г(КГХ Е[ ! 3) !!(МРЗ) = 75жН(МР2) -1- 25е[)(МР[) 72 [Н~[=[)(!)жН(1) [)13/ Н(!) жЕ[3) [133= 1) (3) еС [3) ЗЯ = 84.
ч: 0 1 ! --12, и [) 13 Ч-9. е НЗЗ [)(2) = (16 жН1! —.! ж[)!3+[)33)/(2е([)(1)+ С(3))+5[ч«ЙТ(5©) У(2] =Т[З) ж [В [2]+[)(3) — 2аС[!]] .5[[[) [2) -,[) [3) + [[[!]) Перехоч к 74 73 1Г(КЙА[).МЕ.О) Переход к 89 [)(2) (4 е С [!) Н«(3))[3 '«'[2) =О. Переход к 74 89 [)(2) = В (1) У(2) =- 7 [3) /3. 74 Переход н (75.76.77) КЕХ 75 0(МР2) =1)[МР1)[(1.+2.ФВЕТА) У(МР2) =Т(МРЗ) — (т(МРЗ) — У(МР[) ) еВЕТА/(2+ВЕТА) Переход к 78 76 [Н1=[)(МР[) еН(МР1) [) ! 3 = [) (М Р ! ) е [[( М Р 3) [)33=[)(МРЗ) е[)(МРЗ) 59=64 Ф[)33 — 12 е[)!3+9 «(Н 1 ЩМР2) = (!6.и[[33 — 4.з[)13«-[]!1)[(2.е([)(МР1)+1!(МРЗ)) +ЗЕЙТ(БЯ)] У(МР2)=т(МРЗ) — (У[МРЗ) — У(М[Х1))е([)(МР2)+[)(МР!) 2еУ(МРЗ))е5/ ! [Н[МР2) +[) [ХР!) -,Н«[МРЗ) Переход к 78 77 [)(МР2) =. (4,ж[[[МРЗ) — 1!(МР1))[3 У(МР2) =у[МРЗ) 78 Продолжение 1Г(МЕЯ.ЕО.!) Переход к 45 []О 88 3=1,МРН Начальные проф~««и друга«зависимых переменных ЙЕАО(5 43) Г(3, !), [Г(3,1),1 = 3 МР[),Г(З,МРЗ) Вычисление соответству«о«цнх величин скольжения ОЛМА(3) =.143 Переход к (81,82,83),К1М 81 Г(32) =Г(3;!)-!-(Г[33) Г(3,1)) а(1.+ВЕТА — ОАА[А(3))[(1.+ +ВЕТА«-ОАЗ[А(3)) Переход к 84 82 О= ([](2) + [)(3) — 6 Ф[](!))[(5 е ([)(2) 4-[][3)) + 8 е[)(!)) С«Г=- (1.— РЙЕГ[3) ) /[1.
+РЙЕГ(Л) ) СГ= [О ф ОГ) [[ 1, '- О С«Г) Г(3,2) =Г(3,3) еОГ-~'-1.— ОГ) еГ(3,!) Перехо ! к 84 83 Г(3,2) —.Г[3,!) 1Г(КЙА[) ЕОО) Г [3 2) - [ч,жГ(3,1) — Г(33) )/3. 84 Переход к (85,86,87) КЕХ 85 Г(ЛМР2) =Г[ЗМРЗ)+(Г(3 МР[) — Г(ЛМРЗ))а(1.+ВЕТА — САМА(3))/(1.+ !ВЕТА+ПАХ[А[3) ) Перехоч к 88 86 О=([[(МР2) 4[)[МР!] — 8 жИ(МРЗ))/(5е([)(МР2) +Н«(МР1))+87а[)(МРЗ)) ОГ = (1.— Р ЙЕГ [3) ) /(1 + РЙЕГ (3) ) ОЙ= [О+ОГ) [(1с-ОеОГ) Г(З,МР2) =Г[1,МР1) аОà — , '[1.— ОГ) еГ(Л,МРЗ) Переход к 88 87 Г(З,МР2] = (4.жГ(Л,МРЗ) — Г(АМР!))/3. 88 Продолжение 45 Продолжение СЛЕЕ РЕХЯТУ Вычисление радиусов ЕЛЕЕ КЛР(хи,К(Ц,СЗАЕГЛ) !Г(СЗАЕГЛ.Щ.О.ОК.ККЛР.ЕЯ.О) Переход к 27 РО 28 1=2,ХРЗ 28 К(Ц = К(Ц 4-У(~Ц вСЗАРГЛ 72 Продолжение )Г(ХьЕО.ЕО.Ц Переход к 30 РО 45 3=1, ХРН РО 46 1=2,ХР2 АП(В =А(З,Ц ВР(Ц =В(З,Ц 46 СБ(Ц =С(З,Ц РО 47 1=1,ХРЗ 47 ЗС(Ц =Г(З,Ц СА1.!.
50!Л'1'. (ЛЦВ(/,С!!,ЗС,ХРЗ) РО 48 1=1,ХРЗ 48 1-'(З,Ц =ЗС(Ц !Г(КАЗЕ.Е0.2) Переход к 81 Подстановка !'-величин на стенке (Г(к(хее(3)лхР!хР)(3).е1/2) Г(ЗЦ = (((.—,ВетА24-Олй(А(3) ) вГ 1(32) — (1.Ф ВЕТЛ вЂ” ОЛМА(ц) вГ(1 3)) во/ОЛМА(3) 1Г (КЕХЕ01 АХР(ХРЕ (3) .ЕС) 2) Г(3 Х Р3) = ( (1. + В ЕТЛ+ ОЛА!А (д) ) вГ ! (З,ХР2) (1,+ВЕТА — ОЛМА(!)) вГ(ГЗ,ХРЦ) ж.й/ОЛМА(3) Подстановка Г-величин на,зинни симметрии 8! Пт(К!Х.ХЕ.З) Переход к 82 Г(3,Ц =Г(3,2) 11з(ККЛР.ЕЯ.О)Г(З,Ц =-.754Г(3,2) +.25:нГ(33) 82 1Г(КЕХ.Е6)З)Г(З,ХРЗ)=-.754Г(З,ХР2)-1-.25еГ(З,ХРЦ 45 Продолжение ЗО ХР=Х11 ХР=ХР РЕ1=РЕ1+РХв (К(ЦвЛЛ!! — К(ХРЗ) жАМЕ) Г!редельное условно !Г(ХШТ Х1.) Переход к 15 Переход к 16 Окончание ВЕОГХ Подпрограмма ВЕ01Х СОММОХ /ОЕХ/РЕКАМ(нйй!Е,РРРХ,РКЕГ(2),РК(2),Р(2),РЕХ,АМП, !ХО,ХР,ХР,ХЕ,РХ,! Х ТРОС ЗАЕГА 1/!/Х,ХР1,ХРЗ,ХРЗ,ХЕ(чт,ХРН,КЕХ,К1Х,КЛЗЕ,ККАР 1/В/ВЕТЛ,ОЛМА(2),ТА(/!,ТАЙЕ,ЛЗ!(2),АЗЕ(2),)ХР1(2),1ХРЕ(2) !/Ъ/0(43),Г(2,43),К(43,КНО(43),05)(43),У(43) Характеристики задачи КЕЛР (5,42) ККЛР/Х ЕЯ,КЕХ,К1Х,Х 42 Формат 411, !2 КАЗЕ=2 (Г(К!Х.Е0.1.0К.КЕХ.Е(З.ЦКАЗЕ=! Х(!=О.
ХРН=ХЕΠ— 1 ХР1=Х 61 ХР2=-Х+2 ХРЗ=Х+3 Начальный гзрофиль скорости КЕЛР(543)У(Ц,()(Ц,У(Ц,(/(Ц,1=3,ХРЦ,У(ХРЗ),0(ХРЗ) 43 Формат (ОГ!0.0) Вычисление скоростей схольжсиия и расстояний ВЕТА=.!43. Переход к (71,72,73) К1Х 71 (/(2) — 3(3)/(1.42.жВЕТА) У(2) =У(3) вВЕТА/(2.+ВЕТА) Переход к 74 Перехот к 29 27 РО 30 1=2.ХРЗ 30 К(Ц=К(Ц 29 Продолжение Вычисление ы-величин 101 ом(ц =ю. ОМ(2) =Ю '00 49 1=3,ч(Р2. 49 ОМ(Ц =Ом(! — Ц +5з (РНО(Ц вП(Ц в К(1) ГКНО(1 — Ц:з ГР(! — Ц*К(! — Ц) (У(Ц вЂ” У(1 — Ц) РЕ!=ОМ()СР2) 00 59 1=За(Р! 59 ОМ(Ц =ОМ(Ц/РЕ! ОМ(г(Р2) — — 1.
ОМ()4РЗ) =1. !Г(5(ЕС.Е(З.Ц КЕТ()КК 00 69 3=1, ХР11 1Г(кех.е(З ц ес Ре(1) =-! !Г(КЙ4.ЕО.Ц (ЙР!(3) =1 69 Продочженнс Возврат Окончание СОЕГГ Подпрограмма СОГГГ СОММОзц/ОЕгК/РЕ1,ЛМ1,ЛМЕ,РРОХ,РКГГ(2),РК(2),Р(2),РЕз(,АМО, !ХО,ХР,ХР,Х1.,0Х,!ч(ТОСЯЛ1.ГА 1/1/!4,х(Р),ЧР2зч)РЗ,)(Е©КРН,КЕХ,К(К,КАЯЕ,ККЛР 1/В/ВЕТЛ,ОЛЧЛ(2),ТАЙ,ТЛ()Е,ЛЗ! (2),АЗЕ(2),1140! (2),1К(РЕ(2) 1/Ъ/()(43),Г(2.43).К(43),КНО(43),ОМ(43),У(43) 1/С/ЯС (43),АГ) (43),В(1 (43),С П (43),Л (2 43),В (2,43),С (2 43) СОММОК /Е/ЛК, А!.МО Размерность 01(43) т02(43),03(43),Р(2,43) Я1(43),52(43),ЯЗ(43) Вычисление малых С-везнчин РО 98 1=2с4Р! ЙЛ= 5 з(К(! е Ц -!-К(Ц) РН =.5:з (КНО (14-~Ц + КНО (Ц ) РМ=.Зз (()(1 г Ц ГР(Ц) САП.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
