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

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

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

Численный расчет таких течений наиболее эффективно может быть осуществлен на основе 11ХБ-уравнений. Исторически сложилось так, что вопрос о том, можно ли для данной )тХ8-системы получить устойчивое решение за один проход по маршевой переменной, решается эмпирически. Проводится некоторая аппроксимация и изучается численное решение. Если получается физически реальное решение, система ЙХВ-уравнений считается устойчивой. В работе [Вгйеу, МсРопа16, 1984[ более систематично исследован вопрос о том, является ли )сХ8-система уравнений эллиптической по маршевой переменной. Как отмечено в гл. 2, такой подход требует построения расширенной системы уравнений в частных производных первого порядка, которая может иметь сингулярную характеристическую форму. Это обстоятельство делает рассмотрение неубедительным.

В данной главе предпочитается другой более прямой подход, основанный на анализе Фурье (и. 2.1.6), примененном к системе уравнений. В п. 16.1.2 он будет использован для того, чтобы априори определить, можно ли для данной тсХВ-системы получить устойчивое решение за один проход по маршевой переменной.

Как и анализ на основе характеристик (гл. 2), метод Фурье устанавливает формальный тип уравнений. Однако анализ Фурье имеет значительное преимущество, заключающееся в том, что он точно указывает, какой член ответствен за эллиптическое поведение. Может оказаться, что потребуются дополнительные предположения, часто связанные с давлением, обеспечивающие неэллиптическое поведение (например, п. 16.2.2). Практическое применение тсХ8-подхода рассматривается в п. 16.1.4 на примере задачи о распределении температуры на входе в канал. 16.1.1. Анализ порядков величин В данном разделе тсХБ-уравнения будут выведены из уравнений Навье — Стокса для двумерных стационарных сжимаемых и несжимаемых ламинарных течений.

Для описания турбулентных течений цХВ-уравнения необходимо модифицировать. Однако метод вывода ЙХЯ-уравнений для осредненных параметров турбулентного течения по существу остается тем же, т. е. основан на сравнении порядков величин. Вывод укороченных уравнений Навье — Стокса в основном совпадает с выводом уравнений пограничного слоя [СеЬес1, ВгабзЬаь., 1977]. Эффекты вязкости здесь также считаются су- 5 16.1. Введение щественными лишь в слое, толщина которого 6 (рис. 16.2) мала по сравнению с характерным размером в направлении потока Е. Для ламинарных пограничных слоев величина 6/Е порядка 0(йе-ыа). Для типичных чисел Рейнольдса Ке = 10', это означает, что 6/Е = 0.001.

Одна из причин необходимости вывода КХБ-уравнений связана с рассмотрением вязких слоев, толщина Невязнеа (Ь) Рис. 16ен Типичная толщина вязкого слоя (а) во внешнем течении и (Ь) во внутреннем течении. которых больше толгцины пограничных слоев. Так, для значений 6/Е в диапозоне от 0.1 до 0.01, по-видимому, следует использовать КНБ-уравнения.

В свою очередь это означает, что в полных уравнениях Навье — Стокса можно опустить члены порядка 0((6/Е)'), но необходимо оставить члены порядка 0(6/Е). Вывод уравнений пограничного слоя основан на пренебрежении членами порядка 0(6/Е). Для стационарного ламинарного несжимаемого двумерного течения уравнения Навье — Стокса (11.81) запишутся в безразмерной форме — + — =О, ди до дх ду (!6.1) до до др 1 т дзо дзо т и — +о — + — = — ( —,+ —,), дх ду ду )(е (, дхз дуз ) ' где число Рейнольдса Ке=р(/ Е/р. Безразмерная форма урав- нений (16.1) — (16.3) получена аналогично (!!.42). Для определения относительной величины различных членов в уравнениях (16.1) †(16.3) предполагается, что для и и о д/дх порядка 0(1), д/ду порядка 0(Е/6) и дя/ду' порядка 0((Е/6)а).

Поскольку и — величина порядка 0(1), а о порядка 0(6/Е), уравнение (16.1) упростить нельзя. Все члены в левой части УРавнения (16.2) порядка 0 (!). В правой части уравнения (!6.2) дзи/дх' порядка 0(1), а д'и/дув порядка 0((Е/6) ); поэтому членом дзи/дхе можно пренебречь. (16.3) 298 Гл, !6. Течения, описываемые Й!чо-уравнениями Навье — Стокса В классической теории пограничного слоя величина 1/Ке порядка 0((6/Е)а).

В приближении КМ3-уравнений предполагается, однако, что 1/Ке « 0((6/(.)а); а именно что 1/Ке порядка 0 ( (6/1 ) а), и поэтому (1/Ке) д'и/ду' порядка 0 ( (6/1 ) ). Все члены в левой части уравнения (16.3) порядка 0(6/1,). В правой части (16.3) порядок члена (1/Ке)дао/дх' равен 0 ((6/1.)'), а порядок члена (1/Ке)д'о/дуа равен 0((6/Е)а). Поэтому производной д'о!дх' можно пренебречь. Таким образом, К(ч!Ь-уравнения имеют вид ди до — + — =О, дх ду (! 6.4) ди ди др ! даи и — +и — + — = — —, дх ду дх ке дуа ' (16.5) до до др ! дао и — +о — + — = — —,.

дх ду ду ке ду' ' (16.6) дх (Ри) + д (ро) = О д д ди ди др ! Г дтхх дтхв 1 ри — +ро — + — = — ! + — 1 дх ду дх Пе ь дх ду (16.7) (!6.8) В классической теории пограничного слоя уравнение (16.6) упрощается до др/ду = О. Однако в случае более толстых вязких слоев, для описания которых требуются КИЬ-уравнения, необходимо сохранить все члены в уравнении (16.6). При значительной кривизне линий тока в (16.6) следует также включить дополнительный центробежный член. При использовании обобщенных криволинейных координат (гл. 12), что обычно делается в случае областей расчета сложной формы, этот член возникает естественным образом. Согласно проведенному анализу, порядок члена (1/Ке)д'о/дуа равен 0((6/Л)а), и, как рекомендует Рубин [КцЫп, 1984), им можно пренебречь.

В данной исходной форме КХ8-уравнений этот член сохранен, поскольку он не приводит к эллиптическому взаимодействию (п. 16.1.3). Связанная с ним дополнительная диссипация может быть полезна для численных расчетов. Уравнения (!6.4) — (16.6), если пренебречь вязкостью, сводятся к уравнениям Эйлера, которыми описывается практически невозмущенное течение вдали от тела.

Граничные условия для несжимаемых КМЗ-уравнений зависят от типа системы (16.4) — (16.6); данный вопрос будет рассмотрен в п. 16.1.3. Двумерные стационарные сжимаемые ламинарные уравнения Навье — Стокса в безразмерной форме имеют вид $16!. Введение де де др 1 Г дтих дтвв ! ри — + рп — + — = — [ + — ~, (16.9) дх ду ду Ке ь дх ду 1+ тМ' р = рт. (16.11) Обезразмеривание, использованное при выводе уравнений (!6,7) — (!6.1!), аналогично приведенному в п. 11.2.5, за исклю- вФ И и чением лишь того, что здесь р =(р — р )/р (/ .

В уравнениях (16.8) и (16.9) т, и т. д.— вязкие напряжения; они связаны с градиентами скоростей соотношениями (!1.27), форма которых в безразмерном виде сохраняется. В уравнении (!6.10) вязкая диссипация Ф определяется соотношением (! 1.39). Оценка порядков величин различных членов в системе (16.7)— (16.11) проводится аналогично тому, как зто было сделано для уравнений (16.1) — (16.3). Порядок компонент скорости и и и и производных от и, и и р такой же, как и для несжимаемого течения. Порядок Т, дТ/дх, р и др/дх равен 0(1). Если определена температура стенки Т, то дТ/ду и др/ду порядка 0(ь/б). Однако, если стенка адиабатическая, т.

е. дТ/ду),=е = О, и течение дозвуковое или трансзвуковое, более правильно считать, что средние значения дТ/ду и др/ду поперек слоя порядка 0(1) или даже 0(б/~). Следовательно, при выводе укороченных сжимаемых уравнений Навье — Стокса необходимо рассмотреть два случая: (1) Определена температура )( зт зр стенки или течение сверхзвуковое 1 ду ' ду п р " ('/')' (2) Адиабатическая стенка и ( дт др течение до- или трансзвуковое 1 ау ' ду п Р - (1). Случай (1) соответствует большим изменениям температуры в расчетной области, связанным с большой разницей температур стенки и набегающего потока или с существенным сжатием, обусловленным большими скоростями (большие числа Маха).

Случай (2) соответствует меньшим температурным изменениям в области расчета, обусловленным лишь сжимаемостью. При Уменьшении числа Маха в случае отсутствия внешних источников тепла производные от р и Т стремятся к нулю. В случае (1) все члены в уравнении (16.7) порядка 0(1) н их необходимо сохранить. В уравнении (16.8) все члены в левой части порядка 0(1). В правой части (16.8) член т,„порядка ЗОО Гл. 16. Течения, описываемые й!ЧЯОравиеииями Навье — Стокса —,' (ри)+ д' (р)=6, (16.12) ди ди др 1 д Г ди т ри — + ро — + — = — — ( !х дх ду дх йе ду (, ду ) ' до ди др 1 Г 4 д Г до т и д'и 3 — +р~ — + — — 1ь — (1х )+ — ), (!6 !4 дх ду ду йе 13 ду (, ду ) 3 дхдуз' ( — 1)М д В данных уравнениях сохранены все «невязкие» члены из уравнений (16.7) — (16.11).

Уравнения справедливы при больших изменениях температуры в расчетной области. Во втором случае (адиабатическая стенка в дозвуковом или трансзвуковом потоке) основное отличие заключается в правой части уравнения (16.14). Можно заметить, что одр/ду в (16.7) мало, если др/ду порядка 0(1). Следовательно, если в тии подставить Я из (16.7), та (16.14) заменится следующим уравнением: (16.13) (16.16) 0(1). Безразмерная вязкость ведет себя примерно как безразмерная температура.

Порядок компоненты тензора сдвиговых напряжений т,„= р(ди/ду+до/дх) равен 0(Ь/б) и обусловлен членом ди/ду. Следовательно, членом до/дх, порядок которого 0(б/Е), можно пренебречь. Таким образом, в правой части уравнения (168) остается лишь член (1/Ке)д(рди/ду)/ду. В левой части уравнения (16.9) все члены порядка 0(б/С). В правой части т„- рди/ду. Отдельные части других членов можно сгруппировать, записав их в виде соответствующих производных от скорости. При умножении на (1/Ке) порядок многих членов в правой части уравнения (16.9) оказывается равным 0((б/С)а). В левой части уравнения (16.10) все члены порядка 0(1).

Исключение составляет лишь член одр/ду, порядок которого равен 0((б/Е)а). Однако и его следует сохранить, поскольку в не- вязкой области он может стать порядка 0(1). В правой части уравнения (16.10) членом д(идТ/дх) /дх/Рг Ке можно пренебречь. В вязкой диссипации единственным членом порядка 0((б/Ц) оказывается !т(ди/ду)а/Ке. Таким образом, согласно сказанному, укороченные уравнения Навье — Стокса имеют вид зо! э !б.!. Введение (16.17) За исключением члена, связанного с вязкой диссипацией в уравнении (16.15), все диссипативные члены в уравнениях (16.13), (16.15) и (16.17) имеют одну и ту же форму.

Из рассмотренных выше примеров ясно, что К(чБ-уравнения отличаются от полных уравнений Навье — Стокса только диссипативными членами. 16.1.2. Анализ Фурье качественного поведения решений Чтобы решение укороченных уравнений Навье — Стокса можно было получить за один проход в направлении х, необходимо, чтобы система уравнений, например (16.4) — (16.6), была неэллиптической по отношению к направлению х. Это может быть определено методами, описанными в п. 2.1.4, для чего необходимо построить систему уравнений первого порядка, эквивалентную системе (16.4) — (16.6). Однако больше информации о качественном поведении решения можно получить, если разложить зависимые переменные в комплексный ряд Фурье. Такой подход можно проиллюстрировать на упрошеином уравнении энергии дТ дТ деТ деТ и — +о — — 6 — — е — =О, (16.18) дх ду дхх дух где и и о — постоянные порядка 0(1), но могут быть и положительными, и отрицательными.

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

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

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

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