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

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

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

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

— Ддр(дх и влияние вперел от точки (х, 0 распространяется вдоль характеристики, проходящей через точку (х, !) и имеющей наклон !/д 306 а7. Метода~ решения уравнений для физических перененных Используя схему с разностями против потока, показать, что схемная вят.

кость в случае уравнений для простейшнх фнзнческнх переменных та же, что н в случае уравнения переноса вихря. (См. уравнения (3.)73) — (3.179) нз равд. 3.!.8.) 3,7,6, Сравнительные достоинства систем уравнений для переменных (ф,т.) и для переменных (и,,п, Р) Сравнительные достоинства (тр, Д -системы и (и, и, Р)-системы зависят от решаемой задачи.

Главную роль всегда играет опыт предшествующих расчетов, но при выборе системы уравнений мы увидим, что в большинстве случаев (за исключением задач со свободной поверхностью или других задач о движении жидкостей с поверхностями раздела) целесообразно брать (ф, с) -систему. В качестве эталонной задачи для сравнения этих двух систем рассмотрим сначала задачу о плоском течении жидкости при отсутствии свободной поверхности, предполагая, что прн этом уравнение Пуассона решается при помощи итерационных методов. Если ие требуется находить нестационарное решение для давления, то в (тр, ~)-системе приходится решать одно уравнение переноса вихря параболического типа и одно уравнение для функции тока эллиптического типа т)етр = ~ с условиями Дирихле на некоторых (возможно, па всех) границах. (Стационарное решение эллиптического уравнения для давления находится только на последнем слое по времени, и поэтому выбор метода решения этого уравнения не имеет особого значения,) В (и, о, Р)-системе надо решать два уравнения переноса количества движения, имеющих параболический тип, и одно уравнение эллиптического типа для давления таР = 5р с граничными условиями Неймана на всех границах.

При решении уравнения переноса вихря 1 необходимо дополнительно выполнить две операции дифференцирования функции тока ф для нахождения составляющих скорости, но уравнения переноса количества движения усложняются нз-за членов с дивергенцией бс 7 (в методе МАС эти члены значительно сложнее) и из-за специальных приемов, которые здесь требуются для обеспечения сохранения массы (объема). Решать уравнение переноса вихря ~ можно по неявным схемам, хотя при этом может потребоваться дополнительный итерационный процесс для неявно~о вычисления значений ~"ь' на стенках при условии прилнпания.

В случае же (и, о, Р)-снстемы значения и"ы и и"тт известны точно в тече. ние всего времени, но здесь существует трудность, связанная с неустойчивостью из-за нелинейности (см. равд. 3.7.2). Дости. жение итерационной сходимости прп решении уравнения ЧЯР = = 5р эллиптического типа требует значительно больше времени, Здхб. Достоинства систем уравнений дяя Ц, т) и дяя (и, о, Р) 307 чем при решении эллиптического уравнения >у>>)> = ь. Это объясняется различием граничных условий.

Если же требуется получить также нестациопарное решение для давления, то и в (>)>, с)-системе необходимо решать уравнение Пуассона ииР =5 с граничными условиями Неймана. В случае когда применяются неявные схемы и требуется вычислять поле давления на каждом шаге по времени Л1, результаты можно получить быстрее решением (и, с, Р)-системы. Однако заметам, что при решении (>)>, ~)-системы при помощи явных схем (которые, как сложилось исторически, чаше применяются для решения (и, и, Р)-системы) шаг по времени Ы настолько мал, что значение давления можно не рассчитывать на каждом шаге по времени, а находить только время от времени. (Во всяком случае, обычно оказывается затруднительным разумно использовать все это множество значений давления.) Если поле давления рассчитывается один раз за каждые десять шагов по времени илп резке, то снова рекомендуется применять (ф >)- систему.

Сравним эти две системы и в том случае, когда нужно получить нсстационарную картину линий тока. В случае (>)>, ~)-системы ливии тока >у = сопя| строятся с помощью простой интерполяции. В случае же (и, о, Р)-системы они определяются интегрированием, причем наиболее точный способ состоит в решении уравнения Пуассона Уе>у = би/бу — бо/бх. Теперь рассмотрим случай, когда в исходной тестовой задаче для решения уравнения Пуассона используются прямые методы.

Здесь важную роль играет время разработки программы (оно, конечно, зависит от предшествующего опыта)„и если надо вычислять поле давления, то время разработки программы для решения (ф, ~)-системы будет больше, так как при этом необходимо при помощи прямых методов решать два уравнения Пуассона с различными граничными условиями. Если же поле давления вычислять не надо, то время разработки программы для решения (>(>, ь)-системы будет несколько меньше, поскольку в большинстве прямых методов поставить условия Днрихле проще, чем условия Неймана.

(Метод расчета распространения век~ора ошибки из равд. 3.2.8 является исключением.) В этих случаях для решения уравнения Пуассона требуется меньше времеви, чем для решения параболического уравнения переноса. Поскольку время решения одного уравнения переноса вихря ~ меньше, чем время решения каждого нз двух уравнений количества движения, и в этом случае (>(>, ~) -система оказывается предпочтительнее. Учет внешних сил в уравнениях не влияет на обсуждавшиеся выше сравнительные достоинства двух рассматриваемых систем уравнений.

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

Например, если жидкость пмест переменную вязкость, то уравнение количества двилсения можно переписать так: да ди дн — +и — +о — = дт дх ду = — — д + — (2 — (И дя)+ д (И(д + д ))), (3.607) где Це = Рот'о(4М, И = Й/Иа = И (х, У), (3. 608) а Ио — хаРактеРнаа вЯзкость (скажем, вЯзкость в невозмУщенном потоке). Соответствующее уравнение переноса вихря можно вывести (см, задачу 3.34) в следующем виде '): — = — Ч ° (Уй) + дС дт При помощи (чр, (.)-системы можно исследовать и гораздо более сложные задачи; Шавит и Ливан [!971) решали уравнение переноса вихря в цилиндрической системе координат для течения двух жидкостей с переменвыми свойствами.

Еще раз повторим, что определить положения маркеров в случае (ф, в)-системы ничуть не сложнее, чем в методе МАС для (и, о, Р)-системы, и поэтому возможность получать картину пиний о~меченных частиц в методе МАС не является преимуществом (и, о, Р)-системы перед (чр, В)-системой. Картину же линий тока, несомненно, проще строить при решении (ф,~)-системы. Однако исторически сложилось так, что при решении задач со свободной поверхностью нли задач с поверхностями раздела жидкостей рекомендуется брать (и, о„Р)-систему, поскольку именно таким образом чаще удавалось получить хорошие результаты.

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

следующий раздел). 3.8. Трехмерные течения где уг)= —, + —, + —,, ~=и, о, ит. (3.612) дт( дт! дт) Эти уравнения для физических переменных (составляюших вектора скорости и давления) можно решать теми же методами, что и уравнения в случае плоских течений. Некоторые обобщения на случай пространственных течений уже рассматривались в равд. 3.1. Например, ограничение на числа Куранта для явных схем (при отсутствии расщепления по времени) записывается так: С„+ С„+ С,<1, (3.613) или — + — + — ( 1. иит оЫ вИ Ьх Ьу Ьг (3.614) Можно вывести и уравнение Пуассона для давления (3.615) в правую часть этого уравнения по-прежнему входит член д!9/д! (Уильямс (!969) ).

Современные вычислительные машины дают возможность рассчитывать некоторые трехмерные течения жидкости. В случае пространственного течения несжимаемой вязкой жидкости уравнения Навье — Стокса, соответствующие уравнениям (3.509) — (3.510), имеют вид ди + д(ит) + д(ио) + д(ив) — др + ) трг (3 610 ) д( дх ду дг дх Яе до 1 д (ио) 1 д (о') 1 д (ов) др 1 1 ттго (3 6105) д( дх ду дг ду йе — + дв д (ив) д (ои) д (в') дР 1 д( дх + — + = — — + — чгв, (3.610в) ду дг дг Гче Й вЂ” = Ч ° 'ч( вв — + — + — = О, ди до дв дх ду дг (3.6! 1) 3!О 8.В.

Трехмепиме течении или (в покомпонентной записи через орты 1, 1, !г) в виде Ь=~„!+ в".„1 + 1;,К, (3.6.17 а) див дв х — ду дх' ди двв дх дх ди ди дх ду ' (3.617г) (3.617б) (3.617в) Если течение происходит в плоскости (х, у), то ~, = ье = О, а ~, = до/дх — ди/ду, что отличается знаком от ранее принятого нами определения вихря ~. (Определение ь = — ~, в случае плоского течения наиболее распространено, но не универсально.) Соответствующее уравнение переноса вихря представляет собой векторное уравнение, три его составляющих записываются так; д~" = Ч ' (Чвгх) + ие Ч гвх + Г Чи, (3.618а) дье 1 д, = — Ч (Ч1,)+жЧ21,+1 Чп, — ' = — Ч (Чье) + — Ч=ьх + ь ° Чш (3,618в) (3.6186) На разностной сетке в методе маркеров и ячеек составляющая скорости ш определяется в точке (й 1, й +'/з) и т.

д. Для трехмерного уравнения Пуассона также ставятся граничные условия Неймана; введение дополнительной пространственной переменной существенно снижает скорость сходимости, и в связи с этим желательно рассчитывать решение при помощи прямых методов. Даже в этом случае время, необходимое для решения задачи, очень велико, Здесь большое значение приобретают контрольные и пробные расчеты на грубой сетке. Например, Уильямс (1969] для конвективных членов применял схему Аракавы (равд.

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

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

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

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