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

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

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

В используемых обозначениях однородные граничные условия Неймана для р будут определены в узле (1, 2) вблизи границы АВ на рис. 17.3. Дискретизированное центральными разностями относительно узла (1, 2) уравнение Пуассона может быть представлено в виде — + -' ~= -- + оч! ле! Р!, г + Рг, г Р!, ! Р!. 2 + Р!. 3 (Р!, 2 Ро, ' ) дхг + Луг /Зхз Г Ро Рл ба йо + ~ 3/2 2 п2 2 + ! 5/2 ! 3/2 ~ ( 1 7 . 1 8) Л/ Л» ду 400 Гл. !7. Несжимаемые вязкие течения Граничное условие Неймана для давления с центром в узле (1/2, 2) может быть получено из (!7.8) в виде «т! се и+! (Р!, з Ре, 2) !Ги 2 !Ги а (17.19) ьх а! Если использовать уравнение (17.19) для исключения (р!,ив — рв,а) из (17.18), получится выражение, не зависящее от Г!..

и, следовательно, не зависящее от значений и и о вне области, появляющихся в уравнении (17.9). Поскольку решение не зависит от Г!тз,м в уравнении (17.18) может быть сделана подстановка Г!м, и = и!!и,'и и, согласно (17.9), ре, е = р!, ь Таким образом, вводится однородное граничное условие Неймана для давлениЯ, Член ии!тез! известен из гРаничных Условий.

Описанный метод в основном совпадает с методом, описанным в книге 1Реуге1, Тау!ог, 1983], где отмечены также важные аналогии между методом МАС и методом проекции (17.22) — (17.24). Для многих зависящих от времени задач ограничение на шаг по времени (17.15), связанное с использованием явных формул (17.8) и (17.10), является слишком обременительным. Обобщение метода МАС, позволяющее проводить интегрирование по времени уравнений (17.2) и (!7.3) с использованием неявной приближенной факторизации членов, содержащих скорость, сделано в работе (!)еи!!!е, 1974] для очень малых чисел Рейнольдса и в работе [ОЫа е! а!., 1979] для течений с большими числами Рейнольдса. Общее описание можно найти в книге ]Реуге1, Тау!ог, 1983]. Дискретные уравнения для определения поля скоростей (17.8), (17.!О) могут быть записаны в символическом векторном виде ц"+' =Г" — Ажир ', где Г = (Г, ст)', Чи — разностный оператор градиента.

Уравнение Пуассона для давления (17.13) может быть записано в форме (!7.20) !7йр"" = „!, т Г . (17.21) Уравнения (17.20) и (17.21) являются краткой формой записи метода МАС. Существуют иные методы, подобные методу МАС. С методом МАС тесно связан метод проекции, предложенный в работах ]Сапог!и, 1968; Тешат, !969]. В принятых обозначениях в методе проекции уравнение (17.20) разделяется на два этапа: и'= Г", (17.22) па+! = и Ест!зри+!.

(17. 23) $17.1. Исходные переменные: нестанионарные течения 40! Подстановка (17.23) в уравнение (17.12), записанное в виде 11ап'+' = О, приводит к соотношению 112 рп.~-! 1 111 (17.24) г[1+ АГ (и1.„— — 1,„,) ~ и' = ц", '[ т' + Ы [ п7.„— — 7.„н) ] ц" = и', [т' + АГ (х '7, — ~~ 7ае) 1 и* = и". (17. 25) Уравнения (!7.25) являются приближенной факторизацией, аналогичной рассмотренной в $ 8.2 и 9.5 и позволяющей получить последовательность трехдиагональных систем уравнений, если Е„ 7., и т. д. — трехточечные конечно-разностные операторы.

Года отмечает, что для устойчивости решения необходимо использовать ограничение типа КФЛ на шаг по времени: А! ( :е~ Ах/] йгпах ] Другой вариант реализации метода МАС (17.20) и (17.21) предложен Хиртом и Куком [Н(г(, СооК !972]. В принятых обозначениях предварительные значения скорости определяются из уравнений (17.2) и (!7.3) в виде и* = Р" — Агрдр".

(! 7.26) 26 К. Фчетчер. т я В результате этого п"чч удовлетворяет как уравнению (17.1), так и уравнениям (17.2) и (17.3). В методе проекции из уравнения (17.22) определяется и", из (17.24) определяется р"е' и из (17.23) определяется н"е'. Изначально метод проекции был разработан на неразнесенной сетке. Однако в работе [Реуге(, Тау!ог, 1983] рекомендуется использовать его на разнесенной (МАС) сетке (рис. 17.1). Из уравнений (17.22) и (17.24) видно, что метод проекции совпадает с методом МАС во внутренних точках. Однако граничные условия реализуются несколько по- иному. В работе [Реуге1, Тау!ог, 1983] показано, что в методе проекции решение не зависит от значений и* на границе.

Это по существу эквивалентно исключению анка из уравнения (17.18) после подстановки в него (17.19). Года [Сода, 1979] использовал метод проекции для расчета вязкого течения в двумерной и трехмерной движущихся полостях. Чтобы избежать явного ограничения на шаг по времени (17.15), уравнение (17.22) заменено уравнениями 402 Гл. 17. Несжимаемые аяакие течения Поправка к давлению бр = р"+' — р" определяется из уравнения ( 17. 27) Эта поправка обеспечивает выполнение уравнения неразрывности при подстановке в него ц"+', т. е. и"+' = и' — — 74 б,п. .З1 (17.28) (17.29) Наконец, новое значение давления равно Рч+! рй + бр Постановка граничных условий осуществляется, как в методе МАС (п.

17.1.3). Хирт и Кук использовали такой подход для исследования вязких несжимаемых (ламинарных) течений у трехмерных структур. В работах (5а)сагпо1о, Ма1цо, 1980; Ка1о, Мпга1сапп, 1986) использовался метод Хирта и Кука для исследования нестационарных турбулентных трехмерных течений, возникающих в задачах вентиляции помещений. Применялась (й — е)-модель турбулентности (п. 11.5.2).

Сравнение экспериментальных н расчетных данных для этой задачи приведено на рис. 17.4. Данные результаты получены на сетке 20(х) к', 24(у) Х 15(г). Это самая грубая сетка, на которой удалось удовлетворительно получить основное циркуляционное и вторичное течения. Течение вызывается нагнетанием воздуха через крышу, Вдуваемая струя ударяется об пол и вызывает циркуляционное в плоскости симметрии течение у стен (рис. 17.4(а) и (с()). Хорошо моделируется картина вытекания воздуха у пола и крыши (рис. 17.4(с), (е)); хорошо видны возвратные течения, вызванные нагнетаемой струей. В приведенном выше описании семейства методов МАС предполагалось, что границы расчетной области совпадали с декартовыми координатными плоскостями.

Для расчетных границ произвольной формы можно ввести связанные с границами криволинейные координаты (гл. 12), преобразовать уравнения к криволинейным координатам и применить метод МАС в регулярной расчетной области. В работе (Ра1е1, Вг1йдз, 1983] применен метод МАС в первоначальном явном варианте для расчета в криволинейных координатах задачи о естественной конвекции в нестационарном двумерном ламинарном течении.

Использовалась разнесенная сетка, на которой коитравариантные компоненты скорости (12.65) определялись на границах ячеек, а давление — в их $17.!. Исходные переменные: нестационарные течения 403 центрах. Поскольку в эквивалентном (17.13) уравнении появляются производные от давления Ро и рто для вычисления этих производных по значениям в соседних точках использовались две перекрывающиеся сетки. Это приводит, однако, к тому, что Иллмсч(мрвммеи плоскости (а1 и1 Гариеончеленеи праекни чипе 1 Молелнроеенив (е) (Уе И').лласкааче Эксперимент (е) (О»М).плоское~ Поперечное сечение чипе 1 ° и л и овос) Эксперимент (Ы (оеи(. Молелироовние (ы (о м) Кя (Рв измласкссв Ч".)» н чр р л мл(р тля т нилл о л °,- — в вмн »лаяло о( ' ' ', '; О Молепираеение (Ы (о в и) .плоскость Эксперимент (о] (ив Р 1 "омск (И (о+ ч)1-пласкасче Рис.

17ип Моделирование вентиляции комнаты Яка1о, Мига(сат(, 19861( пе. чатается с разрешения Яарап зос. о1 Совр. Г!оЫ 1)упаш(са). система уравнений Пуассона, эквивалентная (17.!3), в два раза больше, чем в случае применения метода в декартовых коорди- натах. 77.1.5. Разности высокого порядка против потока В первоначальном методе МАС используются центральные конечно-разностные формулы (17.7). Для течений, в которых основную роль играет конвекция, использование центральных 404 Гл. 17.

Неезкииаечые вязкие теяеиив — < 1.0 и — < 1.0. и Л1 ив) Гзх ду Многомерные разностные формулы третьего порядка точности против потока являются обобщением одномерной квадратичной интерполяции против потока Леонарда [1еопагб, !989]. Одномерная схема третьего порядка с разностями против потока может быть продемонстрирована на уравнении переноса (9.56), записанном в консервативном виде дТ д (иТ) даТ вЂ” + — — а — =О, д) дх дха (17. 30) где и известно и изменяется внутри области. Консервативное разностное представление уравнения (17.30) может быть запи- разностей приводит к появлению в решении сильных нефизических осцилляций (п.

9.3.1). Делались попытки стабилизировать решение путем дискретизации конвективных членов двухточечными разностями против потока (п. 9.3.1) или суммами с весами центрально-разностных выражений и разностей против потока. Однако, как правило, при этом получается неточное решение, особенно если локальное направление вектора скорости совпадает с направлением сетки и велики локальные градиенты скорости. Более точное решение получается, если для аппроксимации конвектнвных членов использовать разности высокого порядка точности против потока, подобные четырехточечным формулам, рассмотренным в п.

9.3.2 и 9.4.3. Дэвис и Мур [Оач!з, Мооге, 1982] использовали исходные переменные для решения нестационарных уравнений Навье— Стокса методом, аналогичным методу МАС. Применялась разнесенная сетка; уравнение Пуассона для давления решалось на каждом временнбм шаге прямым методом (п. 6.2.6) [Зтиаг1г!гацЬег, 1974]. Отличительной чертой метода Дэвиса и Мура является использование многомерных разностей третьего порядка точности против потока для аппроксимации конвективных членов.

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

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

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

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