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

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

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

Простые течения в пограничном слое 249 Гхьквке-БЕБИ Боьотзои вата- ло ЗИБх и Оти= .Яо ВГ= 1.оо ИИВХ 19 ОХ= . 100Е+00 Х57= 1.00 ВЕ . 100Е+06 Рис. 16.6. Типичная выдача программы 1АМВ1.. тем интерполяции (1В и среднеквадратичное отклонение между и и ио„. Как видно из рис. 15.5, полученные численные результаты весьма близки к решению Фолкнера — Скан. 15.1.З.

Схема ячеек Келлера Другой метод дискретизации уравнений пограничного слоя (15.1), (15.2) предложен в схеме ячеек Келлера. Основное отличие итого метода состоит в том, что в нем фигурируют лишь первые производные. Уравнение (15.2) заменяется уравнением ди ди ди, дт и — +о — =и — +— дх ду Я дх ду ' (15.23) 1-1 1-1/2. где ди т = 11— ду ' (15.24) и и+112 и+1 Рнс. 16.6.

Сетка в схеме ячеек Кел ет ис- лера. Таким образом, возника куственно введенная дополнительная переменная — сдвнговое напряжение т. Дискретизации проводится внутри «ячейки», как показано на рис. 15.5. Н 1 и г и= 3 И= 4 И 5 И 6 И= 7 Н- В и 9 Н ЭО н= ы И 12 и ы И 14 И 15 Н- 16 И 17 Н 1$ И- 19 ВИВ Х 1.20 КХСГ Х 1.ЭО КХСГ Х-З.ЯО КХСГ= Х 1.50 ЕХСГ Х*З.ЯО КХСГ= Х ЗЛО КХСГ Х=Э.ВО КХСГ= Х 1.90 ЕХСР х Э.оо ысг= х Э.зо ехсг= Х.=г.го кхс. Х=2.30 ЕХСГ-" Х=2.40 КБСГ= Х=2.50 ЕХСГ х 2.6о кхсг Х 2.70 ЕХСГ х Х.во кхсг- Х 2.90 КХСГ Х 3.00 ЕХСГ .674Е-ОЗ .004243 .004022 .ООЭВ2$ .003656 .003502 .003364 .Оозгзв .003123 .00301В .002922 .002В32 .002750 .002673 .002601 .002534 .002471 .002412 .00235$ .оогзоз СГ .004242 СГ .004024 сг= .оозвзг СР .003660 Сг .003506 СГ .003367 СГ= .003241 СГ .003126 сг- .оозогз сг= .оогыЯ СГ= .002В35 СГ .002752 Сг .002675 сг- .Оогбоз СГ .002536 сг= .оо247з СГ= .002414 Сг .00235$ СГ .002305 01БР 01БР= ИБР 01$Р 015Р ИБГ 01$Р ОЛР О15Р= 01БР 015Р ИБР 015Р= 01БР= 01БР ОЫР 015Р= 01$Р 015Р= .003343 ОЕ 1.063 .003430 Ок 1.091 .оозиг Ок= З.ИВ .003591 ОЕ 1.145 .00366$ Ок 1.170 .003742 Ок 1.193 .ООЗВ1Э ОЕ 1.216 .ООЗВ$1 ОЕ 1.239 .ООЭ947 ОЕ 1.260 .004011 Ок 1.2$1 .004072 Ок 1.301 .оовыг Ок З.зго .004191 ОБ 1.339 .00424В ОЕ 1.357 .00430Э ОБ 1.375 .004357 ОЕ= 1.392 .004409 Ок 1.409 .004441 ОЕ= 1.426 .004511 ОЕ 1.442 й 15.2.

Сложные течения в пограничных слоях где ахни+! =вне+' — вне, хн! — — (и1, о1, т!)г, 1г — номер итерации. Таким образом, в начале итераций вне=хне, а после их окон! !' чания вн"+!=хне+'. Обычно на каждом шаге вниз по потоку ! достаточйо трех — четырех итераций для определения значений хнв+'. Можно отметить, что при дискретизации (15.25), (15,26) ! матрица Якоби обычно получается трехдиагональной и система (15.30) может быть эффективно решена при помощи алгоритма„ описанного в п.

6.2.5. Схема ячеек Келлера описана в работе [Ке!!ег, 1978] и более подробно в книге [СеЬес(, ВгадзЬахч, 1977]. й 15.2. Сложные течения в пограничных слоях В пограничных слоях с обратными градиентами давления (давление увеличивается в направлении потока) толщина пограничного слоя быстро увеличивается. Для любого типа пограничных течений компонента скорости и, направленная по потоку, сильно изменяется в направлении, перпендикулярном пограничному слою. По этой причине желательно преобразовать уравнения к новым зависимым и независимым переменным, менее чувствительным к приведенным выше эффектам. Если сетка в преобразованной области близка к однородной, ошибки, связанные с дискретизацией (см. $ 3.1), будут значительно меньше. Общие проблемы, связанные с использованием переменных сеток и возникающими при этом ошибками, рассмотрены в работе [Ь!оуе, 1983].

15.2.1. Замена переменных Как правило, в результате преобразования уравнений пограничного слоя возможно получить определенные преимущества. В результате преобразования Манглера (см. [БсЫ!сЫ1пд, 1968] ) осесимметрнчный пограничный слой сводится к эквивалентному двумерному. Преобразование Хоуарта — Стюартсона [3сЫ1сЫ!пй, 1968] позволяет заменить сжимаемый пограничный слой эквивалентным несжимаемым. Преобразование Блазнуса компенсирует увеличение толщины пограничного слоя и значительно упрощает уравнения в том случае, если задача имеет автомодельное решение. Преобразование Леви — Лиза [В!о(!пег, 1975а], рассматриваемое в п.

15.2.2, сочетает основные свойства преобразований Хоуарта, Манглера и Блазиуса. В преобразовании Дородницына ($ 15.3) скорость и используется в качестве независимой 252 Гл. 15. Течения н пограничном слое переменной. Это позволяет представить уравнения в интегральной форме и делает возможным применение метода Галеркина с конечными элементами и спектральных методов. В преобразовании Дородницына безразмерный градиент нормальной составляющей скорости рассматривается как зависимая переменная, Это позволяет с большой точностью получить сдвиговое напряжение на стенке, а следовательно, и коэффициент поверхностного трения (11.66). 15.2.2. Преобразование Леви — Лиза Отправной точкой преобразования Леви — Лиза являются уравнения, описывающие стационарный ламинарный сжимаемый двумерный (п.

11.6.2) или осесимметричной пограничный слой: уравнение неразрывности (~еРи) + д (гейо) 0 д д (15.31) уравнение импульса в направлении оси х ри — + ро — = — — '+ — ( 1я — ), (15.32) ди ди Нре д г дн х дк ду Нх ду (, ду)' уравнение энергии Для двумерных течений ! =О, для осесимметричиых / = 1, а ге — радиус тела; х в (15.31) — (15.33) измеряется вдоль тела от носка или точки торможения, а у — по нормали к поверхности; р, в (!5.32), (15.33) — известное давление иа внешней поверхности пограничного слоя.

Для сжимаемых течений уравнения (!5.31) — (!5.33) должны быть дополнены уравнением состояния, например (11.1), и зависимостью вязкости от температуры р(Т). Для уравнений (15.31) — (15.33) можно поставить начальные условия и(хо, у) =ио(у) Т(хо у) =То(у) (15 34) и граничные условия .и(х, О) = в(х, О) = 0 (отсутствие потока массы через стенку), Т(х, О) =То(х) или — (х, О) = — 1;1„(х), (15.35) и(х, Ь)=и,(х), Т(х, б)=Т,(х). $ !5.2. Сложные течения в пограничных слоях 253 В преобразовании Леви — Лиза вводятся две независимые переменные х ~( ) Кз(р )тес е о (15.36) т1(х, р) = и г', ( 2 ) ~ рг(р', о где К вЂ” константа, зависящая от рассматриваемого течения. Новые зависимые переменные вводятся следующим образом: и Ь Т а т и уравнения (15.31) — (15.33) принимают вид 2$ — + — +Р=О, дя дп (!5.38) 2ьР д + У д +р(Р— О) — д (1 д ) =О, (15.39) где и., г 2$ '! дн, рн Граничные условия (15.35) записываются в виде Р=г =О и О=О.

при п=О, 7=У=1 при п=п' (15.42) Начальные условия можно получить, если положить $ =0 в (15.35) — (!5.40). Если члены, содержащие производные по $, положить равными нулю, то автоматически получаются урав- нения, описывающие автомодельные течения (например, тече- ния Фолкнера — Скан, п. 15.!.2) [В!о!!пег, 1975а]. Преобразование Леви — Лиза обладает следующими замеча- тельными свойствами: (1) влияние сжимаемости становится несущественным; р явно не входит в уравнения; (2) осесимметричные течения можно рассматривать как эквивалентные двумерные; (3) использование координаты т! компенсирует увеличение толщины пограничного слоя. 254 Гл. !5. Течения в пограничном слое 15.2.3. Связанная схема Дэвиса В данном разделе будет описана связанная схема Дэвиса (РСЬ) решения системы (15.38) — (15.40) для несжимаемого пограничного слоя.

В этом случае 1= 1, 8 = 1 и необходимо рассмотреть лишь уравнения (15.38), (15.39). Применение схемы РСЬ рекомендуется Блотнером (В!оНпег, 1975Ь) после проведенного им сравнения по точности и эффективности ряда Р=п,т'=О Рис. !5.7. Сетка в пограничном слое в пространстве 15, Ч). схем Кранка — Николсона. Термин «связанная» объясняется одновременным неявным рассмотрением уравнений неразрывности и х-компоненты импульса в противоположность последовательному решению уравнений (15.7), (15.8). Если (15.39) записать в виде дга дл дтл э — = ттНЯ = — 1l — — 5 (ге — 1) + —, (15.43) д5 дч дт!е ' то, используя разностную схему Кранка — Николсона в узле (а+ 1/2,1), можно построить маршевый алгоритм (здесь тсНЗ означает «правая часть уравненияв).

Другими словами, (15.43) заменяется выражением (раун- ! (Га)ч 0.5($" + 5" ) ~ ' = 0.5(ЙНЗ" + КНЗ"~ ). (15.44) Производные по т! в (15.43) аппроксимируются трехточечными центральными разностями (рис. 15.7). Нелинейныс неявные члены в (!5.43) и (15.44) лииеаризуются разложением Ньютона — Рафсона, т. е. (ра а+' йрарет~ (рт)е (! 5.45) $15.2. Сложные течения в пограничных слоях 255 где итерации по тт проводятся на каждом шаге по и в направлении потока. После достижения сходимости по л получается решение на слое и+1. Подставив (15.45) в 115.44), можно получить линейную относительно Р'+' и У'+! систему уравнений а! Рте+ !! + ЬтР~"' + с! Р~~+~! + д; )г~+' = с1~т, / = 2, 3, ..., У вЂ” 1, (15.46) где а~т = — 0.5 (1 + 0.5р~тЛт1), б =1+6' 15"+'+ ~ ~ Р~ т= т1! с~! = — 0.5(1 — 0.56т1У!), д! =0.25Лт~(Р~!+! — Р! !), с1~т = — агР! ! — отР~те ! — [1 + 0.5АЧ 5"Р! [ Р! + + ~,е ( т+! ! — !) 1 нетРт)х [ (2лч) ! 1' Уравнение неразрывности (15.38) в разностной форме имеет вид 2 л+цт [( / !)+( ! — ! / — !)! 2Ь1 [(ул+! 1го+!) 1 (1го рл )) 2ЬЧ + О 25[(Рт+' + Р,"" !) + (Р ! + Р! !)[ = 0 (15.47) Уравнение 115.47) можно разрешить относительно У~~~'.

Уравнения (15.46), 115.48) решаются одновременно модифицированным трехдиагональным алгоритмом 1п. 6.2.2). Первый проход осуществляется от внешней границы пограничного слоя к стенке. На границе слоя Е, = О! =О, ет =1.0. При умень- У~!+ ' — — К~!" ! ~— г! (Р~т+ ! ~+ Р~т') + тт, 1 = 2, 3, ..., 7, (15 48 ) где зт — — 26т1(0.25+$"+н'/Л$), йл+ пг ~ л е 1! = — 28т1(0.25 — — ) [Р; + Р! !] — (т'; — У! !). Гл. 15.

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

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

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

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