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

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

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

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

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

Уэбб с соапторамп [1967) проводил расчеты иа ЭВМ с удвоенной точностью для уменьшения ошибок округления и был вынужден сглаживать начальные данные для обе печения устойчивости. Их схему нельзя считать удовлетворительной, так как неясно, чем определяется полученное решение: физикой явления нли сглаживапием начальных дапных. При помощи обратного метода Брнггс [1960) решал задачу об обтекании эллиптяческого конуса под углом атаки. Много вопросов, связанных с обратным методом, затронуто ть работе Пауэрса с соавторамп [1967]. Джонс [1968) разработал обратный метод для конических тел, 5 ДД Рвечет течений дев ударивших волк, вмделекие ударных волн 337 в котором вместо итераций по отклонениям в форме тела (на поверхности которого нормальная составляющая скорости обращается в нуль) рассматриваются повязки этой составляющей скорости в точках желаемого положения тела. При этом Джонс исходил из точного решения Тейлора — Макколла (см.

Овчарек (1964) ) для обтекания конического тела при нулевом угле атаки, вводя малые возмущения угла атаки или формы обтекаемого тела. Обтекание затупленного тела под углом атаки рассчитывали обратным методом Махин и Сягаев (1966). Кайрис 119701 разработал метод, в котором объединены конечно-разностная схема на смешанной лаграпжево-эйлеровой сетке для расчета течения без скачков и обратный метод для расчета скачка. Много алгоритмов обратного метода и метода характеристик используется в промышленности (см, работу Морено (1967], содержащую также библиографию по этой теме).

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

Такие методы носят название методов «размазывания» скачка (з)тос)г-згпеаг(пд) илн методов «улавливаниях скачка (зйос1<-сар1пг(пд) в англоязычной литературе и методов сквозного счета в литературе на русском языке (см., например, Годунов и Семепдяев (1962)). Если нас интересует только стационарное решение задачи, то можно проводить не- стационарные расчеты и выйти на стационарное решение (если оно существует) при больших значениях времени (так же, как в гл. 3).

В литературе на русском языке этот подход называется методом установления (Браиловская с соавторами 11968]) или методом асимптотичсского стационнрованпя. Конечно, если нестациоиарное решение не представляет интереса, то при вышеуказанном подходе появляется дополнительная гибкость, связанная с отсутствием необходимости моделировать реальные физические иестационарные процессы. При таком подходе, нанример, можно использовать схемы первого и меньшего порядка точности по времени. Крокко (1966) предложил схему, в которой конечно-разностпые уравнения фактически б.б Предварительные соображению аппроксимнруют исходную систему уравнений в частных производных лишь в стационарном случае, и эта схема обладает хорошей сходимостью.

Симуни (см. Браиловская с соавторами [1968] ) разработал метод установления с рассмотрением зависящих от времени граничных условий,стремящихся к граничным условиям исходной задачи только при достижении стационарного решения. Фройдигер с соавторами [1967] также рассматривал нестационарный расчет лишь как итерационную процедуру и указал на возможность изменения Ггг по пространству для ускорения сходимости.

Поскольку ни одна из этих идей пе была развита в систематический метод, они не будут обсуждаться в настоящей книге, однако их следует иметь в виду для дальнейшего развития метода установления. 6.1.2. Исследованнв устойчивости Замечания, сделанные в разд. 3.1.5.г по поводу критериев устойчивости и методов ее исследования, остаются в силе и здесь.

Дополнительные сведения, касающиеся, в частности, устойчивости в гиперболических системах, можно найти в работах Куранта, Фридрихса и Леви [1928], Лакса [1954, 1957, 1958, 1961], Лакса и Рихтмайера [1956], Лакса н Вепдроффа [1960], Крейса [1964], Стренга [!964], Лакса и Нирепберга [1966], Парлетта [1966] и Бранловской с соавторами [1970]. Определения и условия, рассматриваемые в этих работах, конечно, не всегда справедливы для нелинейных задач, в особенности прн наличии скачков. Исследовать устойчивость системы уравнений, подобной (4.55), труднее, чем в случае течения несжимаемой жидкости. В системе уравнений с переменными ф, с исследовалось на устойчивость независимо одно параболическое уравнение, а влияние ф устранялось с помощью липеаризацнн.

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

Поэтому модельное уравнение ди дгг дггг — =- — й — + а— дГ дх дх' ' (5 1) па котором часто изучается устойчивость уравнений. описывающих течение несжимаемой жидкости, не моделирует поведение ЗЗ9 б!.2. Исследование дстопчиаости уравнения количества движения сжимаемой жидкости (4.42б) 11аже в случае одного измерения, так как в нсм нет членов с градиентом давления ').

Для исследования устойчивости полной линеаризованной системы уравнений нужно воспользоваться матричными методами исследования устойчивости (см, Эймс (1969], Митчелл (1969) ), что выходит за рамки настоящей книги. Используем, однако, прием, который позволит сделать некоторые выводы об устойчивости уравнений для сзкимаемого газа па основании простого модельного уравнения (5.1). Как мы видели в гл.

3, типичным условием устойчивости для уравнения (5.1) является ограничение на число Куранта С =~ й)Л1/Лх(1, (5,2) Это неравенство следует пз рассмотрения певязкпх членов в уравнении (5.1), т. е. прп сс = О. В уравнении (5.1) прп сс = О информация переносится с копвективпой скоростью непрерывной среды. На каждом расчетном шаге по времени возмущение в Лй точке влияет па новое значение в (с'+ 1)-й точке.

Это означает, что за каждый шаг по времени Л1 информация переносится на расстояние Лх и, таким образом, вьшислигельная скорость распространения информации равна Лх/Лб Неравенство (5.2) означает, что скорость распространения информации в непрерывной среде нс должна превышать вычислительной скорости распространения информации, т. е. )й ~(«Лх)ЛЛ (5.3) В течениях сжимаемой жидкости члены с градиентом давления меняют скорость распространения информации в среде; она уже пе равна с), а несколько больше. Рихтмайер и Мортон [1967) словесно описали путь получения зависимости скорости распространения информации в среде от давления.

Мы же здесь просто ограничимся злементарпымп газодинамическими соотношениями. Малое возмущение давленпя распространяется с местной скоростью звука а относительно газа, который сам движется со скоростью й. Возмущения давления распространяются во всех направлениях, п необходимо рассматривать толысо а ) ) О. Таким образом, скорость распространения информации в сжимаемой жидкости равна 1й1+ а, и тогда ограничение на число Куранта записывается в виде С= ( бх ') К счастью, однако, вязкие члены уравнений для сжнмаемоб зкидкости, содержащие смешанные ироизводные, здесь ие существенны, и в этом отношении уравнение (бд) моделирует уравнение (Ллйб); см.

равд. б.б. З4О б.т Предваригельнь>е соображения 1 1 С =- = плп С(= 1з ч>3 (5. 4б) соответственно для двух или трех простраиствеипых коордииат '). Использование схемы расшеплеиия по времени Марчука [1965! с раздельным расчетом вкладов от каждого измереиия приводит обычно к одномерному условию (5.4а); см.

также работы Гурли и Митчелла [1969б), а также Мак-Кормака [1971~, В литературе приводилось еше одно дополнительное ограиичеиие па шаги расчетной сетки. Гудрич [1969), пользовавшийся иестациоиариой схемой Русакова (равд. 5.4.3) и полпыми урав- ') В более обмяк случаях ямеем ду дкдг с~< иле С ~< — (бые) задке+ ду' ~дльдрх ! духдее ! дехдхх Это и есть условие устойчивости Кураита — Фридрихса — Леви [1928), или КФЛ-условие. Переформулироваипое в других терминах, оио констатирует, что область влияния коиечио-разиостиых уравнений должна бь~ь по меньшей мере столь же велика, как область влияния исходных диффереицпальиых уравнений (см.

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

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

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

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