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

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

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

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

Ларсен [1969) обсуждал возможности применения суперэллнптнческой системы координат. Эта ортогональная система координат основана на кривых Ламе, определяемых уравнением (6.34) При различных значениях параметров а, Ь и п эти кривые образуют прямоугольники, эллипсы, ромбы, звезды, квадраты, круги, а также промежуточные фигуры.

Крайне сложной системой координат являются рассмотренные Копалом [1969) координаты Роша, при построении которых используются эквнпотенциальные поверхности вращаюнгегося гравитационного дпполя. Он предложил использовать эту систему координат при расчете газовых потоков в замкнутых бинарных системах (звездах). Можно настойчиво рекомендовать учебник Берла с соавторами [1960), который содержит уравнения газодинамики с учетом вязкости, записанные в прямоугольной, цилиндрической н сферической системах координат, а также много другой полезной информации по газодипамнке и другим процессам переноса.

Цянь Сюз-сень [!958] приводит уравнения Навье — Стокса для течения сжимаемой вязкой жидкости в ортогональных криволинейных координатах. (Однако пн в книге Берда с соавторамн, ни в работе Цяпя не приводится консерватинная форма уравнений.) В работе Богачевского с соавторами [1965) дана консервативная форма уравнений течения певязкой сжимаемой жидкости в цилиндрических и сферических координатах. (На.

помним отмеченный в гл. 4 факт, что введение коисерватнвнгчх 6.8. Другие ортогонакьные системы координат 445 переменных влияет на вид вязких членов в уравнении энергии для сжимаемой гкидкости.) Применение различного рода систем координат, включающих в качестве координаты радиус (цилиндрической, сферической, параболоидной и др.), требует некоторых пояснений. Что касается свойства консервативности, здесь имеет место некая двусмысленность. В случае осевой симметрии уравнение неразрывности сжимаемой жидкости в цилиндрических координатах имеет впд (Берл с соавторами [!960)) — + — — (гро) + — (ри) = О.

др 1 д д д< г дг д» (6.35) (6.36) где Ч = рг. Аналогичные соотношения справедливы и для уравнений количества движения и энергии (см. Богачевский с соавторамп [!965)). Заметим, что в этой дивергентной форме значение зависимой переменной на оси симметрии (г=О) просто т) =О. Однако для вычисления градиента давления в узлах, отстоящих на один шаг от оси симметрии при < = 2, и для вычисления некоторых вязких членов необходимы значения основных переменных на осн симметрии. В осесимметричпом случае наиболее надежным способом их получения является запись уравнений контрольного объема для конвективных членов на оси симметрии.

(Заметим, что на оси симметрии, очевидно, о = 0,) Уараганение. Показать, кто уравнение неразрывности, ззг1нсзпкос на осн симметрии (1 = <), имеет внд др д<ри) 4 <ро)з<з — + — + =0 дг дк Лг Такой подход обеспечивает консервативность, однако недает возможности вычислить вязкие члены. Вопрос сщс более услож. Применяя любую из схем с разностями против потока или с центральными по пространственным переменным разностями (см.

гл, 5), можно показать (задача 6.2), что это уравнение консервативно в смысле сохранения массы при протекании из одного цилиндрического контрольного объема в другой. Однако эта форма уравнения не является дивергентной формой, которая также имеет отношение к консервативности и важна с теоретической точки зрения (Лаке [1954)).

Дивергентная форма уравнения неразрывности получается при введении консервативной переменной т) = рг, как в следующем упражнении. Укражнение. Показать, что уравнение (6.35) можно записать в дквер. гентной форме — + — (чо) + — (чи) = О, дч д д дг дг д» 446 6 4. Другие системы уравнений няется при использовании цилиндрических или сферических координат при отсутствии осевой симметрии. Довольно грубым, но, очевидно, приемлемым подходом является возврат кдекартовым координатам вблизи центральной линии. В настоящее время удобной методики учета особенности на центральной линии не имеется. Ясно, однако, что прием размещения узловых точек на расстоянии половины шага от центральной линии (г = О при 1 = '/,, например) и использования здесь тех же, что и во внутренних точках, аппроксимаций некорректен. Конечно-разпостные уравнения не могут быть записаны с переходом через линию г = О, потому что при г = О направление возрастания г меняет. ся на обратное.

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

При помощи схемы Русанова (равд. 5.4.3) Итон (личное сообщение) рассчитывал осесимметрнчные вихревые течения и обнаружил, что ошибки можно значительно уменьшить, положив диффузионные члены на центральной линии равными нулю. 6,4. Другие системы уравнений В настоящем разделе мы будем обсуждать следующие вопросы в порядке перечисления: (1) Системы уравнений, являющиеся настолько сильными упрощениями системы уравнений Навье — Стокса, что прп этом меняется тип уравнений. (2) Не столь радикальные упрощения системы уравнений Навье — Стокса. (3) Системы уравнений более сложные, чем уравнения Навье — Стокса.

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

Исзорпчески в этом случае ранее всего был получен обширный класс решений, причем болынинство пз нпх в замкнутой форме (см. любой курс 6.4. другие сисгеггвг уравнений 447 гидродинамики). Для тел сложной формы часто предпочтитель. нее рассчитывать потенциальное обтекание численно (см., например, Дуайер и др. [1971]). Однако для решения осесимметрнчных задач со свободными поверхностями численный подход действительно необходим. Бреннен [1969] рассчитал осесимметрнчное течение с каверной, Джепсон [1969] рассматривал осесимметричные течения со свободными поверхностями, Биззел с соавторами [1970] рассчитал истечение жидкости из осеспмметричпого бака, Константинов [1970] н Уитни [1971] рассматривали двумерные пестационарные течения со свободными поверхностями.

Другим классом задач, для которых необходимы численные методы, являются течения с завихренностью невязкой несжимаемой жидкости (как в классической задаче о подъемной силе профиля, помещенного в свободный сдвиговый слой); такие задачи рассматривал, например, Чау с соавторами [1970]. В работах Хауэлла и Спонга [1969], а также Гельдера [!971] численно решалось уравнение потенциала скорости для дозвукового течения сжимаемой жидкости; в последней работе учитывалось также влияние острых углов, В случае сверхзвукового стационарного течения невязкой жидкости уравнения становятся гиперболическими.

Задача сводится к задаче Коши для пространственных переменных, и здесь становится применимым широко известный мощный метод хирактеристик. Основная теория и численные алгоритмы метода характеристик для двумерных гомоэнтропическнх течений можно найти в различных учебниках; см., например, !11апиро [1953], Липман и Рошко [1957], Аббот [1966], Овчарек [1964], Чепмен н Уолкер [1971]. Представляет исторический интерес оригинальная работа Толлмнна [1949]. В первоначальном варианте метода характеристик положения узловых точек сетки не известны заранее, а являются частью решения, Существуют две разновидности метода характеристик: так называемый метод волн (илн метод ячеек) и метод узлов характеристической сетки, В методе волн несколько упрошается арифметика вычислений прн выполнении их вручную, н он более нагляден физически, в особенности при наличии в задаче границ с постоянным давлением.

Однако метод волн требует перехода в плоскость годографа скорости, приводит к вычислительным трудностям на границе, являюшейся линией симметрии, и не может быть обобщен на осесимметричные течения (поскольку основывается не на строгой математической теории характеристик). Поэтому ему следует предпочесть метод узлов характеристической сетки, Узлы характеристической сетки определяются пересечением характеристик (лпний Маха) различных семейств, т, е, линий. е 4 Другие системы ураеиеииа составляющих с вектором скорости углы ~!л, где 9= = агс з!п(17М), В курсе Липмана и Рошко [!957] изложены также метод характеристик для осесимметричпых течений н усложнения, возникающие при отсутствии гомоэптропичности (т.

е, в том случае, когда энтропия постоянна лишь вдоль линий тока, а по нормали к ним может быть переменной). Однако многие другие проблемы, связанные с получением решений по методу характеристик, пе обсуждаются в имеющихся учебных руководствах. Так, в учебниках отсутствует описание наиболее современного варианта метода характеристик, согласно которому решение продвигается вниз по потоку на заданный шаг сучетом условия устойчивости Куранта — Фридрихса — Леви [1928], а характеристики выпускаются вверх по потоку, причем необходимые величины определяются с помощью интерполяции. На возникновение скачка прн расчете по методу характеристик указывает пересечение характеристик (линий Маха) одного семейства. В этом случае должен быть выделен косой скачок, угол наклона которого определяется совместным решением соотношений Рэнкина — Гюгонио и характеристических соотношений за скачком. Выделение скачков в расчетах плоских задач по методу характеристик изложено в работах Хартри [1958) и Ричардсона [1964] с обсуждением вопросов программирования.

Кеннеди [!956), Вейс с соавторами [!966], Морено [1967], Аббетт [1970] рассматривали выделение скачков в расчетах осесимметрпчных задач по методу характеристик. Метод характеристик удобен для расчета одномерных нестацпонарных певязких течений, Такие расчеты выполняли Шапиро [1953), Хоскин [!964), Хоскпн и Лембурн [1971), причем в работе Хоскина [1964) обсуждаются также вопросы программирования. В методе характеристик пе применяются консервативные схемы, и поэтому проверка консервативности может служить указателем величины ошибок аппроксимации. Пауэрс и О'Нейлл [1963) указали, что ошибки, связанные с нарушением свойства консервативности, возрастают при гиперзвуковых скоростях (при малых углах Маха и больших градиентах), и дали метод определения энтропии в узлах сетки при помощи расчета потока массы.

В случае стационарных трехмерных течений возникает определенная свобода прп разработке метода характеристик, что привело к появлению пяти или шести различных вариантов трехмерного метода характеристик. Чу [1964) приводит простой вывод трехмерных характеристических соотношений. Чушкин [!968) дал обзор четырех вариантов трехмерного метода характеристик и учел влияние неравновесности. Пауэрс с соавторами [!967] разработал трехмерный метод характеристик и включил в свой алгоритм расчет пограничного слоя.

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

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

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

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