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

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

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

Параметры 6 и е — положительные константы, 6, е «1. Уравнение (16.18) — стационарное двумерное уравнение переноса Я 9.5). Из сравнения его с (2.1) и (2.2) следует, что уравнение (16.18) эллиптическое. Типичные граничные условия для этого уравнения приведены на рис. 16.3. Во втором случае члены одТ/ду и одр/ду в (16.!0) малы. Однако их следует сохранить, чтобы в невязкой области получить уравнения Эйлера.

Если пренебречь изменением плотности и зависимостью вязкости от температуры, вид уравнения (16.13) для х-компоненты импульса в случае сжимаемого течения совпадает с уравнением (16.5) для несжимаемого течения. Члены в правой части уравнений (16.14) и (16.16) порядка 0((6/С)е), если предполагается, что 1/Ке порядка 0((6/Е)е). Следовательно, как и для уравнения (16.6), правая часть уравнений (16.14) и (16.16) может быть отброшена.

Однако с точностью до 0((6/Е)е) уравнения (16.14) и (16.16) можно привести к виду, эквивалентному (16.6), т. е. до до др Г ! 1 дхо ри — +ро — + — = ~ — ) —, дх ду ду ~ Ке) дух ' 302 Гл. 16. Течения, описываемые й)Ч5-уравнениями Навье — Стокса Предполагается, что решение уравнения (16.18) можно представить комплексным рядом Фурье, т. е. Т = —, ~ ~~~ Та;ехр(1 (о,)ух) ехр(1(а„)ау), (16.19) а- — / — — о Однако, поскольку уравнение (16.18) линейное, для определения качественного поведения решения достаточно рассмотреть лишь 'у=! т(х,1) =0 т(о,у)и =т (у) дт/дх(оо,у) =О у=п Х Х=О д~/ду(х,п)=О Рис. 16.3. Типичные граничные условия для уравнения (16.18).

одну компоненту разложения (16.19): Т = —, ехр (! (а„) х) ехр (1 (а„) у), т (16.20) где Т можно рассматривать как преобразование Фурье Т. Подстановка этого выражения в (!6.18) позволяет получить следующий полинам относительно а и а„: 6а' + уио' + во' + инт = О. (16.21) Уравнение (16.21) иногда называется символом (16.18). Уравнение (16.21) будет использовано для определения зависимости а от произвольного вещественного а„.

Подстановка этой зависимости в (16.20) определит соответствующее поведение решения. Будет ли конкретная мода (16.20) фигурировать в решении, зависит от граничных условий. В результате решения уравнения (16.21) можно получить и и / ео„+ гоо„т ~~ а = — 1 — ~ 1 — )ч1+ 46 26 26 ч и' ) Поскольку 6 « 1, это выражение можно упростить: и и ( ео„ + 1ооа 2 а = — ! — -~-ю' — 1+26 " "). 26 — 26 иа 5 !6.1. Введение Таким образом, еи-„вин в /и ео~! и а =1 —" — —" или — 11х — + — в(+ о —.

(16.22) и и Хй и ( ии' При заданном вещественном значении ои первый корень соответствует решению Т, которое, как следует из (16.20), согласно члену — о„о/и, осциллирует в направлении х и, согласно (во~/и, экспоиенциально убывает по х, если и больше нуля. Если и меньше нуля, имеет место экспоненциальный рост решения. Второй корень соответствует осциллирующему экспоненциально возрастающему при положительном и решению. Таким образом, из рассмотрения двух корней следует, что Т осциллирует по х и экспоненциально увеличивается при любом знаке и. Однако, поскольку уравнение (16.18) эллиптическое, соответствующее граничное условие при х = ео, изображенное на рис.

16.3, не допускает появления в решении экспоненциально нарастающей моды. Приближение, приводящее к укороченным уравнениям Навье — Стокса (п. 16.1.1), эквивалентно отбрасыванию члена бдеТ/дх' в уравнении (16.18). В результате получится параболическое по х уравнение (5 2.3), для которого не требуется граничного условия при х = оо. После подстановки (16.20) в укороченное уравнение вместо (16.21) получится уравнение ви„ и„в (ио + ео„'+ !ой„= 0 или о„=( —" — — ", (16.23) что совпадает с первым корнем в (16.22). Соответствующее решение для Т, согласно (16.20), будет осциллирующим и затухающим по х, если только и и е одного знака.

Если это не выполняется, имеет место осциллирующее нарастающее решение, что исключает возможность использования для его определения одного маршевого прохода в направлении х. Параметр е эквивалентен коэффициентам вязкости или теплопроводности, которые всегда положительны. Таким образом, для устойчивости решения необходимо, чтобы и везде было больше нуля. Поскольку при б = 0 уравнение (16.18) параболическое, положительное значение и соответствует переносу информации в положительном подобном времени направлении. Из приведенного выше примера видно, что анализ Фурье позволяет точно определить, какого типа решение можно ожидать и какие члены ответственны за данный тип решения. В частности, применительно к укороченным уравнениям Навье— Стокса этот анализ позволяет определить возможность появления экспоненциально нарастающего решения, при котором 304 Гл.

!6. Течения, описываемые ЙХБ-уравнениями Навье — Стокса невозможно получить устойчивое решение за один маршевый проход вниз по потоку. Существует очевидная параллель данного метода с применением анализа Фурье для определения характера решения дискретных уравнений (п. 9.2.1). Можно ожидать, что решение дискретных уравнений будет сходиться к решению задачи, описываемой уравнениями в частных производных, если в результате анализа Фурье исходных уравнений получится поведение решения, сравнимое с поведением, определенным на основе анализа Фурье разностных уравнений.

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

В этом состоит суть подхода, используемого для аналитического исследования различных течений (см., например, [81паг1, 1963] или [Огах(п, Рте!о, 1981]). Однако столь исчерпывающий метод не нужен для определения необходимой формы укороченных уравнений Навье — Стокса. Настоящий анализ Фурье можно сравнить с традиционным характеристическим анализом уравнений в частных производных (гл. 2). При традиционном методе характеристик сохраняются лишь высшие производные и уравнения приводятся к характеристической форме, т. е. получается характеристический полипом, например (2.36).

Можно заметить, что, если в проведенном анализе Фурье оставить лишь высшие производные, полипом по о, а„, например (16.21), совпадет с характеристической формой (см. также п. 2.1.5). Таким образом, характеристическая форма уравнения (16.18) может быть получена из (16.21) и имеет вид б,в+„'=6. к в Это уравнение имеет мнимые корни, поэтому (16.18) классифицируется как эллиптическое уравнение в частных производных.

Однако для рассматриваемых задач метод Фурье определения возможности появления точек экспоненциального нарастания решений уравнения предпочтительнее характеристического анализа по следующим причинам: $16.!. Введение и й ехр (1о,х) ехр (Ыру), о-5схр(Ест„х) ехр(1аеу), р р ехр (1олх) екр ((ару), где знак — означает, что подразумевается такая форма решения. После подстановки (16.24) в замороженные уравнения (16.1) — (16.3) можно получить | (а„ 1ар 0 1Л + (а,'+ а„')/Йе 0 ЕЛ + (о'„+ а'„)/)ке 1о„ б = О, (16.25) р где Л вЂ” = ио, + оа„.

20 К. Флетчер, т е 1. Учитываются все члены уравнений, а не только высшие производные. 2. Можно непосредственно определить вклад различных членов системы уравнений в возможно появляющееся экспоненциально растущее решение. 3. Метод более работоспособен в том смысле, что возможен анализ и вырожденных систем (п. 2.1.4). 4. Решение задачи на собственные значения, например (16.21), имеет больший физический смысл, чем решение характеристического полинома, например (2.36). 16.1.3. Качественное поведение решений укороченных уравнений Навье — Стокса Исследование различных приближенных уравнений Навье— Стокса может быть проведено, как это сделано для уравнения (!6.!8), после их локальной линеаризации.

То есть, значения и и о в конвективной части уравнений (16.5) и (16.6) предполагаются замороженными. Таким образом, настоящий анализ не учитывает явлений, связанных с нелинейными взаимодействиями. Так как уравнения (16.5) и (16.6) аналогичны уравнению (16.8) при 5=0, можно ожидать, что решение укороченных уравнений Навье — Стокса подобно решению уравнения (16.18) при 5 =0 будет нарастать (затухать) вдоль линий тока. Имеет .пи это место на самом деле, будет показано ниже.

При применении анализа Фурье к КХ5-уравнениям (и. 16.1.1), подобно тому, как это было сделано в п. !6.1.2, необходимо рассмотреть не одно скалярное уравнение, например (16.8), а систему уравнений. Обобщение анализа Фурье на системы уравнений будет сначала продемонстрировано на уравнениях Навье — Стокса (16.1) — (16.3), описывающих несжимаемое стационарное двумерное течение. Вместо (16.20) предполагается, что 306 Гл. 16.

Течения, оиисываемые ЯХ5-уравнениями Навье — Стокса Чтобы однородная система уравнений типа (16.25) имела решение, необходимо, чтобы де![ ]=О. Для (16.25) это позволяет получить полинам относительно а„: (ав + а'„')(г(ио„ + оо„) + — (а'; + а~~)) = О. (16.26) Вид второго множителя совпадает с (!6.21). Его корни равны о, = (а„'/(и Ке) — а„о/и, — с [и Ке + о„'/(и КеЦ + ст„о/и. Экспоненциальный рост связан с первым корнем, если и меньше нуля, и со вторым, если и положительно.

У первого множителя корни о, = -~-1а„. Мнимый корень со знаком минус после подстановки в (16.24) приводит к энспоненциальному росту по х. Система уравнений (16.1) — (!6.3) эллиптическая. Это можно установить путем введения дополнительных переменных для вторых производных и построением на основе этого эквивалентной системы уравнений в частных производных первого порядка. Для анализа этой системы пригоден метод, описанный в п.2.1.4, приводящий к характеристическому полиному (2.39). Настоящий анализ Фурье дает идентичный полинам, если в уравнении (16.26) отбросить член более низкого порядка !(ио + оа„). Это согласуется с классификацией уравнений в частных производных, основанной на наиболее высоких производных по каждой независимой переменной. Из (16.26) следует существование мнимых корней и при отбрасывании члена 1(иа„+ оа„).

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

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

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

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