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

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

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

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

Робертс и Вейс [1966[). Объем вычислений при этом существенно увеличивается. Моленкамп [1968[ отмечает, что при использовании схемы с расположением узлов в шахматном порядке требуется в 45 раз больше машинного времени и в 4 раза больший объем памяти, чем при использовании схемы с разностями против потока. Формальная ошибка аппроксимации Е = 0(Л1е, Лх') не выдерживается глобально во всех точках, если граничные условия тоже нс будут иметь ошибку порядка 0(Лх4), чзо, как правило, не выполняется (см.

равд, 3.3.2). Неустойчивость, связанная с расчленением решения по временным шагам в схеме «чехарда» (равд. 3.1.6), приводит к появлению двух расчлененных решений; данная схема допускает появление четырех расчлененных решений. Для объединения этих решении, очевидно, требуется наличие диффузионных членов и (если считать, что опыт применения схемы «чехарда» может служить некоторым руководством) требуются малые числа 1(е при условии вероятного достижения стационарного решения. Шахматная сетка приводит также к некоторым затруднениям при постановке граничных условий; постановка граничных условий, предложенная Робертсом н Вейсом, приводит к тому, что интерпретация значений в узлах границы прн помощи метода контрольного объема оказывается несогласованной с интерпрстацпей значений во внутренних узлах сетки, а это приводит к снижению точности вблизи границ.

Другая схема четвертого порядка точности с меньшей фазовой ошибкой предложена У. Кроули [1967[, Здесь на первом шаге по схеме Лейта (3.224) с М/2 вычисляются предварительные значения 7«+'м в точках й 1~ 1. На втором шаге используется схема «чехарда»: ~о+! ~а ~ (~а-';Ш ~«+17>) (3.346) Фазовая ошибка данной схемы обсуждается в работах У.

Кроули [1967[ н Фромма [1968[. 158 Ю.д Методы решеная уравнения аереноеа вихря 3.1.20. Схема Фромма с нулевой средней фааовой ошибкой Большинство явных схем обладает только запаздывающей фазовой ошибкой, т. е. рассчитанная скорость конвекции 0-компоненты составляет иг(0), где г(0) ( 1. Фромм (!968) построил комбинацию схем с опережающей и запаздыва1ощей фазовыми ошибками для получения схемы, которая 1) была бы безусловно устойчива и 2) имела бы нулевую среднюю фазовую ошибку. Идея такой схемы заключается в определении точного решения (для линеаризованного уравнения с постоянным и) в новый момент времени при С = 1.

Из этого то(ного решения прн помощи схемы Лейта с разностями назад по времени получается искомое решение при С ( 1. Для простоты продемонстрируем эту идею сначала на схеме с разностями вперед по времени п с центральными разностями по пространственной переменной, а не на схеме Лейта. Для точного решения при числе Куранта, равном единице, имеем б(' = лх/и. Тогда перенос сеточных значений от точки к точке в случае точного решения происходит таким образом: (3.347а) (3.347б) (3.347 в) Искомое решение находится переходом по времени в обратном направлении от 1+Я' до 1+ ог, где Лг < Лг', тогда 41н-лп 4г-наг л4г+лн лн — ле лх лх (~ьн ~е-~ ) (3 348) Разрешая это уравнение относительно ц+л' и замечая, что " (л' '1 1 С (3.349) получаем Члены в правой части определяются точным решением (3.347), и, учитывая обозначение Ц~-" = — Цн,+', имеем , + (С вЂ” 1) (ь", — с",,)).

(3.351) Уравнение (3,351) показывает построение одномерной схемы с опережающей фазовой ошибкой. Вместо цеитральио-разностного представления второго порядка точности для б~'+л'~бх, принятого в уравнении (3.348), Фромм берет схему Лейта а.!.20. Схема Фромма 159 дробных шагов по времени (3.254) для двух пространственных переменных. Схема с опережающей фазовой ошибкой комбинируется (в смысле осреднения) со схемой Лейта (3.254), которая имеет запаздывающую фазовую ошибку.

В результате получается схема Фромма, названная «схемой с нулевой средней фазовой ошибкой». Отметим значения, вычисленные на полушаге и+ '!з и не имеющие физического смысла, знаком гильда и положим С„= иб!!тзх и С„= об!!Лу. Тогда схема Фромма с нулевой средней фазовой ошибкой будет иметь следующий вид: гс,!=!с,!+ 4" (г"-и! гн" !,!+а!-и! '", )+ + 4 хь!-, ! — 2ьс,!+ь!+ь!)+ тСз 2С +~, 4 )(1! е ! — 2~! . т+ь! !). (3.352а) 4 хьс, !-! ьт, /+!+и!, у-2 ~с,!)+ ! Сз — 2С„'з + 4 (ьс,1-< — 2»с, !+»с,!+!) + ~, 4 у (»с, 1-з — 2»с,! ! +»! !).

(3.352б) Эта схема устойчива при С + Се ( 1. Фромм ]1968] построил изолинии модуля ]6] и фазовой ошибки в зависимости от параметров С, С, и 6. Несмотря на то что схема формально имеет второй порядок точности, ее фазовые свойства существенно лучше соответствующих свойств схем четвертого порядка точности Робертса — Вейса [1966] и Кроули ]1967], рассмотренных в предыдущем разделе. Как и для этих схем, затраты машинного времени для схемы Фромма значительно больше затрат для более простых схем. Как н в схеме Лейта и во всех схемах дробных шагов здесь имеется трудность, связанная с постановкой граничных условий на первом полушаге (3.352а).

Эти трудности можно преодолеть, выбирая в качестве значений !. на стенке значения ~ с первого полушага или получая их итерационным путем (см. равд. 3.1.16). Фромм ') рекомендует вблизи границы переходить к болсе простым разностным схемам с центральными разностями или с разностями против потока. Разностные схемы типа (3.352) с учетом диффу.

зионных членов пока еще не появлялись в открытой литературе. Упражнение, Построить схему с нулевой средней фазовой ошибкой иа основе разностей против потока. ') Личное сообщение. в /. Методы ретиения уравнения иереноеа вихря /60 В одной из последующих работ Фромм [1971[ использовал разности против потока и центральные разпости на чередующихся шагах по времени, а не в двухшаговой схеме, уменьшая тем самым время вычислений, по разности против потока брались при числах Куранта С - !/,, С, ( !/м Расчет точек вблизи границ здесь также требует специального рассмотрения (см. Фромм [!97!) ). 3.1.21. Схема Аранавы Схема Аракавы [1966[ часто применяется для решения метеорологических задач, в которых рассматривается уравнение переноса вихря в невязкой жидкости, Это существенно двумерная девятиточечная схема с узлами типа (1+ 1,! — 1) и т.

д. В пей составляющие скорости непосредственно выражены через функцию тока, т. е. принято и = д!р/д//, о = — дф/дх. Не приводя вывода схемы, мы просто выпишем ее не в обозначениях автора, а в наших обозначениях: н-! а/ 6ххлу ллт!л!, / 1!-!,/) (~!, !л! ь!, / !) (ф, !)!, . )(7,»! 7. ! )+л!,. ! (7. ! / ! 7...) — Ф-!,/К!-!,/е! — 1!-!,! !) — ф!,/»!(41~!,!л! — ."; !,/7!)+ +ф/,/ !Ы!+!,! ! — 1! !./ !)+16/л!(!р/~-!,/л! — !р! !,!л!)— — 1!./ !Й!+!,/-! — ф!-!,/ !) — 1! 7!,/Й! л!,/+! — ф!л!,/ !)+ +ь!-ь/(ф!-!./+! "ч!!-!,/-!))" (3 353) Несмотря на очевидную сложность, данная схема обладает некоторымп преиму!цествами.

У нее формальная ошибка аппроксимации составляет Е = 0(Л/Я,/Лх!, тл//4). Это одношаговая схема, и поэтому здесь пе возникает проблем, связанных с граничными условиями. Для этой схемы тождественно выполняется равенство [6[= 1 и тождественно сохраняются величины ~, ~' и кинетическая энергия и'+ о', эти свойства схемы делают ее особенно удобной для решения задач гидродинамической устойчивости. Поскольку схема сохраняет величину ~я, она не подвержена нелинейной неустойчивости Филлипса [1959[, возникающей из-за обусловленных неразличимостью ошибок (такие ошибки имеют место, но остаются ограниченными, так как ~Я остается ограниченным).

Хорошие свойства этой схемы, относящиеся к фазовой ошибке и обобщающие ее на случай метео. рологическпх уравнений в приближении «б-плоскости», рассмотрены Граммельтведтом [!969[. Используя подход Дюфорта — Франкела (равд. 3 1.7), Феста [1970[ включил в данную схему диффузионные члены. 8.! 22. Схемы дан расчета стаеианарных течений К сожалению, эта схема обладает недостатками, присущими другим схемам, используюшим разности по времени типа «чехарда» (см. равд. 3.!.6), которые чувствительны к неустойчивости, связанной с расчленением решения по временным шагам (см. Уильямс [1969[).

Для достижения устойчивого стационарного решения Феста [19?0] время от времени проводил усреднение по временным слоям. 3.1,22, Замечания о схемах дян расчета стационарных течений Многие свойства схем для уравнения переноса вихря, описанные выше, имеют смысл только для нестацпонарпых решений (напрпмер, фазовая ошибка).

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

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

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

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