Anderson-et-al-1 (1185923), страница 30

Файл №1185923 Anderson-et-al-1 (Андерсон Д., Таннехилл Дж., Плетчер Р. Вычислительная гидромеханика и теплообмен) 30 страницаAnderson-et-al-1 (1185923) страница 302020-08-25СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

Ключевым моментом, определяющим эффективность применения неявного метода переменных направлений для решения эллиптических уравнений, является выбор значений итерационных параметров р». Этот выбор можно осуществлять либо по методике, предложенной Писменом и Ракфордом [Реасешап, цасЬ1огд, 1955[, либо по методике, предложенной Вахспрессом [%аспзргезз, 1966[. Хотя имеющийся опыт применения этих методик не позволяет узнать, какая из них лучше, результаты некоторых исследований указывают на то, что вторая из них предпочтительней. Читателю, который захочет воспользоваться 168 Гл. 4.

Метод конечных разностей длн модельных уравнений неявным методом переменных направлений, мы советуем сначала ознакомиться с работами, посвященными выбору итерационных параметров рь Сопоставить времена счета при использовании точечно- н блочно-итерационных методов с последовательной верхней релаксацией не просто из-за сложности подбора оптимального значения параметра верхней релаксации. Кроме того, результат такого сопоставления во многом определяется конкретными особенностями рассчитываемой задачи, граничными условиями и числом узлов разностной сетки.

Блочно-итерационные методы требуют обычно меньше итераций, чем точечно-итерационныс методы, но, как мы уже указывали, на каждой итерации приходится проводить больший объем вычислений. Имеющийся опыт показывает, что для большинства задач время счета для достижения одинаковой точности при применении методов последовательной верхней релаксации по строкам и Гаусса — Зайделя с последовательной верхней релаксацией практически одинаково. Неявный метод переменных направлений с последовательной верхней релаксацией (при постоянном параметре итераций) позволяет часто сократить время счета на 20 — 40 $ по сравнению с методом Гаусса — Зайделя с последовательной верхней релаксацией. При переменном параметре итераций можно обычно достичь еще большего сокращения времени счета. Сильно неявные методы.

В последние годы разработан новый тип блочно-итерационной процедуры решения системы алгебраических уравнений, возникающей при дискретизации эллиптических уравнений в частных производных. Для иллюстрации этого подхода рассмотрим систему алгебраических уравнений, получающуюся при использовании пятиточечной схемы для решения уравнения Лапласа. Запишем эту систему уравнений в виде [А]и=С, где [А] — сильно разреженная матрица коэффициентов, н — вектор-столбец неизвестных, а С вЂ” вектор-столбец известных величин.

Стоун [Ропе, 1968] предложил решать эту систему уравнений методом факторизации с использованием сильно неявной процедуры. Суть этого метода состоит в замене разреженной матрицы [А] матрицей [А + Р], которая представляется в виде произведения верхней [Ц и нижней [Е] разреженных треугольных матриц. Если матрицы [Е] и [Ц неразреженные, то эффективность рассматриваемого метода оказывается близкой к эффективности метода исключения Гаусса; следовательно, успех применения метода факторизации определяется выбором мат- э 4.3.

Уравнение Лапласа 169 рицы [Р). Важно, чтобы элементы этой матрицы были малы по абсолютной величине, а получающаяся в результате система уравнений должна быть более неявной, чем при использовании метода переменных направлений. Для построения итерационного алгоритма перепишем систему уравнений [А]ц = С в виде [А+ Р]ц"+'=С+ [Р) п".

Представим матрицу [В] = [А + Р) как произведение верхней [0] и нижней [Ц треугольных матриц. Тогда [Ц [и] ц"" = С+ [Р] ц". Если ввести промежуточный вектор Ч"+' =[(7]ци+', то получим следующий двухшаговый алгоритм: Шаг 1 [ь] Ч"+ = С + [Р] ц". [и] и"" = Ч"". (4.!27а) (4.127Ь) Шаг 2 По этому алгоритму расчет проводится на каждой итерации. На шаге 1 осуществляется прямая подстановка, а на шаге 2 — обратная. Стоун [8(опе, 1968) выбрал матрицу [Р) таким образом, чтобы в матрицах [Е] и [(7] было только по три ненулевые диагонали, а главная диагональ матрицы [(7] состояла бы из одних единиц.

Кроме того, элементы матриц [Ц и [Ц были подобраны так, что ненулевые элементы матрицы [А] совпали с расположенными на их месте элементами матрицы [В), т. е. матрица [В] отличалась от матрицы [А] появлением двух новых ненулевых диагоналей. Элементы матриц [ь], [Ц и [Р] можно определить по заданным уравнениям, основываясь на произведении [Ц [Ц.

Подробно этот метод описан Стоуном [8(опе, 1968). Рассмотренный метод неявный как по х, так и по у. Проведенные численные расчеты показали, что время решения уравнения Лапласа этим методом составляет лишь 50 — 607а времени решения неявным методом переменных направлений. Шнейдер и Зедан [ВсЬпе(бег, Еес(ап, 1981] предложили новый способ построения матриц [Ц и [(7], позволяющий сократить время расчета в два — четыре раза по сравнению с методом Стоуна (Яопе, 1968). Получившийся в результате метод авторы назвали модифис(ироеанным сильно неявным методом.

Алгоритм решения уравнений этим методом также описывается соотношениями (4.127а) и (4.127Ь); улучшение же достигнуто благодаря распространению метода Стоуна на девятиточечную схему 170 Гл. 4. Метод нонечных разностей длн модельных уравнений решения уравнения Лапласа. Модифицированный сильно неявный метод можно использовать и для решения разностиых уравнений, получающихся при применении пятиточечной схемы, при этом, как уже отмечалось, время расчета сокращается в два— четыре раза. По-видимому, новый метод является эффективным и довольно общим методом решения алгебраических уравнений.

Подробно он описан в приложении С. (4.128) Первое и второе слагаемые в левой части этого уравнения являются соответственно нестационарным и конвективным членами, а в правой части стоит вязкий член. Если вязкий член не равен нулю, то уравнение (4.128) параболическое; если же он равен нулю, то в уравнении остаются лишь нестационарный и нелинейный конвективный члены. Такое уравнение гиперболическое и имеет вид — + и — =О. ди ди д~ дх (4.129) Его можно рассматривать как модельное для уравнений Эйлера, описывающих движение идеального газа. Уравнение (4.129) есть нелинейное уравнение конвекции и обладает некоторыми математическими особенностями, к рассмотрению которых мы сейчас перейдем.

После этого мы опишем различные разностные схемы, используемые для решения невязкого уравнения Бюргерса. При этом будут приведены типичные результаты, полу- 9 4.4. Уравнение Бюргерса (невязкое течение) Мы изучили различные конечно-разностные методы и применили их к решению простых линейных задач. Это позволило нам лучше понять эти методы и познакомиться с их основными специфическими особенностями.

К сожалению, в гидромеханике обычно приходится решать нелинейные задачи, так как давление, плотность, температура и скорость должны быть определены из решения нелинейной системы уравнений в частных производных. Полезно сначала изучить какое-то одно простое нелинейное уравнение, аналогичное уравнениям гидромеханики. Это уравнение должно включать в себя члены, описывающие те же физические процессы, что и члены, входящие в уравнение гидромеханики, т. е. конвективный, диффузионный или диссипативный и нестационарный члены.

Такое простое нелинейное уравнение было предложено Бюргерсом [Впгйегз, 1948). Оно имеет вид ди- Ь ди дзи — -+ и — =И вЂ” ° дт ' дк дкв ' й 4.4. Уравнсние Ьюргерса (невинное течснис] В общем случае и неизвестная и, и функция Р(и) — векторы. Перепишем уравнение (4.130) в виде — +А — =О, ди дн д1 дх (4. 131) где в общем случае А = А(и) — матрица Якоби дР;/диь а в нашем простом примере А = с(Р/с(и. Так как наше уравнение (или система уравнений) в частных производных гиперболическое, то все собственные значения матрицы А вещественные.

Гладким называется такое решение уравнения (4.131), когда функция и непрерывна внутри области, а ее производная может иметь разрыв на границе (т. е. решение уравнения является непрерывным по Лнпшицу). Слабым называется решение уравнения (4.131) гладкое всюду, кроме некоторой поверхности в пространстве (х, 1), на которой функция и может иметь разрыв. На величину скачка функции и при переходе через поверхность разрыва накладываются определенные ограничения.

Пусть се — произвольная непрерывная векторная функция с непрерывной первой производной, равная нулю вне некоторой ограниченной области, тогда и называется слабым решением уравнения (4.130), если ~~ (ао,и+ се,Р) с(хе(1+ ~ св(х, 0)ф(х) с(х= О, (4.132) о причем ф(х)=и(х,О). Гладкое решение всегда является одновременно слабым решением, а всякое непрерывное слабое чаемые при расчетах по многим широко используемым разностным схемам, и выяснена роль нелинейных членов. Уравнение Бюргерса будет рассмотрено в $4.5. Уравнение (4.129) можно интерпретировать и как нелинейное волновое уравнение, при этом скорость распространения волны в различных точках будет разной.

В противоположность этому в ранее изученном линейном одномерном уравнении конвекции (линейном волновом уравнении) (4.2) скорость распространения любых возмущений была постоянна. Так как скорость распространения возмущений меняется, то характеристики начинают пересекаться и в решении возникают разрывы, аналогичные ударным волнам в газовой динамике. Следовательно, рассматриваемое одномерное модельное уравнение позволяет изучать свойства разрывных решений. Нелинейные гиперболические уравнения в частных производных обладают решениями двух типов согласно Лаксу [1ах, 1954].

Поясним это на примере простого скалярного уравнения (4.130) 172 Гл. 4. Метод конечных разностей для модельных уравнений решение является гладким. Подробно теория слабых решений изложена в прекрасных книгах ['(1ьг)11(Ьагп, !974) и 1Лейгеу, Тая)ц(1, 1964). Теория слабых решений гиперболических уравнений в частных производных — относительно недавно созданная математическая теория. Одним из примеров слабого решения являются ударные волны, возникающие в сверхзвуковых течениях невязкой жидкости. Интересно заметить, что решения уравнений газовой динамики с ударными волнами были известны за 50— 100 лет до того, как была построена теория слабых решений систем гиперболических уравнений в частных производных.

рис. 4.23. Типичная задача о распространении разрыва для уравнений Бюр- герса. Вернемся к невязкому уравнению Вюргерса и найдем условия существования слабого решения этого уравнения, т. е. необходимые условия существования решения с разрывом, как изображено на рис. 4.23. Пусть ш(х,() — произвольная непрерывная функция, имеющая непрерывную первую производную. Кроме того, пусть она обращается в нуль на границе В области Р и вне Р (на дополнении к Р); Р— произвольная прямоугольная область в плоскости (х,1).

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

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

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