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

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

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

Метод Дородницына описания пограничного слоя 233 Типичные результаты расчетов коэффициента поверхностного трения и изменения толщины пограничного слоя в направлении потока при нулевом градиенте давления приведены на рис. !5.10 и !5.1!. На рисунках приведены также результаты, полученные по методу Дородницына в сочетании со спектральным методом (ООКОО-5РЕС) (п.

15.3.3) и по типичной программе метода конечного объема 5ТАЬ!5, основанной на алгоритме ОЕ!ч[М!Х [Ра1ап!саг, Бра!б!пд, !970]. Очевидно, что точность всех трех методов примерно одинакова. В алгоритмах 8ТАЬ)5 и ОЕ!ч[М1Х Таблица 13.2. Сравнение алгоритмов !гОЙО!л-ГЕМ и ЯТА!Ч б Нулевой граднент давленвн Обратный градиент давленнн зтАыб ионов-рем зтАив Роноо-рвм используется преобразование Мизеса [БсЬ!!сЬ1!пд, !968[ уравнений, в которых функция тока тр используется в качестве независимой переменной вместо нормальной координаты у.

Однако из табл. !5.2 видно, что ООКОО-РЕМ примерно на порядок экономичнее 8ТА!ч!5. Преимущество получается частично за счет возможности использовать в три — четыре раза меньшего количества точек поперек пограничного слоя и частично за счет уменьшения требуемого числа шагов в направлении потока. Описанный выше метод Дородницына в сочетании с методом конечных элементов применялся для расчета ламинарных [Р!е1сЬег, Р!ее1, 1984а] и турбулентных [Р!е1сЬег, Р!ее1, 1984Ь] несжимаемых пограничных слоев.

Данный метод применялся также для расчета турбулентного сжимаемого пограничного слоя [Р!ее1, Р!е1сЬег, 1983], течений в пограничных слоях с поверхностным переносом массы [Р!е1сЬег, Р!ее1, 1987] и внутренних пограничных слоев с вихрями [Р!е1сЬег, 1985], 264 Гл. 1З. Течения н пограничном слое 15.3.2. РОКОР: течение в турбулентном пограничном слое В данном разделе описана программа РОКОР, в которой реализуется метод Дородницына в сочетании с методом конечных элементов (п.

15.3.1). Блок-схема программы РОКОР приведена на рис. 15.12, а распечатка — на рис. !5.13. Параметры, используемые в программе РОКОР, описаны в табл. 15.3. Рнс. 15.!2. Блок-схема программы 1ЗОйОР. Для начала необходимо определить профиль т(хм и;), соответствующий иь Он получается из расчета ди/дп (15.58) на равных интервалах по у (и, следовательно, по г!) и интерполяции т на равные интервалы по ир До начала интегрирования вниз по потоку в подпрограмме Т()К ч'5 (рис. 15.14) определяется распределение турбулентной вязкости тг(у)7ч. Параметры, используемые в подпрограмме ТОКУ5, описаны в табл, 15.4.

Турбулентная вязкость рассчитывается в подпрограмме Т()КЧ3 на основе двухуровневой алгебраической модели, описанной в и. !1.4.2. Таким образом, вблизи стенки вводится длина смешения, а во внешней области пограничного слоя используется формула Клаузера (11.79). Для определения 0; (в подпрограмме АВСР) в подпрограмме Т13КЧЯ рассчитывается также д(чгЯ7д9.

Коэффициенты СС, ЕЕ и АА, определяемые выражениями (15.63), вычисляются в подпрограмме ЕЕРАК (рис. 15.15) по )5.3. Метод Дородницына описания пограничного слоя 265 ?С ХС' 3 С 4 5 б 7 8 9 С 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 ЗО 31 С Э2 С ЗЗ С 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 ВОВОЮ УБЕ5 ТНЕ ?1М1ТЕ ЕЬЕИЕИТ НЕТИОО ТО БОЬТЕ ТИЕ ООЯООИ1ТБТИ ВООМОАЯТ ЬАТЕВ ГОЯНУЬАТ?ОН 01КЕИБ10И АВС(5,65),О(65),ТАУ(41),ОТНЕТЫбй О1нем810и тАОО(41),тВТ(41),атяч(41) 01НЕМБ10И СС(41.5),АА(41,51.ЕР(41,5) 01НЕИБ10И ХУЕ(24),УЕЕ(24),УЕХ(24),Т1ТЬВ(8) соннои сс,АА,ХР,АЙс,а,тху,тят,агат ОРЕМ(1,?1ЬЕ='ООЕОО.ОАТ'1 ОРЕМ(б,г1ЬЕ='ООВОО.ОУТ') ОРЕМ(2,?1ЬЕ 'ООЯОО.БТА') ВЕАО(1,1)ИНАХ,1НАХ,ВЕТА,АКР,АРЕ,РСОМ,АТВ,ЯЕЬ 1 ГОЮ(АТ(215,495.2,2Е?0.3) ЯЕАО(1,2)ОХ,ОХЙ1,ОЮ(А,О?СИ,ХО,ХНАХ,МА?Си 2 РОЮ(АТ(76?О.И УЯ1ТЕ(6,3)МНАХ,1НАХ 3 гоанат( иихх,?нхх = ,юи УЯ1ТЕ(6,4)ВЕЬ,ХО,ХМАХ,ОХСН,ЯАТСН 4 ГОЮ(АТ(' ЯЕ =',Е10.3,' ХО ',ГБ.З 1,' ХНАХ '.ГБ.З.' ОХСЙ ',Е10.3,' ЯАТСН ',Е10.3! НЯ1ТЕ(6,5)ОХ,ОХН1,0?ИА,ВЕТЯ 5 ГОВНАТ(' ОХ ',Е10.3,' ОХМ? ',Е10.3,' ОХНА ',Е10.3,' ВЕТА ',В10.3) МЯ1ТЕ(6,6)АКР,АРЕ,АТЯ,РСОИ б ГОЯКАТ(' АГР '.Е10.3,' АРХ=',Е10.3,' АТВ ',810.3,' РСОН ',Е?0.3,/) 1ИАР ?Ю(Х - 1 КСТ 1 1ВВГ 1 ВБО БОЯТ(ВЕЫ ВЕО = БОВТ(ВБО) ВБАО 1Н Бтаат1ис теьос1тт Ряо?1ье и(0 ехтеВию 7еьос1тг РВОР1ьк ВЕАО(2.71Т1ТЬЕ ВХАО(2,8)с?Бт,Вхь,акьт?,тигмон,яРО ВЕАО(2,9)ТАМО ЯЕАО(2,10)(ХУЕ(И),И=1,НРС) ВЕАО(2, 10)(УЕЕ(М),М 1,НРС1 ЯЕАО(2, 101(УЕХ(И),И=1,ИРС) 7 РОВНЯ?(БА41 Б ГОЯИАТ[4Е10.3,15) 9 ГОЯНАТ(10?8.5) 10 РОЮ(АТ(ЗЗРБ.З) УЕ УЕК(1) ОУКОХУ УЕХ(И РОККИ) 10ЕЬ 1 1ГИИАХ .ЕО.

6)1ОЕЬ 8 И'ЮЮХ .КО. Ы)?ОБА=А 1?(1НАХ .ЕО. 21)10ЕЬ = 2 ОО Ы 1 1,?МАХ ?А "- 1 + (1-11*?ОЕЬ 11 ТАУ(1) ТАУО(1А) УЯ1ТЕ(6,12) 12 ГОЯНАТ(! 1И1Т1АЬ ТАУ РРОГ1ЬЕ') МЯ1ТЕ(6,13)(ТАУ(1).1 1,1ИАХ! 13 ГОЯИАТ(' ТАУ=',БЕ12.5) УЯ1ТЕ(6,14)Т1ТЬЕ 14 РОЯНАТ(' ',ВА41 УЯ?ТЕ(6,15)1МАР 15 ГОЯНАТ(( ',12,' Ь?НЕ?Я Г1Н?ТЕ ЕЬЕНЕИТБ, ТУВВУЬЕИСЕ НОВЕ?0 Ь 1' Н1Х1НС ЬЕМСТЯ г ЧИ! ОЯ?ЕБТ ОИ(Р1ИС') Рис. 15.)3. Распечатка программы ))ОПОР (начало).

Гл. 15. Течения н пограничном слое 62 С 63 С 64 С 65 66 67 С 68 69 70 71 С 72 С 73 С 74 75 С 76 77 78 79 во с В) С вг с Вг 84 85 86 87 ° ав 89 90 С 91 С 92 С 93 94 95 96 97 98 99 1ОО 1О) гог 101 104 105 106 107 !Ов 109 жо 111 112 113 114 С !15 С 116 С 317 жв с 119 120 )гг с ЕЕТ 1Н1Т185 ТКЧ САЬЬ ТРАТЕ(АЕО,ОКОТА,)НАХ,ОЕ.ООЕОХО,ТАО,ТКЧ,ОТНЧ,АКР, 1АРЕ,АТН,РСОН) ОО 16 1 = 1,1НАХ РТЯЧ(1) О. СОНТ1НОЕ 16 СОНРОТЕ СС,АА АНО ЕР САОО РЕРАН()НАХ) 00 18 1 = 1,1НАХ ОО 17 д = 1,5 17 АВС(3,1) = О. 18 СОНТ1НОЕ ВЕС1НН1КС ОЧ ООЧИ5ТЛЕЛН ЬООР ХСИ = ХО Х = ХО А)Н = 1МАР ОУ = 1./й1Н КН=2 ОО 27 Н = 1,ННАХ КВТ = КН СОНРОТЕ РКЕ99ОКЕ СЯАО1ЕНТ РАВАНЕТЕКЕ ОО 19 КА = К5Т.НРС К = КА 11(Х .СТ. ХОЕ(К))СОТО 19 КН = К-1 17(К+1 .СТ.

ИРС)К = НРС"1 ОШ( = ХОЕ(К+1) — ХОЕ(К) Х1 = (ХОЕ(К-1)-ХОЕ(К))/ООН Х4 = (Х - ХОЕ(К))/ООН Хб (Х + ОХ - ХОЕ(К))/ООИ ООН = ОЕЕ(К-1)-ОЕЕ(К)-Х1*Х1*(ОЕЕ(К+1)-ОЕЕ(К)) А! = ВОН/Хг/().-Х)) 82 = ОЕЕ(К~1)-ОЕЕ(К)-А) ОЕ = ОЕЕ(К) + А1 Х4 + А2*Х4*Х4 ОЕР "- ОЕЕ(К) + А1*Х5 ~ А2"Х5*Х5 РОН = ОЕХ (К" 1) "ОЕХ (К) -Х1*Х1* (ОЕХ (КАИ -ОЕХ (К)) 61 = ООН/Х1/!1.-Х1) А2 ОЕХ(К+1)-ОЕХ(К)-А1 ООЕОХО (ОЕХ(К)+А1 Х4+А2*Х4"ХЯ)/ОЕ РОЕОХР = (ОЕХ(К) + А1'Х5 + А2"Х5*Х5)/ОЕ СОТО 20 19 СОНТ1НОЕ САОСОЬАТЕ ТН1Р1АСОНАЬ ТЕКНВ 20 САЬО АВСО(1НАХ.ОХ,ООЕОХО,ООЕОХР,ОЕ,ОЕР,ВЕТА) САЬО ВАНРАС(АВС,1НАХ.!) САЬЬ ВАНВОО(О,ОТВЕТА,АВС,!НАХ, М Рис.

15.13 (продолжение). Р (5.3. Метод Дороднипынв описания пограничного слоя 207 122 С 12З С 124 С 125 126 127 128 129 ЭЗО ЭЗ) 1Э2 ЭЗЗ 134 135 С 136 С 137 С ЭЗ8 139 ЭАО 141 142 С 14З С 144 С 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 С 162 С 163 С 164 165 166 С 16Т 168 169 170 171 172 173 С 174 С 175 С 176 177 178 179 180 181 КАТ 1$ ТНЕ ТЧР1САЬ РВАСТ10ИАЬ СНАИОЕ 1М ТВЕТА ййт = АОБ(тяп(Эйке! ОТВЕТА(зйкг!) 1Г(0.5"ОХ .ЬТ.

РХН1)СОТО 21 1Г(ААТ .СТ. ЯАТСН)РХ = 0.5 РХ ОЕР = 0.5*(ОЕеОЕР) РОЕРХР -" 0.5*(ООЕОХОБРОЕОХР) 1Г(ААТ .ОТ. ЕАТСВ)СОТО 20 21 Х = Х + РХ 1Е(1.5*ОК .ОТ. РХНА)СОТО 22 1Г(ВАТ .ЬТ. 0.1*ААТСН)ОХ 1.5*ОХ 22 СОИТ1НОЕ ЕЧАЬОАТЕ ИЕМ ТАО АВААЧ РО 23 1-"1,1НАХ 23 ТАЮ[1) 1./(1./Тйа(1)+ОТНЕТА(1)) ст * 2.*Тяп(1)/Бойт(йеь! 1Г(ТАО(1) .ЬТ. 0.)ООТО 25 СЯЬСОЬАТЕ 01$РЬАСЕНЕИТ АИО ИОНЕИТОН ТН1СКИЕББ РЕЬТА О. тнхмон=О. 1И 1НАХ/2 РО 24 1А 1,1Н 1 = 2"1А - 1 81 = 1 - 1 О1 - "ОааА1 РЕЬТЯ РЕЬТА + 2./ТАО(1) + 4./ТЯО(1+1! такнон = тнкиОИ + 2. 01/тАО(1) + 4.*(01ГОО)/ЕАО(1+1! 24 СОИТ1ИОЕ РЕЬТА = (ОЕЬТАЕЭ./ТАО(1НАХ! "1./ТАО(1))*ОО/3.

тнкнон = (тнкнон+1./тяп(1НАК)!"РО/3. БНРАТй = РЕЬТА/ТНКМОН ОЕЬТА - "ОКАТА/Ок/БОВЕ(ВЕЬ! тнкнон - такиои/Ок/БОМТ(йкь! АТН = КЕЬ ОЕР*ТНКНОИ САЬСОЬАТЕ ЕОРЧ Ч1$С0$1ТЧ САЬЬ ТОАЧБ(ВЕЬ,РЕЬТА,1НАХ,ОЕР,РОЕРХР,ТАО,ТАЧ,ОТЕЧ,ЯКР, 1АРХ,АТВ,РСОИ) 1Г(Х .ЬТ. ХСН)ООТО 27 ХСН = ХСН + ОХСН 25 че1те(6,26)И,Ох,х,тАО(1),сР,РеьтА,тнкион,бнРйтй,йта,тАЕ(4) 26 ГОВМАТ(' И='.1Э,' ОХ=',Г5.3,' Х=',Г4.2.' ТЯО(1)=',Г5.3, 1' СЕ-"',Е).5,' О"ТН ',Г6.4,' Н-ТН ',Е6.4.' Бн=',Г5.3, 2' АТН ',Е10.3,' ТВЧ(4)=',Гб.з! ТЕБТ Гой БЕРАААТ10И АИР Х .СТ.

ХИАХ 1Г(ТАО(1) .ЬТ. 0.)ООТО 28 1Г(Х+0.001 .СЕ. ХНАХ)ООТО 28 27 СОИТЭМОЕ 28 СОИТ1ИОЕ БТОР ЕИО Рис. !б.!3 (ОВоичание). 268 15.3. Таблица Описание Параметр следующим (внутренним) формулам: и/ — Пз и/ сс...—, сс,,— —, 1 — и', 1+ Ли ЕЕ/ 1, — — 0.5 / + 0.5и/— (15.68) ЕР/ / — — — 2/3, 1МАХ (т(МАХ 11, 011 Х, 0Х 0ХМ!, 0ХМА 0ХСН ХО, ХМАХ ХУЕ НЕЕ, 13ЕХ 13Е, ПЕР 0УЕОХБ 0БЕОХР ТАО 0ТНЕТА йАТ йАТСН АВС 0 ВЕТА йЕЕ СР 0Е1.ТА ТНКМОМ ЗНРйТй йТН Гл.

15. Течения в пограничном слое Параметры, используемые в программе 00й00 Число 1пагов 6$ Число точек сетки поперек пограничного слоя и, Ли Е 6$ ~!4в и, .'4вПриращение 6$ при печати решения аа Ива» Точки Е в которых заданы и и и и. и и,1 на границе пограничного слон (задаются) н н+1 е е ((~иа/66)/ие)" ((йиещ)/ие)м ~ т 66 66п+ 3/6п у ССС в (15.66) Р в (15.66) 6, контролирует степень неявности, (15.64) йеы относительное число Рейнольдса сь коэффициент поверхностного трения 6*, толщина вытеснения 6 ', толщина потери импульса Н = 6"/бв ", фактор формы и,б ' йес, число Рейнольдса по толщине потери им- пульса 1 — и/ 1 — Ьи ЕР/,/+~ = — 0.5 + 0.5и/— Ьи 6 й !6.3. Метод зсородницына описания пограничного слоя 26Э 1 — и! ! 1 АА! ! ! — — — (1 — я!) д (1 — и,.) 2 АА =2 ' + —, Ь! аиа 3' 1 — и+, 1 АА! 1+' (1 и!) '~ 3 ' (15.69) Уравнения (15.67) — (15.69) получены в результате аппроксимации (15.63) одномерными линейными элементами при разбиении на Ли.

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

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

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

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