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

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

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

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

При первом направлении обхода по схеме (3.316а) должно быть известно значение сп,-' с (и+ 1)-го временнбго слоя; при втором направлении обхода по схеме (3.3166) должно быть известно значение ~п,' ', где ! = шах й Это обстоятельство не вызывает осложнений в случае задач теплопроводностп, где температуры или градиенты температуры на границах. как правило, известны для всех моментов времени. Но значения вихря на стенке не известны и, как уже было отмечено при обсуждении неявнтих схем метода чередующихся направлений, это вызывает затруднения. Во-вторых (и это гораздо важнее), если данная схема комбинируется с другими схемами и в ней для конвективных членов используются какие-либо варианты аппроксимации из схемы с разностями против потока, схемы с разностями вперед по времени и центральными разностями по пространственным переменным, схемы «чехарда» и явной схемы метода чередующихся направлений, то полученная комбинированная схема либо оказывается безусловно неустойчивой, либо для нее опять появляются ограничения вида С ( 1 и с) ( '/я '), характерные для явных схем.

Единственной сравнительно успешной комбинацией является комбинация схемы, в которой по обоим чередующимся направлениям обхода точек используются разности против потока для конвективных членов и явной схемы метода чередующихся направлений с осреднением для диффузионных членов (см. Ларкин [1964[): ь'= ь", — С(ь", — ьт' ь) + г[(~",~, — ~", — ь, + ~,,), (3.328а) ь*'=т",— С(ьт' — ьт' т) + с((~,~, — ~,.

— ~", + ~",,), (3.328б) ~"," =Ц~',+~",*). (3.328 в) ') Неопубликованный результат автора. ДД17. Явные схемы метода чередующихся направление 149 Для такой схемы анализ устойчивости при помощи метода фон Неймана дает Алгебраическое исследование уравнений (3.329) весьма затруднительно. Численное исследование их для полного интервала изменения параметров С, с! и О приводит к условию устойчивости С ( 1; ограничения на с( здесь нет. Для решения задач гидродинамнки с ненулевыми конвекгивными членами ни схема (3.328), ни какая-либо другая явная схема метода чередующихся направлений на практике не использовалась.

Единственное исключение составляет рассчитанное Сакураи и Ивасаки !1970) решение задачи о структуре одномерной ударной волны, в которой не ставятся краевые условия. Робертс и Вейс [196б) предложили удачный вариант явной схемы метода чередующихся направлений для решения уравнения переноса для невязкой жидкости, который они назвали «схемой с разностями по диагонали». Эта схема основана на центральных разностях как для производной но времени, так и для производной по пространственной переменной, которые вычисляются в точках, расположенных на полушагах сетки; тн+!д рн, !д й!з.!д, ъ! — !1л и ьн~ ! ьн (3. 330) Члены, стоящие в правой части, вычисляются как средние значения вдоль диагонали ячеек пространственно-временнбй сетки. Если расчет ведется в направлении увеличения ! (см. рис.

3.13), то будем иметь ь'ч!1а = 2 (~т + ~ з-!) пРн ' '!' ~, „, = 2 (~; ! +ь,) приеТ ° (3.331б) Подставляя эти выражения в уравнение (3.330), получаем =ь", — э(~тч! — ь",'!), (3.332а) где 1+ С/2 ' (3. 332б) ! — С вЂ” д+ де в+ Се 1.1- л — ле тв ! — С+ Се тв — д+ де 1+ д — де 2( + (3.329а) (3.329б) (3.329в) )бй 8./. Методы решения уравнения переноса вихря Упритгнение. Показать, что в случае проведения расчета в направлении убывания ! схема с разностями по диагонали Робертса и Вейса имеет вид тьп- ! тьп З' (Ьтп+! ьтп ) (3.333а) где С/2 1 — С/2 ' (3.333 б) Поскольку в рассматриваемой схеме используются центральные разности как по времени, так и по пространственной переменной, схема имеет второй порядок точности; однако из-за Рис. 3.13.

Схема с рааностями по диагонали. Кружки соответствуют известным значениям при обходе в направлении роста !', ромбик — неизвестному значению в точке (!', п+ 1), треугольники — значениям, определяемым иа середине диагоналей ячеек на полушагах. Стрелка справа показывает направление обхода. ~,п+! )т (~ и те „,и,-! — та) 1 — $е те 1 — »е (3.334) (3.335) откуда ) 6)Я = 1 тождественно. Значит, эта явная схема чередующихся направлений безусловно устойчива по фон Нейману и для нее )ст) = 1, как это имеет место и для схемы «чехарда» наличия смешанной производной (по времени и по пространственной переменной) в разложениях в ряды Тейлора формально это трудно доказать.

Как установили Пначек и Уильямс [1970), эта схема обладает достаточной точностью в практических расчетах. Подобно схеме Саульева, рассматриваемая схема сказывается явной во внутренних точках и неявной при расчете граничных условий. Для достижения симметрии при расчете можно чередовать направления обхода точек !т и Ц.

Исследование устойчивости схемы (3.332) при помощи метода фон Неймана дает Зддв, Схема «классики» 151 (см. равд. 3.1.6), но в отличие от последней не обладает неустойчивостью, связанной с расчленением решения по временным шагам. Эта схема двухслойная, но при ее программировании на Фортране достаточно одного массива для хранения 1ч И 1л+! Однако любая ошибка, допущенная при вычислении значений на границе в начале обхода узловых точек, может расти не с увеличением времени, а при изменении пространственной координаты. Пусть все ~л = О, и пусть на граничное значение наложена некоторая ошибка е. Тогда из уравнения (3.332) следует, что "',)~т', = — се.

Для того чтобы избежать увеличения ошибок (скажем, машинных ошибок округления) вдоль пространственной координаты, необходимо, чтобы при возрастании К (17) было /$/:~ 1, а при уменьшении г (Ц) было )й'1«= 1. Для схемы чередующихся направлений оба этп ограничения приводят к условию [С/ = 1 (см. задачу 3.!9). Это один из примеров, когда обычный метод фон Пеймана не дает ответа относительно действительной вычислительной неустойчивости. Подобно неявным схемам, явные схемы метода чередующихся направлений в применении к уравнению конвекции для невязкой жидкости приводят к появлению бесконечной скорости распространения возмущения, что нс является свойством дифференциального уравнения. Рассмотренная явная схема метода чередующихся направлений не комбинировалась с явной схемой метода чередующихся направлений для уравнения диффузии с целью получения безусловно устойчивой явной схемы для полного уравнения, включающего конвективные и диффузионные члены, и не использовалась для решения реальных задач гидродинамики.

Улразгнеиие. Показать, что явная схема метода черелующнхся нвпрввленнй Свульевв (З.З!5) для уравнения днффузнн не приводит к росту ошибок прн переходе от одной прос~рвнственной точки к другой. 3.1.18. Схема «ияасоини» Эта схема была впервые использована для итерационного решения уравнения Пуассона Шслдоном [19621 (см, разд.

3.2.7), а для решения уравнений, описывающих течение сжимаемой жидкости, Скалой и Гордоном [1966, 1967). Гурли [!970а] обобщил эту схему, доказал се устойчивость для общего случая уравнений многомерных течении, применил схему к нестационарному уравнению теплопроводности п к эллиптическим уравнениям и дал ей точное образное название «классики», ад Методы ретаения уравиения переноса вихря 152 Это название связано со способом обхода расчетных точек, показанным на рис. 3.14. На каждом слое по времени п обход точек пространственной сетки совершается дважды.

На первом и последующих слоях по времени с нечетнымн п вычисляют Рис. 3.!4. Схема «классикиж Ромбики соотаетстауат первому обходу, четному числу !+1+ и; кружки — второму обходу, нечетному числу 1+ 1+ а, значения в узлах с нечетной суммой 1+1 (на рисунке они отмечены ромбиками). Для простого уравнения диффузии + дЬ д»Ь д»Ь д! дх' ду' (3.336) получаем и«! и (»!«!.! 2с~,1+ »"-! 1 (',1«! 2»г,1+ »от-! ) =1. + б! '! "! ! Ьх» дуе ! + 1 нечетное.

(3. 337) Этот первый обход основан на схеме с разностями вперед по времени и с центральными разностями по пространственным переменным в узлах с нечетными !'+1. При втором обходе на том же слое по времени выполняются те же вычисления в узлах с четной суммой !+1 (на рисунке онн отмечены кружками), причем в соседних точках берутся значения на (и + !)-м слое, полученные при первом обходе; и«! их! е+! и-,-~ Еч! е«1 ~'."'=~"..+Л1( "' ' ' '.' *-"+ '"' ' " ь1-'), Ь«,1 =Ьсы+ Х, Дхх оуа Е+ / четное. (3.338) Этот расчет при втором обходе полностью неявен в том смысле, что значения в точках (!,1), (! ч-1,1) и (!,/~ 1) бе'рутся на (и + 1)-м слое, но зта «неявность» не влечет за собой Вах!В.

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

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

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

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