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

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

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

Если 1сД1/Дхе ограничено при стремлении Дс и Дх к нулю, неявная схема Мак-Кармана, как и явная (18.50), (18.51), имеет второй порядок точности по времени и пространству. Сохранение второго порядка точности по времени [МасСогтаск, 1982) обеспечивается третьим порядком точности по времени дополнительных поправок в (18.53) и (18.54). Двухдиагональность алгоритма (18.53), (18.54) означает, что «неявные» поправки Дс)' ' и Ддл+'с могут быть определены явным образом. На шаге предиктор (18.53) представляется в виде ьч'; е+(Л ас(ах) ьч',,с, дссс С + Л аСслх Гл.

18. Сжимаемые вязкие течения 488 Шаг корректор Лй",",' '= — Лг(/.„-Г'ь, + /.„С,'.,), (! — ЛГУ.„А )(! — Л1/.„В )Ы~",+,~' =ЬП +ь~", (18.59) <~л+! = 0.5 (<!л ! Я ! Л<~л+~, ь) ьм Односторонние разностные операторы, применяемые к Г и Сз, совпадают с описанными в 9 18.2. Для смешанных производных используются центральные разности. Модифицированные матрицы Якоби Ам и Вм связаны с дей- ствительными матрицами Якоби для невязкого течения, но их собственные числа всегда положительны. Матрицы Якоби не- вязкого течения могут быть модифицированы, как в (14.107), т.

е, А = Тл 'ЛлТл, В = Тв'ЛвТв, (!8.60) При этом предполагается, что из граничного условия Дирихле на правой границе /=УХ можно определить 7вх. Вычисления по формуле (18.56) проводятся в сторону уменьшающихся значений / до тех пор, пока не будет достигнута левая граница. Полный неявный алгоритм численного решения уравнения (18.49) состоит из следующих шагов. Сначала проводятся вычисления по формулам (18.50), затем — по формулам (18.53). После этого следует шаг корректор, состоящий из вычислений по формуле (18.51) и затем — по (18.54).

При обобщении метода на сжимаемые уравнения Навье— Стокса (18,6) неявный алгоритм следует представить в виде ~!+ ЛГ ( — А+ —,В)1ла"" =Ла""'е (18.57) где якобианы А=дГ/дй и В=дб/дй и Ьй"+„ь'= — Л/1дГ/дх+ + дкз/ду]; м Пространственные производные в левой части (18.57) действуют на произведения АЛе! и ВЛе). Из сопоставления с (14.!03) следует, что у = О, 6 = 1, и в (18.57) требуется ввести еще пространственную дискретизацию. Полный алгоритм, эквивалентный (18.50) — (18.54), может быть представлен в виде Шаг предиктор (! — Л(/.+А ) (! — /зт/„+В ) Лг(,",' = Лг(,",, (18 58) е)ь я Чг, я+Лт(г,'я.

4 18.3. Неявные схемы 489 где Лл и Лв — диагональные матрицы, составленные из собственных чисел матриц А и В, т. е. с((анЛл — — (и, и+а, и, и — а), с(1апЛв=(о, о, о+а, о — а). (18. 61) Матрицы Тл и Тв приведены в работе [МасСогшасК 1982]. Модифицированные матрицы Якоби имеют вид А =Тл0Тл, В =Тв0Тв. (18.62) Здесь 0' и 0в — диагональные матрицы, диагональные элементы которых имеют вид Р~~с =щах~~ З,ли~1+ — — ', О~, (18.63) Рс,1=гпах((З,вь с(+ — „— ' ~, О~, (18.64) где т = шах(414/3, 714/Рг). Значения 14 и Рг в зависимости от рассматриваемого случая можно считать ламинарными или турбулентными.

Правые части уравнений (18.63) и (18.64) могут быть лишь неотрицательными. Данные уравнения являются обобщением условия (18.55) для двумерных систем уравнений. То есть, если соответствующие диагональные элементы Рь1 и Рь 1 больше нуля, в (18.58) и (18.59) включаются неявные в шаги. Если Рь ~ и (или) Рь г равны нулю, явная схема будет л в устойчива и в неявных шагах нет необходимости. Неявные шаги осуществляются в два этапа. На шаге предиктор сначала в результате решения уравнения определяется промежуточная неявная поправка Л41" '. Начиная с линии сетки Л = Л1У и правой границы 1= )хл, уравнение (18.65) используют для расчетов при уменьшающихся значениях 1' до тех пор, пока не будет достигнута левая граница.

Далее процесс повторяется в сторону уменьшающихся л, пока не будет достигнута нижняя граница. Из сопоставления уравнений (18.65) и (18.58) следует, что (1 — Ла.„'в ) Лп",", = Лй,','. (18.66) Следовательно, Лй'*„' определяется уравнением (!+ — В ~Лй'"'е =Лс)"'е+ — „В,Л41'*'л',,(=%). (18.67) Гл. 18. Сжимаемые вязкие течения Уравнение (18.67) используется начиная с линии сетки 1= Л~Х и верхней границы А=Фу в сторону уменьшающихся значений й. Далее процесс повторяется при уменьшающихся значениях 1 до тех пор, пока не будет достигнута левая граница.

Матрицы А~ и Вж в уравнениях (18.65), (18.66) имеют блочную структуру. Поэтому на каждом шаге требуется решать систему уравнений размерности 4зт',4. Однако этого можно избежать, если использовать формулы (18.62). Уравнение (18.67) можно представить в виде (!+ — (Т ) 0 Т ) Лп*'е =%. (18.68) Здесь % означает правую часть уравнения (18.67). Умножая обе части уравнения (18.68) на Тв и выполняя несложные преобразования, получаем Т Ле)*"е~=(!+ — !з ] Т %=У, (18.69) Ье)""! =(Та) 1У. (18.70) Следует отметить, что матрица (Тв) — ' может быть определена аналитически.

Мак-Кормак [МасСогшаск, 1982] отметил, что при программной реализации вычислений по формулам (18.68)— (18.70), заменяющим (18.67), требуется около 28 фортрановских операторов. Аналогичный алгоритм может быть построен для решения (18.65). Описанный алгоритм при применении его к сжимаемым уравнениям Навье — Стокса имеет, как и в случае модельного уравнения (18.49), второй порядок точности по пространственной переменной. Мак-Кормак [МасСогшас(е, 1982] использовал описанный метод для расчета отрывного течения, обусловленного взаимодействием скачка с пограничным слоем. При образовании модифицированных матриц Якоби Ам и В~ вязкие члены учитываются лишь приближенно.

Поэтому метод работает лучше при определении стационарных, чем нестациоиарных решений. При расчете трансзвуковых течений около агродинамических профилей под углом атаки [Когг(ц!!а, МасСогшас)т, 1982] обнаружено, что если использовать большие шаги по времени, необходимо в явную и неявную части (18.58) и (18.59) включать численную диссипацию. Из сравнения известных значений времени счета следует, что неявный алгоритм (18.58), (18.59) требует для расчета одной точки примерно столько же времени, сколько и блочно-трехдиагональный алгоритм, который будет рассмотрен в п. 18.3.2. Однако преимуществом схемы Мак-Кормака является то, что для 491 $18.3. Неявные схемы многих точек сетки явный алгоритм, лежащий в основе схемы, является устойчивым.

Поэтому для многих задач общая эффективность неявного алгоритма Мак-Кормака оказывается более высокой. Ханг и Кордулла (Нппд, КогдиПа, 1984] использовали неявную схему Мак-Кормака для решения уравнений тонкого слоя 7 8 сь 1 -4 -3 -2 -1 О л7'и Рис. 18.4. Давление на плоскости вдоль линия симметрии ((Нппи, Когдп!!а, 1984]; печатается с разрешения А!АА). (18.31).

Они применили пространственную дискретизацию по методу конечных объемов и одномерное расщепление, эквивалентное описанному в п. 18.2.1. Рассчитывалось сверхзвуковое течение около затупленного вертикального стабилизатора, расположенного на пластине при М = 2.95. Число Рейнольдса, рассчитанное по скорости набегающего потока и диаметру стабилизатора, равно 0.8рс, 10а. Для учета турбулентных эффектов использовалась алгебраическая модель турбулентной вязкости Болдуина — Ломакса (18.17) — (18.21). Уравнения решались в обобщенных координатах на неравномерной С-сетке с 40, 32 и Гл.

18. Сжимаемые вязкие течения 492 32 точками в окружном, радиальном и вертикальном направлениях соответственно. Рассчитанное распределение давления (рис. 18.4) очень хорошо совпадает с экспериментальными данными [Ро!1!пд, Воддопо(1, 1982]. Первое увеличение давления связано с вихрем, ось которого расположена в точке х/О ж — 0.75. Второе увеличение давления связано с прохождением искривленной ударной волны и торможением в точке торможения х/О = О. Неявная схема Мак-Кормака [МасСоггпасК 1982] весьма эффективна, если на границе определены условия Дирихле, но двухдиагональный метод решения непригоден при других типах граничных условий.

Поэтому Мак-Кормак [МасСогшасК 1985] предложил заменить двухдиагональный метод итерационным методом прямых Гаусса — Зейделя ($ 6.3) или методом Ньютона ($ 6.1). При этом, чтобы сделать систему (18.57) с диагональным преобладанием, используется расщепление потока [5!едег, Юаггп!пд, 1981]. 18.3.2. Схема Бима — Уорминга Данная схема предшествовала схеме приближенной факторизации (14.104), (14.105), использованной в п.

14.2.8 для решения уравнений Эйлера. Схемы Бима — Уорминга [Веат, вагш!пд, 1978] и Брили — МакДональда [ВН!еу, МсРопа16, 1977] являются схемами приближенной факторизации и тесно связаны со схемами, описанными в 9 8.2 и п. 10.4.2. Для применения схемы Бима — Уорминга к сжимаемым уравнениям Навье — Стокса уравнения (18.6) записываются в виде (18. 71) де)/д1 = ГхНЗ. В КНЯ входят все пространственные производные, т. е. ~НБ= — —,„— —,„= — —,„~à — Г,(9, 9,) — Г,(й, й„Д— дн дб д ~ е е — — [б' — б;(а, й„) — б,"(Ч, йи)1, (18.72) В (18.72) невязкие потоки Гг и б' совпадают с Г и б, определяемыми уравнениями (14.95).

Другие слагаемые Г1, Гз, б~ и бт' связаны с вязкими членами в (18.7). Эти члены получаются в результате подстановки выражений (18.9) и (18.10) в (18.7) и группировки их таким образом, что Г", и б", содержат производные по х, а Г," и бе — по у. 4 18.3. Неявные схемы 493 В результате применения трехслойной схемы (п. 8.2.3) к уравнению (18.71) можно получить (1+ а) Лал — ада" = Ж (й РСНВл + (! — 8) КН8"), (!8 73) где дал-~-! — ал-~-1 ал н дал ал ал — ! Параметры а и 8 выбираются так, чтобы обеспечить требуемую точность и устойчивость; и эквивалентно у в (8.26); КН8л+' — НЕЛИНЕйНая фуНКцИя а, а, И а„.

В рЕЗуЛЬтатЕ ЛИНЕаризации, аналогичной (14.101), (14.102), относительно п-го временнбго слоя имеем мН8"" = мН8" + ( — ',"' ) Ла"" + ( — ', ' ) Ла. + + ( ) Дан+ ... (18.74) или — — (В Ла — чда — Лб', "' — 8 да„"' ) . (18.75) В (!8.75) А и  — невязкие якобианы (14.99) и дн", дк~~ доле до~я Р= — ', К= — ~, Я= — ~, 8= — ~. дя ' дях' дя ' дяя Уравнение (18.75) можно упростить, если испольэовать равенство 1(Д л+1 (1(л ле~) (х Дал.н и аналогично для 8да„~~. Члены ЛРе' ~~ и Лбр "~~ приводят к появлению смешанных производных при подстановке (18.75) в (18.73).

Этими членами легче оперировать, если заметить, что ЛРе'"+'=ЛРе'л+ 0(ДГ ), Дбб" =Дбб" + 0(лг ). (18.76) Таким образом, с точностью до второго порядка неявные поправки Ре и б1 могут быть заменены соответствующими известными поправками с предыдущего временибго шага. В работе [Веат, багги(пд, 1978] отмечается, что применение (18.76) к модельному уравнению (18.49) не нарушает безусловной линейной устойчивости рассматриваемого алгоритма.

494 Гл. !8. Сжимаемые вязкие течения Подстановка (18.75) и (18.76) в уравнение (18.73) позволяет получить линейную относительно Лг!""' систему уравнений Ы ПН$в (!8.77) Аналитические выражения для — Р+ ʄ— !1+ 8„, К и 8 получены Бимом и Уормингом. Значение р' в (18.77) при решении нестационарных задач полагается равным р, поскольку при р= =44+0.5 схема имеет второй порядок точности по времени. Однако при решении методом установления (9 6.4) стационарных задач для сокращения числа операций целесообразно использовать схему первого порядка по времени и положить 8' = О. Точность схемы (18.77) по пространству зависит от используемой дискретизации; Бим и Уорминг рекомендуют использовать центральные разности второго порядка.

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

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

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

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