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

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

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

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

Итерационные методы, описанные ранее, применимы и в случае рассмотренных разностных формул для лапласиана, хотя скорость сходимостп может при этом измениться. Здесь можно приспособить и некоторые прямые методы (Лебейль [1969], Роуч [1971а, 19716] ), но в этом отношении они обычно оказываются менее гибкими, чем итерационные методы. Наконец, остановимся на вопросе согласованной аппроксимация при дискретизации уравнения Пуассона и при определении скоростей. Решение уравнения Пуассона для з(т используется только для определения скоростей конвекции, входящих в уравнение переноса вихря.

Уравнение Пуассона Чззр = ~ представляет собой пе что иное, как определение вихря н в дискретизи. рованной форме будет записынаться так: би бо бу бх ' (3.431) где и = бтр/бу и и = — бар/бу. Решение уравнения Чаз]з = с следует рассматривать не как решение для хр, а как решение для бар/бх и бтр/бу. Только пятиточечная схема для оператора Лапласа, по которой записано уравнение (3.365), соответствует разностным представлениям для и и и со вторым порядком точности ') и,;=(зйт 1„— зРт ( !)/(2 АУ), (тут+! ( фт-! 1)/(х Ах).

Если лапласиан представляется по схемам более высокого порядка, то может оказаться, что точность решения полной задачи фактически снизится, если и и и не аппрокснмируются согласованно с тем же порядком точности, что и лапласиан. Аналогичные замечания относятся и к методам наименьших квадратов для решения уравнений эллиптического типа (например, Дэвис и Рабинович [1960]). В этих методах величины хР аппроксимируются полиномами высокого порядка или некото- ') Точнее, пятиточечная схема для лапласнаиа соответствует определению скоростей на границах ячейки между узлами (см. равд, 3.!.2), как в уравнениях "т, (н!з=(трт, 1+! трт, 1)/ у (3.433а) ит ! ИЗ = (зрт 1 трт — ~)/ду (зАЗзб) причем значения в узловых точках находятся следующим образом: ит;='/ (и, (, + т т,та).

8.2.!!. Об оценке рассмотренных методов 21! рыми другими т-параметрическими функциональными формами, Функциональные формы обычно выбираются таким образом, чтобы удовлетворить дифференциальным уравнениям во внутренних точках (что не трудно для уравнения Лапласа), а т параметров (нли коэффициентов полиномов) находятся при помощи минимизации ошибок иа границах методом наименьших квадратов '). Как правило, от этих методов нс следует ожидать успеха, особенно когда ф меняется не монотонно, и их не рекомендуется использовать в вычислительной гидродинамике. Но если эти методы применяются, то определение и и и должно быть согласовано с выбранной аппроксимирующей функциональной формой.

3.2.11. Об оценке рассмотренных методов Замечания, сделанные в равд. 3.1.23 относительно оценки методов решения уравнения переноса вихря, применимы также к оценке методов нахождения решения и конечно-разностным представлениям уравнения Пуассона. Следует также учесть замечания, сделанньге в предыдущем разделе относительно согласованности уравнения Пуассона для функции тока н уравнения переноса вихря как в отношении порядка ошибки аппроксимации, так и в отношении вычисления скоростей. Наконец, при оценке быстроты проведения серийных расчетов в зависимости от сложности метода и времени разработки программы необходимо рассматривать оба эти уравнения вместе. Если при решении уравнения Пуассона для функции тока татр = !, используются итсрационныс методы, а в уравнении переноса вихря для дц!д! берется простейшая одношаговая явная схема, то при пестационарном подходе обычно около 90а!с машинного времени затрачивается на решение уравнения хтзтр = ~.

Поэтому, если при представлении д~/д! перейти к двухшаговой явной схеме (например, к схеме Аллена — Чена из равд. 3.1.15), то машинное время прн решении всей системы уравнений для ф и ~ не удвоится, а только увеличится приблизительно на 1Оо!ю Тогда отношение скоростей расчета по схеме с разностями против потока и по схеме Робертса — Вейса четвертого порядка точности (равд. 3.1.!9), равное 45, при решении всей системы уравнений для ф н ~ намного уменьшится (хотя и останется все еще значительным) н станет равным примерно 6. ') Либо функциональные формы выбираются таким образом, что граничные условия точно удовлетворяются, а невязка дифференциального уравнения в частных производных минимизируется. Возможно также, что не удовлетворяются ви граничные условия, ни двффереипиальные уравнения, а минимизируется комбинация ошибок иа границах и невязкн уравнений (личное сообщение А.

Руссо), 2!г З.Д Методы реитения уравнений дяя функции тока Аналогично, число итераций, необходимых для сходимостн решения уравнения й'ф = ~, отчасти зависит от начального приближения для фи р В качестве начального приближения для фн+' при решении уравнения невун+' = ьннн естественно брать ранее полученное решение уравнения уатрн = ~н. При больших б! величина ьй"н' существенно отличается от ь", и тогда ф" будет уже недостаточно хорошим начальным приближением для врн+'.

В пределе при Лг- 0 вихрь аманн-э-ьн, а начальное приближение трн-э-фн+'. (Это наблюдается и при приближении к стационарному состоянию.) Таким образом, выигрыш в машинном времени за счет увеличения допустимых шагов А! при использовании неявных схем для д",/д! по меньшей мере частично теряется из-за увеличения времени, требуемого для каждого шага итерационного решения уравнения йети = Ь, а также из-за дополнительного времени вычислений за счет неявности самой схемы. Иллюстрируя это, Фромм (!964) привел ряд примеров расчетов, в которых в определенных пределах машинное время практически не зависело от А!! При использовании прямых методов решения уравнения 7инй = ~ возникает противоположная ситуация.

Прямой метод расчета вектора распространения ошибки (равд. 3.2.8) приводит к уменьшению машинного времени, необходимого для ре- шениЯ УРавнениЯ йтиф = ы пРиблизительно в 100 Раз; это значит, что на решение этого уравнения теперь потребуется только 10в(в машинного времени, а на расчет д~/д! по одношаговой явной схеме — около 90отв. Если же для расчета дь/д! использовать двухшаговую схему Чена — Аллена, то полное машинное время почти удвоится, а схема Робертса — Вейса четвертого порядка точности окаэкется приблизительно в 40 раз медленнее схемы с разностями против потока.

С другой стороны, увеличение допустимои величины шага Л! (если не учитывать дополнительного усложнения самого уравнения для дьтд!) непосредственно приводит к сокрашению машинного времени, так как время решения уравнения уявр = ~ при помощи прямого метода не зависит от выбора начального приближения для ф. Выбор между итерационными и прямыми методами определяется также видом граничных условий, Если на всех границах ставятся условия типа Неймана, как, например, при решении уравнения для давления (разд. 3.5), то итерационный процесс сходится очень медленно, Если же уравнение для давления желательно решать на каждом шаге по времени, то прямые методы оказываются более эффективными. Рассмотренные здесь методы, конечно, применимы не только для уравнений для функции тока и давления, но и для других гидродинамических задач, описываемых уравнением Пуассона.

Например, в работе Розенбаума [!968) приведен расчет невяз- 8.3 Граиичвые условия 2!3 кого потенциального течения. Итерационные методы для нели. нейных уравнений эллиптического типа рассматривали Эймс [1965, 1969) и Лик [!969]. 3.3. Граничные условия для уравнения переноса вихря и уравнения для фуннции тона Нетрудно представить себе некоторые типы правдоподобных граничных условий для ф и ь, но попытки определить реальные точные условия, приводящие к устойчивым решениям, могут оказаться весьма неудачными.

Было обнаружено, что адекватность любого граничного условия, определенного при помощи численных экспериментов, может зависеть от числа Рейнольдса, разностных схем, используемых во внутренних точках, других граничных условий, а иногда н от начальных условий. Многочисленность этих факторов затрудняет аналитические исследования и ограничивает их применимость. Тем не менее в этом направлении известен целый ряд работ; см., например, Эдди [1949), Кист и Митчелл [1967], Кемпбелл и Кист [1968), П. Дж.

Тейлор [1968, 1969, 1970) и Чен [1968, 1979). Для многих задач отсутствуют математически строгие решения. Наши выводы в основном будут основываться на интуиции, на экспериментах и аэродинамических трубах и на численных экспериментах. Большинство численных экспериментов по исследованию граничных условий осушесгвлялось при помощи простых двухслойных явных схем для уравнений переноса вихря. Заметим, что известно несколько случаев, когда те же граничные условия, взятые в иных схемах, приводят к неустойчивости. (Термин «неустойчивость» используется здесь в смысле отсутствия сходимостп итераций, а не обязательно в смысле экспоненциального роста ошибки.) Эти примеры могут служить предостережением от применения таких существенно частных методов.

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

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

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

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