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

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

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

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

(3) Если для сходимости ья проводится количественная про. верка, то должна быть учтена возможность неустойчивости, связанной с расчленением решения во времени; такую неустойчивость нельзя обнаружить при проверке сходимости по любому 268 3.4 Критерии сходилсосги и наналвнвсе услоеиа критерию, основанному на величине вида 1(,"ва — ьо~, где р— некоторое четное число.

(4) Скорость сходимости различна как для разных функций, описывающих течение, так и для одной п той же функции в разных точках области течения. Если известна функция с наименьшей скоростью сходимостп итераций, то проверку можно проводить только по этой функции; в противном случае следует проверять все переменные.

(Обычно скорость сходимостп для вихря Ч меньше, чем для составляющих вектора скорости, для которых в свою очередь скорость сходимостн меньше, чем для функции тока ар, из-за операции дифференцирования.) Рискованно также проводить проверку по какой-либо переменной только в одной точке: параметры течения в разных областях задачи могут сходиться с существенно различпымп скоростями.

Например, при использовании стационарных методов с «комби1гированными итерациями> (см. равд. 3.1.22) было обнаружено (Текстор (1968), Текстор с соавторами (1969) ), что скорость сходимости для функции ф"4' меньше, чем для вихря ~а+г1 противоположное явление наблюдается при использовании нестационарных методов. В задаче об отошедшей ударной волне давление в точке торможения сходится гораздо медленнее плотности в той же точке, а это означает, что в качестве переменной для проверки сходимости пе следует брать плотность. (5) Величина в должна зависеть от величины шага сетки; если шаг сетки уменьшается, то соответственно должно уменьшаться е.

(6) В формулировке любого утверждения, относящегося к критерию достижения стационарного значения ь (как долго решение должно устанавливаться илн какой критерия по а должен выполняться), лучше использовать соответствующее физическое время г, а не помер итерации п. Например, критерий (3.501) может выполняться на любой стадии расчета п при любом е просто за счет выбора достаточно малого М. Более осмыслсн критерий гпах ~ с.сг+гг — ~' ~ < е, (3.504а) где Г = 7((уо/ь), для задач с преобладающей копвекцией (Ке» 1) или критерий гпах) "вг — ~,",( < е, где г' = 7(ч!ьа) = Г/Ве, для задач с преобладаюгпей диффузией (Ке « 1); см.

равд. 2.4. Точно так же сообщение, что решение получено при и = 1000, несет мсныне информации, чем указание для этого решения максимальных значений г или г', которые, конечно, зависят от принятой в методе величины И. 84. Кригерни сходнмосги и нональньае «словил 269 Сон и Ханратти [1969[ оценивали стационарное значение коэффициента сопротивления для сферы, графически представляя зависимость Со от 1!4ьааа и экстраполируя иа 1/4ьаа» = О. Никаких подробностей, касающихся экстраполяции онн не приводили, но в действительности проведение такой экстраполяции по лекалу или на глазок», вероятно, так же справедливо, как и любая другая процедура, хотя ни одна из этих процедур неповторима.

(7) Наконец, сходимость численного метода в целом, включая расчет граничных точек, особенно рекомендуется проверять на крайне грубой сетке, но с такой точностью решения, которую допускает длина машинного слова. Использование грубой сетки с несколькими стандартными внутренними точками обычно позволяет получить решение с «машинной» точностью за приемлемое время. Затем тест может быть повторен с существенно отличными начальными условиями и желательно с «противополо>кным» первому случаю направлением счета в итерационном процессе. Если два таких решения согласуются, то можно утверждать (хотя и не строго, но по крайней мере с чистой совестью), что данный метод сходится.

Имеются и некоторые другие способы проверки сходимости, которые целесообразно рассматривать для каждой задачи. Том [1933[, а также Том и Апельт [1961[ предложили критерий сходимостн, основанный на величине «невязки» (см. равд. 3.2.3 и 3.2,4). Вообще говоря, даже для линейных уравнений такой тип проверки может оказаться ненадежным; см. Форсайт [1970]. Браун [1967) отметил, что в задаче тепловой копвекции температуры и скорости переноса тепла (которые представляют наибольший интерес) сходились задолго до того, как сходились скорости течения; убедившись в сходимости по скоростям, в дальнейшем он мог прекращать итерационный процесс раньше, как только устанавливалось поле температур.

Брили [19?О[ проверял сходимость итераций для (а[сне')»+ при решении уравнения Пуассона, определяя значение вихря на стенке ь„ на каждой итерации ф и проверяя выполнение условия Л'~ ( аь Этот чувствительный и рациональный тест для ар обладает также и тем преимуществом, что алгоритмически вычисления имеют такой же вид, как в тесте для расчета неявного граничного условия на стенке Лаь ( ег. Для того чтооы этот второй тест имел смысл, очевидно, должно выполняться условие е1 ( ег, 'Брили [1970) положил е1 = 'бег.

Роуч и Мюллер [1970) наряду с другими авторами обнаружили, что существенный выигрыш машинного времени может быть достигнут за счет применения гр бой величины критерия, скажем е = 10-4, при решении уравнения Пуассона на 270 34 Критерии еходияоети и начальные уставил существенно нестационарном этапе и последующего его уменьшения до 10 — ' при подходе к стационарному решению.

Другой общепринятый прием заключается в том, чтобы итерационный процесс для уравнения Пуассона заканчивать либо тогда, когда шах (три' ь.' — хрх, 1( ( е, либо при А ) А ... полагая, .1 возможно, й,,х = 50, Торранс [1968) применил этот прием и отметил, что при приближении к стационарному решению критерий итерационной сходимости для уравнения Пуассона легко удовлетворялся при й ( й .. Далее он вновь решал уравнение Пуассона только на (и+ р)-м слое по времени, где ! + (йыхх + Ьсх) /4, ПРИЧЕМ )Ехх х1ИСЛО ИтсрацИй нх, трсбО вавшнхся для сходнмости решения уравнения Пуассона при предыдущем его итерировании на слое и !).

Следовательно, поле скорости при решении уравнения переноса вихря остается «замороженным» в течение р шагов по времени, что приводит к существенной экономии машинного времени. Многие из сделанных выше замечаний о итерационной сходпмости приложимы также и к понятию аппроксижаг)ионной гкодимости. После того как на разностной сетке с величиной шага Л! получено решение конечно-разностного уравнения, можно рассчитать другое решение на сетке с Лз = Л1/2, где Л, может быть любой из величин Лх, Лу и Л! (если эти величины равны, то пм равна и Л!). Затем сходимость проверяется по равенству (3.505) Г (Лз) = ~ (Л!) + е и т. п. Конечно, о сходимости можно судить по графику, изображающему решение в зависимости от Л или от !/Л. Далее вычисляется разность между ~(Лз) и Т(Л1), и она служит той величиной, по которой можно судить о сходимости, но это суждение остается субъективным.

Выбор величин в обычно определяется имеющимся машинным временем и объемом памяти машины. Как правило, при решении дифференциальных уравнений в частных производных нерационально уменьшать величину шага сетки даже до величины Лз = Л!/4; в действительности обычно рассчитывается только одно решение; оно и публикуется — хорошее ли оно или плохое.

Было бы прекрасно, если бы мы могли с уверенностью хотя бы качественно представить поведение параметров течения исходя из решения, полученного ') Предположим, например, что й,„= 80, н допустим, что нтераднонный процесс прн решении уравнения Пуассона на слое по времени и = 1000 сходился прн й,„= 18 итераций. Тогда р = 1+ (80+!8)!4 = 18, поэтому уравненне Пуассона для хр можно не решать заново вплоть до слоя и = 1018. 84 Критерии стодииости и накальные условия 271 на грубой сетке; иной раз это возможно, но такая практика опасна. Например, Бургграф [1968[, решая задачу о течении жидкости внутри замкнутой прямоугольной области с подвижной стенкой, обнаружил противоположные тенденции в движении центра вихря в зависимости от числа гсе для решений, полученных на грубой и мелкон сетках.

Стоит отметить, что скорость сходимости решения в случае схем второго порядка точности с ошибкой аппроксимации О (Л') подчиняется квадратичному закону, а в случае схем первого порядка точности с О (Л) — линейному закону, поэтому о реальной сходимости легче судить в случае схем второго порядка точности. Здесь стоит вспомнить о известном способе, называемом экстраполяционным методом Ричардсона ') (Ричардсон [19101; Шортлн и Уэллер [1938[; Сальвадори и Барон [!961)) и служащем для оценки окончательной аппроксимационной сходнмостн разностного решения в схемах второго порядка точности.

Пусть ~ — истинное решение дифференциального уравнения в частных производных с вычислительными граничными условиями, т. е, ~= !ип ~(Л). Тогда ошибка, допущенная в решении л-»е ~(Л~) в некоторой точке (хмуь11), может быть записана в виде ряда Тейлора Е1 = ~ — ~ (Л~) = аЛг + ЬЛ[+ (3.506) Если Лг = (1/р)Ль где р — некоторое целое число, то можно получить разностное решение Р(Лг) в той же точке хь уь 7~ н записать ошибку в виде Ег = » — «(Л ) = иЛг+ 5Лг+ (3,507) Исключая из этих двух разложений а, получаем аппроксима- цию более высокого порядка; » (Лг) ! (гг /д )г + О [Л!Лг, Лгг (3.508а) и при Лг = Л1/2 имеем 1=~/,1(Л,) — ~/,1(Л,)+ О(Л«), (3. 508 б) Неудобство, связанное с этим методом, заключается в том, что практически нет ни путей для нахождения коэффициентов прн Л» и при членах более высокого порядка, ни способов ') Этот прием известен под названием «экстраполяция к пределу», «продолжепима позкод к пределу» и «итерационная экстрапояяиия».

д4. Критерии сходалости и начальные условия установить, является лн сходимость монотонной по Л'). Если сходимость не монотонна, то экстраполяпия Ричардсона может привести к оценке й, худшей чем ь(Л!). Экстраполяционные методы могут привести также к увеличению как ошибок округления (Берджес [1971]), так н ошибок в итерационном процессе. Для того чтобы оценить аппроксимациопную сходнмость решения по шагу сетки, не меняя этого шага, можно провести пересчет задачи по схеме другого порядка точности, Том и Апельт [1961] предложили при Лх = Лу пересчитывать оператор Лапласа (ьвф в уравнении Пуассона и ьгз(/гхе в уравнении переноса вихря) при помощи оператора, построенного на пятиточечном двагональном шаблоне (см.

равд. 3.2.10), который имеет порядок точности О ( Хгг2Л), или при помощи других шаблонов для лапласиана. Пересчет с помогцью схем первого, второго и четвертого порядка точности, рассмотренных в равд. 3.1, предполагает то же самое. Заметим, что соответствующим образом должен быть изменен п порядок точности граничных условий.

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

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

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

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