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

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

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

Следовательно, для стационарных течений уравнения (18.6) можно упростить, заменив уравнение энергии алгебраическим соотношением. Это можно сделать следующим образом. Суммарная энтальпия Н =(Е + р)/р. Поэтому стационарное уравнение энергии (18.6), (18.7) может быть записано в виде — „(риН)+ — (ров)= д (и „,+ „„— д,)+ д д д + д (от„„+ ит„„— Яу).

(18.24) д Используя (18.9), (18.10) и соотношение для идеального газа и=с,т+ — ,'(и+ ), (18.25) уравнение (18.24) можно преобразовать к виду + д ~'роН вЂ” (р + р ) д 1=и:НЗ, (18.26) где д т ди 2 дв т д Г дв 2 ди т ЙН8 = — ( 1хо — — — 1хи — ) + — ри — — — 1хо — ) + дх ч ду 3 ду ) ду [, дх 3 дх ) + д И И р )+„~(4 )1„д 1+ + "(-'- — ',)1" В++И ('-+)' + ф Ц1х (1 — — ) + 1хт(1 — 1, )~ и ф 11 . (18.27) $ !В.1. Фнзнческне упрощения 477 Для течений около тел и в однородном потоке при больших числах Рейнольдса вязкие и турбулентные эффекты существенны лишь в тонком слое вблизи тела и в следе.

Из анализа порядков величин различных членов в уравнении (18.27) следует, что только член ау ~~Р(1 Рг1'+!гт[1 Рг )~и д сравним по величине с членами, входящими в уравнения им- пульса и неразрывности, и с левой частью уравнения (18.26). Для чисто ламинарных течений предположение Рг = 1.0 приво- дит к обрашению в нуль выражения (18.28). Для турбулентных течений !гг » р и предположение Рг,=1.0 обращает (!8.28) в нуль. Можно отметить, что для воздуха Рг = 0.7 и Ргг = 0.90, так что малый вклад все же остается.

В областях отрыва член (18.28) является доминирующим и приведенные выше замеча- ния применимы и в этом случае. В невязкой области, т. е. вдали от твердых поверхностей, можно пренебречь всеми членами в правой части уравнения (18.27). Следовательно, для стационарных сверхзвуковых течений (18.26) можно заменить уравнением + д [РРН вЂ” ( Р + Р ) д ~ =О. (18.29) Очевидно, что уравнение (18.29) выполняется при Н =сонэ! или, согласно (18.25), при с 7'+ 0.5(и'+ о') =сопя(, что можно представить в виде т Р— + 0.5 (ие + ое) = У Р + 0.5(уе . (18.30) т Р Уравнение (18.30) является алгебраическим уравнением, связы- ваюшим значения р, р, и и о.

Данная связь заменяет уравнение энергии в системе (!8.6), (18.7). Такой способ описания рас- сматривался в работе [Вг(1еу, МсРопа16, 1977] и использовался Флетчером и Сринивасом [Р!е(серег, Ьг(п)чаэ, 1985[. 18.1.3. Приближение тонкого слоя Как отмечено в п. 18.1.2, при больших числах Рейнольдса вязкие и турбулентные эффекты сушественны лишь вблизи твердой поверхности и в следе. Если в течении нет больших Гл. 18. Сжимаемые вязкие течения 478 отрывных зон, на основе сравнения порядков величии в направлении, примерно совпадающем с направлением течения, можно пренебречь многими диссипативными членами, т. е.

составляющими т и ('„1 в формулах (18.9) и (18.10), Для трехмерных течений из-за ограничений, связанных с объемом памяти компьюте- (а) м Рис. 18.2. Измельчение сетки у твердой поверхности, (а) Декартовы коорди- наты; (Ь) обобщенные координаты. ров, обычно невозможно построить достаточно мелкую во всех направлениях сетку, на которой бы правильно аппроксимировались все члены трехмерных уравнений, аналогичных уравнениям (18.9) и (18.10).

Эти два свойства (одно — физическое и второе — расчетное) сочетаются в приближении тонкого слоя !Ва!бту(п, 1 оптах, 1978]. Мелкая сетка используется лишь в направлении нормали к твердой поверхности (рис. 18.2). На грубой сетке в направлении, параллельном стенке (направлении к), невозможно правильно представить производные по х, входящие в т и С) в (18.7). Однако на основе сравнения порядков величин эти члены могут быть отброшены. В приближении тонкого слоя уравнения (18.6) и (18.10) принимают вид — + — + — + — =О, дч др~ дб~ дС д( дх ду ду (18.31) где рт (рц рви ! р рыв (Е ! о) и)г б'=(рв, рив, рв'+ р, (Е+ р) )), б'= — (О, т„„, с„„, (тхаи+ т„„о — (,)„")Г. (18.32) 479 4 18.2.

Явные схемы В уравнениях (18.32) "р — (!к+!кг) д 'рр — З (!к+!кг) д З р~ дв 2 Рт 'к дт Яр = — (е+ ср— Ргг ) дд Можно заметить, что при подстановке различных членов б' в (18.31) получаются лишь производные по у, а смешанные производные исключаются. Это облегчает построение неявных схем (п. 18.4.1). Приближение тонкого слоя обычно используется в обобшенных координатах (гл. 12 н п.

18.4.1). Обобщенные координаты для большинства геометрий выбираются таким образом, чтобы физическая диссипация была существенна лишь в одном направлении. Однако для течения вблизи точки пересечения двух стенок следует сохранить диссипативные члены в направлении нормали к обеим стенкам. В приближении тонкого слоя при этом отбрасываются некоторые смешанные производные с порядком величины, равным оставленным членам.

На практике это мало влияет на общую точность решения [Нипц, Котс(и!! а, 1984) . Приближение тонкого слоя для стационарного течения можно ' рассматривать как укороченную форму уравнений Навье— Стокса (гл. 16). По существу то же самое приближение используется при построении пространственного маршевого алгоритма для расчета вязкого сверхзвукового течения (п. 16.3.1). Однако в приближении тонкого слоя не требуется, чтобы и было больше нуля. Вместо этого в приближении тонкого слоя для определения стационарного решения используется метод установления ($16.4).

Поэтому данным методом могут быть рассчитаны течения с небольшими отрывными зонами в направлении потока. Приближение тонкого слоя широко используется в расчетах [СЬацззее, 1984) и позволяет получать результаты, хорошо совпадающие с экспериментальными данными. $18.2. Явные схемы В п. 18.2.1 будет рассмотрено применение явной схемы МакКормака для решения сжимаемых уравнений Навье — Стокса. Хотя данная схема и экономична при расчете существенно не- стационарных течений, ограничение КФЛ с числом Куранта, равным единице, делает ее менее удобной при расчете стационарных течений. Схемы Рунге — Кутты, позволяющие использовать большие шаги по времени, описаны в п.

18.2.2. 480 Гл. !8. Сжимаемые вязкие течения !8.2.1. Явная схема Мак-Кормака Наиболее широко используемой явной схемой для решения сжимаемых уравнений Навье — Стокса является схема Мак-Кормака (МасСоппас!4, 1969]. Для одномерного невязкого течения эта схема описана в п.

14.2.2, а для многомерного иевязкого течения с пространственной переменной, играющей роль времени,— в п. 14.2.4. Применение схемы Мак-Кормака для решения уравнений со вторыми пространственными производными будет продемонстрировано здесь на примере двумерных уравнений Бюргерса (10.57) и (10.58), При обобщении на сжимаемые уравнения Навье — Стокса (18.6) — (18.10) возникает дополнительная трудность, связанная со смешанными вторыми производными, которая будет рассмотрена позже. Двумерные уравнения Бюргерса записываются, как и в (10.52), в виде — + — + — — 8=0, дч дР дС д1 дх ду (18.33) где ит — и т(и(т(х ч Г ив — и ди(с(у Г= ). С=~' -(' ) =( ио — т йо!т(х ) и' — и Йо!т(у 0.5и (из+ о')/и 8= 0.5о (и'+ оз)/и =(:) (! 8.34) Применительно к системе (18.33) схема Мак-Кормака иа однородной сетке имеет вид Шаг предиктор (1 8.35) Шаг корректор 4!0+ 0 5(4)и + а ) 0 5 (Г" Г' — 0.5 —," [С",.

„— С;,,1+ 0.5518",, (18.36) Можно заметить, что каждая пространственная группа Г или С на шагах предиктор и корректор аппроксимируется односторонними конечно-разностными операторами. Вся схема имеет второй порядок точности по времени и пространству, если для аппроксимации производных, фигурирующих в Г н С (см. (18.34)), используются разности, обратные используемым для Г и С в уравнениях (18.35) и (18.36).

Например, на шаге $18ха Явные схемы предиктор (Р )о (Н2)о Ч !Ч! В Ь й) ах Ьх (18.37) Применение разностей вперед на шаге предиктор и разностей назад на шаге корректор можно поменять местами. Порядок использования разностей может быть различен для разных пространственных направлений.

Однако для сохранения второго порядка точности необходимо сохранять симметрию разностных формул на шагах предиктор и корректор. Вопросы устойчивости схемы Мак-Кормака для решения уравнений (!8.33) и (18.6) рассматривались в работе )Реуге1, Тау!ог, 1983); точных результатов ие получено. Поскольку схема явная, можно ожидать, что невязкая часть уравнений должна приводить к ограничению типа КФЛ, аналогичному (9.11), а вязкая часть — к ограничению диффузионного типа (п. 7.1.1). Для скалярного эквивалента уравнения (!8.33) в работе [Реуге1, Тау!ог, 1983) рекомендуется при бх = Лд использовать следующее необходимое условие устойчивости: Ахт 42+ ( ! и )+ ! о 1) ах для ламинарных сжимаемых уравнений Навье — Стокса (18.6) рекомендуется при Лх = бу использовать следующее условие (18.40) 31 К. Флетчер, т. 2 Для сжимаемых уравнений Навье — Стокса (18.6) можно непосредственно применять схему предиктор — корректор (18.35), (18.36), положив 8 =0.

Однако в выражениях (!8.7) — (18.10), определяющих значения Г и С, фигурируют производные по разным направлениям. Например, чл2 Рмо (1+122)! д + л ) (18.38) Дискретизация ди/ду проводится так же, как в (18.37). Для до/дх используются центральные разности. Таким образом, иа шаге корректор 0.5(о1„~ 2 ~ — о! ~ 2 ~)) ах Гл. !8. Сжимаемме вязкие течения устойчивости: Л/ (~ о ь о ь, (18.41) (2Гь/Пе Р) [22/Рг+ (2/3)о'ь) + [! а!+ ! о 1+ (2) ' а) ах где а — скорость звука. Из вида условий устойчивости (18.40) и (18.41) следует, что в трехмерном случае ограничение на шаг по времени более жесткое, чем в двумерном и одномерном случаях.

Чтобы избежать этого ограничения, Мак-Корман [МасСоггпасК 197!] в описанную выше схему ввел расщепление по времени. Расщепление по времени, аналогично описанной в 9 8.5 процедуре, приводит к последовательности одномерных пространственных операторов. Для уравнения (!8.6) одномерные операторы могут быть представлены в виде ь(/,'е = Р„(/)/„'/а/ „ ь)/'ь = Ро (Ид) и/' Здесь уравнение (18.42) эквивалентно уравнениям Аналогичные выражения для (18.43) могут быть получены из уравнений (18.35) и (18.36). Весь алгоритм определения решения на новом временнбм слое вместо (18.35), (18.36) принимает вид е+~ Р '[ о)Р [ х)Р [ ")Р ( — ")йе (18.44) Симметричная картина повторяющихся пространственных операторов от А//2 необходима, чтобы получить алгоритм второго порядка по времени.

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

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

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

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