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

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

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

Самая мелкая сетка состояла из 65 Х 17 узлов. Для расчета трансзвукового обтекания препятствия в канале потребовалось 900 шагов по времени, чтобы достичь стационарного состояния, т. е. (14.117) и (14.122) использовались только на самой мелкой сетке. При использовании описанного выше многосеточного алгоритма требуется 130 многосеточных циклов с сокращением времени счета примерно в 4 раза. Тот же алгоритм был использован для расчета течений около осесимметричиой гондолы и каскада лопаток ротора турбины.

Гл. 14. Неиязкие течения л1б В работе [Зо)1пзоп, 1983] описано обобщение алгоритма Ни на двухшаговую явную схему типа Лакса — Вендроффа (п. 14.2.2); в работах [Рач!з е! а1., !984; СЫта, Ло!1пзоп, 1985] применен метод типа Ни для решения сжимаемых уравнений Навье — Стокса. Использована также дискретизация по методу контрольного объема уравнений Эйлера с расположением узлов сетки в центрах контрольных объемов [Дашезоп, 1983]. Вместо описанной выше схемы Лакса — Вендроффа использовалась четырехшаговая схема Рунге — Кутты ($ 7.4). Многосеточный алгоритм Джеймсона аналогичен алгоритму Ни, за исключением того что использовался лишь один цикл типа зуба пилы при переходе на самую грубую сетку и интерполяция на самую мелкую сетку проводилась через все промежуточные сетки. В работе [Мц!бег, 1985] скомбинирована неявная схема маршевого интегрирования по времени с многосеточным методом и методом расщепления потоков [5!едег, вагш!пй, 1981].

Чтобы получить более крутые фронты скачков, использовался метод ограничения потоков. Применялась также многосеточная схема РАБ (п. 6.3.5) с релаксацией Гаусса — Зейделя к стационарным уравнениям Эйлера при помощи дискретизации методом конечного объема второго порядка с расположением центров контрольных объемов в узлах сетки [Нешкег, !986]. Быструю сходимость многосеточного метода можно использовать в схемах ТЧГ! первого порядка (п. 14.2.6) до момента, когда сходимость будет почти достигнута. Ограничения потоков второго порядка можно вводить лишь для окончательного определения стационарного решения. 3 14.3.

Трансзвуковые невязкие течения Выделение трансзвуковых невязких течений в отдельную категорию связано с тем, что такие течения можно рассматривать специальным образом на основе уравнений потенциала, если образующиеся в иих скачки малы. Трансзвуковые невязкие течения характеризуются наличием зон до- и сверхзвукового течения (рис. 11.15). Основное внимание будет уделено методам решения стационарного уравнения потенциала. Дополнения, необходимые для решения нестационарного уравнения, будут приведены в краткой форме. ' 14.3.1.

Общие замечания Если течение безвихревое, можно ввести потенциал скорости (11.102). Уравнения Эйлера тогда можно свести к одному уравнению в частных производных и вспомогательному алгебраиче- $ !4.3. Трансзвуковые невязкие течения (14. 128) р'= Р =11+0.5(у — 1)Мз„[1 — М вЂ” Жт (14.! 29) где у — отношение удельных теплоемкостей, а М вЂ” число Маха набегающего потока.

Если для (14.128) используется консервативная дискретизация, то решение будет соответствовать слабой форме (14.128) (например, (5.6)). Следовательно, «ударные волны» будут скому уравнению. Для двумерного стационарного течения уравнение для потенциала скорости имеет вид (а' — из) Ф„„— 2ипФ„в + (а' — о') Фвв —— О, (! 4.126) скорость звука а определяется выражением (11.104). Если (14.126) записать в естественных координатах (з, п), т. е.

з параллельно локальному направлению потока, а и перпендикулярно ему, то уравнение примет вид (1 — М') Фаз + Ф„„= О, (14.127) где локальное число Маха М = д/а и дз = и~+ оз. Очевидно, уравнение (14.127) превращается из эллиптического в гиперболическое уравнение в частных производных, если число М становится локально больше единицы, т. е. течение становится локально сверхзвуковым. Как следует из рис.

11.15, сверхзвуковая область в направлении течения обычно заканчивается ударной волной. Поскольку вывод (14.126) основан на предположении о том, что течение безвихревое, то по теореме Крокко [Е!ертапп, Коз%о, 1957] течение при обтекании тела однородным потоком должно быть изэнтропическим. Однако, хотя условия Репнина — Гюгонио (!1.110) гарантируют сохранение массы, импульса и энергии при переходе через скачок, энтропия увеличивается пропорционально третьей степени от интенсивности скачка (М~ — Мв)з). Если число Маха для нормальной компоненты скорости перед скачком (М~) меньше 1.1, изменение состояния течения при переходе через скачок можно считать приблизительно изэнтропическим и течение можно описывать уравнением потенциала.

Уравнение (14.!26) обычно рассматривается в эквивалентной безразмерной консервативной форме (т. е. уравнение неразрывности) (р'Ф'„) + (р'Ф„')„= 0, где Ф'=Ф/((/ Л), 7.— характерная длина. Безразмерная плотность получается из (11.104): Гл. 14. Невязкие течения 2!8 улавливаться решением. Однако в разрывных решениях будут сохраняться энтропия, энергия и масса, но не импульс. Поэтому условия Ренкина — Гюгонио будут выполняться лишь приближенно. В работе 1.)ашезоп, 1978) отмечается, что скачок импульса на рассчитываемом разрыве является мерой волнового сопротивления. Численное решение уравнений (14.128) и (14.129) можно получить весьма эффективным образом.

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

14.2.8, 14.2.9) решать полные уравнения Эйлера (11.22) — (1!.24). Эти методы обычно значительно медленнее методов, предназначенных для решения уравнений (14.128) и (!4.129). Экономичность последних особенно заметна при расчете трехмерных или нестационарных течений. Поскольку уравнение (14.126) описывает изэнтропическое не- вязкое течение, его решением могут быть как ударные волны сжатия, так н ударные волны разрежения (рис. !4.28). В численных схемах должны быть предусмотрены специальные процедуры, делающие невозможным появление (физически невозможных) скачков разрежения. Это может быть сделано путем использования разностей против потока или введением искусственной вязкости в сверхзвуковых областях. Оба эти метода создают диссипативный механизм, который делает невозможным образование ударных волн разрежения.

14.3.2. Трансзвуновое уравнение малых возмущений В случае обтекания тонких тел возможно дальнейшее упрощение уравнения (14.126), в результате чего можно получить уравнение (11.107) для возмущения потенциала чР. Если про- $ !4.3. Траясзяукояые яевязкяе течеяяя 219 филь тела задается уравнением у=т/(х), где т предполагается малой величиной, уравнение (11.107) можно представить в виде [Со!е, 1975] — [Кфк — 05(У+1)(фх)']+ д з =О, (14130) где К= (1 — Мз )/тмз, ф„=- дф/дх. В уравнении (!4.!30) величины у и ф приведены соответственно к масштабам тц' н т-зм. Уравнение (14.130) называется гранезвуновым уравнением малых возмущений. Совместное с (!4.130) условие на ударной волне имеет вид [фз] —,У [Кф„— 0.5(У+!) [ф„]з] = О, (14.131) где [ ] означает скачок величины, а ду/дх — наклон ударной волны.

Можно также заменить граничное условие дф/да=0 на поверхности тела условием на оси х: — — при у=О. дФ д/ ду ах (14. 132) Уравнения (14.130) — (14.132) правильно описывают трансзвуковые течения около тонких тел, за исключением области, расположенной непосредственно перед передней затупленной кромкой, поскольку в этой области возмущений скорости и имеют тот же порядок, что и скорость набегающего потока (/ . Для сверхзвуковых скоростей, например М ) 1.3, возможно дальнейшее упрощение (14.130), поскольку при этом 0.5 (у + + 1) (фх]з « Кф„. Уравнение (14.!30) сводится к линейному уравнению (11.109), для решения которого можно использовать панельный метод (например, РАЙ-А1Рт), описанный в п. 14.1.6.

Консервативную форму (14.130), необходимую для правильного расчета скачка, можно компактно представить в виде д +д ду дО (14. 133) где Р = Кф, — 0.5 (у + 1) (ф.)', 6 = дф/ду. Центрально-разностиые операторы, определенные через значения в полуцелых точках сетки, имеют вид ["1+ца з "1-цз, я] [Рь яз цз Ог з — цз] Риз- — ' — . Яня= ах ' ьу (14. 134) Если значение 6, определенное вновь при помощи центральноразностных операторов, подставить в (14.134), то (91, ят-! ~~Ь Й + ти Ф вЂ” 1) Яг,з= стук Гл.

14. Невязкие течения Аналогичное выражение можно записать для Рь «. При помощи (14.134) алгоритм, удобный для дискретизации (14.133) при описании трансзвуковых течений, можно представить в виде Р1 «+ Я! « — )«1 «Р1 «+ )«1, «Рг, „—— О, (14.135) где ~ = О в дозвуковых точках, т. е. К ) (у + 1) ф„, )«1, « т =! в сверхзвуковых точках, т. е. К ((у+ 1) ф„. Из (14.135) следует, что в сверхзвуковых областях член дР/дх в (14.133) аппроксимируется по трехточечной направленной против потока схеме (рис. 14.29), а в дозвуковой — центральной (а (Ь к+1 Н-! ззис. 14.29, Точки, используемые в (14.13о). (а) Лозвуковая точка (й «); (Ь) сверхзвуковая точка (й «). разностью. Отличие формул в сверхзвуковой области правильно отражает гиперболическую природу уравнений в этой области, заключающейся в отсутствии влияния вверх по потоку.

На звуковой линии или на ударной волне (см. рис. 11.15), расположенными между точками (1 — 1, )з) и (1, )з), уравнение (14.135) принимает вид На звуковой линии: 1;)1 « = О, (14. 136) На ударной волне: Рг «+ Р1, «+(41 «= О. :Эти специальные представления необходимы для правильного описания ударных волн и для построения эффективного алгоритма решения уравнений.

Преобразуя выражения (14.135) через потенциал скорости, получаем систему уравнений, которая может быть решена методом итераций. Специальные процедуры описаны в п. 14.3.5. Изложение данного метода следует в основном работе [Мцгптап, 1973). Разложением в ряд Тейлора Р; ь» в окрестности узла (1, й) можно показать, что этот оператор аппроксимирует дР/дх в данной точке, но приводит к появлению искусственной вязкости, 5 14.3.

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

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

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

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