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

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

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

Исхпнные переменные: стапионарные течения 427 где с( = Е Лу/((1 + Е) а," ), Е = Л!а,". е,/Лх Лу. (17,76) Можно получить аналогичные выражения, связывающие о' е и (бр, — бр,,). Подстановка н".''=н*. +н'. „в (17.67) и использование выражений (17.75) и т. д. позволяет построить следующий явный алгоритм для определения бр! е1 а» бр = ~ а» бр„+ Ь», (17.77) где Ь» = — (и' и — и", ) Лу — (о,'. е — о',) Лх.

Уравнение (17.77) есть замаскированная дискретная форма уравнения Пуассона, которое в символическом виде можно записать как Ч'„бр= — „',Ч, и*, (17.78) что эквивалентно (17.59). Можно отметить, что (17.75) эквивалентно уравнению и — — — !74 бр. бт (17.79) Из сравнения уравнений (17.79) и (17.62) следует, что бр— эффективный потенциал скорости, а поправка ц' безвихревая. Полностью алгоритм 81МР).Е можно представить следующим образом: 1) н* находится из (17.72) и (17.73), 2) бр находится из (17.77), 3) и' определяется из (17.75) и эквивалентной формы для о', 4) р"+' определяется из соотношения р"ч' = р" + се»бр, где а — релаксационный параметр.

В алгоритме 81МР(.Е два релаксационных параметра: а» и Е(— = бт). Решение стационарного уравнения для импульса соответствует значению Е = оп. В этом случае для устойчивости сходимости рекомендуется значение а» =0.075. Эмпирически обнаружено, что более высокая скорость сходимости получается при Е = 1 и а = 0.8 [Ра1ап1саг, 1980]. В работе [Ка)(ЬЬу, ЯсЬпе!бег, 1979] проведено систематическое исследование алгоритмов типа 5!МР(.Е и сделан вывод, что наиболее эффективному алгоритму соответствуют значения Е 4и а» = 1/(1 + Е) (17.80) В работе [Чап Рооппаа!1, Ка!1ЬЬу, 1984] выражение (17.80) названо согласованным $1МР).Е-алгоритмом или сокращенно 81МР).ЕС. Однако здесь дана и другая интерпретация алго- 428 Гл 11. Несжимаемые вязкие течения ритма 51МРЕЕС.

Остается неясным вопрос, увеличивает ли число итераций необходимое для сходимости приближение, имеющееся при переходе от (17.74) к (!7.75). Несомненно, что такое приближение повышает экономичность каждой отдельной итерации. Более точное приближение уравнения (!7.74) получается в результате вычитания ~, а'„'„и'„от обеих частей уравнения и отбрасывания члена ~ аи (и„' — и' ) в правой части. Вместо (17.75) тогда получится уравнение (17.81) где с(!, а = Е ЛГД(! + Е) Пч и — Е ~ ааа). Если поправка и' медленно меняется по пространству, отбрасывание члена ~ ая (и'„в — и' ) приводит к небольшой ошибке. Вместе с тем (17.81), будучи явным выражением, сохраняет экономичность алгоритма.

При выводе уравнения Пуассона для бр в этом случае используется уравнение (17.81), а не (!7.75); точно так же и при расчете и' л. Однако, если в алгоритме 5!МРЕЕ используется уравнение (17.81) при определении р"+' на четвертом шаге алгоритма, нет необходимости вводить релаксационный параметр св„т. е. ал = 1. Аналогичная модифи. кация алгоритма 51МРЕЕ рассмотрена в работе [Соппе!1, 5!очи, 1986). Применение алгоритма 51МРЕЕ в оригинальной формулировке к широкому кругу задач позволило сделать вывод о том, что введение бр эффективно подстраивает поле скоростей, но не позволяет получить быструю сходимость для давления.

Для исправления этого недостатка Патанкар [Ра1ап(саг, 1980) предложил алгоритм 51МРЕЕК, который осуществляется следующим образом: 1. Поле скоростей и определяется из решения уравнений (17.72) и (17.73), в которых члены с давлением исключены из правой части. 2. Уравнение (17.77) становится уравнением Пуассона для определения р"+', а не бр после замены и* на и в членах Ь. 3. Значение р"чл (найденное на шаге 2) подставляется вместо р" в (17.72) и (17.73).

Полученные уравнения решаются методом 5!МРЕЕ, в результате чего получаются значения ц*. 4. Уравнение (17.77) решается для определения бр. В результате определяется иачл = н* + и'. Уточнение значения р"чл, полученного на шаге 2, не производится. 4 17оь Исходные переменные: стационарные течения 429 Очевидно, что в методе ЫМРЕЕК приходится дважды решать уравнение Пуассона и уравнение импульса на каждой итерации.

Хотя число операций на каждой итерации больше, чем в методе 81МРЕЕ, для сходимости достаточно нескольких итераций. Таким образом, алгоритм ЫМРЕЕК оказывается примерно на 50 о(о более эффективным. Можно отметить, что шаги 1 и 2 в 81МРЕЕК соответствуют методу проекций (17.22) и (17.24). В работе 1Чап Роогптаа!1, Ка!!ЬЬу, 1984! проведено сравнение применения методов ЫМРЕЕ, ЫМРЕЕС и ЫМРЕЕК к рас" 400 200 150 300 в 200 100 в !00 50 чету циркуляционного течения и течения за уступом. При решении уравнение (17.77) повторялось т раз до тех пор, пока не выполнялось условие !!г,!~' ~ тдг„!!о, где !!гп~! — среднеквадратичный остаток уравнения (17.77), т.

е. аа брав+ Ьп — аа бр! и !!г (~=[~' ~„г-'~п~, Оптимальное значение у лежит в диапазоне от 0.05 до 0.25. Сравнение вычислительных затрат (время СР(1 в секундах), необходимое для достижения сходимости, приведено на рис. 17.!О, Очевидно, что методы и ЫМРЕЕС, и 81МРЕЕК эффективнее ЫМРЕЕ. Несколько предпочтительней метод 51МРЕЕС. Однако оптимальный выбор Е и в меньшей степени у, зависит от задачи. Алгоритмы типа ЫМРЕЕ на различных сетках применялись также в обобщенных координатах (связанных с телом) (гл.

12). В работе [Ка!!ЬЬу е! а!., 1986) использовался алгоритм ЫМРЕЕС в ортогональных обобщенных координатах. Оказа- О 0 0.5 ! 2 5 1О 20 0.5 1 2 5 10 20 Е Е (а! (Ъ) Рис. 17.10. Сравнение методов 51МРЕЕ, 8!МРЬЕС, 81МРЕЕЙ ((Роогшаа1, йаЦЬЬу, 1984); печатается с разрешения Неш!врьеге РцЫ(аЫпп Со.). 430 Гл. Л.

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

В разностные формулы более высокого порядка входит большее число точек сетки, и алгоритм становится менее эффективным. В работе [5Ьуу е1 а1., 1985] метод $1МР1Е использовался на разнесенной сетке в неортогональных обобщенных координатах. Проведено сравнение использования разностной схемы ЯШСК (п. 17.1.5) и трехточечной схемы второго порядка с разностями против потока (д= 1.5 в (9.53)) для аппроксимации конвективных членов. В качестве тестовой рассматривалась задача о двумерном турбулентном течении в почкообразном канале на сетках 31Х 26 и 56зс', 36. Хотя данная задача и не имеет точного решения, можно сделать вывод о том, что схема второго порядка с разностями против потока в целом более предпочтительна.

Эта схема оказалась более работоспособной и не приводила к очевидной потере точности решения. Схема Я1ЛСК (4=0.375 в (9.53)) расходилась на сильно деформированных сетках, а там, где она сходилась, требовалось большее число итераций. Неработоспособность схемы ЯШСК отмечалась также в работах [Ро1!агб, Ьш, 1982; Ра!е!, Маг!са1оз, 1986], где алгоритм 51МР1Е использовался на декартовой сетке. Из соотношения (9.53) видно, что уменьшение д соответствует приближению к трехточечной центрально-разностной формуле (д = О). Поэтому меньшая по сравнению с трехточечной с разностями против потока схемой (д = 1.5) работоспособность схемы ОШСК (и = 0.375) неудивительна, Филлипс и Шмидт [РЫ!11рз, ЯсЫп!с(1, 1985] использовали алгоритм 51МР1Е в сочетании со схемой ЯШСК для аппрокси. мации конвективных членов на разнесенной сетке.

Для ускорения сходимости к стационарному состоянию использовался многосеточный подход (п. 6.3.5). Филлипс и Шмидт рассматривали задачу о движущейся полости при )се = 400 и задачу о естественной конвекции в вертикальной полости [де ЧаЫ Оач!и, Яопез, 1983] при 1(е=10. Использовалась многосеточная процедура с различным измельчением сетки в различных подобластях. Обычно наиболее мелкие сетки (й = 1/32) вводились вблизи стенок, а менее мелкие (й = 1/16) — во внутренних областях. Самая грубая сетка (й = 1/4) в многосеточной процедуре использовалась во всей области. Гл. 17.

Несжимаемые вязкие течения 432 Ае! = К', 8®й+Вр=В", (17.86а) (17.86Ь) где с) — вектор, содержащий все неизвестные узловые значения обеих компонент скорости, а К и К' — известные из граничных условий Дирихле векторы. Уравнение (17.86) является глобальной нелинейной системой. Решение обычно получается нтерационно методом Ньютона ($ 6.1). На каждой итерации приходится решать разреженную линейную систему уравнений. Обычно это делается методом исключения Гаусса для разреженных систем (9 6.2), в основе которого лежит фронтальный метод [Ноос), 1976[.

Основная проблема такого применения метода конечных элементов связана с выбором аппроксимирующих и весовых функций ф", ф' и фа в (17.82) — (17.85). Можно было бы ожидать, что для и, и и р в (17.82) достаточно интерполяции одинакового порядка.

Однако это может привести к сингулярности, связанной с применением дискретного уравнения неразрывности (!7.83) в слишком большом числе узлов. Это было обнаружено эмпирически Худом и Тейлором [Ноог), Тау!ог, 1974[. Даже если уравнение (17.86) с математической точки зрения ведет себя хорошо, в получающемся решении в поле давления имеются сильные осцилляции; поле скоростей обычно получается гладким. Можно напомнить, что такая же ситуация возникает при центрально-разностной аппроксимации производных от давления в уравнениях импульса, если компоненты скорости и давления определяются в одних и тех же точках сетки (п.

!7.1.1). Это послужило основной причиной введения разнесенных сеток. Тейлор и Худ преодолели проблему осцилляции давления путем введения смешанной интерполяции, биквадратичной для компонент скорости и и о и билинейной для давления р. Это позволило получить гладкое решение и широко использовалось в дальнейшем. Однако такой подход малоэффективен в том смысле, что глобальная точность второго порядка определяется билинейной интерполяцией давления, в то время как экономичность метода определяется биквадратичной интерполяцией ком.

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

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

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

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