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

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

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

Точное решение этой задачи Ф„= ехр(0.5пх)соз(0.5пр). 14.14. Примените приближенную факторизацию в сочетании с многосеточным подходом на сетие 65Х33 к задаче 14.11. Сравните приблизительно эффективность метода с методом, использующим просто приближенную факторизацию на самой мелкой сетке. Глава 15 Течения в пограничном слое Течения в пограничном слое традиционно выделяют в отдельную категорию течений (см. табл. 11.4 и 2 !!.4).

Применительно к численным расчетам течение в пограничном слое удобно определить как поток, для которого диффузия, связанная с вязкостью, существенна лишь в направлении, нормальном к поверхности, на которой возникает пограничный слой и=не(м) и=о=в Рис. !зл. Течение в пограничном свое. (рис. 15.1). Уравнение для нормальной составляющей импульса может быть заменено условием постоянства давления в этом направлении. Если распределение давления известно, то уравнения, описывающие такое течение, перестают быть эллиптическими, что позволяет разработать весьма эффективные однопроходовые маршевые алгоритмы (в направлении х на рис.

15.1) для их решения. В пограничных слоях имеют место большие градиенты скорости в направлении, нормальном к поверхности. Поэтому желательно использовать преобразования координат, для которых эти градиенты стали бы меньше в новой системе координат. Наиболее эффективные преобразования описаны в Ч !5.2 и 15.3. Кроме того, для получения хорошего разрешения вблизи стенки используются сетки, шаг которых геометрически возрастает в направлении нормали к поверхности (п. 15.1.2). Уравнения, описывающие движение в трехмерном пограничном слое, являются гиперболическими в плоскостях, нормальных к поверхности, на которой развивается пограничный слой. Это усложняет области зависимости и влияния в этих пло- й 15.!.

Простые течения н пограничном слое 239 скостях (п. 2.2.1). Область влияния в явных маршевых алгоритмах определяет возможную величину шага (9 15.4). Для расчета турбулентных пограничных слоев используются те же численные методы, что и для ламинарных. Однако увеличение градиента нормальной составляющей скорости вблизи поверхности может привести к необходимости введения более сильного сгущения сетки в направлении нормали.

Эта проблема может быть решена применением метода Дородницына ($ 15.3), в котором компонента скорости и, направленная по потоку, рассматривается как независимая переменная. Другой путь избежать сильного сгущения сетки вблизи стенки состоит в введении пристенных функций, которые позволяют локально получить аналитический профиль скорости вблизи стенки. Пристенные функции рассматриваются в и. 18.1.1. Пренебрежение диффузией в направлении потока и отбрасывание уравнения нормальной составляюгцей импульса может быть также использовано при расчетах струй, течений в следе и трубах. Подобные течения с тонкими сдвнговыми слоями могут быть весьма эффективно рассчитаны при помощи методов, применяемых для расчетов течений в пограничном слое.

9 15.1. Простые течения в пограничном слое Уравнения, описывающие ламинарный несжимаемый двумерный пограничный слой, могут быть представлены в виде (15. 1) и — +о — =и — +о —, ди ди дие дан дх ду а дх дуя ' (15.2) где известное распределение скорости и,(х) на внешней границе пограничного слоя (рис. 15.1) получается из уравнения Бернулли (11.49).

Поскольку система уравнений (!5.1), (15.2) относится к смешанному параболическо-гиперболическому типу с переменной х, играющей роль времени, необходимо задать начальные условия и(хо у) =ио(у) (15.3) и граничные условия и (х, О) = О, о (х, О) = О, и (х, 6) = и, (х). ( 15.4) Уравнение импульса (!5.2) можно сопоставить с уравнением одномерной диффузии, рассмотренным в гл. 7, и с одномерным уравнением переноса из 9 9.4. Основные отличия заключаются в нелинейности конвективных членов и в связи с уравнением Гл. !З. Течения в пограничном слое неразрывности через нормальную компоненту скорости о.

Поскольку и,г(и,/г(х известно, то этот член действует как источник и мало влияет на выбор численного метода. Все схемы, описанные в ф 7.2 или 9.4, в принципе пригодны для численного решения уравнения (15.2). Явные схемы ($ 7.1) исключаются, поскольку их применение привело бы к неприемлемому ограничению на размер шага Лх по маршевой переменной, следующему из условия устойчивости. Схема Кранка — Николсона (п. 7.2.2) и полностью неявная трехслойная схема (п. 7.2.3) являются безусловно устойчивыми схемами второго порядка точности (по Л1 и Лх) для уравнения диффузии.

Чтобы получить второй порядок точности по Лх при решении (15.2), необходимо аппроксимировать нелинейные члены иди/дх и оди/ду со вторым порядком точности. Для схемы Кранка — Николсона это приводит к необходимости введения итераций для расчета всех точек, расположенных вниз по потоку. Для полностью неявной трехслойной схемы (п.!5.1.1) удается избежать итераций путем использования значений и и в в точках, расположенных вверх по потоку (см. (!5.6)). 15.1.1. Неявная схема Для разработки численного алгоритма конечно-разностные выражения, аппроксимирующие различные члены в (15.1) и (15.2), на равномерной сетке вводятся следующим образом: дц (! Зци+~ — 2и" + О Вч ~) + 0 (Лха), ди (и" + ' — и" +1) — г+~ г ~ -1-0 (Луа) (15.5) ди 2Ьу ""' + 0(Лу').

Назначение индексов в этих выражениях иллюстрируется на рис. 15.1. Верхние и нижние индексы введены таким образом, чтобы подчеркнуть роль координаты х, подобную времени а Для получения линейной относительно иичч системы уравне- ний недифференцируемые компоненты скорости и и о, входя- щие в левую часть (15.2), экстраполируются следующим об- разом: ице~ 2ип ии — 1+ 0(Лха) ои 1 2ои ои — ~+ 0(Лха) (15 6) ! ! ! ! Подставив эти выражения в (15.2) и сделав необходимые преобразования, можно получить трехдиагональную систему $15.1.

Простые течения н пограничном слое 241 уравнений, связывающую значения функций на слое и + 1 по- перек пограничного слоя: а ип+' + Ь ип " + с ип" = с(. ! / — 1 ! / / !+1 !' (15.7) где (2 Оп Пп — 1) т ах ах злу ! ! аух ' 1,5(2ил — ип-') + 2т —, ах / ауе ' ах ах — (2сл — Ол 1) — т— хау l / Аут ' а.= ! Ь.= 1 с = / с(. =(2и". — и" — ') (2и,". — 0.5и",.

') + Ах~и,— „"'~ 1В К Флетчер, т. т Уравнение (15.7) неприменимо при / = 1 (у = 0) или = ЯМАХ(у = учлех); и! лчех = ие при / =ЯМАХ вЂ” 1, следовательно, /1/ в (15.7) следует заменить на /1 — с/и,"+', а затем положить с! равным нулю; и/ =0 при / =2. Наиболее эффективно система (15.7) может быть решена при помощи алгоритма Томаса (п. 6.2.2). После нахождения ип+' значения оп"' могут быть получены / / из уравнения (15.1), которое представляется в следующем дискретном виде: оле1 — оле1 0 5 У ~(! 5ил+1 2ип 1 0 5ил-1) + а / /-' ' ах ' / 1 ' / + (1 бил+' — 2и". + 0.5ип ')1 (15.8) где оп!+1=0. Схема (15.7), (15.8) имеет второй порядок точности по Лх и Лу, является безусловно устойчивой (по Нейману), работоспособной и эффективной.

Однако она должна быть либо дополнена однослойным алгоритмом для начала маршевого расчета вниз по потоку, т. е. при и=1, либо должны быть заданы два слоя (и — 1 и и) начальных данных (15.3). Если для решения уравнения (15.2) использовать схему Кранка — Николсона, значения ип-' не понадобятся.

В этом случае уменьшается требуемый объем памяти и начальные данные должны быть заданы лишь на одном слое. При этом вместо экстраполяции (15.6) используются соотношения ип+' = / =ил+0(Лх) и оп"'=оп+0(Лх). Для получениявторого по/ / / рядка точности по Лх требуется проведение итераций для всех точек, расположенных вниз по потоку. Полученные в результате решения уравнений (15.7) и (15.8) текущие итерации их+' н о/х+1/ используются вместо (15.6), после чего снова решаются уравнения (15.7), (15.8). На первом шаге итераций и" = ип и 242 Гл. !5.

Течения в пограничном слое 15.1.2. 1АМВ).: ламинарный пограничньгй слой Неявная схема, описанная в п. 15.1.1, использовалась при расчете пограничного слоя, возникающего при плоском обтекании двумерного клина однородным потоком (рис. 15.2).

Рис. !5.2. Течение около клина. Эта задача относится к классу течений в пограничных слоях Фолкнера — Скан, обладающих автомодельным профилем скорости (см. (ЬсЫ1сЫ(пд, 1968)). В данном случае компоненты скорости являются функцией одной переменной У[(2 — (!) хч1 (!5.9) ое = о". Итерация продолжается до тех пор, пока с некоторой допустимой точностью ие будут выполняться равенства иле,' = =и,",. После этого значения и"+', о"+' полагаются равными ил-н оае!. На практике более эффективным, несмотря на формальное уменьшение скорости сходимости, оказывается не проведение таких итераций, а уменьшение шага Лх до величины, обеспечивающей необходимую точность. Основной трудностью при использовании однородных по х и у сеток является необходимость введения специальных процедур для учета роста толшииы пограничного слоя.

Кроме того, для правильного описания распределения скорости вблизи стенки должны использоваться весьма мелкие в направлении у сетки, что становится особенно существенным при рассмотрении турбулентных пограничных слоев. 4 !5.1. Простые течения а пограничном слое 2243 (15.10) (! 5.! 2) Точное численное решение для 1(ч) при различных углах раствора клика 6 приведено Розенхедом [тсозепЬеаб, 1964) . Здесь эти табличные данные будут использованы для задания начальных значений и и о и сравнения полученных результатов дальше вниз по потоку с «точным» решением.

Удобно ввести следующее обезразмеривание уравнений (15.1), (15.2): х' = —, у' = —" Ке~~~, и' = —, о' = — Ке~'~, (!5.13) где число Рейнольдса Ке= (1,Ь/ч, а Ь и (,1,— характерные длина и скорость. Для данной задачи (/, соответствует значе- нию и, (15.12) при х = Е. Преимущество соотношений (15.13) состоит в том, что без- размерная координата у и нормальная составляющая скорости оказываются нормированными множителем Веце таким образом, что они становятся величинами одного порядка с х' и и' соот- ветственно.

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

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

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

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