отчет Югов и др Сеточные методы расчена НИР (1047878), страница 13
Текст из файла (страница 13)
'.-'»»Е — Г! 1 ."-', 4 4Š— Г!'=, '-'.:-'6Е+ !'.!;:.' Й, !.лГлЕ-»»1 .л, Г!Г»Е+Г!"л ,",. ЙЙЕ+О.! й. »лС»Е-й1 1. !.»1»Е+ОО ".1-,— й ~ О„Г!цЕ-Й1 1 „: '!Е- »!, ,;1- »Г! „'л, ~~ -'Š— О1 ».1ОЕ-О1 Э. ЙГ!Е-Г!'- ,л. '"" „'"Е+О 1г 1. 41Е-О!» '.4. 44Е+Гл~~ «.;,гЕ+Г!-, 9. 7 г Е-О!1 ГЕ-О1 4.О6Е-О" .. 4':-'Е+ГлЯ 1О! 1. л»»лЕ-1-!: 9.
1; Е+»1л5 ~: 3., 1 -,ЙЕ-Г!1 «»,.Е -»!~, '... 1-+~ !,.' 1.:..'6Е- О,=' !",', 1»!-!Е»О! ! О. 4'.'Е+й,-' 9.'1~ Е-О1 „''. О»„1 -Й1 л, ОЗЕ-Г!'.".' л. 35Е+!.'! 1, ' ЙЕ-Й.' "«.,«! ! »Йлс, Я', 'ййЕ+Гл:- :л. 46Е--О1 ='. 4ОЕ-О1 ;. ",1Е»»л. 1. л5Е-Г»л:й '~. '-",Е+ Г!, ' 1Е+О'=' '=:, 71Е-Г»! 1. 5'=:Е--Г!1 4.О2Е-О ,,4ЙЕ+О=' 1Л 1. ":О!Е-Г!'й 9.
51Е+»л~ '„-3, !.!!"!Е+Й'.' О 4Е-й1 1. Й1Е-О1 4. 11Е-О л "'. 44Е+О ' 14 й 1. О..'Е+й6 9. !'3Е+Г! ' 1. ООЕ+Слй й.ООЕ-й1 .;. «О!Е+1!. ' !!»;И1 ! л. "ОЕ+О':» "л. 4 "5Е»Г!::.' , 4! »Ел О„,' 3. "'5Е»»:!.=' .. л!»лЕ»6! ' ~Е+1! -' !.!Е+!.! ' 1' лЕ»! !„* 1»Р1»Х 1. 1йЕ+й-. 1. !.!' Е'+!" ' 1, ЙОЕ+!, ", 1, »»; 1-+»! ', 1. !.":Е+~.'! '.'" 1, !»-',Е»»», 1. »»4Е"-»!'"' 1. »::! 'Е»1!.-: Т<.И1 ! «». ~;«лЕ+Г»О 'л.
»',ОЕ+» '-!. ~-, -'Е+»! й 6!, »',4 1. +! »,.! "! „~.'-',6Т +!'! ' Э. 6 'Е+!".1;=' '» «,.~Е»1»," '3. 7Г»Е+»,'! -' ;::!! :-'. »ГлЕ»Й1 ;', „~~Е+О! 1 »-"!'! - ф й ! »»л ! О!» С: 1' '-'. О!.:!" » !.'1 '»»! И1 '! й 1, »6».!Е'-!!.' 1. » »- Е-»»!.:.' 1, » ':Е-О! ' 1. О»::1"-О '-' '-:. »л4! -О ' . 44Е-!.'. " .:-'6!Е-»!.» !'-'!11ЙС'Н ".!.
64Е-О»1 ~!. " «ТŠ— Й1 г '"'Е-Й1 »!.: ~.'лЕ-О1 «,.' 'à — й1 ':, ! "'Е-О1 '. О !Е-Й1 »»» !!!' !.!"»»»»Ней! ие!"'е !!йн!» ."!НГ»чен!1е 611»!Р!»!. т!! !Гтйлг' Г»тр!1цйтельн1!11»»» О -'Е-О ' ,'= й Е»Г!'! О» Е+О.' »л»!Е-!1! ! !"»! 1Е+ Г!» ! -'~-,Е -Г! ' ='»1»О»1 4. Г!1Е-Г12 !='.:-'.»Е+Г»5 ».,=!/Е+Г! ' ~р - «»Е »-.! 1 «,Š— О! 1 ;". „:"4Е-с!» :-'. =,:6»Е+Г!'» 4. 61Е-Й» : . 1 ГЕ л О!"' « .:=:»л Е+»! .= '.-.'. 76Е-Г»1 '.;„Г!»лЕ-!»1 , 41Е-»л5 "Е+О!'.' !О!. ЙЙЕ-!.!1 ,' „.'»ОЕ+»й!";. ".
':-"=!Е+О:- »!, О!ОЕ-»31-. -' 1, »»! !Е+Ой 6Т.'-й', '!. !'!»'4»Е- »,»; -!.= 4, 14Е-!!:..' 4. "1Ел !.!'" 4. 5:.,!Е+»! ' »= .:"Е- '~.»1 '='. '31Е-!'.'1 «1Е !)" 1. л»,Е»!."»; 4. »ГЕ-~".2 . ' Г»1-"-- »! 4, О'='Е-й -' О.:37Е+! !'. ~-',, 4;Е» »л,-' '='. '5'.:»Е-!" ! .: 1Е-! 1 »Е-Г!',.' -'„»;, »=, Е+ О! ':Е+О! ,»:.Г-»»1 1 й , Г!;-'Е-й '»!. О»!Ел й"~ '.!»-,'-,! -;" ! 1» :. 14Е-Гл2 '3 .:-':=, Е+ !.!" ;.,:=, . Е+Г!; '3. 'Е!Е-Ол1 1.. »Е-й1 4, О!..'Е-О~: «»~",Е-Г! ' '„-',, ',Е+»!' ~ := .
""4!'='+!"'' : '. "." Е-О1 л. ' '1Е-!" ~ 1. О'='Е-й:=' 'л, ! 'Е+»»,'» ''-', ' ' Е+О" 4. О1Е-»!5 НЙ ! ГИ~ !!, !..!Е ТЕТЙ .»1ЗС !.!':; 1ЕЙР ! !.«! !Е ТЕ 1! й !! ! Н':- тейг ! !,'! !Е ТЕ Т'Й !! ~ 3 Н'1 ТЕ!1Г ц ЦЕ тЕТЙ ! !1'-,!' ,! 1 ) Н':-' ! ЕИР !!:!!1- !Етм ц '1 *-",!' 3.69Е+Ой "л. г 1Е+ОО .л.
» »Е»Гл-' "'. 74Е+Г»! =' .л. »6Е»й=' 3. 7 Е+О«' 3.79Е+ОО т ..1:! 3. ОГлЕ+Гл'-' .л. Е1Е+Г»= :. ОЗЕ+Г! =' .л. 34Е+ОО .л. О~»Е+Г»2 :. 6!6Е+О! й , . Г!'»Е+О 3. ЕОЕ+О! «' ТКЧКНМК В 0СКСИ!ИКТРМЧНСМ НЛНЛЛЕ На рис.П,4.1 представлена схема т~чения в рассматриваемом КанаЛ6В а На РИС.О, 4за: ПОКаааНО ПРЕДПОЛаГаЕМОЕ ИЗМЕНЕНИЕ СКорсоти пс длине кан6ла ° Геометрия канала определяется величинами УМ1 = У~М1)- иирина ~8~8~8 по н~р~~~~ к нентральнсму телу и Ц1) — радиус пентральнсгс теле, В задаче мсжнО выделить дв6 части Часть 1 : МР/дХ - неизвестно ~ 127РРФ = .РЛ1$Е.
), озона 2 : оР/ох — задано а Ванно 0 о хоРЕ = .тВое. ). Ба Хннонин 25 Р Х> 20 ооононоаан нанн нана: Тцн3к= 5 температура вдуваемого Газа, ЙМ1 = М~И = 0,1 Й~1) расход вдуваемого газа, О „„ = М С Т„„ тепловой поток, обусловленный вдуваемым Газсмт Сз = 1 теплоемкость вдуваемого и основного газа, Число Нуссельта Ии= ЕЛ3./Л , где 1. - характерный размер,~в программе ЦМ1)-Ц1) ) ,Л вЂ” козКипиент теплопроводности, который можно определить из числа Прандтля Рк= ~РСГ/Л ", т.к.
~х =1 и Ср=1тто цЛ = 1/Рх, к - козФфипиент теплоотдачи на поверхности стенки, ~ = ФДТ,~- Т~ ). Т.о. ~исло Нуссельта определяется по следух~ией формуле аач = Р Где Т~= ТВ - среднесбъем- ная температура. ~МВ СГТду ТдФ ТЕ2= 0н= нонна Т = + = ~Т дзн- ~ТВВ 0РЕТ такое простое выражение для Тн, является одним из преимуществ Х Ы КОординат. В р~с~матрнв~ем~й ~~даче ПЛстясотва д~~а~~ч~~~ий Ксаффнпнеят вязкости, число Прандтля и число Ымидта: = 1, ~4, Рх =0 Ьо= 5.
ЙЕ~36МЕЛНЫМ ЛРИСНОИННВтСЛ СннлДУВЩИЕл ЗНОЧОНИЛ". ЧЬИХ = 10 ИУ = 11 Нйиб = 12 И9М = 15 ЫТОР = .ГАЫК. Логическии переиенным ~ЗОНЕ , соотнетстнувщии ИР = 1-11, приснаизавтся сдед~вщие значения ЫО1МЕ~1) "- .ТИМЕ, Р) = .РАМЕ. ~10) ~11) = .ТВ)Е. т.е ° б1ед~т рещатьсн ~ранненип ДДЯ Т~ Й Р-"1)в 0 ~ МУ=11)и Аогическии переченньм 1.РйМТ состьетстнувщин ИГ=1-11 приснаимавтся одедувщие Оначения: ЬЛж'~1) = .ТВОЕ. Р) = .РАЫЕ, ~5) = .ТЙУЕ. ~ч) = .ыы. ~л5) = .ТВЕ, ~б) е — РА1 ~Е ~10) ~11) = .ТРК. теи сан,:и подпрограмие РйОК1. предписнняетсн нынод на печать величин т',э), т~э~ - т~и1~, иле~, ичв) ТБ - Б„'41) МАХ 15ТКР = О АН1~11) ЫЬТЕРии 1000 ЫБХ1~11) ХО =О АГХК~11) лЬАЬТ = 10 ВУХЦ11) ОИЫАли 1 К301~117~ Р511 = О :Ь' К~11)) ,40'ев ~ТаХ = 10 А3!лЦ 11)~ РОЖЕЙ = 1 Й::~~1 = О Й,.1Е = О 1ЮД11 РО1'РАМИА И1Р 1.
Число уплоных точек М1 = 12 . 2, Ра11иус нентрального тела 1~1) = О,5 5. Осесимметричнап геометрия МОРАМ = 2 . 11, <~Ру'дХ вЂ” неианестен, но ианестна геометрия копала ЮР0=,ТА1.5Е. 5» ~Максимальное число итерапий 1НЬ11.'1 '= бО б, Ра:30инка Области на контролънме ооъемм Осуйестнляется подпрограммо11 'мУЗк111, для этого аадается аанисимость у от 3 н11бираетсй ли11ейпая аанисймОстн у =: — ). 3-1 й1- 1 7. Задается исходное распределение величин: плотнОсти ~ , скоростиУ, температуры 1, коннент11апии ндунаемого гана О 1 рис. П.4.3 ) ЧОД1) РОРРАММА ВТАЙНЕ 1. '1емперагура ндунаемого Гапс 11Ю =-5. г. 11опраака на 110ток массы чорез Гранину Е 5. 1'ронины 1 й Е 11эл11ются стенками К1И =1, ХЕК=1. 5- 1 5-1 2 2 5 5 5 б б 7 6 6 9 Ч 1Π— о 1 1О 11 о 1 11 12 12 О 1О ао11 о 5 5 б б 7 7 6 6 з з 1О 1О 11 Р— ~1 11 Π— С» 12 б — Π— -Ф Т О '~ ~ б Рис.
11,4,6 1 Ьякисляетсл максимальиая скОрость н канало М,,ЬИ = еах У~д), 2 ° .:!мкислнется среЛнеобьемноя темнература Тй = 1~ = 14ьз ° 5. '1ереа к'.анно Леснть йавов вмамааетсн Яолиротрамма РВОГ1). и Го *. Оьктсл к н' Кати следуюлис воли Б~н 4 'И~'Д '~д )/Е11), М/О:,~ЛХ. ГН Г~ „,11) на каином оаво амволятся на некать величина: текутее анамени"' Х ЬМХ дР/дХ 1:1) температура Контрольного тела, ЛИО Кисло Р~~- с4ъ тын или длн среднеобъемная температура~ Й1) конментраиия здуааемоио оааа, К~,;~.1 ) внутренний ралиус канала или для ХЕХ=2 ралиус свооолной транииь~, максимальная скорость В канале, 1~ОЬ11РЫ'РА,ИА 6ь 15СЙ 1.
Если реп аетсн .1-оо уравнение',май=1 ),:~=Т,то ГД) = 64ЬЭ) = 1/Гх, 2 ' ~ '~и ран~а'т~'и 2-ое '"равнение ~мЬ= ) '~=6 то Г1Э) = ЙЬ'2,3) - "1/Яо. 5. Если ре~иетсн 11-ое уравнение ~1%=11 ), Ф=М,то ГД) = ВЛАД) =~4, ~А = 1. РаооматРинаемое течение мойнО 1разделитй на 4 Нонн, показан- нме па рио, П.4.4. 1е Зона 1 1.1,Ь зоне 1а аХ = 0,5, 1идХ = 0,5, 1о ьХ = 1,2 .
1е2ейирина кана~а по нормали к централнному телу УМ1 = 1. 1.5вДлл уразненил ЙР = 1 ~4=7)р на гранино 'й задана температура ЦЯ1)=1 а на границе 1 задан Обцийа поток .г' =0 ~адиабатичеокал отенка), 2. Зона р! Н.р.наг оо носрннааге Х высврвегон равенн аа=р,р, 2е2аЫИРИНа КНННЛа ПО НорМЯЛИ К ЦантРНЛЬНОМУ ТЕЛУ УМ1 = 1 + 0 1~ ХО +,4Д 20 ) 1.5.Поток маооы через границу 1 принимаетол разным ЯЧ) = 0,1 К~1). 5.5.В отои зоне будет рецатьол ураннение длл М~ = 2 ~~=С) и БЫВОДитЬОЯ На ПЕЧатЬ ВЕЛИЧИЯМ С И С/01е 5.~.Длл уранненил МР =,2 ~4р=С) на границак 1 и Е заданы обцие потоки АК'1~2) = 3 = -' ~- р ЛУХЦ2) =3 = О.
Ц1) 5.5.Тепловой поток на границе 1 ранен ЛН1~1) =3, = -~-)- Т„„4, И~1) 5» 5.1. 5.2. «. Зон» 1У «,1. Нон оо ноорннонн» 1 »НО»ОНО»ОН ООНННН Нх = 1. ~~.2. Границей 1 Явлнетск ось симметрии К1М = 5 и $~1) = О. 5. ~'.ели ~ра~ица Е пвлнетск ~~обод~ой границей, т ° е ° зоны 111 и тс поток массы чсрез границ~ Е считается слел~юции Образом' 5.1. Расчитываетск теоретический поток массы ~ З н нн М~М1)-У(М2) 5.2.













