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

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

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

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

Коэффициент при дв1/дкв равен нулю, и в схеме Лейта нет схемной искусственной вязкости'). Таким образом, схема Лейта, фактически представляя уравнение для невязкой жидкости в разностной форме с ошибкой порядка 0(Л!и, Лх'), приводит к той же форме (3.224), что и схема с разностями вперед по времени и центральными разностями по пространственной переменной, дающая в применении к полному уравнению для вязкой жидкости ') Диссипаиия появляется лишь в члене четвертого порядка. 3.! ! '!.

Схема Лейта 121 ошибку порядка 0(31, Ах) при а, = 'Йи'А! (сравните с формулой (3.18) ). Однако при исследовании только стационарных уравнений для этой схемы снова получается я, = '/аи'А1, откуда следует, что стационарное решение зависит от А! н имеет только первый порядок точности. Такая аномалия связана с необходимостью постановки дополнительного условия на выходной границе потока при использовании центральных разностей для производных по пространственным переменным.

На практике высокая точность обеспечивается за счет постановки на выходной границе потока условия градиентного типа (равд. 3.3.7). Эту аномалию можно рассматривать только совместно с граничными условиямп; подробности можно найти в статье Роуча !!97!в]. Отметим один важный момент, на который обычно не обращают внимания. Если в уравнении, включающем конвективный и вязкий члены, для конвективного члена используется схема Лейта, то схсмная искусственная вязкость имеет вполне определенный вид а, = !),иаА1, за исключением единственного случая, когда С =! (см. приложение Б). Анализ устойчивости прн помощи метода фон Неймана можно провести очень просто, используя пример в равд, 3.1.5.

б, относящийся к схеме с разностями вперед по времени и центральными разностями по пространственной переменной, и замечая, что теперь в уравнении (3.!05) надо заменить е( на Сз/2. Тогда в силу формулы (3.108) для множителя перехода в схеме Лейта будем иметь а = 1 — Са (1 — соз 0) — УС з!п 0, (3.231) откуда следует, что 16!а= 1+ Сз(Са — 1) (! — сов 9)а, (3.232) и для устойчивости должны выполняться неравенства — 1(1+ Са(С' — 1)(1 — соз0)з( + 1. (3.233) Легко проверить, что левое неравенство выполняется всегда, а правое, как обычно, только при С (!.

При помощи формулы (3.232) можно также исследовать ошибку, обусловленную затуханием, и фазовую ошибку рассматриваемой схемы. Лейт [1965) провел такой анализ для случая, когда С « 1 и 0 мало, что представляет интерес при метеорологических расчетах. Первое условие имеет место потому, что метеорологические расчеты проводятся при малых А(, 3.1. Методы реаявнал уравнения переноса вохра а второе связано с тем, что наибольший интерес здесь представляют длинноволновые (по сравнению с Лх) возмущения, когда 0 = й,Лх = 2иЛх/А « 1. Полагая С4 « Св и 1 — соз О ж жОа/2, формулу (3.232) можно переписать в следующем виде: ~ 6 1в 1 — СвО'/4; (3.

234) так как (1 — е)'1 — 1 — е/2 при е « 1, отсюда следует, что ! 61= 1 — С 0'/8. (3.235) Поскольку для исходного дифференциального уравнения ~ 6~ = = 1, член Св04/8 представляет собой приближенное выражение для ошибки, обусловленной затуханием. Заметим, что ошибка, обусловленная затуханием, имеет четвертый порядок по Лх. Аналогично можно оценить фазовую ошибку. Как уже обсуждалось выше (см, равд, 3.1.6), точная фазовая скорость для дифференциального уравнения при всех 0 будет и. Решение дифференциального уравнения можно (пренебрегая влиянием границ) записать в виде (3.160) илн в следующем эквивалентном виде: (3.236) ~(х, 1) =~(х — ит, 1 — т), где т — произвольный сдвиг по времени.

Это решение молгно записать через фурье-компоненты с волновым числом й и соответствующей амплитудой У '): св = 1' ехр [И, (х — и/)), (3,237) или (х, 1) = 1'ехр [/(Π— /е,и1)]. (3.238) Значит, точный сдвиг по фазе за время т для решения дифференциального уравнения равен ЛО = — й,ит. Сравним теперь этот результат со сдвигом по фазе для численного решения. Вычислим сдвиг по фазе для дифференциального уравнения (ДУ) в частных производных за время т = Л1: (ЛО) = — /ахи Л1 = — й,С Лх, (3. 239) т. е (ЛО)„= — СО. (3.

240) Фактический сдвиг по фазе для конечно-разностного уравнения (КРУ) находится из геометрических соображений (см., например, рис. 3.8) и определяется равенством "" РЛО)к.Л (3.241) ') Проаифференнировав (3.237), получим дй/д1 = — /а/е ь и дь/дх = ИяЬ, т.е. уравнение дь/д/ = — ад(/дх выполняется точно. 8. Л ! 3.

Схема Лейта Из уравнения (3.23! ) имеем 1гп 6 = — С айп 9. Полагая айп(ЛО) — ЛО и используя равенство (3.235), получаем — Смяв (ЛО)крэ = ! севу~а (3.242) или, поскольку (1 — е)-' 1+ а при а « 1, (ЛО), р, = — С з!п 0(1+ СтОе/8). (3.243) Для удобства сравнения с (3.240) перепишем (3.243) в виде (ЛО)арр —— — СОг, (3.244) где (3.245) Поскольку 9 мало, раскладывая з(п 9 в ряд, получаем  — В'/31 + О (В') (3.246) В или тж! — О'/6 < 1.

(3.247) Сравнивая равенства (3.244) и (3.240), мы видим, что в конечно-разностном решении каждая фурье-компонента переносится вдоль оси х медленнее нз-за наличия множителя т(0) ( 1. В точном решении дифференциального уравнения в частных производных (ДУ) все компоненты переносятся за счет конвекции со скоростью и; в решении же конечно-разностных уравнений остаются все фурье-компоненты точного решения, но различные компоненты переносятся с различными скоростями. Эта ошибка больше для больших О, т. е. для более коротких длин волн Л.

Таким образом, в процессе численного решения различные фурье-компоненты будут отклоняться одна от другой или диспергнровать; это явление часто называется дисперсионной ошибкой. (Одно из первых исследований дисперсионной ошибки было дано в работе Стоуна и Брайена (!963).) Фазовая ошибка за один шаг по времени будет Ее= =(ЛО)крр — (ЛО)дю Из формул (3.240) и (3.243) находим, что Еа — С з!п Π— ( — СО) = — С(0 — Ое/3! + ...

— 01, или Еа СОз/6. (3.248) Таким образом, фазовая ошибка имеет третий порядок по О, а следовательно, третий порядок и по Лх. Предельный случай коротких длин волн можно рассмотреть без дополнительных приближенных допушений, При С 1 из уравнения (3.232) следует, что !6!'( 1, за исключением случая, когда 0 = О. Из равенства (3.23! ) имеем ! п| 6 = — С з!п О ЗЛ Методе~ решения уравнения переноси вихря 124 В пределе при В- и 1т(6)- 0 и из равенства (3.241) получаем, что з(п(ЛВ)- 0 или (М)кри-пО.

Таким образом, фазовая ошибка будет полной, причем трурье-компонента с наименьшей длиной волньв становится полностью стаиионарной. Этот эффект также имеет место в схеме «чехарда со средней точкой» (см. равд. 3.1.6) и типичен для всех схем, использующих центральные разности для члена бь/бх (Фромм [1968]). Прн С = 1 из равенства (3.232) следует, что ) 0[Я = 1. Здесь фазовая ошибка также исчезает, поскольку равенство (3.241) принимает вид [з(п (ЛВ)] 1 (3.

249) или (58)к„ж — В, (3.250) что соответствует равенству (3.240) для точного решения при С = 1. Это легко проверить, полагая С = 1 в конечно-разност- ном уравнении (3.224): ~и+~ ~п (~п ~п ) + (~п + ~п 2~п) (3 251) Упратинение. Определите фввовую и дпсперспопиую ошибки для схемы с ревпостямя против потопе. Другой тип ошибки, которая оказывает влияние на точность нестационарных и стационарных расчетов при использовании как схемы Лейта, так и любых других схем, называется ошибкой, обусловленной неразличимостью. Ошибка, обусловленная неразличимостью, впервые была описана и проанализирована Филлипсом [1959]. В качестве других важных работ на эту откуда Ц'+'=Ц", и что соответствует точному решению.

В работе Фромма [1968] приведены изолинии фазовой и амплитудной ошибок схемы Лейта в зависимости от параметров С и В. Отмечается хорошее повеление этих линий в предельных случаях при С-и1 и В-пО, 2п. В качестве фундаментальной работы по фазовым ошибкам можно рекомендовать статью Кроули [1968а]. Если интересуются только стационарным решением, то фазовые ошибки не представляют интереса.

Однако в нестационарном решении они могут играть важнейшую роль. Основываясь на формуле (3.248), Лейт [1965] нашел, что число й1 шагов по времени, при котором получается фазовая ошибка в один радиан, равно й( = 6/(СВв). Лейт считает, что при численных расчетах фазовые ошибки оказываются важнее амплитудных ошибок. З.Л18. Схема Лепта тему можно указать статьи Граммельтведта [1969[ и Роберта с соавторами [1970]. Здесь мы только опишем это явление. Неразличимость связана с обменом энергией между фурье- компонентами (термин «энергия» здесь используется в общем смысле первого момента переносимой функции сх). Как уже было указано выше, наименьшая длина волн компонент, которые можно различить при данной разностной сетке, равна Л = = 2Лх.

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

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

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

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