Anderson-et-al-1 (1185923), страница 64

Файл №1185923 Anderson-et-al-1 (Андерсон Д., Таннехилл Дж., Плетчер Р. Вычислительная гидромеханика и теплообмен) 64 страницаAnderson-et-al-1 (1185923) страница 642020-08-25СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

В случае сверхзвуковой маршевой задачи требуется решать уравнения (6.132)— (6.134) с граничным условием непротекания на поверхности дфг'дп = 0 (6.135) и заданным числом Маха в свободном потоке М . Такая постановка задачи существенно проще, чем решение полных уравнений Эйлера. Уравнениеполного потенциала в недивергентной форме было использовано в работах [члгоззшап, 1979; бгоззшап, Б)с!аг1, 1980) для расчета обтекания конических и закрученных с несимметричным профилем дельтовидных крыльев.

Для получения решений была использована недивергентная формулировка и релаксационная схема для трансзвуковых течений. Сейчас мы опишем маршевую процедуру «8)тап)саг, 1981; 8)тап)саг, С)та)сгачаг1Ьу, 1981] решения уравнения потенциала, записанного в дивергентной форме, чтобы пояснить идею линеаризации плотности, использованной в этих работах.

Прежде чем обрисовать в общих чертах всю процедуру, необходимо обсудить конечно-разностную аппроксимацию уравнения потенциала. Так как, считая течение безвихревым, мы вводим потенциал, то в нашей системе уравнений отсутствует механизм диссипации. Вследствие постоянства энтропии уравнения потенциала допускают решения в виде скачков разрежения и уплотнения. В расчете дозвуковых течений нет особенностей. В сверхзвуковых потоках невозможно появление скачков разрежения, и их следует исключить из рассмотрения. Эту трудность можно обойти явно или неявно путем введения искусственной вязкости. где звездочка, обозначающая обезразмеривание функций, опущена, и и = дф/дх, о = дф/ду. (6.133) 5 6.5.

Методы решения уравнения потенциала Н Н Н ° Н Н Направление (( .)) Н Н (1,! ) Н Н. (а) (ь) Рис. 6.21. Аппроксимация, зависящая от типа уравнения; (а) точки эллиптич- ности; (Ь) точки гиперболичности. оси х. Если течение дозвуковое, то уравнение эллиптическое и производные аппроксимируются центральными разностями. Если течение сверхзвуковое, то уравнение гиперболическое в рассматриваемой точке и вторые производные в направлении течения (в продольном направлении) берутся со смещением шаблона вверх по потоку.

Имеем следующие конечно-разностные выражения для вторых производных в точке (1,1): д! 2ов — ь 1+ Ею — 2, У (Ьл)я ав 1+~ еьг — ~ ч'в-ьн1+'в'-ь) — 1 т ху 2 Ьлау Фв 1+~ — ив 1+ Ф~ 1 (ау)я Сеточный шаблон для точек, в которых скорость потока меньше (точки эллиптичности) или больше (точки гиперболичности) скорости звука, показан на рис.

6.21. Изображенные на рнс. 6.21 узлы сетки иллюстрируют корректную зависимость от типа уравнения для до- или сверхзвукового течения. Расположение точек шаблона для стационарного уравнения потенциала говорит в пользу применения неявной (6.136) 12е Мерман и Коул [Мцгшап, Со!е, 1971] в своей знаменательной работе о трансзвуковых течениях указали, что производные в каждой точке области расчета должны быть корректно аппроксимированы в соответствии с типом уравнения. Они рассматривали трансзвуковое уравнение малых возмущений, но эта же идея применима н для уравнения потенциала.

Чтобы проиллюстрировать конечно-разностную аппроксимацию, зависящую от типа уравнения, рассмотрим уравнение в недивергентной форме (6.129). Тип этого уравнения гиперболический в точках, где (иэ+ пэ)/аа — 1) О, и эллиптический в точках, где (пэ+ пя)/ая — 1(О. Пусть поток направлен вдоль 666 Гл. 6. Численные методы решения уравнений течения схемы. Если мы имеем дело только со сверхзвуковым течением, т. е. в расчетной области нет точек эллиптичности, то решение можно получить и по явной схеме. Но этого не рекомендуют делать, если в некоторых точках поля течения скорость потока ненамного превышает звуковую, так как критерий устойчивости Куранта — Фридрихса — Леви запрещает иметь разумные размеры шага. В этом случае применение явных схем, даже для чисто сверхзвуковых течений, становится непрактичным.

Линия Маха ° а ° ° Направление швченыя ° ° ° ° Ланыя Маха Рис. 6.22. Случай, когда направление сетки не совпадает с направлением течения. Если построить дифференциальное приближение для конечно-разностного представления рах в точках гиперболичности, то окажется, что старшие отброшенные члены имеют вид йх (и' — й) 16„„„.

(6.137) Это эквивалентно введению положительной искусственной вязкости в точках, где иа ) иа. Если конечно-разностная аппроксимация (6.136) применяется в точке эллиптичности, то искусственная вязкость становится отрицательной и возникает проблема устойчивости. Джеймсон [,)ашезоп, 1974) указал, что эта трудность возникает в тех случаях, когда поток сверхзвуковой и компонента и скорости в направлении х меньше скорости звука.

Проблему можно понять, рассматривая случай, когда поток не строго параллелен направлению х, как это показано на рис. 6.22. Шаблон разностной схемы неправильно учитывает область зависимости точки (1,1). Одна из точек шаблона (1, 1 + 1) в этом случае расположена перед характеристикой, проходящей через точку (1,1). Чтобы преодолеть эту трудность, Джеймсон предложил схему, учитывающую разворот потока. $ 6.5.

Методы решения уравнения потенциала где з н и — расстояния, измеряемые вдоль линий тока и по нормали к ннм. По правилу дифференцирования сложных функций вторые производные ф„и ф,„выражаются через х н у следующим образом: 1 = — „„(и'ф„„+ 2иоф,„+ ьафаа), (6.139) 1 ф„„= уа (оаф„„— 2ивф„а+ итф „). В разностной аппроксимации для ф„производные по х н у берутся со смещением назад (с запаздыванием), тогда как для ф„ используются центральные разности. Когда поток направлен вдоль одной из осей координат, схема Джеймсона сводится к уравнению (6.!36) н вводит искусственную вязкость с главным членом вида (1 — —,) (Лзиаф„, + ...).

(6.140) Позтому мы имеем положительную искусственную вязкость во всех точках, где поток сверхзвуковой, н можно надеяться, что скачки будут скачками сжатия. Так понятие искусственной вязкости используется для объяснения свойств решения уравнения полного потенциала. К тем же самым выводам можно прийти путем тщательного анализа конечно-разностных уравнений, рассматривая разные его члены. В работе [На1ег е! а1., 1979) при рассмотрении трансзвуковых течений применена ндея искусственной сжнмаемости для введения искусственной вязкости в сверхзвуковых областях течения.

Эту идею впервые предложнл Хартен [Наг1еп, 19781, пытаясь улучшить методы сквозного счета для сверхзвуковых течений. Для получения требуемой искусственной вязкости в некоторых работах [НоЫ, Ва!!Ьапз, 1979; Но!з1, 1979] использовалось смещение шаблона плотности вверх по потоку. Описанный ниже метод включает в себя этн идеи и очень полезен прн решении уравнения полного потенциала.

Чтобы понять, как через смещение шаблона прн аппроксимации плотности или искусственную сжнмаемость вводится искусственная вязкость, полезно рассмотреть одномерное уравнение потенциала (6.141) Идея состоит в том, что уравнение потенциала записывается в потоковых координатах в виде (а' — !уа) ф„+ а'ф„„= О, (6.138) Зоз Гл. 6. Численные методы решения уравнений течения Его аппроксимация имеет второй порядок и выражается в следующем виде: Р(Р„п,йф,) =О, (6.142) где обозначения имеют тот же смысл, что и ранее. В точках эллиптичности уравнение (6.142) корректно. В точках гипербо- личности следует добавить искусственную вязкость так, как это делает Джеймсон [.)атезоп, 1975[: — Лх (1аф„,)„, (6.143) где О, н = 1п)п ( е„') (6.144) Как говорилось выше, такое явное добавление искусственной вязкости эквивалентно введенной Мерманом и Коулом [Мпгтпап, Со!е, 1971[ разностной аппроксимации, зависящей от типа уравнения. Джеймсон [Яатезоп, 1975[ показал, что член (6.143) эквивалентен члену вида (6.145) — Лх (тр„ф,)„, где О, и=шах ~~х (6.146) Такое представление получают дискретизацией одномерного уравнения энергии.

Если искусственную вязкость в таком виде ввести в уравнение потенциала, то конечно-разностная аппроксимация уравнения (6.141) будет выглядеть так: — (р — ) = Ч [р,, айф,~ — Ч ~и,(р...— р,, а)Аф,~ =О. (6.147) Было показано [НоЫ, Ва!!Бааз, 1979[, что она имеет второй порядок точности и в дозвуковых зонах эквивалентна аппроксимации уравнения центральными разностями. В сверхзвуковых зонах из-за добавления искусственной вязкости конечно-разностное представление (6.147) является схемой первого порядка с разностями вверх по потоку. По мере увеличения числа Маха шаблон все более смещается вверх по потоку. В дозвуковых зонах шаблон производных плотности не смещается.

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

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

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