Главная » Просмотр файлов » Н.Н. Калиткин - Численные методы

Н.Н. Калиткин - Численные методы (1133437), страница 85

Файл №1133437 Н.Н. Калиткин - Численные методы (Н.Н. Калиткин - Численные методы) 85 страницаН.Н. Калиткин - Численные методы (1133437) страница 852019-05-12СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

Разностная запись первого краевого условия сводится к заданию решения ул в граничных узлах сетки, т. е. при п=О, и=/[/, я[=0 и т=М. МНОГОМЕРНОЕ УРАВНЕНИЕ (59а) (59б) где разностные операторы Л„определены формулами (57). Как обычно, под у подразумевается значение на полуцелом слое Исследуем продольно-поперечную схему. Вычисление р аз ноет ного решен н я. Переходна полу- целый слой делается при помощи уравнений (59а). Согласно определению оператора Л, каждое такое уравнение содержит три неизвестных значения: ул,, ул и ул„т; остальные значения у берутся с исходного слоя.

Иными словами, при таком переходе схема неявна по направлению х, и явна по направлению х,. При любом фиксированном индексе т уравнения (59а) образуют относительно неизвестных ул линейную систему с трехдиагональной Если вести расчет по абсолютно устойчивому варианту схемы (О '(,), то можно брать т 6. Но тогда на каждом слое надо решать линейную систему Л1» уравнений. Даже с учетом того, что ее матрица ленточная с шириной ленты 2)У'»-', решение этой системы методом Гаусса требует М'»-' действий.

Поскольку для расчета до момента Т теперь надо делать )у шагов по времени, то полный расчет требует )у'л»-' действий. Значит, для двумерной задачи (У=2) неявная схема (55) и явная схема приводят примерно к одинаковому объему вычислений, а в 9 1 мы видели, что явная схема обладает плохими свойствами и невыгодна для расчетов. При р~2 неявная схема (56) даже невыгодней явной. Однако для многомерного параболического уравнения построены абсолютно устойчивые схемы, позволяющие вести расчет шагом т Й и требующие только У» действий для перехода со слоя на слой (т. е. число действий в расчете на одну точку сетки не зависит от шагов Л,). Такие схемы называются экономичными.

Подавляющее большинство многомерных расчетов проводится по таким схемам. В следующих пунктах те+У мы рассмотрим два основных вида экономичных схем для параболического уравнения — продольно-поперечную и локально-одномерную схемы. 2. Продольно-поперечная схема, называемая также схемой переменных Рис. 83. направлений, является одной из лучших двумерных экономичных схем. Выберем изображенный на рис, 83 шаблон, содержащий полуцелый слой (=1-)-т~2, и составим на нем эту схему: 1 О Олт (Улт Улт) = Л1У»т+~»Улт+1лт~ 1 (Упт Улт) =Л1рпт+Л»Упт+~пт~ ПАРАБОЛИЧЕСКИЕ УРАВНЕНИЯ 1гл.

х1 матрицей. Поэтому значения у легко вычисляются одномерной прогонкой по индексу н, т. е. по направлению х,. Наоборот, при переходе при помощи уравнений (59б) с полу- целого слоя на целый схема явна по направлению х, и неявна по х,. Поэтому решение у„,„на целом слое вычисляется тоже одномерной прогонкой, но в поперечном направлении х,.

Нетрудно подсчитать, что для перехода с целого на целый слой нужно всего 20 — 30 действий на каждую точку сетки независимо от величин шагов, так что схема экономична. Как и в одномерной схеме (6), диагональные матричные элементы в уравнениях (59) преобладают; следовательно, прогонка устойчива, а разностное решение существует и единственно. У с т о й ч и в о с т ь продольно-поперечной схемы исследуем методом разделения переменных. Множители роста гармоники на первом и втором полушаге по времени могут быть различными. Поэтому положим учш=ехр (щх1в+1гхва), у =рч~у, у =уды (60) Подставляя соотношения (60) в схему (59), получим множители роста р~,=(1 — — „., Е(п" — '$1+ — „, з(п'- — "), (61а) (61б) Нетрудно заметить, что для всех гармоник при любых шагах выполняется неравенство ~р,',р,",~==1.

Таким образом, при переходе с одного целого слоя на следующий целый слой ошибки начальных данных не нарастают, и схема (59) равномерно и безус'ловно устойчива по начальным данным. Нетрудно проверить, что дополнительный признак устойчивости по правой части (9.54) выполняется на каждом полушаге по времени. Следовательно, схема (59) устойчива по правой части. Замечание 1. Если йт)й',, то существуют такие гармоники, которые усиливаются при переходе с пелого слоя на полу- целый; например, ~р',„~) 1. Зато при переходе с полуцелого на следующий целый слой эти гармоники настолько затухают, что в целом усиления не происходит.

Аналогично, при нт ) 61 есть гармоники, усиливающиеся при переходе с полуцелого слоя на целый. 3 а меча н ие 2. Суммарный множитель роста р,=р,',р,", таков, что ~ ре,(=1 только при о=с=0; для всех остальных гармоник ~ р, ~ (1. Следовательно, продольно-поперечная схема обладает аппроксимациоиной вязкостью и расчет по ней должен приводить к сглаживанию разрывов. МНОГОМЕРНОЕ УРАВНЕНИЕ 393 Аппроксимация.

При переходе с целого на полуцелый слой каждая пространственная разность вычисляется несимметрично но времени и погрешность равна 0(т+й'). Но ошибка на второй половине слоя компенсирует первую, и в итоге при переходе с целого слоя на целый погрешность локальной аппроксимации на равномерных сетках есть 0(та+А',+й!).

В этом легко убедиться при помощи следующего преобразования. Вычитая уравнение (596) из (59а), получим ! т уппг 2 (ул'в + улм) 4 Лз (уюл ул~л) (62) Сложим уравнения (59) и подставим в них полученное значение у„„, исключив тем самым полуцелый слой: ! 1 т (Улпг Улм) = 2 (АХ+Аз) (Улм+Упм) 4 ~ГЛа (Улм Упт)+пгпт. (63) Предпоследний член справа есть тзи, „.г/4=0(т'), а остальные члены в (63) совпадают с симметричным вариантом схемы (56), который соответствует а='(з и имеет аппроксимацию 0(т'+й', + + Ц). Поскольку продольно-поперечная схема отличается от этого варианта на член 0 (т'), она также имеет второй порядок аппроксимации по всем переменным. Остааовимся на аппроксимации краевых условий (55б).

На целом слое в уравнения продольно-поперечной схемы (59) входят значения решения рл на сторонах прямоугольника х,=о и ха=5; очевидно, надо положить у~о !"з ( тл' )' "«и ро ( т»' )' Гьаи полУцелого слоЯ тРебУютсЯ значениа Уп на стоРонах х, О и х,=а. Полагать (у — )ьо)трпа О невыгодно, ибо полуцелый слой не вполне соответствует моменту 7 и такая аппроксимация внесла бы погрешность 0 (т). Следует воспользоваться уравнением (62), отнесенным к стороне хт=О: 1 уот о (Ргпг+ Шм) 4 Ла (Шм Шм) 1 ( гл ( гн ! (546) аналогичное условие записывается для стороны хт=а.

Граничные условия (54) обеспечивают погрешность аппроксимации О (то). Сходи мое ть. Проведенное исследование аппроксимации и устойчивости показывает, что схема (59) безусловно сходится в (( ((!ы причем в прямоугольной области на равномерной сетке и при краевом условии (64) она имеет точность 0(тз+й'!+)т',) на решениях с непрерывными пятыми производными. Более сложными методами можно доказать равномерную сходимость со вторым порядком точности. Отметим некоторые усложнения исходной задачи (55). ПАРАБОЛИЧЕСКИЕ УРАВНЕНИЯ [гл.

х! Произвольная область. Пусть для уравнения (55а) в области произвольной формы заданы краевые условия первого рода [и — р (х, ())г = О. Тогда разностное краевое условие (64б) не удается применить, ибо неясно, как вычислять Л,р. Приходится ограничиться усло- виями [У вЂ” (А(х, й)1т — — О, (65) [д — (А(х, г)1„=0, где у — множество граничных узлов, Погрешность аппроксимации условия (65) на полуцелом слое равна 0(т). Поэтому в произвольной области схема (59) сходится с точностью 0(т+й1+Ье). Если область ступенчатая, т. е. составлена из прямоугольников со сторонами, параллельными осям координат, то в ней можно написать краевое условие повышенной точности (64б). В этом случае схема (59) имеет второй порядок точности. Переменный коэффициент теплопроводности.

Для уравнения теплопроводности с переменным коэффициентом можно составить два варианта продольно-поперечной схемы, являющихся обобщением наилучшей схемы (34). В первом варианте на всех слоях — исходном, полуцелом и новом целом — разностный коэффициент теплопроводности к приписывают полуцелому слою Г; во втором варианте на этих слоях берут соответственно х (Г), к(у) и х(1). Оба варианта успешно применяются на практике. Второй вариант лучше исследован теоретически; для него доказана безусловная сходимость в прямоугольной области с точностью 0 (тз+ + Й1+й;.), если коэффициент й(х, г) непрерывен со своими вторыми производными. Анизотропная теплопроводность в простейшем случае приводит к тому, что по каждому направлению имеется свой коэффициент й,(х, 1). В этом случае уравнение теплопроводности принимает вид дГ ~~~ адх [ а( ' ) д ~+[( ~ ) а=1 (66) Продольно-поперечная схема, ее обобщения и все теоретические обоснования переносятся на этот случай практически без изменений.

3. Локально-одномерный метод. Продольно-поперечная схема на задачи с числом измерений р=З непосредственно не обобщается. В самом деле, введем р — 1 промежуточный слой и на МНОГОМЕРНОЕ УРАВНЕНИŠР— А„и, Д'2 А„=йа д; ~ха ка = сопз1, а=! х = (х„х„..., х ). Аппроксимируем это уравнение симметричной неявной схемой, которую назовем исходной: Р ! .

! — (у-у)= —,,~ Л (у+у) где Л,— разностные операторы, аппроксимирующие А, с погрешностью 0(й'); обычно для них используют формулы (57), соответствующие г = 2. Благодаря симметричной форме исходная схема имеет погрешность 0 ~т2+ Х й„'!. а Однако эта схема неэкономична, потому что не найдено хорошего алгоритма вычисления у. Наряду с исходной схемой построим локально-одномерную схему. Введем промежуточные слои и на каждом слое в правой части (68) вместо,У',Ла возьмем рЛ„; в левой части поставим а шаг т/р. Обозначим решение на промежуточных шагах через ш„ (я=1, 2, ..., р).

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

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

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

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