Главная » Просмотр файлов » Роуч П. Вычислительная гидродинамика

Роуч П. Вычислительная гидродинамика (1185913), страница 19

Файл №1185913 Роуч П. Вычислительная гидродинамика (Роуч П. Вычислительная гидродинамика.djvu) 19 страницаРоуч П. Вычислительная гидродинамика (1185913) страница 192020-08-25СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

Используя метод дискретных возмущений, исследовать устойчивость схемы с разностями против потока (см. раэд. 3.1,7 — 3.1.9) (л+1 4» 4» 4» — и)о, Ы Ьх (3.92) соответствующей модельному уравнению течения левизной жидкости. Показать, что условие устойчивости накладывает на число Куранта ограничение С = ибиах ( 1 и что оно включает критерий отсутствия осцилляций, обус.

ловленных чрезмерно большим шагом по времени. 3.1.6. б. Анализ устойчивости по фон Нойману Наиболее распространенный метод анализа устойчивости был предложен Дж. фон Нейманом в Лос-Аламосе в 1944 г. В то время с этим методом был частным образом ознакомлен сравнительно узкий круг заинтересованных в нем сотрудников (Эдди [1949]). Краткое описание метода впервые появилось в работе Кранка и Николсона ]!947], а затем в работе Чарни, Фьертофта и фон Неймана ]1950].

Наиболее раннее полное обсуждение метода было дано в работе ОчБрайеиа, Хаймена и ') Это ограничение на гсе, не следует путать с ограничением устойчи. ности в работе Тома и Апсльта (1961, с. 1Зб), В этом случае нтернрование стационарных уравнений выполнялось без «нижвей релаксапиим Это аналогично интегрированию нестацнонарных уравнений с фиксированным д( (см, равд, 3.1.2), и ограничение на Ве, у указанных авторов фактически соог. ветствует ограничению, накладываемому на йб 3.!.5. Исследование устойчивости Каплана [1950).

Как мы покажем ниже, в этом методе решение модельного уравнения представляется рядом Фурье с конечным числом членов и устойчивость (или неустойчивость) определяется тем, что каждое отдельное колебание затухает (нли нарастает). Рассмотрим сначала линейное модельное уравнение с одним только диффузионным членом, снова используя схему (3.18) с разностями вперед по времени и центральными разностями по пространственной переменной: или глч! йл + С( (~л +»Гл 9Г»л) (3.93) где с( = аЛ(/Лхз. Каждая фурье-компонента решения записывается в виде й," = Ь'" ехР1И„(! Лх)1, (3.94) ') При таком подходе возникают некоторые вопросы, которые лучше игнорировать при первом чтении и которые выясняются после ознакомления с соответствующей литературой, Формула (3.94) представляет собой решение конечно-разностного уран. пения (3.93) при нулевых граничных условиях.

Общее решение получается заменой в (3.94) Ул на Аь", где л интерпретируется как показатель степени. Это конечно-разностное решение можно использовать для того, чтобы наглядно продемонстрировать некоторые свойства сходимости конечно-разностиой схемы (3.93) (см. Рихтмайер и Мортон (1967)). Хотя (3.94) не являешься решением уравнения с конвектпвиым и диффузионным членами, мы хотим вскоре обратиться к такому полному уравнению, избегая, однако, обсуждения вопроса о влиянии граничных условий. Это легко сделать (после дополнительной аппроксимации), анализируя устойчивость для бесконечной области согласно фои Нейману. Для рассматриваемой здесь бесконечаой пространствевной области й, принимает все целые значения, й, = 1, 2, ....

Если желательно исследовать влияние граничных условий, то максимальное й, должно быть конечным и должно зависеть от шах ! = !. Некоторан путаница может возникнуть из-эа того, что в литературе используются две различные «нормирующие» системы для длин. В одной систе»1е шах й, = 1 — 1 выбирается как п)Лю При этом требуется, чтобы Х( = шах х сетки) было задано как Х = и.

Если фазовый угол определен как 0 = й»Лх и й еэ (1, à — 1), то это дает ппп 0 = Лх = = и)(! — 1) и шах 0 = (! — 1)Лх = и. Во второй системе шэт х сетки задается более естественной нормированной величиной Х = 1. Тогда Лх = = 1)(! — 1) и фазовый угол определяется как 0 = пй„Лх. Отсюда снова получается соответственно ппп 8 = пЛх = п((1 — 1) и шах 0 = п(!— — 1)Лх = п Понятно» где Ул — амплитуда отдельной компоненты с волновым числом й„ (длина волны Л = йя/а,) на и-м временном слое и 1 = ~/ — 1. Пространственная область считается бесконечной '). 8.П Методы решения уравнения переноса вихря 70 Если ввести фазовый угол 0 = )е„бх, то (3.94) примет вид ~, = 'тт"е (3.95) Аналогично, +1 )т»+~ ыеьив (3.96) ся! Подставляя в уравнение (3.93) выражения (3.95) и (3.96), получаем 17»+' пв )т» пв+ н()т» си+ив+ )т» ти-ов 2)т» пв( (397) или, после деления на общий множитель впе, )г"~' = )т"(1 + с(( ' + е гв — 2)) (3.98) Используем тождество е" +е 'в=2созО (3.99) и определим множитель перехода 6 равенством (3.100) (3.101) 16)(1.

(3. 102) Это условие является критерием устойчивости для уравнения (3.93) с диффузионным членом. Из (3.101) и (3.102) получаем условия — ! (1 — 2с((1 — созО)-='1, (3.103) которые должны выполняться для всех возможных О, т. е. всех возможных фурье-компонент. Правое неравенство выполняется для всех О. Левое неравенство становится критическим при шах(1 — сов О) = 2, что накладывает на с( условие устойчивости д ~~ '/я илн А(» — —" е а (3. 104) Это условие совпадает с критерием (3.73), полученным при по' Мощи метода дискретных возмущений.

Из (3.98) для 6 имеем 6 = 1 — 2сХ (1 — соз О). Заметим, что 6 = 6(0), т. е. в этом случае множители перехода для различных фурье-компонент различны. Равенство (3.100) ясно показывает, что для того, чтобы решение оставалось ограниченным, для всех О должно выполняться условие 3.!.д. Исследование устойчивости Используя тождество (3.99) и тождество е'е — и — 'в=21 з!пО, получаем (3, 10?) 6 = 1 — 2с((1 — соз О) — 1С з1и О. (3.108) В отличие от предыдущего случая уравнение, включающее коивективный и диффузионный члены, приводит к комплексному множителю перехода (3.108). Этот комплексный множитель 6 сводится к действительному множителю 6, определенному равенством (3.100), при С -н О, т. е.

когда уравнение, йеа1 иеа! Рис. 3.8. Голограф множителя перехода О, записанного в виде 13.110). При С < 1, и < '/з и С' < 2й эллипс лезкит внутри единичного круга, что соответствует устойчивости. включающее конвеитивный и диффузионный члены, сводится к уравнению, содержащему только диффузионный член. Условие устойчивости в рассматриваемом случае имеет вид 16 ~<1, (ЗЛ 09) где теперь 16) — модуль комплексного множителя перехода 6, годограф которого построен на рис.

3.8. Уравнение (3.108) мож- Теперь рассмотрим схему с разностями вперед по времени и центральными разностями по пространственной переменной для уравнения (ЗЛ8), включающего конвективиый и диффузионный члены; зто даст С глы ~л 1 ~ (~л ~л )+су(~л 1 ~л 2~л)1 (3 105) Подставляя (3.95) и (3.96) и сокращая на епе, снова получаем (3.100), но с множителем перехода 6, имеющим вид 6 = 1 — — (е'е — е са) + сК(все+ е-та 2) (3 108) С 2 72 а!. Методы решения уравнения переноса вихря но переписать в виде 6= 1 — 2г/+ 2Ысоз0 — /С з!и О, (3.110) что соответствует уравнению эллипса с центром в точке 1 — 24! на действительной оси и с полуосями С н 2с(.

Устойчивость ([6[( 1) имеет место в том случае, когда этот эллипс целиком лежит внутри единичного круга. Для устойчивости, очевидно, необходимо, чтобы С ( 1 и й ( '/,. Более общее условие можно найти, используя для модуля 6 следующее выражение; [6 !т = 66 = — (1+ Ы (соз 0 — 1)]э+ С'(1 — соз'О) (3.111) (здесь через 6 обозначена величина, комплексно-сопряженная 6). Используя элементарные методы определения максимума !6[а в зависимости от соз9, можно убедиться') в том, что при С' ~ ~24! (3. 112) внутри интервала — 1 ( соз О ( 1 максимума [6[ не существует. Этот максимум достигается при сон О = — 1 и просто дает условие г(( Ча, которое было получено для уравнения с одним диффузионным членом.

При [се, ) 2 максимум имеет место в интервале — 1 ( сов 0 ( 1 и всегда 6'г',„> 1. Следова. тельно, двойное неравенство С' ( 2с( ( 1 является необходимым и достаточным условием для устойчивости. Условие (3.112) можно записать в видех) с(/ ~( 2а/из, (3. 113) откуда сразу следует, что при отсутствии вязкости (а = 0) схема с разностями вперед по времени и центральными разностями по пространственной переменной неустойчива при всех М.

Два условия с( ( '/я и тсе, ( 2 являются достаточными для устойчивости в случае линейного уравнения в бесконечной области при постоянном и. Случай, когда и является функцией пространственной переменной, также можно исследовать при помощи данного метода, но это трудно. В случае более общих конечно-разностных схем, использующих не менее трех временных слоев, уравнение, соответствующее ') Этп методы элементарны, но выкладки несколько громоздка. Резуль.

таты для двумерного случая, приведенные Фроммом [!984! н в первом англяйском язааянн настоящей книги, были неверны в том отношении, что в качестве необходнмого условня накладывалось ограннченне на сеточное число Рейнольдса; см по этому поводу равд. 3,3.8. ') Это ограанченне может быть получено также нэ условия, что крнвнзна !дту/дхт! эллипса 6 (8) а точке (1,0) должна быть больше кривизны еднннчной окружности (У.

Д. Сандберг, лнчное сообщение). 73 Дг.б. Исследование устойчивости (3.101), становится матричным уравнением. Для устойчивости при этом требуется, чтобы ~)с ~ ( 1, где Х вЂ” все собственные значения') матрицы сл. Когда Π— просто число, это условие эквивалентно условию (3.102). Пример такого случая будет приведен в равд. 3.1.6.

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

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

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

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