Fletcher-2-rus (1185919), страница 46

Файл №1185919 Fletcher-2-rus (Флетчер К. Вычислительные методы в динамике жидкостей) 46 страницаFletcher-2-rus (1185919) страница 462020-08-25СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 46)

На стенке и внешней границе пограничного слоя Таблица 13.4. Параметры, используемые в программе Тг!КНБ Параметр Описание »г!» дй (»гу») д (»гЯ „! „ определяется моделью Клаузера Постоянная Клаузера, п. 11.4.2 н, постоянная Кармана Ае+, основной демпфирующий параметр ван Дриста, п. ! !.4.2 Ае, демпфирующий параметр ван Лриста, связанный с давлением Параметр контроля давления ре, и+, п. !8.1.1 ит, скорость трения 1, длина перемешивания д!/да ТВН ПТКН ТЯНО АТК АКЯ АР2 АР РСО!ч УР, ОР 1!Т Е1. 01ТН имеется лишь один элемент и и! = 0 и 1 соответственно. Соответствующие формулы используются в подпрограмме гЕ1-Агт (рис. 15.15). В начале каждого шага вниз по потоку величины и,", и,"+', йй/ии и й+!/ип+! интерполируются по вводимым значениям ие(э) и г(и,Щ В подпрограмме АВСР (рис.

15.16) определяются коэффициенты трехдиагональной системы (15.66), а в подпрограммах ВА)чгАС и ВАЫБ01. (п. 6.2.3) осуществляется ее решение. В результате определяются значения Лба+!. Гл. 15. Течения в пограничном слое -$70 1 5ОВАООТ1ИЕ ТОЕЧ5(ЕЕЬ,ОЕЬТА,1НАХ,ОЕ,ООЕОХ 2 1АРХ,АТК,РСОИ) 3 С 4 С САЬСЦЬАТЕ ЕООУ У15С0$1ТТ, ТКЧ АИО ОТЕЧ 5 С 4 7 С 8 С 9 С 10 550 $0ЕТ(КЕЫ 11 еец $0ат(8$0) 12 тауО = Ата*це аеь ОеьтА 13 1НАР - "ПМХ " 1 14 11Н 1НАР 15 ОЦ = 1./11Н 16 УР(1) = О. 17 ОР (1) О.

18 цт 50ат(тАО(ы)/кео 19 РР = -ОЦЕОХ/ЦЕ/АЕЫЦТ/ОТ/ЦТ 2О ЕИ = 1. + РСОИ РР 21 АР = АР2/ЕИ 22 ацт - цт 850 23 ЦР(2) ОШЦТ 24 УР(2! ЦР(2) 25 теу(1) = О. 26 С 27 28 29 30 31 32 33 34 35 36 37 3$ 39 40 41 42 43 44 45 46 47 4$ ТАЦ,ТЕЧ,ОТЕЧ,АКР, 01НЕИ$10И ТАЧ(41),УР(41),ТАЦ(41),ОТЕЧ(41),ЦР(41) тачц 1$ тне ецву у15с051ту 1Н тне Оцтее Аец10Н(сьацхеа РоенцьАт10и) 00 2 К = 2,1МАР АК "- К - 1 Ц АК*ОЦ 1Г(К .ЕО.

2)00 ТО 1 ЦР(К) = Ц/ЦТ цх 1./(1. - ц + г. Оц) цт = 1./(1. — ц + Оц) Цх 1./(1. - Ю Окн - Оц (цх/ТАО(к-г) ч 4.*цыган(к-ы + цитхц(к))/3. УР (К) = УР (К-2) + ОЕН*АЦТ 1 ОЦН = 1. — ЕХР(-УР(К)/АР) ЕЬ = АКР*УР(К)*ОЦН ТАУ(К) ЕЬ*ЕЬ*ТАЦ(К)*(1.-Ц)/ТАЦ П) 11(теч(к) .От. теуо)сото Отач(к) = -(1.-0)*кь кь тац(к) тхц(к)/тхц(ы 2 СОНТ1ИЦЕ 3 ОО 4 Ь - К,ХНАХ ОТЕ7(Ы = О. тау(ы = тачо 4 соит1ицк ЯЕТОЕИ ЕИО Рнс. !5.14. Распечатка программы Т()мч'5. Прежде чем получить новое решение 0"+' = 0" + ЛО ~), рассматривается необходимость изменения шага маршевой переменной /Ц в зависимости от значения ЛО"+ /О". Если ЛО"+'/(О" ) ) у, шаг Л$ уменьшается вдвое при условии, что полученный шаг будет больше минимального значения (30.

Если Л$ уменьшается наполовину, й н т. д. пересчитываются и рассчитывается новое решение ЛОце(. Обычное значение Т порядка 0.02. % !5.3. Метод Дородиицыиа описания пограничного слоя 27! 1 20 3 С 4 С 5 6 7 в 9 3О 11 12 13 14 35 С 16 С 17 С 1$ 19 20 21 гг 23 24 25 26 27 С 2В С 29 с зо ы Эг зз 34 35 36 Зт ЗЗ 39 С ЗО С 41 С 42 43 44 45 46 47 4В 49 5О 51 52 53 54 55 56 57 5В 59 60 01НЕИ$10Н СС[41,5),АА(41,5),ЕГ(41,5),АВС(5,65),Р(65) 1, ТАО (41), ТЯЧ (41), ОТЕЧ (41) СОНИОИ СС,АА,БГ,АВС,О,ТАО,ТЯЧ,ОТЕЧ $1 1./б. Т1 = 1./3.

ттг 2."тг 13ВГ = 1НАХ - 1 ЯИИ = 1ИАГ ОО = 1./)ИИ ООЗ Оп*ВО ССИ,11 ССИ,2) ССИ,З) ААП,1! ААИ,2) Я)И1,3) БГИ,1) ЕГИ,21 ЕГИ,З) АТ ООТЕЯ ЕООЕ ОГ В.Ь СС ПЗИХ, 1) СС ( ПВХ. 2) СС(13ИХ,З) АА НИАХ, 11 ААПВП,21 АА ИВАХ,Э) ЕГ ПНЯХ, 1) ЕГ ИВАХ. '2) ЕГ ПНЯХ, Ы РО 1 д 2.1НАГ ° Ад 0-1 Ц = Ад*Во ССПьы = $1*(О - 0.5чОО) СС(д,21 ТТ1*О СС(4,3) $1*(О + 0.5*ВО) ООН И. - 01/ОО ООЗ = ООН ООН ААП,11 = - ООЗ - ООН - Т1 ААП,г) 2.*006 + ттх А!Палы = - ооэ + Оаи - тх ООН 0.5" И. - О*О) /ОО РО$ = 0.5*О - $1 егы,и Оои + ооз - (и*оп ег(д,2) = — тт1 ЕГ(0.31 - ООИ + ООЗ + $1"РО 1 СОНТ1НОЕ Яхтоан ЕИВ Рис.

!5.!5. Распечатка программы ГЕРАЙ. ЗОВЯООТ1ИЕ ГЕГАЯ ПНАХ) СОНГОТЕЗ ЧАЬ(ЗЕЗОГ СС,ЗВ, АНО ЕГ АСВОЗЗ ВООН(ИВТ ЬЛТЕЯ (15 Э 13) Ат ВАьь о. = ооиг. - ОО/12. о. - э ПРОЗ + тх - -ААИ,2) + 1./РО - о. = -ол/оо - тх = -ОЛ/ОО - 61~И.-ОО) 61*И. - 0.5"РО) Т1*И. " 0.25"ОО) - о. -Т1 тз О. = тэ" и. " ел*оп) -тэ О. АТ 1ИТЕЯ10Я ИОРЕЗ Гл. )6. Течения в пограничном слае 272 Если ЛО"+'/О" ( 0.1у, шаг (Ц увеличивается на 50 О/О до тех пор, пока не будет достигнуто максимальное значение шага. Изменение шага позволяет вести расчет величины О с малым шагом в области ее быстрого изменения и с большим шагом там, где О изменяется слабее.

В программе РОКОВ коэффициент поверхностного трения с! определяется непосредственно через значение т по формуле с! = 2т (1)/Ке(с, (15. 70) где число Гейнольдса Кес = (/ /./ч. Толщины вытеснения б* и потери импульса б ' получаются в результате численного интегрирования методом Симпсона с использованием выра- жений б'=,, ~ (1/с) с(и, и Кеи' 3 б ' = , ~ (и/т)(Аи.

и ИЕА (15.71) (15.72) Программу РОКОО можно использовать и для расчета пограничных слоев с обратным градиентом давления. Экспери- 1 2 С 2 С 4 С 5 6 7 а с 5 19 11 12 11 14 15 16 17 1В 15 29 21 22 21 24 25 26 27 26 ЕОВЯО(п19Е АВСР(1ВАЯ,РХ,РОЕРХО,РОЕРХР,ОЕ,ОЕТ,ВЕТЪ) СОМРОТЕЕ АВС,Р АСЯО$$ ВОВМРЯЯЧ ЬАТЕЯ [15.1.16) 913ЮВ$10И СС(41,5),АА(41,5),ЕТ(41,5),АВС(5,65) 1,Р(65),ТАО(41),ТВЧ(41) .РТЯЧ(41) спм(ОИ СС,АА,ЕЧ,Авс,в,тАв,тяч,втяч 90 2 Х 1,1ИАХ двТ 1 дти 1 11(Х .ЕО. 1)дЕТ 2 11(к .ев. 1нях)дтм 2 91Н О.

РО 1 д дЕТ,дтИ дт = д + 1 ЕР К-2+д ООН РОЕРХРвкт(к,д) Рки АА(к,д)* ип + тяч(кв)) РМ( ТАВ(кв)*(РЕН"ТЯО(КР)" АА(К,д) РТЯЧ(ХР))*ОЕР АВсЫР,Ю сс(х,д) + ВетА*3п" (РАИ-Рвн) РОИ = РОЕРХО"Ет(к,д)*(1.-ВЕТА) + ОРИ*ВЕТА 91Н Р1Н + РОИ/ТЯО(кв) + ВЕН*ЧАР(кв)*(ОЕ'(1.-ВЕТА)+ОХР*ВБТЯ) 1 спит(мве Р(Ю 91Н*РХ 2 СОИГ(М)Е ВЕТОИ( ЕИО Рис. )6.)6.

Распечатка программы АВС0. 9 15.3. Метод Дороднииына описания пограничного слоя 273 1.63548 .381$1 .02981 .28916 2.4ЗЯ .813 -.081 Е-О! .33270Е.О1 ментальные данные [Вгаг)з(7аэг, 1967) для таких течений приведены в трудах Стэнфордской конференции по турбулентным течениям 1968 г. [Со1ез, Н1гз1, 19681. Эти результаты использовались для определения входных данных (рис. 15.17) в программе 130)чОО.

Типичная выдача для этого случая приведена на рис. 15.18; видно, что шаг (Лх — = 635) увеличивается по мере продвижения вдоль х. Это согласуется с уменьшением скорости изменения решения. Очевидно также, что алгоритм не допускает неприемлемо большого увеличения шага Лз. Изменение коэффициента поверхностного трения и толшины пограничного слоя при продвижении вниз по потоку приведены 18 К Флетчер. т. 2 500 11 0.60 0.4125.0021.00 0.168Е-01 1.601Е+06 0.0018-00 0.001Е-ОО О. 1008-00 0.2008-00 0.918Е-ОО 3.810Е+00 0.0205-00 ОААО5НАЫ РЬОН С) АОЧЕЫЕЕ Р.б.

.265Е-02 .160Е+07 .73)Е-02 .532Е-02 22 1.67653 1.7153) 1.75605 1.79910 1.82957 1.85014 1.83343 1.80281 1.72612 1.5Э116 1.39436 1.26574 1.10231 .95764 .81205 .67268 .55842 .45748 .31169 .23790 .18126 .09885 .07391 .05997 .04830 .03961 .03327 .02949 .03254 .03581 .04051 .04739 .05840 .07624 .10474 .15709 .44260 .938 1.063 1.188 1.313 1.438 1.563 1.688 1.813 1.938 2.063 2лаа 2.313 2.563 2.688 2.938 3.188 3.313 3.438 3.563 3.688 3.813 1.011 .989 ..968 .949 .925 .904 .890 .872 .859 .847 .835 .822 0.800 0.790 0.774 0.760 0.750 0.742 0.736 0.729 0.726 -.119 †.196 †.181 -.167 -.156=.141 -.126 -.113 †.105 -.098 -.093 -.086 -0.076-0.072-0.067-0.064-0.061-0.059-0.057-0.056-0.055 Рис.

!5.17. Начальные данные для программы 00й00. ОКАЯ, пиа ьоо ы ай ° .)ьоаыт йо ° .338 мост* з.8)о оаса. .Эоой+оо аатса .ыоа-о) ы .1оеа-О! Оан) .Эаеа-оз они ° .)Оее+оо Ыата. .Боое+оо АКР .4!Ой 00 ЯРЕ .)30Е+02 АТй .168$ 01 РСОЯ .2)ей+О) 1О1Т18$ ТЯО РЯОПЕЕ ТАО .16765Е+01 .)8296Е О! .)Т)6)Е+О! .!)657Е+О! .6)2$8Е+00 .З))49Е+ОО .73910 .)$$)ой-83 .)МБОЕ.О),44ЫОЕ+Оо оааньоао ГБОН с; яочензй Р.с. 30 ыаейй Г1Я)те еоекент$, тнаЫОСеасе нерео: н)21нс мнстн + чяа ояпьт ОАЫРЭЫ6 й ! Оа ло! а ль ТАО(!) !.ьм ст .оогьз о-та .Оозз н-то.

.Ноьз ьа-!.зьо ага* н 6! Оа= .0)! а !.)4 тйе(!) 1.433 сГ .08227 О.ты .0090 и-тн .0064 $$ !.409 ать .Оы 3"!.)4 тм(м !.)34 «Г .Оо)оя о тк .О)оь н-тк .0874 зя Э.мт йтк и 97 Рх .О)! а !.и тю(!) 1.2)9 ст .00394 О.та .0)23 н.та .0046 $$ !.44! йта ы=))4 ох .О)! х"!.74 ты(!) !.!)9 сГ' .00)86 о-та .0)39 н тя .0096 Зн !.448 атя к !)3 Оа .0)6 а*1.94 тяо(!) !.)38 сГ .08)80 О.тн .0)54 н.тн .0)0$3$ !.453 атн= а )3! Оа .026 х 2.)5 ты(!) !.)07 сг .00!)$ О.тн .0)70 н-то ,0)17 за !.4$6 йтн н з)9 оа .0)6 3*2.эь тйо(!) !.082 сГ .00!)! О.та .О)47 н.тн .0)28 $$ !.439 атн а )47 ох .026 3 2.м тЯО(!) !.062 сГ .00)$$ О-та .0204 и-та .0)39 за=!.460 ата н )$5 Оа= .0)6 х 3.76 т)л(!) !.043 ст ° .00)65 О.та .02)0 а.тн .0)5! ьы !.46) атн Н 362 оа .026 Л 2.94 ТЯО(!) 1.026 СГ .00362 О.ТН .0)35 Н.ТН .0)ЬО 3Н=).463 ЯТН н !)О ох .0)Б х 3.)$ тйе(!) !.003 сг .00)$9 о-то .0)м и-та .0)72 За !.466 ата я пь Ох .026 3 3.3$ тзл(!)= .985 сг .00)м О-тн .0)70 а.тн .0)84 зн !.469 ата н )86 Оа" .024 х 3.$6 тйо(!) .96) сГ .00)$3 О.тк .0288 и"тн= .0)96 зн=!.47! Ятй Н')94 Оа .0)6 Х 3.76 ТАО(!) .948 СГ ,ОО)$0 о-ТН .0307 Н-ТН .0308 $$ !.475 ЯТН Ы 201 Ра .026 А*Э.94 ТЮ(!)= .930 СГ .0034) о-ТЫ .0320 Н.ТН= .03)7 Бо !.478 ЯТН НБОЕ+04 996Е+04 3 !)Е+05 324Е+0$ )35Е 0$ )46Е 05 )$7Е 0$ )ьае+0$ Ы9Е+05 ЭЭОЕ 0$ )99Е+0$ ))ой+05 22ОЕ+0$ 23)а+05 24)Я+05 2$2Е+05 Рис.

15.18. Типичная выдача программы 00й00. ТАЧ(4) .919 тач(4) !.3$4 тят(4)"!.63) тат(4) !.896 Тм(4)'2.0)) ТРЛ(4) 2.24) тйт(4) 2.388 ТЯЧ(4)*2.$35 Т$ Ч (4) 2. 614 Тат(4) 2.743 Тат(4)=2.8$7 ТАЧ(4) 3.0)8 Чат(4)*3.!)) тат(4)=3.33) Тат(4) Э.323 ТАЧ(4) 3.708 Гл. 15. Течения в пограничном слое 0.003 0.002 ст 0 001 !.0 2.0 3.0 ж/ь Рис. 15.19. Сравнение поверхностного трения при обратном градиенте дав- ления. 0.03 ф Г 6 Т 0.02 0.01 !.0 2.0 З.О дусь Рис. 15.20. Сравнение толщины вытеснения и потери импульса при обратном градиенте давления. на рис. 15.19 и 15.20. Видно, что РОгсОР-РЕМ (рассматриваемый метод) позволяет получить решение, сравнимое по точности с решениями, полученными по программе РОТОР-ПРЕС $15.3. 7"ьетод Лороднипына описания пограничного слоя 275 (п.

Характеристики

Тип файла
DJVU-файл
Размер
6,97 Mb
Тип материала
Высшее учебное заведение

Список файлов книги

Свежие статьи
Популярно сейчас
Как Вы думаете, сколько людей до Вас делали точно такое же задание? 99% студентов выполняют точно такие же задания, как и их предшественники год назад. Найдите нужный учебный материал на СтудИзбе!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6418
Авторов
на СтудИзбе
307
Средний доход
с одного платного файла
Обучение Подробнее