Главная » Просмотр файлов » Fedorenko-RP-Vvedenie-v-vychislitelnuyu-fiziku

Fedorenko-RP-Vvedenie-v-vychislitelnuyu-fiziku (810773), страница 21

Файл №810773 Fedorenko-RP-Vvedenie-v-vychislitelnuyu-fiziku (Fedorenko-RP-Vvedenie-v-vychislitelnuyu-fiziku) 21 страницаFedorenko-RP-Vvedenie-v-vychislitelnuyu-fiziku (810773) страница 212020-08-18СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

Но мы не знаем з)„знаем только оценку для нее и то, что это есть величина сколь угодно малая, если шаги сеткя достаточно малы. Мы не имеем права утверждать, что из совпадения (с точностью до погрешности аппроксимации) уравнений следует, что и их решения совпадают с точностью до величин порядка погрешности аппроксимации. Это право мы получим, если докажем устойчивость разност ной задачи. Разностная задача называется устойчивой, если из Ьзи, = Г, + в',, следует, что Зи, — иД а С(яв'З+ Зв"З), причем постоянная С не зависит от сетки з (т.е.

от т, Л). Установить устойчивость обычно бывает очень трудно. Но для линейного уравнения теплопроводиости это можно сделать. Устойчивость явной схемы. Рассмотрим уравнение с теми же краевыми условиями и начальными даннымн: —,= — х(г, х) — +у(г, х). ви вГ ди1 9111 интеп'НРОЕАние УРАВнениЯ с члсзиь$ми ИРОизеодными 109 Это есть линейное уравнение с переменнымн козффициентами. Прежде всего заметим, что нужно конкретизировать вид норм. Устойчивость можно устанавливать в разных нормах. Здесь мы используем самую наглядную и надежную: ««2«« =шах «Г(С, х)«, ««и«« =шах «и(С, х)«и тд.

С,к С,к Введем аналогичные нормы для сеточных функций: ««и,«« =шах «и,"„«, ««и"««шах «и" «н т.д. ,л Сл Основой для установления устойчивости является следующая лемма. Лемма 1. Нормы сеточных функций и" и и"+' связаны между собой неравенством Ци" Ц ай!пах (Ци" Ц + т««1'««, ««ср««), ««ср««штах (««ср1ср1««, ««~рзсрз««), если выполнено условие Куранта ««х««тсссзн1/2 (мы ограничимся случаем «31 > О; в случае Ц = О, а = 1 оценки проще). Доказательство. Имеем соотношение +! л лл — «х" лп2(и" +, — и" ) — хл и (и" — и",) «+ ~'", и = — '2х -Сгзи,+ ~ — —,х -!сз — '2н +,сз и + СС Сс Л сл+ 2 лс+ 1 Отсюда следует (так как «и" «и ««и" ««, х > 0) « "+'« — ', х",с «« "««+ — ', 21" +,сз««и"««+ к! лс И + 1 — —,(х",С2+ х" +па) Ци"««+ т«Я.

СС2 Условие Куранта дает важное соотношение ! к л л 1 2 (Хсл-1/2+ Нсл+!/2)1 = 1 — 2 (Хли !я+ ХИЧП2), и мы получаем «ц"„+' «и ««и" ««+ 2«Я, т = 1, 2, ..., М вЂ” 1. Заметим, что ««и"+'««=(либо шах «и" '«, либо «ив+!«, либо «им~1«). сл 1, .... М - ! основы вычнслнткньной млтвмлтнкн ыо В первом случае имеем Ц и"т'Ц И Ци" Ц + тЦу'Ц, Во втором случае используем краевое условие: (а, + лр,)и" ~' = аи",~'+ Ь р,"+', аы Ц, > О. Вели Цил+1Ц [иат![ то [и"т'[ И [ит~'[, (а, + Й[3~) [и$+'[ И а,[и"+'[ + Ь[ р"+'[ т.е. Цй"Ц=[<"'[ [р,""~В[ ЦрЦ. Такую же оценку мы получим и в случае Ц и" "'Ц = [и„"+'[.

Итак, либо Цй+'Ц к Цф), либо Ци"+'Ц и ЦйЦ + тЦД. В любом случае Ци"+1Ц ишак (Ци"Ц + тЦуЦ ЦфЦ) Теперь докажем следующую теорему. Теорема 1. Явная линейная разностная схема при выполнении условия Куранта ЦнЦт/Ьз < 1/2 устойчива по начальным данным, краевым условиям н правым частям. Докаэатлельства Используем лемму 1 рекуррентно, обозначая ради простоты и" ш Ци" Ц, У вЂ” ЦУЦ и т.дл и» ~$ шах (ип — 1+ ту 1р) к и шах (шах [и" з+ ту р[+ ту ~р) = =шах (и" т+2т/, 1р+ ту) и к шах (шах [й '+ т/, та[ + 2ту, р + ту) = =шах (и" з+Эт/, ~у+ 2тУ) < к 1пах (ио+ нт~, ~у+ (л — 1)ту) Так как лт и Т, получаем результат Ци" Ц И Ци'Ц+ ТЦУЦ'+ Ц РЦ. Обозначая ЦттД = ЦиоЦ + ТЦД+ Цф~, запишем оценку в форме ЦиД и ЦттД, т.е.

ЦХ,, 'Ц к 1. Мы установили оценку нормы решения разностной задачи через нормы начальных данных, правых частей и краевых условий. Это еще й !11 интеггиговАние л«внений с чАстными нгоизводными 111 не совсем то, что нужно. Нам нужно установить, что при малых возмущениях начальных данных, правых частей и краевых условий решение изменится соответственно мало. Но это следует из линейности задачи (как известно, ограниченность и непрерывность для линейных операторов — это одно и то же).

Воспроизведем это рассуждение. Если г.,и, = Г„г,,й; = .и, + т!„ то в силу линейности г.,((1, — и,) = т!, из ограниченности г,, ' получаем Цс!', — и«Ц а Цг!,Ц. Таким образом, для линейных разностных задач устойчивость есть равномерная (по всем сеткам) ограниченность обратного оператора. Устойчивость неявной схемы. Покажем, что неявнав схема дает разностную задачу безусловно-устойчивую, т.е. для ее устойчивости не требуется выполнения условия Куранта.

Ограничимся доказательством следующей леммы, Лемма 2. При любых шагах Л, т, нормы сеточных функций и" и и" +' связаны неравенством Ци"+!Ц и шах (Ци"Ц + тЦЯ, ЦгрЦ). Доказательство, Имеем альтернативу: Ци«+!Ц = (либо 1и«е+1~, либо !и"„+!1, либо п!ах (и"+!!). В двух первых случаях, как было установлено выше, Ц и"+1Ц ж Ц!рЦ. Нужно исследовать третий случай.

Для неявной схемы имеем « ~ + г (и — !!г+ и +!гг) Л -тг«.! и«з «««! з «« «+! — +" + г'! -!!га -!+ г" +итие+!. Л Л Пустын — внутренняя точка, для которой Ци"+'Ц = 3 и" «1~. Тогда с +лг(и.". не+и."+иг) Ци""Ц и ж Ци"Ц + тЦУЦ + — '(и«нгЦи"+'Ц + и«! Ци"+'Ц) нли, после сокращения, Ци«+!Ц а Ци"Ц + тЦЯ, На этом мы закончим исследование устойчивости неявной схемы. Еще раз подчеркнем, что она носит безусловный характер: схема всегда устойчива. В этом ее отличие и решающее преимущество перед явной схемой, счет по которой возможен лишь при 1ч.! ОСНОВЫ ВЫЧИСЛИТЕЛЬНОЙ МАТЕМАТИКИ 112 ! с < 0.5Л9~~иЙ. В других задачах, как мы увидим, тоже появляется. зто характерное условие: явные схемы устойчивы лишь пря некоторых ограничениях на шаг по времени ж он должен быть достаточно малым относительно шага по пространству Л.

Переход к неявным схемам, как правило, либо снимает условие устойчивости, либо су1цесгвенно его ослабляет. Возникает естественный вопрос: действительно ли условие Куранта существенно для явной схемы (ведь оно было необходимо для проведения достаточно простых оценок) нлн, может быть, оно связано с грубостью оценок, а не с существом дела? Оказывается, условие Куранта носит принципиальный характер, его нарушение делает результаты расчета совершенно бессмысленными.

Покажем зто на простом примере (который, кстати, иллюстрирует возможный экспериментальный прием исследования устойчивости разностной схемы: он состоит в фактическом вычислении последствий «единнчной погрешности в начальных данныхэ). Если в качестве начальных данных взять иа = О, то решение будет нулевым (мы не учитываем здесь краевых условий, считая, что задача решается на бесконечном интервале: — м < т «). Таблица 9 Теперь возьмем начальные данные с изолированной погрешностью: иа = 1, остальные ио = О, и станем решать задачу по явной схеме, нарушив условие Куранта.

Для иллюстрации удобно взять т = Лз. Тогда и«+' = и", — и" + и,"„„, В табл. 9 представлены результаты, полученные для п = 1, 2, ..., 6. Из таблицы видно, что решение возрастает почти в три раза за шаг и примерно через 20 шагов достигает катастрофического значения (порядка 1О«). А ведь в расчетах делаются сотни шагов по времени! Обратите внимание на характерный (по т) профиль функции и".

Он носит «пилообрааный» характер: и" ж (-3)"( — 1)"О", где еа — «гладкая» сеточная функция. Это характерный признак вычислительной неустойчивости. Обмчно вычислители просматривают 5 111 интзгеиговАниз г»лвнкний с чАсгными пгоизводнымн 1Ш полученные результаты, строят графики сеточных функций. Часто уже внешний вид таких функций содержит «намек» на какое-то неблагополучие, иа сомнительность результата.

На рнс. 11 показаны два примерных графика и", л! = 1, 2, 3, ... Первый, естественно, воспринимается как сеточная проекция «хорошей» функции, второй — типичный пример «подозрительного» решения. Общее качественное соображение носит простой характер: в методе конечнмх и разностей каждое «событие» должно быть , .'ьО разрешено несколькнмн точками. «Сабы- 'ь с''ч Р тием» мы называем колебание функции, в переход с одного уровня на другой и т.п. ° ' ' ° - ° ° ' Если такое «событие» происходит на одном ' и ° ° ' счетном интервале — зто явно подозрительно, настораживает, делает результаты »! сомнительными.

Вычислители очень не любят «пилооб- Рис. 1 1 разных» графиков. Однако не следует все абсолютизировать. Не следует думать, что если получены точки в", легко укладываююциеся на гладкую функцию, то имеется гарантия правильности расчета. «Пила» на решении — тоже не 100 7,'-ная гарантия ошибочности расчета, хотя ничего хорошего в этом нет, Появление «пилы» на графике сеточной функции часто является признаком вычислительной неустойчивости разностной схемы. Но настоящая вычислительная неустойчивость сопровождается еще и очень быстрым нарастанием амплитуды «пнлы», настолько быстрым, что за несколько шагов решение может вообще выйти за пределы машинной бесконечности.

Сходимость равностных схем. (Точнее, следует говорить о сходимости приближенного решения к точному при с, Л- О.) Установив устойчивость схемы н оценив погрешность аппроксимации, воспользуемся теоремой Рябенького-Филиппова и получим оценку 110,— и,11 = О(т+ Л). В явной схеме т О(Лз). Такое же соотношение во многих случаях приходится выдерживать по соображениям точности расчета и в неявных схемах (см. 5 21). Было бы желательно иметь в оценке О(т+ Лз), тем более что почти во всех узлах сетки невязка есть О(Лз); мешает только аппроксимация краевых условий с погрешностью О(Л). Улучшим ее, используя характерный прием. Выпишем погрешность аппроксимации (2) более аккуратно, используя ряд Тейлора: У~~ ~ У(1„, Л) = ЯВ„, О) + ЛУ„+ -»»" У „+ О(Лз). основы вычислитгльной млтвмлтики сс4 !ч.

! Тогда !со ! — а,, + Рсзо — Ч>," = — а,(с', + РУ вЂ” 1с! — — асЬУ„„+ 0(Ьа). Перенося главную часть погрешности — а,Ь!С„из правой части ! в левую н заменяя 4У„„разностной аппроксимацией, получаем аппроксимацию краевого условия второго порядка: — ' „'+ ВР;+ — „пс(и, «вЂ” Щ+ и«) — р!- = -а,!с„+ В,У вЂ” ~р, + О(Ь'). В реализации явной схемы никаких осложнений не возникает. В неявной схеме это приводит к нарушению трехдиагональной структуры уравнений на верхнем слое. Предоставим читателю внести необходимые дополнения в алгоритм решения уравнений на верхнем слое прогонкой.

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

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

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

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