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

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

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

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

(5.170) Этн два способа не эквивалентны. Опишем вкратце другие способы постановки граничных условий на выходе (здесь важно подчеркнуть, что к постановке граничных условий на выходе нельзя подходить легкомысленно). Как было установлено в равд. 3,3.9 по поводу парадокса на выходе для квазиодномерных уравнений, граничные условия иа выходе играют большую роль в случае сверхзвукового течения на входе, чем в случае дозвукового течения на входе. Вычислительные эксперименты Крокко [1965], а также проведенные автором расчеты двумерных задач показали, что именно граничные условия на выходе определяют положение скачка ').

Еше более важен следуюший факт. Рассмотрим часто повторяемое утверждение: в случае сверхзвуковых скоростей на выходе граничные условия на выходе не сушествениы. Это утверждение неверно. Как мы видели в равд. 3.3.9, если бы не было влияния этой границы, то было бы невозможно «выключить» аэродинамическую трубу со сверхзвуковым течением на входе. Более того, это утверждение неверно даже в смысле малости ошибок. Аллен 11968] обнаружил, что при расчете течения вязного газа около обратного уступа использование простого граничного условия на выходе (7,=(7,, (5.171 а) илн, что эквивалентно, (5.171б) приводит к монотонной расходимоств решения даже при полностью сверхзвуковом потоке на выходе.

Руо ]1967], наоборот, ') Для задачи об обтекании обратного уступа приемлема линейная эксграполяция (б.!70) Однако на очень грубой сетке н том случае, когда рас. стояние выходной границы от уступа было вчетверо болыие его высоты, а за уступом продолжалась прямзя стенка с условиями првлппания, решение «раз. наливалось». Хотя граница расчетной области нахолилась вне области вторичного сжатия потока, ошибка от грани шаго условия на выходе вызывала появление сильной ударной волны с осцилляциямн, распространявшейся от выходной границы п разрушавшей замкнутую застойную область.

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

Описанное выше явление дает еще одни пример весдинственности решения задач вычислнтельной газо. динамики. 5.7.6, Вьыойтгч ероннци 4!5 нашел, что для задач о дозвуковых течениях этот способ дает сходяшиеся решения, и даже утверждал, что он ведет к ускорению сходнмости по сравнению с другими способами. Итон и Цумвальт (1967], а также Кссслер [1968] обнаружили, что этот способ дает хорошие результаты в том случае, когда скачок проходит через границу. (В последних трех работах применялась схема Русанова, описанная в равд. 5.4.3.) Причиггы различных отзывов о способе постановки на выходе этих условий неизвестны. Итон и Цумвальт [1967] нашли, что в случае пересечения ударной волной выходной границы применение линейной экстраполяции (5.170) приводит к неустойчивости даже в том случае, когда первоначально поток был сверхзвуковым.

Ясно также, что в случае пересечения хвостовой волной выходной границы В 6 использование на ней линейной экстраполяции может привести к серьезным ошибкам, если ие к возникновению неустойчивости. Если отвлечься от случая выхода ударной волны на границу В 6, то способ линейной экстраполяции (5.170) оказывается, вообше говоря, приемлемым; его использовали Сондерс [1966], Лапидус [1967], Итон и Цумвальт [1967], Скоглунд и Коул [1967], Руо [1967], Аллеи [!968], Роуч и Мюллер (1968], Скоглунд и Гей (1968], Аллен и Чен (1970].

Аллен (1968] экспериментировал также с квадратичными и кубическими экстраполяциями вида )г гу-! (7-2+ ~г-ъ' ~У-1 7г-2 + ~г-в ~1-4 В своих последних расчетах (Аллен и Чен [1970]) он пользовался соотношением (5.172), но обнаружил, что разница между способамп (5.172), (5.173) и (5.170) незначительна. Аллен [1968] также менял положение выходной границы и не обнаружил значительной разницы в двух полученных решениях — одном частично дозвуковом (около границы В! на рнс. 3.22) и другом полностью сверхзвуковом на выходной границе,— разрешив таким образом вопрос о необходимости сверхзвукового потока на выходе (см.

также численные эксперименты с выходной границей Росса и Чена (1971]). Выяснилось, что условие сверхзвукового потока на выходе не является ни необходимым, ни достаточным условием для обеспечения устойчивости и точности. Руо (1967] экспериментировал с осредиснием граничных условий (5.171б) и (5.170) и с введением сдвига по времени '): (5.! 74) ') Для модельного уравнения (5.1) способ [5.174) при числе Куранта С = 1, конечно, дает точное решение. 4!6 б 7 Граничные условия для течения гясил~авлгой пчдкагти он обнаружил, что лучше использовать условие (5.17!б), как это было указано ранее.

Эрдош и Заккаи [!969[, а также Шихна с соавторами [1970) применяли квадратичную экстраполяцию, причем полиномы строились по значениям в последних пяти узлах сетки по методу наименьших квадратов. Представляется, что этот способ дает не ббльшую точность, чем менее сложные экстраполяции '). При пересечении границы ударной волной здесь также развивалась неустойчивость. В расчетах невязких течений по схеме Русанова (равд. 5.4.3) Итон и Цумвальт [1967) вернулись к схеме с разностями против потока в направлении х, причем вторые производные, входящие в члены с искусственной днффузией, на выходной границе полагались равными соответствующим производным в узлах, лежащих на расстоянии одного шага от границы: д„[п,„д ) ~ = б [п,„б ) ~ ° (5175) Производные по у в точках выходной границы вычислялись при помощи стандартных аппроксимаций, принятых во внутренних точках расчетной области.

Как и экстраполяция (5.171б), этот способ давал хорошие результаты даже в случае пересечения выходной границы ударной волной. Роуч и Мюллер [19681 рассмотрели сходный способ постановки граничных условий на выходе, аналогичный их способу для течений несжимаемой жидкости (равд. 3.3.7). Конвективные члены уравнения количества движения в направлении х аппроксимировались по схеме с разностями против потока. Диффузионные члены с пропзводными по х, члены со смешанными производнымп н составляющая градиента давления по х вычислялись в точках ! — 1; сама по себе эта процедура порождает тенденцию к дестабилизации расчета, которая подавляется за счет сдвига по времени. Как и в случае расчета течений несжимаемой жидкости, члены с производными по у могут вычисляться на входной границе при помощи стандартных аппроксимаций, принятых во внутренних узлах.

1-1апример, уравнение количества движения в направлении х (4.426) может иметь следующий конечно-разностный аналог: ') Как отметил '!сп [!970), любые зкстраполяпии основаны на значе. пнях, полученных в результате расчета, а пе на точных значениях и, сяе. Новательно, их точность ограничена. Б 7 7 Верхняя гранииа где б/бх п б/бу означают центральные конечные разности, а Л(рпо)/Лу — вклад в д(ри)/д1 от потока в направлении у, вы шслепный по схеме, принятой во внутренних узлах расчетной области (этот способ может применяться только в сочетании со схемой с разностямн против потока во внутренних точках).

В цитированной выше работе во избежание усложнения программирования все вязкие члены, вкл!очая б(я(ба/бу))/бу, вычисля- лись прп 7 — 1 и и — 1. Представляется, что этот способ должен давать хорошую точность, однако численные эксперименты Аллена !1968) продемонстрировали незначительную разницу между различными способами экстраполяпии, обеспечивающими устойчивость.

В соответствии с этим для расчета течений сжимаемой жидкости при больших числах Рейпольдса можно рекомендовать более простой способ (5.170) при условии, что выходную границу В 6 не пересекает сильная ударная волна. Прн больших значениях составляющей скорости о на выходной границе точность расче~а можно было улучшить, применяя в сверхзвуковой области приближение простой волны (см. следующий раздел), однако до настоящего времени такие расчеты ие проводились. 5,7.7. Верхняя граница Как н в случае течепия несжимаемой жидкости, верхнюю границу (В 3 на рис. 3,22) можно трактовать как простую стенку с прилипанием плн (еше лучше) как стенку аэродинамической грубы без трения. Однако в этих случаях можно ожидать нежелательного влияния отраженных от границы волн (см., например, Унлкипс !19691), Гораздо лучшим способом является воссоздание условий «свободного полета», что возможно в случае сверхзвуковых течений.

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

С Физической точки зрения мы рассматриваем поток через верхнюю границу как не- вязкое стационарное гомоэнтропическое течение, сформировавшееся прп расширении первоначально однородного потока около искривленной поверхности, т. е. течение Прандтля — Майера (см., например, Овчарек 1!964]). Соответствующая геометрия показана па рис. 5.7. Точка р и узловая точка верхней границы (й 7) лежат на одной и той же характеристике, или линии Маха. Эта линия определяется углом 418 5.7. Граничные условия длл течений слсимавмтд жидкости вперссвсвки Рис.

5.7. Постановка условия на верхней границе при помощи приближения простой волны. а — приближение простой волны для течения вблиан верхней границы; (иг 1„, 1„= й... + фбх) Вас-1 — й-се,); б — действнтелы нос течение типа простой волны. ()ам+ 0), с направлением оси х, причем рм = агс з|п(1/М)— угол Маха ') и 0 =агс10(о/и) — угол наклона вектора скорости к оси х. Величины углов рм н 0 в точке р, а также значения всех других параметров течения в этой точке )р — — р„ис, ос, Т, вычисляются линейной интерполяцией по их значениям в узлах (1 — 1,1 — 1) и (й 7 — 1). В простой волне параметры течения постоянны вдоль характеристики.

Поэтому, если положение точки р определено, то условие на верхней границе запишется в виде )с, 1 = )а + — (1, — 1 ), где ) = р и о Т. (5.177) Положение точки р можно определить следуюшнм образом В узле Ь вЂ” (с — 1, У вЂ” 1) вычисляется величина со=1д()тат+ О). ') Очевидно, что рм не существует при а4 ч.!. Д7.7. Верхняя граница 4!9 Если ььь ) Ауггбх, то точка р находится между точками Ь и с, как это показано на рис. 5.7. В точках Ь и с вычисляется вели- чина га' = !п[п/2 — (!ьм + 6)], а затем расстояние 1 находится из соотношения дх/Лд — мь (гьс ььь)/~» + 1/~у (5. 178а) Если гьь < Ау/Ах, то точка р лежит выше точки (1 — 1, У вЂ” 1) и отстоит от нее на расстояние 1, причем др/ах — мь (м — мь)/ар+ !/ах ' (5.

1786) /ьг= /ь-ьг-~ /г-яг-г (5, 179) что равнозначно приближению простой волны в случае совпадения характеристик с диагоналями расчетной сетки. На ранних стадиях вычислений метод (5.179) приводил к неустойчивости и был заменен следующим: /,,=/. ь., +г(/. ь., — /...,). (5.!80) Для обеспечения устойчивости Лапидус в течение первых 500 шагов по времени медленно менял г от г = 0.5 до г = 1 (что дает формулу (5.179)). Итон и Пумвальт [1967] обнаружили, что квадратичная экстраполяция типа (5.172) приводит к неустойчивости при пересечении границы В 3 ударной волной ц Здесь величина гаа вычисляется в предыдущей точке верхней границы, причем процесс построения граничного условия осуществляется последовательными шагами слева направо.

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

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

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

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