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

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

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

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

Фромм также опробовал линейную экстраполяцию для постановки нь на выходной границе и обнаружил, что такое условие обладает дестабилизирующим свойством в случае явных схем для уравнения переноса вихря. Автор настоящей монографии, применяя явные схемы, получил аналогичный результат. В дальнейшем, применяя метод последовательной верхней релаксации для решения уравнения Пуассона, Кемпбелл и Мюллер [1968] и Роуч и Мюллер [1970] с успехом ставили условия Томана и Шевчика (3.478) в целом ряде задач. В работах Хына и Макано (Хын и Макано [1966], Макано и Хын [1967]) задавались следующие условия на выходной границе: 3 3 Граничные условия 240 условие вполне удовлетворительно, хотя условия (3.479), конечно, лучше. Джакуинта и Хын [1968] брали условия (3.479) при исследовании течения неньютоновской жидкости (жидкости Рейнера — Ривлина).

Гринспэн [!969а] на выходной границе потока ставил следующие условия, предложенные Р. Мейером: дф/дх = О, д,'/дх — и (те Я + ди/ду) = О. (3.480а) (3.4806) ') Уравнение движения (2.2) в направлении у в безразмерных переменных имеет вид до до дв дР 1 /д'о д'о т — + и — + о — = — — + — !ч — + —,). д/ дх ду ду ке (, дхз ду' ) В случае стапнонарного течения при о(х, у) = О для всех х, у и при дР/ду=о это уравнение упрощается: до дзо Йе и — = — ~-.

дх дх Из определения ь = до/дх — ди/ду получаем до/дх = й+ ди/ду, д'о/дх' = = дйдх + д(ди/дх)/ду, В силу уравнения неразрывности ди/дх =- †/ду, поэтому д(ди/дх)/ду = — д'и/ду' = О, Таким образом, йе и(С + ди/ду) = = дг/дх. Конечно-разностная форма этого условия такова: Г =2и,йе(~ + — ~ )Лх+й Первое условие означает, что о = 0; с учетом этого равенства и при предположении стационарности течения второе условие означает, что дР/ду = О, т.

е, отвечает обычному условию, принятому в приближении пограничного слоя '). Экстраполяционное условие (3.478б) для зр, может оказаться ошибочным. Такая экстраполяция должна проводиться на каждой итерации при решении уравнения Пуассона, Роуч [1970] исследовал достаточность этого условия для нахождения решения. Пусть условия на входной границе при / = 1 фиксированы. Для одномерной задачи, описываемой обыкновенным дифференциальным уравнением с/зф/с/ха = ь, экстраполяциояное условие с/еф/с/хз = 0 на выходной границе либо противоречит указанному дифференциальному уравнению, если Ь Ф 0 на выходной границе, либо просто совпадает с этим дифференциальным уравнением, если Ь = О. Значит, в одномерном случае такое граничное условие ставить вообще нельзя.

В двумерном случае условие даф/дхз = 0 сводит на границе уравнение Пуассона к обыкновенному дифференциальному уравнению д'ф/ду' = 1. (3 А81) З 3 7 условие на выходной еранице потока Если условия иа границах В 1 и В 3 (рис. 3.22) таковы, что дают досгаточные граничные условия для уравнения (3.481) вдоль В 6, то экстраполяционная процедура может быть достаточной. В частности, если ф = 0 вдоль В 1 и ф или дф/ду фиксировано вдоль В 3 (как было рассмотрено ранее), то экстраполяции вдоль границы В 6 будет достаточно, ио если условие на ВЗ также получаегся с помощью экстраполяции, то экстраполяции вдоль Вб недостаточно.

Постановка же условия даф/ду' = 0 па В 3 либо противоречит (3.481), если ~(В 3) Ф О, либо просто совпадает с (3.481), если ь(В 3) = О. Таким образом, достаточность условия дтф/ду' = 0 на границе, расположенной вниз по потоку, зависит от условий, заданных на смежных границах, причем существенную роль играет размерность задачи. Отметим, что решение задачи методом последовательной верхней релаксации с экстраполяционными условиями как на В 3, так и на В 6 могло бы «сходиться» в пределах некоторой заданной точности; значит, дискретизация могла бы, вероятно, привести к единственному решению, т.

е. к решению, не зависящему от начального приближения. Но полученное таким образом единственное решение зависит от Лх и Лд, и при Лх-э 0 и Лд — 0 задача становится неопределенной. Проведенные выпье рассуждения подсказывают эффективный способ удовлетворения условия батг/бха = 0 на границе, расположенной вниз по потоку (Роуч [!970), Брили (!970)). Вместо того чтобы применять это условие при 1 =! — 1 с помощью линейной экстраполяции, его можно поставить непосредственно при /= /, сводя условие (3.48!) к разиостному аналогу обыкновенного дифференциального уравнения б'ф/бут = ~, (3.482) для которого имеет мес~о краевая задача с условием ф(/, 1) = = 0 в точке на В! и с условием Дирихле нли с условием Неймана в точке на В 3. Величина «(/, /) на Вб может быть найдена при помощи любого из уже рассмотренных способов.

Это одномерное разностиое уравнение можно быстро решить методом прогонки (см. приложение Л), не прибегая к итерационным методам. Если на границе, расположенной вниз по потоку, величина ф определена таким образом, то можно с уверенностью проводить решение разностного аналога уравнения Пуассона с условиями Дирихле на В 6. Скорость сходимости при этом также увеличивается (неопубликоваиная работа автора). Такой способ реализации условия даф/дх' = О, сводящий на границе уравнение Пуассона к обыкновенному диффереппнальному уравнению, необходим, когда для решения уравнения Пуассона л42 8.8.

Грпничные условия применяются неявные схемы метода чередующихся направлений, поскольку расщепление в этих схемах эффективно сводит итерационное решение на полушаге по времени к решению обыкновенного дифференциального уравнения в направлении х методом прогонки, а для такого уравнения условие д'ф/дк' = 0 не является достаточным.

Обыкновенное дифференциальное уравнение (3.482) можно решать более простым одномерным методом, описанным в разд. 3.2.8, вместо более общего метода прогонки (приложение А). Как уже было указано выше, Фромм [1967] и автор настоящей монографии, используя явные схемы для уравнения переноса вихря, опробовали экстраполяцню как для функции так и для к и обнаружили, что такие условия обладают дестабилизирующим свойством. Применяя неявные схемы метода чередующихся направлений, Брили [1970] и Феннинг и Моллер [1971] получили устойчивое решение. Вводя только прп расчете К дополнительную «фпктивную» узловую точку при = /+ 1, после каждого вычисления дь/д/ вплоть до точки 1= / на каждом слое по времени Брили полагал (3. 483) Чтобы удовлетворить условию деф/дх' = 0 в точке 1 = /, решалось обыкновенное дифференциальное уравнение с(пф/с/у' = ~.

Даже при сравнительно мягких условиях (3.478), предложенных Томаном и Шевчиком [!966], можно получить нереально резкое изменение функции ", в окрестности границы В 6 для течений при малых Ке =- 0(!0). Для течений при таких малых Ке па выходной границе Роуч и Мюллер [1970] брали «самые мягкие» граничные условия для ~, которые получаются из уравнения переноса вихря.

Предполагая, что иь; ) 0 (т. е. что В 6 действительно является выходной границей потока), конвектпвный член для и можно аппроксимировать с помощью разностей против потока прн 1 = /, не прибегая к дальнейшим аппроксимациям. Конвсктивный член для о также можно аппроксимировать при помощи разностей против потока (в зависимости от знака оп с) и,чи при помощи какой-либо другой схемы, используемой во внутренних точках; аналогично, для диффузионного члена в направлении у при 1 = / не требуется аппроксимации.

Член, описывающий диффузию в направлении х, мог бы быть вычислен при 1 / — 1. Но само по себе такое вычисление является статически неустойчивым (см. равд. 3.1.4) для члена (д'ь/дх')/Ке, особенно в течениях при малых Ке. В этом легко убедиться, если вернуться к рпс. 3.6; корректирующее смещение, обусловленное членом д'Цдх' для точки 1 = ! — 1, йуйу.

условия на аыходяой границе потока в действительности накладывается в точке г = /, вызывая тем самым монотонный рост ошибки нлп статическую неустойчивость '). Исследование устойчивости при помощи метода дискретных возмущений показывает также, как разрешить эту проблему. устойчивость члена дзс/дхз достигается дополнительным сдвигом этого члена на один слой по времени. В результате получается ~о+1 п + й 1 (пь)г,! (пь)г-ь7 й(оь) (3. 484) Наличие дополнительного временнбго слоя п — 1 в уравнении (3.484) отнюдь не означает, что необходимо хранить в памяти полную матрицу для г",-'. Поскольку это уравнение используется лишь на границе, необходимо помнить только вектор 1'(/) = бзй/бхз 1",, Здесь производная вычисляется в конце каждого расчета нового значения вихря перед получением окончательных значений, а затем осуществляется переход к новому временнбму слою.

Рассчитывая все внутренние точки по схеме с разностями против потока, этот способ с успехом применяли О'Лири и Мюллер [1969), а также Роуч и Мюллер (1970], Независимо от схемы, принятой для расчета внутренних точек, на выходной границе потока рекомендуется использовать разности против потока, хотя бы для представления конвектнвного члена для и. В пределе при Ке — оо это означает, что граничное условие для с на выходной границе ие является необходимым; это аналогично случаю одномерного дифференциального уравнения д~/д( = — д(и~)/дх, где для полной определенности задачи необходимо только условие на входной границе.

Если, например, для конвективных членов выбрана схема «чехарда» (равд. 3.1.6), а для членов, описывающих диффузию, схема с разностями вперед по времени и центральными разностями по пространственным переменным, как в уравнении (3.166), то ') Слсдуя методу днскретных возмущеянй (равд. 33.5 а) н препебрегая конвектнвнымн членамн н членамн, описывающими диффузию в направленнн осп у, получаем й~ ьу~ =ьго 7+(ы/ке)(ь" з г+ й" 7 — зй", у)/аяя. Предположив, что решение стапнонарно, в в Ьс г введено возмущение й, получим (ь)' 7 — ьг 7)/й=й(/(Кеах ), прячем член, стояпшй в пра. вой части, положителен. Следовательно, йь имеет тот же знак, что н й, т.е, вязкий член праводнт к статической неустончнвостн.

244 8.8. Граничные условия при (= с' можно применять следуюШее уравнение: '+ 2Л! ~ бя4л ! ! беел — 1 + — —,, ~ + — ', ~ 1. (3.485) Ле буе !Л! Ле бх' ~П Ь1д' Эти способы были опробованы автором настоящей моногра. фии, и при их помощи удалось добиться плавного измспеппя ~ у выходной границы до тех пор, пока не появлялся срыв вихрей. Когда появляется срыв вихрей, составлявшая скорости и может стать отрицательной при с' = I, поэтому приведенные выше разпостные формулы станут формулами с разностями по потоку, что приводит к неустойчивости.

В этом случае для устойчивости следует ставить условие (3.478) Томана и Шдечика (!966]. Поскольку срыв вихрей происходит только при больших Ке, схемы с разностями против потока в этом случае дадут несущественное улучшение. Устойчивость схемы в целом теперь может определяться ограничением по числу Куранта, соответствуюшему применению разностей против потока, на границе В 6. Отметим, например, что в комбинации со схемой «чехарда», как это имеет место в уравнении (3.486), величина эффективного шага по времени для конвективпого члена и будет 2ЛЛ Критерий устойчивости в одномерном случае будет С, = и ЛИх ( Ъ Однако, как показали Бао и Догерти [!969), разности против потока на границе В 6 можно применять в комбинации с неявной схемой метода чередующихся направлений без каких-либо ограничений па устойчивость.

Другие комбинации должны быть проверены индивидуально. Достаточность рассмотренного условия на выходной границе не была строго доказана даже для линейных дифференциальных уравнений в частных производных, однако некоторые выводы можно сделать из рассмотрения одномерной модельной стационарной задачи — и д + а д, = О. дй деь (3. 486) Применение разностей против потока для конвективного члена не накладывает каких-либо условий ца выходной границе, а из рассмотрения члена, соответствующего диффузии в направлении оси х, в уравнении (3.484) следует, что при Лх- О 6 ~(с), в=О. (3.

487) При и ) О и а ) О этого граничного условия вместе с фиксированным значением величины с на выходной границе достаточно для нахождения решения (см. задачу 3.29). Если же дд7 Уаяааия на выходной ераяиие патака 24й а = 0 (т. е. Ке = аа), то это условие на выходной границе, поставленное для решения исходного дифференциального уравнения, делает задачу переопределенной и создает особенность в решении при сс- 0 (см. равд.

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

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

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

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