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

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

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

Следовательно, источник [и в уравнении (16.99) может быть выражен явно. Расчетной границей при решении (16.99) является ряд точек, лежащих за стенкой. На этой границе для р' ставятся граничные условия Неймана др' дра т — = Рд или — = гд, ду (16. 101) ~д — — 1д+ —,! [ ) Тс(А — ~ д сЬ[ (16.102) ЕА где Рда и Рд определяются выражениями (16.!00). Для выполнения теоремы Грина значение 1« в (16.99) заменяется выражением $16.2. Внутренние течения 335 Для численного решения уравнения (16.99), как и для решения уравнения (16.96), можно использовать прямые и итерационные методы. Весь алгоритм расчета можно представить следующим образом.

Решение уравнений (16.80) — (16.83) получается за один 0.5 0.4 03 В 0.2 0.1 0 0 О , 1.0 2.0 0.5 0.4 0.3 мъ 0.2 ОЛ 0 0 0 гц ЬО 2.0 Рнс. 16.13. Распределение осевой скорости н канале с удлинением УУ/!т' =2.0 ((Вг!!еу, 1914]; печатается с разрешения Асаг!еш!с Ргезз). проход вниз по потоку. В каждой плоскости к"+' значения ие+г, о"+', иге+', р',*+' и р' "+' определяются следующим образом: 1) решение уравнений (16.88), (16.89) и (16.92) позволяет определить и"+' р"+' оя н гпя; 2) решение уравнений (16.96), (16.94) и (16.98) позволяет определить о"+' и ю"+'! 8) решение уравнения (16.99) позволяет определить рь "+'. 366 Гл.

16. Течения, описываемые й!Ч5-уравиеииями Навье — Стокса -0.5 -О. 5 Общая погрешность всего алгоритма 0(Лх, Луе, Ляг). Несмотря на то что некоторые шаги алгоритма требуют проведения итераций, он в целом является весьма аффективным. По сравнению с задачей о закрученном течении (п. 16.2.!) наличие двух поперечных координат 0 2.0 приводит к более сложному )-+--ч-! Ипш)ув """ (игйп)йв способу определения о, ш и по- 0.5 правки давления р' при помощи уравнения Пуассона. В задаче о закрученном потоке радиальная скорость и поперечная поправка давления получались за один мар-' шевый проход в радиальном направлении.

В работе (ВП!еу, 1974) получено решение для ламинарного течения в канале с от;вг ношением сторон 1: 1 и 2: 1. Типичные результаты для канала с отношением сторон 2: 1 представлены на рис. 16.13 и 16.14. Число Рейнольдса йе, вычисленное по средней скорости на оси У и гидравлическому диаметру (рис. 16.7), равно 1333. Решения получены на сетке 21 и, 21 в поперечной плоскости. Для прохода вниз по потоку обычъ/1!у но требовалось 75 шагов. Как видно из рис. 16.13, рассчи- Рис. 16.14. Профили втоРичных сио танные значения осевой соростей в канале с удлинением 2: 1 ((пг!!еу, 19741; печатается с рааре- ставлЯющей скорости в Разшеиия Асадегп!с Ргеаа).

личных поперечных сечениях хорошо согласуются с экспериментальными данными из работы [Ьрагготгг е1. а!., !967). Можно заметить (рис. 16.14), что вторичные компоненты скорости довольно малы. Большие значения поперечных скоростей получаются при неодинаковом нагреве стенок канала и учете плавучести в уравнении вертикальной поперечной составляющей импульса. Этот случай также рассмотрен в работе (ВП!еу, 1974). Описанный выше алгоритм в несколько измененном виде использовался для расчета ламинарного течения в прямом по- $16.2.

Внутренние течения 337 лярном канале [ОЫа е! а1., 1977] и в искривленном прямоугольном канале [ОЫа, Зо)гйеу, 1977]. Обобщение на искривленный полярный канал описано в работе [ОЫа е1 а!., 1979]. Задача о расчете течения в канале за один маршевый проход рассматривалась в работах [Ра1апЕаг, Яра!й!пд, 1972; КпЫп е! а1., 1977; Кгезкочз)гу, ЬЬатго!)г, 1978; Апоегзоп, 1980; Соо1се, Р~чоуег, 1983] и ряде других. В работе [ОЫа е! а1., 198!] показано, что если внутреннее устройство канала приводит к отрыву потока, то укороченные уравнения Навье — Стокса дают правильное решение.

Однако при этом один маршевый проход заменяется повторяющимися маршевыми итерациями, в которых сохраняется и используется на следующей итерации все поле давления. Такой итерационный подход к решению КХБ-уравнений весьма близок к методу, описываемому в п. 16.3.3. Р ! о = 1l + пе + п,р, гп = йг + гве + гве, или в векторной форме ч=т! +че+че, (16. 103) (! 6.

104) 22 к Флетчер, т. е 76.2.3. Течение в искривленном канале прямоугольного сечения Для расчета течения в прямоугольном канале со слабо искривленной осью применим метод, описанный в п. 16.2.2. При малой кривизне поперечные составляющие скорости и и ш малы по сравнению с продольной составляющей и. Это соответствует малым изменениям поправки давления р'' в формуле (16.84), особенно в направлении течения, что позволяет пренебречь членом др" е/дх в уравнении составляющей импульса, направленной по потоку. При большой кривизне оси канала поперечные составляющие скорости и и ге могут стать одного порядка с компонентой и, направленной по потоку.

Вследствие этого возникают сушественные поперечные изменения давления. Расщепление давления, введенное в п. 16.2.2, в этом случае не приводит систему уравнений к неэллиптическому типу и не может обеспечить возможность определения решения за один маршевый проход. Другой способ получения неэллиптических уравнений, пригодный для течений большой кривизны, предложен в работе [ВН1еу, МсРопа10, 1984]. В этом методе для заданной геометрии канала в качестве первого приближения используется не- вязкое решение И, )и, Ю' и р', которое далее модифицируется путем решения К!ч)Я-уравнений, что позволяет учесть вязкие эффекты.

Поперечные компоненты скорости расщепляются следующим об азам: 333 Гн. 16. Течении, описываемые )((ЧБ-уравнениями Навье — Стокса (! 6.105) где каждый член ч состоит из двух компонент (и, ш). Потен- циальная поправка скорости уф= — (пф, вф) создается градиен- том продольной составляющей скорости ди/дх. Ее введение,. как и в (16.95), необходимо для выполнения уравнения нераз- рывности (16.80). Вихревая поправка у = (и„, пгф) создается направленной по потоку компонентой завихрейности ь), кото- рая в декартовых координатах равна двуг дигф д)гз доф ззХ + ду ду ди да Если невязкое решение потенциальное, то дКг/ду — д)г'/дх = О.

Неэллиптические уравнения могут быть получены в связи с тем, что потенциальная поправка скорости уф мала по сравне- нию с невязкой Чг и вихревой у . Это, как будет показано ниже, следует из сравнения порядков различных членов в уравнениях. Ниже описывается подход для расчета несжимаемых вязких течений в ортогональных координатах; обобщение на сжимае- мые вязкие течения в ортогональных координатах описано в работе [ВН!еу, Мс()опа!б, 1984). Ортогональные координаты обозначаются через ($, т), Ь); и, и, ш — локальные составляю- щие скорости вдоль этих координат. Метрические параметры »т, Ьг, »з определяются соотношениями (12.20), т.е, йя г+„г+ г, ),г г+„г+ г йа г+„г+ г(16106) и вычисляются, как описано в $12.2, один раз после построе- ния сетки.

Эквивалентная декартова система координат полу- чается, если положить х1 = у„= хс = 1, а все остальные па- раметры приняты равными нулю. При этом»г =Из —— »з — — !. Поперечное поле скоростей уф и т связано с потенциалом ч) и поперечной функцией тока з(з соотношениями ! дф 1 дф Рф = — —, )гз дч ' гпф = — —, Из дь' д(» ) ! д(» ) пф — —— Ы И,».

дй ф »,», дп Для сравнения порядков величин различных членов предпола- гается, что направленная по потоку координата $ совпадает с линиями тока невязкого течения, поэтому Р = )Р" = О. Под- становка выражений (16.107) в уравнении (16.103) дает и= — — +— 1 дф 1 д (»!гР) 1 дф ! д (»~гР) пт = — — — — ' . (16.108) Иг дп Иг»з дь ' »з дй Иг»г дп При сравнении порядков величин, как и в п. 16.1.1, предпола- гается, что вязкие эффекты ограничены слоем безразмерной толщины б, малой по сравнению с безразмерным продольным $ !6.2. Внутренние течения 339 размером порядка 0(1).

Для обезразмеривания используется характерная длина 1. в осевом направлении. На входе в камал (рис. 16.7) 6 — толщина пограничного слоя; в областях и,о,ой! -О(!) ит,гор,ийр - 016) 0(нг) и, оф,о,р 0(о] игф 0(В~] Рис. !б.!5. Летали поперечной геометрии. вязкого и сформировавшегося течений величина 6 равна половине гидравлического диаметра Ра. Метрические параметры й!, Ья, Ьа в уравнениях (16.108) порядка 0(1).

Вблизи стенки, соответствующей постоянному значению т! (рис. 16.15), и, и, Зе, тй 0(1), о 0(б), ~ — 0 (б). (16.109) Вблизи стенки, соответствующей постоянному значению ь, и, о, д, д 0(1), ш 0(6), д 0(б ]. (16.110) Из сопоставления оценок (16,109), (!6.110) с выражениями (16.107) следует, что !й имеет порядок 0(бт), а тр — порядок 0(6). Тогда в соответствии с оценкой (16.109) -О(6), ш О(6),, О(!).

В соответствии с (16.110) ге~, ше-О(о), ое-0(6'), ое-О(1). Поскольку сильные градиенты возникают у стенок, на которых Ь или т! постоянны, то че 0(6), и, тге- 0(1). (!6.11!) Зап Гл. )б. Течения, описываемые й)ЧЗ-уравнениями Навье — Стокса Следовательно, поперечная потенциальная поправка скорости и мала по сравнению с компонентой скорости и, направленной по потоку. В то же время поперечная вихревая скорость и,з того же порядка, что и скорость в напоавлении потока. Уравнения, эквивалентные (16.80) — (16.83), могут быть записаны в следующем виде: (16.112) (16.113) у п=0, М=(ц у)и+ар — — Г=О, ! Ке где и = (и, и, ш), М вЂ” вектор трех уравнений импульса, Г— сила, обусловленная вязкими напряжениями. После подстановки соотношений (16.104) и (16.107) в (16.112) получается уравнение Пуассона для ф.

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

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

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

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