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

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

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

Для стационарных задач итерации (17.123), (17.124) достаточно проводить три-четыре раза на каждом физическом шаге по времени Лй Сходимость к решению уравнения (17.92) достигается прн приближении решения (17.93) к стационарному состоянию. На каждом шаге по времени Л! для решения уравнений (17.!23) н (17.124) можно непосредственно применять 29 К. Флетчер, т, в 460 Гл.

!7. Несжимаемые вязкие течения на втором шаге (17.128) Дополнительные члены )(о) и д(и) возникают при применении теоремы Грина к деф/дхе и дезр/дуе для узла Галеркина, расположенного на границе (п. 8.4.2). Для поверхности РЕ на рис. 17.!4 7(о) =0 и д(и)= и,и/Ау. Если расчетная граница совпадает с поверхностью уступа, то или = О. Однако, чтобы избежать трудностей, связанных с сингулярностью завихренности в точке Е„Флетчер и Сринивас ввели поверхностный слой. В результате или не равно нулю. При применении метода Галеркина с конечными элементами вблизи границы используются два четырехугольных элемента, а не четыре, как во внутренней области. Следовательно, вид массового и разностного операторов по нормали к границе отличен от операторов (17.115), используемых внутри.

Однако массовый и разностный операторы по касательной к границе имеют тот же вид, как и внутри. Например, на поверхности РЕ (рис. 17.14) — = й Ф 'Г (17. 129) алгоритм Томаса (п. 6.2.2) вдоль линии сетки в направлениях / и й соответственно. Уравнения (17.91) устанавливают связь между полем скорости и функцией тока. Применение одномерного метода Галер- кина к уравнениям (17.91) позволяет получить следующую зависимость решения для скорости от решения для функции тока: Мии = Еизр, М„о = — Е„зр.

(17.! 25) Эти уравнения трехдиагональные вдоль направлений линий сетки й и ! соответственно. Следовательно, после того как найдено решение для функции тока, эти уравнения могут быть легко решены, в результате чего определяется поле скоростей (и, о). Постановка граничных условий требует введения дополнительных процедур. На твердой поверхности для определения граничных значений ~ используется псевдонестационарная форма (!7.92): (17. 126) Применение метода Галеркина с конечными элементами и приближенной факторизации позволяет получить двухшаговый алгоритм, аналогичный (17.120) и (17.121).

На первом шаге (у+ в(зЖ) М„Лэ =аЛ7(Ми 9 5„„+ М„Э Ьии) ф" + + аЛ7(Ми! (и) + М„й(и)), (17.127) $17.3. Переменные завикренность — функция тока 45! Для решения уравнений (17.125) требуется определить гранич- ные условия Дирихле для о на АР. Это условие при помощи уравнения (17.91Ь) определяется через значения ф внутри об- ласти где точки с 1 = 1 совпадают с границей АР. На выходной границе (ВС на рис.

17.14) функция тока рассчитывается из уравнений (17.123) и (17.!24) с добавлением — ЛтМковс/бх. Вид операторов М„и Л„в в (17.!23), (17.!24) приведен в (!7.115). Операторы М, и Е„„определяются выра- жениями 5»»= ~„т (1, — 1, О), (!7,132) -=й~Ч После того как значения ф на ВС определены, овс определяется через значения ф внутри области при помощи выражения, аналогичного (17.13!). На границе ВС ставится граничное условие: да~/дха = 0; следовательно, 1.„„Ь = 0 на ВС в (17.120). Кроме того, операторы М, и Л„имеют вид ̻— = ( г, з, 0~, 5»= д ( — 1, 1, О).

(17.133) Соответствующие операторы в направлении у приведены в (17,115), 29» Уравнения (17.127) и (17.128) трехдиагоиальные и совместно с уравнениями (17.120) и (17.121) образуют трехдиагональную систему, решение которой позволяет определить новое значение завихренности ~»+'. Для устойчивости параметр релаксации сс должен быть ограничен условием: са (0.1. Для других зависимых переменных и, и и ф на твердой поверхности ставятся условия Дирихле. На входной границе (АР на рис.

17.!4) значения и соответствуют профилю скорости в пограничном слое вблизи Р. От верхней границы пограничного слоя до точки А значение и = 1, что соответствует скорости набегающего потока. Функция тока на АР при этом определяется из уравнения (17.91а). Уравнения (17.127) и (17.128) при )(о) = пав/Лх и ат(и)=0 используются для определения завихренности иа АР. Операторы М„и в (17.127) и (!7.128) имеют такой же вид, как в уравнениях (17.115). Операторы М, и /.„имеют вид М„=~о, —, — ~, 7.»» — —,(0, — 1, 1). (17.130) Гл. 11. Несжимаемые вязкие течеиия На удаленной границе (АВ на рис.

17.14) и=1 и ~=0. Функция тока получается из уравнений (17.123), (!7.!24) с добавлением члена ЛтМиилв/Лу к правой части (17.!23). Операторы М, и Ьии имеют вид Ми — (О' з' ~ ' Е л ' [О' 1 1)г (17134) 1 зг 1 ии дуя Как показано на рис. 17.14, расчетная граница смещается на тонкий слой от твердой поверхности. При помощи анализа порядков величин [Р1е!с)зег, Ьг!п!ттаз, 1983] значения и, и и тр на границе слоя выражаются через завихренность ~ на границе слоя.

Введение поверхностного слоя позволяет изолировать угловые точки (с! и Е на рис. 17.14) от расчетной области. Это важно, поскольку завихренность ведет себя сингулярным образом на острой кромке. Для локального определения и, и и тр вблизи острой кромки используется аналитическое решение для завихренности [1 цп(, Зс!зчт1с[егзЫ, 1965[. Детали приведены в работе [Р!е!с!зег, Ьг!п!ттаз, 1983). Вся процедура определения решения на новом временнбм слое состоит из следующих шагов. В результате решения уравнений (!7.120) и (!7.12!) находится Ь"+'. После этого для определения зр"+' три-четыре раза решаются уравнения (17.123), (17.!24).

Значения ил.~' и ия.~' определяются из (!7.125). Процедура повторяется до тех пор, пока решение не перестает изменяться. В результате получается стационарное решение. Описанный метод использовался для расчета ламинарного течения за двумерным уступом, изображенным на рис. 17.14. Характерные результаты, касающиеся длины отрывной зоны за уступом, приведены на рис.

17.!6, Очевидно, имеется хорошее совпадение с экспериментальными данными [3!п)за е! а!., 1981; Оо1бз!е(п е! а1., 1970[. Поток отрывается от острой кромки (Е на рис. 1?.14), и за поверхностью ЕР образуется отрывная зона с медленно вращающейся жидкостью. Граница отрывной зоны образуется разветвляющейся линией тока, которая сходит с кромки Е и присоединяется к поверхности 7!С на расстоянии, зависящем в основном от числа Рейнольдса, рассчитанного по высоте ступеньки (рис. 17.!6), Расстояние от обратной поверхности уступа до точки присоединения х, слабо зависит от толщины пограничного слоя набегающего потока вблизи точки Е на рис.

17.14. В этом причина различия двух наборов экспериментальных данных. Наиболее мелкая сетка используется вблизи точки Е на рис. 17.14. Размер шагов возрастает в положительном и отрицательном направлениях х до тех пор, пока шаг Лх не станет равным Лх,„; после этого г„= 1. В направлении у Лу = Лувия $17.3. Переменные аавихренность — функция тока 453 от 0С до Р'С' (граница набегающего пограничного слоя). Между Р'С' и АВ Лу нарастает с геометрическим отношением г„. Для более мелких сеток г < 1.04, гн < 1.07, Для более грубых сеток г„г„< 1.15.

В зависимости от числа Рейнольдса сходимость к стационарному решению требуется от 500 до1000шагов 24Я О случай А ) сема 111 Х67 Ф слууак В 20.0 слу ай С а слу»аа О сеуке 5 9 Х 5 9 14.00 йч 12.00 15 лла сс ак ева!1 ):.=:: л Юсыас л ее ам енуО! Ваап О.ОО 50.00 100.00 150.00 200.00 250.00 300.00 не„ Рис. 17.16. Положение точки присоединения прн обтекании уступа, по времени. Состояние считается стационарным, если остаток, т. е. КНЯ в (1?.117), становится меньше 10-'.

В результате применения массовых операторов в уравнениях, подобных (17.114), точность метода конечных элементов выше точности эквивалентных конечно-разностных методов. Эквивалентные конечно-разностные методы могут быть получены в результате сборки массовых операторов. Например, М„ в (17.115) заменяется оператором М„~о, —, 0~. Таким образом, путем сквозного деления на 0.25(1+гл) (!+г„) можно избежать явного появления массовых операторов в конечно-разностной форме уравнений, подобных (17.93). Влияние сборки массовых операторов на точность, устойчивость и вычислительную эффективность алгоритма исследовалось в работе [Р1е1с)1ег, ЬНп)паз, 1984] . На грубой сетке Гл.

17. Несжимаемые вязкие течения (29;к', 18) для получения устойчивого решения массовые операторы необходимо сохранять вблизи расчетных границ. Сборка во внутренних точках не влияет на устойчивость и не приводит к значительному снижению точности. В двумерном случае массовая сборка во внутренней области дает небольшой (18 $) выигрыш в экономичности. Если соответствующая формулировка в переменных скорость — завихренность (п.

17.4.2) используется в трехмерном случае, из оценки числа операторов Рис. 17Д7, Область расчета при обтекании полости, обращенной вниз по потоку. следует, что внутренняя массовая сборка позволяет получить экономию порядка 40 с)о. Рассмотренный метод может быть использован для расчета течения около обращенной по потоку полости (рис. 17.17). Полость образуется путем добавления кромки ЕО к уступу, Наличие полости сдвигает первичную область вращающейся жидкости вниз по потоку, а внутри полости могут возникнуть вторичные циркуляционные зоны. Для изменения картины циркуляционного течения по нормали к 0Е вводится вдув и отсос жидкости. Уравнения и граничные условия для этой задачи имеют тот же вид, что и в случае обтекания уступа. Функция тока тр полагается равной нулю на ЕСтЕ.

Распределение функции тока фпл получается в результате интегрирования известной скорости вдува и отсоса жидкости мол. На отрезке 0С значение функции тока постоянно, т(зпс = зрп. $17.3. Переменные завихренность — функция тока 455 Завихрениость в точке 6 (рис. 17.17) неоднозначна. Для преодоления этой трудности вводится поверхностный слой, локальное аналитическое решение внутри которого сшивается с численным решением вне слоя.

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

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

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

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