Главная » Просмотр файлов » Роуч П. Вычислительная гидродинамика

Роуч П. Вычислительная гидродинамика (1185913), страница 95

Файл №1185913 Роуч П. Вычислительная гидродинамика (Роуч П. Вычислительная гидродинамика.djvu) 95 страницаРоуч П. Вычислительная гидродинамика (1185913) страница 952020-08-25СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

Первая точка верхней границы (1, 7) берется на входной границе. В случае пересечения границы В 3 ударной волной или менее сильной волной сжатия в точке (ь,/) могут пересечься две характеристики одного семейства; тогда можно предпочесть методы выделения скачков, однако и более простая процедура, описанная выше, обеспечивает однозначное задание граничного условия. Эту процедуру применяли Аллен [!968], Аллен и Чен [1970], Роуч и Мюллер [1968], а также Гудрич [1969].

В сочетании с тремя различными конечно-разностными схемами во внутренних точках она давала устойчивые н реалистические результаты. Существенно, однако, чтобы граница В 3 не лежала внутри пограничного слоя, иначе можно получить абсурдные результаты. При расчете обтекания затупленного тела Лапидус [1967] применил линейные экстраполяции вдоль диагоналей сетки, по- ложив 420 о.8, Критерии схооимости и иаиальиые условия с успехом использовали на этой границе свои условия па выходе (см, предыдущий раздел). Эрдош и Заккаи [1969[ также применяли на верхней границе свои условия для выходной границы, предполагая, что па пей нет скачков.

5.8. Нритерии сходимости и начальные условия Многие из замечаний разд. 3.4 и в особенности то обстоятельство, что не существует объективных удовлетворительных критериев ни итерационной, ни аппроксимациопной сходимостн, относятся и к течениям сжимаемой жидкости. Вопрос об итерационной сходпмости (об установлении решения по времени) в случае течения сжимаемой жидкости дополнительно усложняется наличием большего числа искомых функций (например, давление устанавливается медленнее, чем плотность) и появлением нового характерного времени — времени прохождения волны давления через расчетную область.

Росс и Чен [1970[ отметили, что в сверхзвуковых течениях вязкого газа можно ожидать очень долгого затухания нестационарных процессов, а это делает суждение об установлении еще более затруднительным. В этой связи можно было бы рекомендовать сравнение окончательных стационарных решений, полученных при раз.тичных начальных условиях (хотя бы для некоторого контрольного варианта задачи).

К несчастью, представляется, что для сверхзвуковых течений задача с начальными условиями будет более критической в смысле устойчивости, чем для течений несжимаемой жидкости. В литературе приводится много примеров неустойчивости при одних начальных условиях и устойчивости прн других. Эта неустойчивость по определению обусловлена нелинейностью, однако, как отметил Моретти [1968а, !968б[, в некоторых случаях источником такой неустойчивости может являться п неправнльная постановка граничных условий.

Верно, однако, и то, что подобные неустойчивости по крайней мере усугубляются, а быть может, н полностью порождаются распространением ложных ударных волн, связанных с неправильными начальными условиями. В том случае, когда интерес представляет только стационарное решение, существуют три способа уменьшить влияние начальных условий. (1) Расчет можно начинать с малыми значениями Лс. Этот способ часто приводит к большим затратам машинного времени. (2) Расчет можно начинать с искусственно заниженным значением числа Рейнольдса, постепенно доводя его в процессе расчета до желаемого значения, (3) Расчет можно начинать по другой схеме, обладающей большим искусственным земпфнровапием. Это легко осуществимо для схем с явной 421 ду, Расчет дозвуковых и еверхввуковых течений искусственной вязкостью.

Аналогично, для схемы Рихтмайера (т. е. для двухшаговой схемы Лакса — Вендроффа, см, равд. 55.6), в которой первый шаг делается по схеме Лакса, а вто! ой — по схеме «чехарда», первые несколько сот шагов можно реализовать только по схеме Лакса и лишь затем пользоваться волной двухшаговой схемой (Лапидус [1967[). При С «! схема Лакса обладает большой диффузией и поэтому даст сглаженное решение. Поскольку при проведении отладочных расчетов нелинейные неустойчивости могут разрушить решение, иа ранних стадиях следует вести решение короткими этапами, попеременно записывая результаты чередующихся этапов на две магнитные ленты. При возникновении неустойчивости можно вернуться к результатам расчета ва предыдущем этапе с меньшей величиной Лй с большиье демпфированием и т, д.

Тайлер и Цумвальт [1965), а также Тайлер и Эллис [1970[ показали, что в задаче об одномерном распространении скачка можно получить более плавные профили скачков, распространяя начальный разрыв иа две расчетные ячейки. Вместо скачка величин от значений перед ударной волной (а) до значений за ударной волной (5) в пределах одной расчетной ячейки Р, = Р„Р,+~ — — Рь (5.18! ) они определяли в промежуточной узловой точке среднее значение Рь = Р Рь+ ~ = (Ро + РьУ 2 Ртч.х = Рь (5 182) (подробности см, в указанных выше работах). Уоткипс [1970[ обнаружил особую важность совместности начальных условий в связи с использованным им преобразованием координат, По опыту автора, а также Л. Д. Тайлера (частное сообщение) обычные двухшаговые методы дают незатухающие колебания примерно в трех узловых точках около начального положения скачка. До настоящего времени ие найдено удовлетворительного объяснения этого явления и способов его устранения.

Гурди и Моррис [1968в[, Вернер [1968], а также Смит и Мак-Колл [1970) для улучшения аппроксимационной сходимогти гиперболических систем применяли методы экстраполяции. 5.9, Замечания о расчете дозвуковых и сверхзвуковых течений Часто возникает мысль о гом, что не стоит пользоваться уравнениями течения несжимаемой жидкости, а достаточно разработать лишь программу для расчета течений сжимаемой жидкости, и тогда соответствующее течение несжимаемой жидкости 59 Рос мт дозвуковых о сзерхзвукоеых течеяаа 422 можно рассчитать по этой программе, считая число М, малым, скажем положив Мо = 0.1. Поэтому, казалось бы, программа расчета течений сжимаемой жидкости обладает болыпей гибкостью.

В общем случае, однако, подобный метод будет и весьма неэффективным, и весьма неточным. Уменьшение эффективности обусловлено очевидным усложнением вида вязких членов (см. гл. 4) н усилением условия на размер шага по времени Лй Значение Л! при расче~е течения сжимаемого газа ограничено условием по числу Куранта О и!+ а) ат ах (5. 184) Выражая при помощи этого равенства коэффициент давления через отношение давлений, получаем Ср (т/2) мт (5.185) что после дифференцирования дает "(Орв) д (Р,) тМа„ (5. 185) ') Заметим, что пра М- О уравнение анергап отделяетса от ураааеаай !теразрыапостп и количества ааажеаая, Прп М, — 0.1 мы имеем а 10)и), и по условию (5,183) Л! ограничено в основном скоростью звука, что уменьшает максимально допустимые значения Л! в десять раз по сравнению со случаем применения уравнений для несжимаемой жидкости.

Более того, паразитные звуковые волны на сетке приведут к возрастанию ошибки, связанной с неразличимостью, и (что, быть может, наиболее важно) ухудшат итерационную сходимость. Черни с соавторами 119501 указал на желательность эффективного отфильтровывання этих паразитных воли и применения поэтому уравнений течения несжимаемой жидкости, Ухудшение точности проистекает из неопределенности прп М- О. Даже при поверхностном знакомстве с газовой динамикой ясно, что для сверхзвуковых течений используются отношения давлений, например отношение давлений в донной области Р, = РвуР, в то время как для течений несжимаемой жидкости используются разности давлений '), например коэффициент давления в данной области Сев = (Рв — Р,) Я , где д — скоростной напор в невозмущенном потоке, т, е.

В!О Схемы высокого ао!ги)ка ляпропсяликаа 423 При у = 1.4 и М = 0.1 получаем отношение относительных ошибок г = 0.007. Это означает, что при вычислении Р, по уравнениям течения сжимаемой жидкости с погрешностью менее 1$ получаем значение Сгв с погрешностью 100%. Для смешанных течений, т.е. для сверхзвуковых течений с областями дозвукового течения, отношения давлений еще имеют смысл п можно пользоваться уравнениями сжимаемой жидкости. Однако задачи с большим интервалом изменения числа Маха по времени очень трудны для расчета.

Примерами таких задач являются разгон тела пз состояния покоя до сверхзвуковой скорости и расчет взрыва от его начала до поздних стадий. Харлоу и Амсден [!968] разработали неявный эйлеров метод расчета движений сплошной среды (метод 1СЕ), дающий хорошие результаты от М = 0 до М » 1. Метод !СЕ основан на расчете уравнения неразрывности по неявной схеме, что придает системе уравнений эллиптический характер [см. Фромм [1963] н Руо [1967]).

В нем нет ограничения на размер шага по времени, связанно~о со скоростью звука. Дальнейшие усовершенствования метода 1СЕ н его приложения указаны в работах Харлоу и Амсдепа [1970], а также Харлоу с соавторами [1971]. 5.10. Схемы высоного порядка аппронсимации Замечания равд. 3.1.10 по поводу ограничений на схемы высокого порядка аппроксимации еще в большей мере относятся к течениям сжимаемой жидкости.

Как уже было указано в равд. 5.5.2, в сверхзвуковых течениях при больших числах Рейнольдса искомые функции нс обязательно непрерывны по пространственным переменным и в этих случаях ряды Тейлора, применяемые для оценки ошибок аппроксимации, непригодны. Берстейн и Мприн [1970, !971] разработали схему расщепления для уравнений течения невязкого газа третьего порядка точности по пространственным переменным и времени.

Схема третьего порядка точности в приложении к задаче о расчете обтекания затупленного тела с отошедшей ударной волной дала более точные значения давления, по менее точные значения плотности в точке торможения и потребовала втрое больше машинного времени, чем схема второго порядка точности. Другая схема третьего порядка точности приводится в работе Русанова [1970]. Глава 6 ДРУГИЕ РАСЧЕТНЫЕ СЕТКИ, СИСТЕМЫ КООРДИНАТ И СИСТЕМЫ УРАВНЕНИЙ Выше мы познакомились с основными понятиями и методами вычислительной гидродннамики иа примере простейших задач, используя некоторые формы уравнений Навье — Стокса в декартовых координатах и предо~валяя их в виде уравнений в конечных разностях на равномерных расчетных сетках с постоянными Лх и Лу.

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

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

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

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