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

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

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

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

й 12. Сцеитраиьный признак устойчивости Рассмотрим основной аналитический аппарат исследования устойчивости разностных схем, который имеет дело не с реальной вычислительной схемой, а с некоторой ее моделью. Он связан с более нли менее обозримой и выполнимой аналитической работой, благодаря чему и получил самое широкое распространение.

Хотя этот метод исследования не дает точного ответа на вопрос аб устойчивости, он позволяет отбраковать подавляющее большинство заведомо неустойчивых схем, а схемы, признанные на основе спектрального признака устойчивыми, как правило, на самом деле являются таковыми. Начнем с двух упрощений, которые приходится произвести, чтобы можно было применять аппарат спектральной устойчивости. Имеются в виду: а) линейные, однородные с посюянными коэффициентами схемы; б) задача Коши на всем пространстве, без краевых условий (нх место занимают условия тина «ограниченности на бесконечности»). Итак, если нас интересует разностная схема для общего уравнения теплопроводностн — — и(с, х, и) — "„+ у(с, х, и), то исследование проводится для уравнения ди ди — = и — + аи, и, а сопз!.

ас ах' СПЕКТРАЛЬНЫЙ ПРИЗНАК УСТОЙЧИВОСТИ Й !21 Рассмотрим явную разностную схему (1.!0): и" +1 — и" и" — 2и" + и" «1, ии х 2 +аи, Т ЛЗ т = О, + 1, ..., + м, и = О, 1, ..., М = Т/т. Исследование основано на следующем общем факте: все линейные однородные разностнне уравнения с постоянными коэффициентами, заданные на всем пространстве (с условием ограниченности на бесконечности), имеют универсальное полное семейство частных решений ии ЛВЕ™Т 0 И ~р И 2я Здесь р — параметр семейства, Л(т, Ь, «схема»; р) — функция, зависящая от шагов т, й, параметра р и вида схемы.

Каждая схема характеризуется своей функцией Л( р) (остальиые аргументы (т, й, «схема») мы будем всегда иметь в виду, не выписывай их явно). Функция Л( р) называется спектральной функцией схемы. Совокупность значений, пробегаемых точкой Л(р) (в комплекснои плоскости), когда р пробегает )О, 2я), называют спектром разностной скемы. Введем формальные определения.

Разностную схему называют спектрально-устойчивой, если ) Л( р) ) и 1+ Ст, Ч р Е )О, 2л), (2) где С вЂ” не зависящая от т постоянная. Другими словами, спектр устойчивой (по спектральному признаку) схемы должен лежать в Ст-расширении единичного круга. Разностную схему называют спектрально-неустойчивой, если существуют а> 1 (а не зависит от т) и р ~ )О, 2я), такие, что (3) Это пока чисто формальные определения. Сейчас мы научимся вычислять спектр разностных схем, а затем выясним содержательный смысл введенных понятий.

Он будет простым: спектрально-неустойчивые схемм не годятся для вычислений, расчеты по хаким схемам юпровождаются катастрофическим нарастанием последствий погрешностей вычислений (т.е, погрешностей машинного представления чисел, округления и т.п.). Примеры вычисления спектра. Рассмотрим примеры вычисления спектра для различных разностных схем с учетом вышеприведенных определений. 1ч. г основы вычислительной млтемлтнкн Явная схема. Вычисление Х( р) проводится просто: нужно решение и" = Х"е'"т подставить в разностные уравнения. для явной схемы имеем гыыеее« вЂ” 1«е""т Ге~'" нт — ЗРен'т+ Ге" +Не Сокращая на Х"е'"'т, получаем Л вЂ ! е ч-3+ее» Используем соотношение е 'т — 2+ е'г = — 4з!пг(р/2), В результате спектральная функции явной схемы принимает вид Л(р) =1 — 4 — 'з1п х.

Л Легко видеть, что 2( р) вещественна и Ци) = 1 — 4т//гз < Х(,р) < 1 Л(0) Итак, спектр есть отрезок 11 — 4т//гз, 1]. Условие устойчивости: ! 1 — 4т//гз и — 1, и т % /гз/2. Это есть условие Куранта, которое нам уже знакомо. Таким образом, явная схема для уравнения теплопроводности устойчива при выполнении условия Куранта (условно-устойчива). Неявная схема. Проведем те же вычисления: Гее«т(1 — О ° ~ ° е'« — а+е ч ~ л"е мтй т * После очевидных преобразований получаем Х( р) = 1 + 4 — ' з!пз х Л~ Очевидно, Х( р) Е 10, 1), 'т' т, Ь. Неявная схема безусловно-устойчива (по спектральному признаку).

Этот факт (правда, с другим пока смыслом термина «устойчивость») нам уже известен. Схе.ча «кресел» для волнового уравнения ин = и,, Схема имеет вид спектРАльный пРизнАк хстойчизости ф 2г! !27 ПОдетаВЛяя и«!!"Е1"'Т И СОКращая На Л" 'Е'вт, ПОЛуЧаЕМ 22 — гг+! Х . 2« = -4 — з!и 2 Л2 г' т.е. )! есть решение квадратного характеристического уравнения 3Р— 2 1 — 2 — 'з!ПЗМЛ+! О. — л2 г/ Исследовать спектр можно, не решая уравнения.

Заметим, что свободный член есть единица, т.е. )11 1 = 1. Здесь имеются две возможности: а) если корни вещественны, то один корень меньше единицы, второй больше единицы, т.е. схема неустойчива; б) если корни комплексно-сопряженные, то «3!1« = ! 3 ! = 1, т.е. схема устойчива. Итак, схема устойчива, если корни комплексные, т.е.

если отрицателен (при всех р) цискриминант 1 — 2 — 'з!Пгх — 1 1 — 4 — 'з«пэ т -«-4 — 'з!П«х — 1 = „2 г! г «л г =4 — 'з!пг~(т з1пгй- ! Л2 г ! Л2 г Очевидно, что при всех г а !О, 2п) это выражение отрицательно только для т/Л и 1, Это н есть условие Куранта для схемы «крест».

Шахматпиал схема. Рассмотрим систему уравнений, описывающих распространение звука и,+22„=0, 221+и„=О. Поясним некоторые новые обьекты. Прежде всего удобно ввести так называемую шахматную сетку, т.е. определить сеточные Функции и и 21 в разных точках. Итак, введем «целые» точки, илн и-точки: Т„° лт, х ° п2Л. В этих точках определим и" . Введем «полуцелые» точки, нли 22-точки: /„+н — — (п+ 1/2)т, х +„— — (п2 + 1/2)Л.

В этих точках опРеделим 22"+ 11!22. На такой сетке удобно аппроксимировать систему: л+1 ил,в+1П,А+И »л+1П »л-1П 1» » +1П » -1П » +1П » +1П " +1 Обобщим конструкцию стандартного решения ОСНОВЫ ВЫЧИСЛИТЕЛЬНОЙ МЛТЕМАТИКИ !!в !ч.! где у, )е — некоторые постоянные. Подставляя это решение в разностные уравнения, после сокрашения на )РСЬВТ и й" 'еп"'Тпг>Т получаем и —,+Р „=О, Р=,+и2=„— '=О, 2-1 ечгг — е 'Т~г Х-1 енп — е ЧЕг Система (относительно !1, !е) имеет нетривиальное решение при Это и есть уравнение, определяющее й( р): — ) +4 — з)п х=О, 1 г 2 ) ьг 2 г ),г — 22 ! — 2 — 'гйпгх~! + ! =О г 21 Такое уравнение мы уже исследовали в связи со схемой «крест» для волнового уравнения, и ответ нам известен: схема устойчива при условии Куранта т < л. Схема «рогеб» б!Вя уравнения нгеплолроводности.

Этот пример интересен тем, что он связан с поиском явных безусловно-устойчивых схем. Схема имеет вид Эта схема трехслойная, расчет требует задания двух начальных слоев иь и и' (и' можно вычислнтгч например, по явной двухслойной схеме).

В приведенном уравнении предполагаются известными значения и" ' и и", и"+! явно выражается через известные величины на двух предыдущих слоях. Стандартное исследование устойчивости приводит к характеристическому уравнению В ре-Ф 2т Ьг Обозначая г = йг/2Т, представим уравнение в другой форме: Хг — 2 еее Х Л + 1 " = О.

1+г !+е Решение выписывается просто: 1+е Ф 121 СПВКТРАИЬНЫЗ ПРИЗНАК УСРОЙЧИВОСТИ Рассмотрим два случая. а) созз р — 1 + гз < О. Корни комплексно-сопряженные, нх произведение ~ Л,1, ~ = ~ (1 — г)l(! + г) ! < 1, т.е. !Л, ! = ! Лз! < 1. б) рзш гз — (1 — созз р) > О. Очевидно, что р< г и корни Л< з = (соз 'р~ р)/(1+ г). Несложный анализ, основанный на том, что ~соз ф ж1, ! ~ р~ < г, показывает, что в этом случае корни ~ Л,! <1, ~Лз~ < 1.

Итак, схема «ромб» безусловно-устойчива. К сожалению, она не годится для решения уравнения теплопроводности, так как не аппроксимирует его. Причина этого состоит в замене значения и" на 0.5(и" +'+ ил '). Погрешность такой замены есть О(тз) и была бы допустимой в других ситуациях, но это среднее используется при аппроксимации второй производной, т.е.

в выражении, делящемся на ЬЗ. В результате в погрешности аппроксимации появляется член О(тз/ЛЗ), что делает такую схему допустимой лишь при очень малых т, например при т = 0(йз), т.е. мы не получаем серьезных преимуществ от безусловной устойчивости. Однако она представляет определенный интерес, особенно в задаче «на установление», когда решение уравнения теплопроводности (с не зависящими от времени правыми частями н краевымн условиямн) в пределе при /- м переходит в решение уравнения Пуассона. В этом случае детали процесса выхода решения на предел игнорируются.

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

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

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

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