Fletcher-2-rus (1185919), страница 85

Файл №1185919 Fletcher-2-rus (Флетчер К. Вычислительные методы в динамике жидкостей) 85 страницаFletcher-2-rus (1185919) страница 852020-08-25СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

Устойчивость алгоритма (18.44) определяется устойчивостью отделвных операторов. Таким образом, аналогичное (18.41) ограничение на Р прн произвольных Лх и Лу имеет вид А/, ~( о ь . (18.45) (ьь/йеа)[2т/Рг+(2/3) ' ( '/аи))+1!а)+а)ах Очевидно, что возможно применять ббльшие шаги по времени. Дополнительное преимущество состоит в том, что для различных координатных направлений могут быть использованы различные шаги по времени. Для обтекания тонких тел, параллельных оси и, необходимы малые в направлении у шаги сетки для разрешения больших пространственных градиентов скорости и тем- 483 й 18.2. Явные схемы пературы, связанных с поверхностным пограничным слоем. Применение схемы расщепления по времени (18.44) позволяет избежать влияния на Л1 ограничения на шаг Л1„.

Схема расщепления по времени эффективна, хотя и требует введения вспомогательных процедур вблизи границ [Реуге1, Тау1ог, 1983]. Схема Мак-Кормака в сочетании с методом установления (п. 6.2.4) малоэффективна для определения стационарного решения, если только не применяется многосеточный подход (п. 6.3.5) [СЫгпа, доЬпзоп, 1985].

18.2.2. Схема Рунге — Кутты В явной схеме Мак-Кормака, описанной в п. 18.2.1, второй порядок точности по пространству достигается весьма экономным образом. Однако ограничение на шаг по времени (18.45) существенно снижает общую эффективность алгоритма при нахождении стационарного решения методом установления для течений с большими числами Рейнольдса. Из условия (18.45) следует, что при больших числах Рейнольдса условие устойчивости схемы расщепления Мак-Кормака эквивалентно тому, что число Куранта С=([и[+а)бг„/сх1 не превышает единицу.

Конечно, влияние больших чисел Рейнольдса непосредственно сказывается и в ограничении на шаг по пространственной переменной для правильного разрешения тонкого пограничного слоя. В результате применения метода прямых (9 7.4) из исходной системы уравнений можно получить систему обыкновенных дифференциальных по времени уравнений.

Для ее решения удобно использовать маршевые по времени схемы Рунге — Кутты, поскольку для этих схем можно использовать большие, чем в схеме Мак-Кормака, значения числа Куранта. Например, для не- стационарных течений можно использовать схему Рунге — Кутты четвертого порядка (7.53). Эта схема при достаточно малых диссипативных членах устойчива, если С (2 1/2. Для нахождения стационарных решений методом установления имеет смысл использовать рациональную схему Рунге— Кутты (РРК) [ЮатЬесй, 1978] первого или второго порядков, поскольку она позволяет использовать еще большие числа Куранта.

Рациональная схема Рунге — Кутты будет продемонстрирована на уравнении общего вида 8д/Ж = йт (д), (18.46) Это уравнение можно рассматривать как одну из компонент системы (18.6) после проведения пространственной дискретизации. Схема РРК использовалась в работе [Ба1о1и(са е1 а1., 31* Гл. 18. Сжимаемые вязкие течения 484 1986). Для дискретизации первых и вторых пространственных производных в уравнениях (18.6) — (18.10) рекомендуется использовать обычные трехточечные центральные разности.

-1.2 Ср -0.8 -0.4 0.0 0.4 0.8 1.8 0.0 1.0 Харда Рис. 18.3. Сравнение с экспериментальным распределением давления в кас- каде 65(12)10 ((за!о1п)та, 1986); печатается с разрешения А!АА). Двухшаговая схема РРК применительно к уравнению (18.46) может быть построена следующим образом. Промежуточные поправки определяются выражениями (ы ! лт(((т(, ) А г А!(()т( и Лс)з =(1 — Ь) М+ (т Лдг, и решение на новом временном слое — по формуле „+ 2Ля' (АЧ',ЛЧ') — Лдз(ЬЧ', ЛЧ') ( 4 ) !7 с) (Ачз А з) В уравнении (18.47) (е,() означает скалярное произведение, т.

е. (е, $)=- ~ еА, $18.3. Неявные схемы где ! пробегает все точки сетки. Скалярные произведения, вычисленные один раз на каждом временнбм шаге, определяют весовые множители, с которыми берутся поправки Лд' и Лде в каждой точке сетки. Поэтому метод явный и экономичный. Схема РРК (18.47) имеет первый порядок точности по времени, если Ьс~ — 0.5. При Ьс = — 0.5 порядок точности равен двум. Схема А(се) устойчива (рис. 7.10), если (18,48) Ьс( 2 сое а (2 — сое а) В работе [За!о!ц(са е! а1., 1986) получено устойчивое решение с числом Куранта до четырех при решении полных уравнений Навье — Стокса с моделью турбулентности Болдуина — Ломакса (п, 18.1.1). Рассматривалась задача о трансзвуковом обтекании лопатки компрессора при М =0.76 и числе Рейнольдса ЗХ10'. Как видно из рис. 18.3, в случае отсоса на стенке имеется хорошее совпадение с экспериментальными данными.

Решение получено в обобщенных координатах на сетке 129 Х ЗЗ (гл. 12). Для решения уравнений Эйлера, где не требуется использовать столь мелкие сетки вблизи стенки, как в случае уравнейий Навье — Стокса, схемы Рунге — Кутты эффективно комбинировались с многосеточным подходом [Затезоп, 1983]. 2 18.3. Неявные схемы Для определения лишь стационарного решения более предпочтительными являются неявные схемы, несмотря на успешное использование явных схем в сочетании с методом Рунге — Кутты интегрирования по времени и многосеточным подходом. Это связано с тем, что в неявных схемах нет ограничений на шаг по времени, следующих из линейного анализа устойчивости (9 4.3). На практике ограничение на шаг по времени существует, но оно значительно более слабое, чем в явных схемах. Это ограничение в неявных схемах может быть связано с нелинейными эффектами, с требованием к точности при рассмотрении нестационарных задач или с медленной сходимостью при использовании метода установления для нахождения стационарного решения.

В данном параграфе будут описаны четыре неявных алгоритма. Сначала будет рассмотрена двухдиагональная схема Мак-Кормака, являющаяся прямым обобщением соответствующей явной схемы (п. 18.2.1). Затем будет рассмотрено применение неявной схемы Бима — Уорминга [Веагп, вагш(пй, 1978) (п. 18.3.2) для решения сжимаемых уравнений Навье — Стокса. Гл. 18. Сжимаемые вязкие течения 18.3.1. Неявная схема Мак-Кормака Данный метод будет продемонстрирован на примере обобщения явной схемы Мак-Кормака (п.

!8.2.1) для решения одномерного уравнения переноса (9.56). Для удобства одномерное уравнение переноса записывается в виде де де дзд — +а — — /х — =О. д1 дх дхз Явная схема Мак-Кормака (18.35), (18.36) применительно к уравнению (18.49) может быть записана в следующем виде. Шаг лредиктор Д/!" ' = — ад/1.„+Н" + /з Д/1.„„//й, (/.е 7Я ! Д/*,е / l Шаг корректор д//"+' = — ад/1.„ //' ' + /х Д/1.х„/7 ', ,7и+ =0.5(/" +д" +д/и+и ') 1 ' 1 / / где Е, и Ь вЂ” односторонние разностные операторы: е/ е/- ~ 1 е е/+~ — е/ кд/= дх ° к/7/= дх (18.50) (18.5 !) а 1.„„— центрально-разиостный оператор второго порядка: д/, — 2е/+ д/,, 1 ""Ч/ дхз Для устойчивости шаг по времени в схеме (!8.50), (18.51) должен выбираться так, чтобы выполнялось условие 2И Дх а+ — — — ((О.

Дх Д1 Данное условие является одномерным аналогом условия (!8.40). Мак-Кормак [МасСогшаск, 1982) предложил неявный алгоритм, эквивалентный (18.50), (18.51). Этот алгоритм может быть представлен следующим образом. Данный алгоритм аналогичен схеме приближенной факторизации, описанной в $ 8.2, 8.3, 9.5 и и. 10.4.2 и 14.2.8. В п. 18.3.3 будет рассмотрено обобщение приближенной факторизации группового метода конечных элементов (п. 17.3.3) на сжимаемые течения. Развитие алгоритма приближенной факторизации путем введения расщепления П) описано в п. 18.3.4. 487 й 18.3. Неявные схемы Шаг предиктор (1+ Л ~' ) До, = Д4с",'+ (Л ~' ) Дд,,'о й*, с йл + д,)«.

с с (18.53) Шаг корректор ~1+Л аС)д л..с,с д,„лсс,е+(Л ~')д, нкс ах с схх с- ' (18 54) ,)л+с 0 5( л+,7* с 1 д,)лес. с) При расчете по формулам (18.53) и (18.54) поправки Дс)" л и Дс)л.+с «получаются по формулам (18.50) и (18.51). При этом с в расчете по формуле (18.51) сс" ,заменяется на с)" ', рассчитанное по формуле (18.53). ная линейная устойчивость (э 4.3). Для этого необходимо, чтобы выполнялось условие Параметр Л выбирается так, чтобы обеспечивалась безуслов- Л ) гпах ) (а + ах — —,1), 01.

(18.55) (18.56) Из сравнения с условием (18.52) следует, что если Дс такое, что явный алгоритм будет неустойчивым, то можно выбрать положительное значение Л, при котором схема (18.50) — (18.54) будет устойчивой. Если Д( таково, что чисто явный алгоритм устойчив, т. е. выполнено условие (18.52), значение Л полагается равным нулю и не требуется проведения неявной стадии алгоритма (18.53), (18.54). Условие (18.55) проверяется в каждой точке. Поэтому для многих задач дополнительные неявные шаги не требуются в частях расчетной области. Это обстоятельство делает весь алгоритм более экономичным.

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

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

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

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