1626435584-7c6402f545ecf856225d6cf8d21519c9 (844233), страница 69
Текст из файла (страница 69)
Получим следующее выражение: АППРОКСИМАЦИЯ 4 2! что на равномерной сетке х„41л — х„112=6, получим разностиую схему !г(Ул уп) = 2 ~ (сл+ 122 (Улл1 Ул) яп — 112 (Ул Уп 1)). (27) Если й=сопз(, то схема совпадает с неявной схемой (!6). Интегро-интерполяциоиный метод особенно полезен для уравнений с негладкими или разрывными коэффициентами, поскольку именно интегральная запись законов сохранения выделяет из всех математически допустимых решений таких уравнений физически правильное обобщенное решение. Метод неопределенных коэффициентов заключается в том, что в качестве разностной схемы берут линейную комбинацию значений разностного решения в узлах шаблона. Коэффициенгы этой линейной комбинации определяют из условия, чтобы невязка схемы имела как можно более высокий порядок малости относительно т и !2.
Например, для уравнения и, = йил на шаблоне рис. 52 будем искать разностную схему в следующем виде: аул, + руп+ уул.„+ буп = О. (28) Подставим сюда разложения (24), ограничиваясь для простоты членами О (т) и О (й'), и вычтем схему (28) из исходного уравнения. Получим невязку (индекс и всюду опускаем) 2Р= и,— й脄— (и+р+у+6) и+тб и,+ +(и — у) йи — — (а+ у) й'илл+6 О (с')+ (а — у) О (62)+ Чтобы сократились выписанные здесь члены, надо положить а+()+у+6=0, сб= — ), и — у=О, — (а+у) й'= — й.
2 Отсюда находим коэффициенты: 1= — „, +-, 6= — —. 2А ! ! Подставляя их в (28), получим разностную схему ()8). Метод неопределенных коэффициентов применйм на косоугольных сетках. Например, при его помощи нетрудно составить пяти- точечную схему для уравнения и1=йи„на треугольной сетке с шаблоном рис. 53: ! - ° Л л1 !3 2!2'1 2т (Уп+Уп-2) (Зт+ Ап/ (Ул-1+Упп2)+(42 Лп) Ул. (29) УРАВНЕНИЯ В ЧАСТНЫХ ПРОИЗВОДНЫХ [гл. 1х в этих узлах нетрудно написать разностные уравнения (17)1 и лог Уо = К ((т), УА1 = Рг (Ге) которые являются точными (их невязка равна нулю).
Более сложен случай второй краевой задачи для того же уравнения (дальше будем рассматривать только левое условие): ти и-т и„(0, ()=р,((), и,(а,()=р,((). (30) Можно аппроксимировать производную односторонней разностью: Л (У1 Уо) = Р1(гто1). (31) Однако невязка этого разностного уравнения равна 1 . Ь оро = (й„) о — — „(й, — й,) = — — и„= 0 (й), (32) т. е. имеет меньший порядок малости, чем невязка (25) в регулярных узлах.
Это приводит к понижению общей точности расчета. Рассмотрим способы написания разностного краевого условия нормальной точности 0 (Ло). Сделаем это на примере явной схемы (18). С п о с о б ф и к т и в н ы х т о ч е к очень нагляден. Введем вне отрезка 0 =хк-а фиктивную точку х,=х,— й и будем считать исходное уравнение справедливым прн х 1~х. Тогда разностное Возможны случаи, когда часть коэффициентов схемы типа (28) определяют из условия наивысшего порядка малости невязки, а часть коэффициентов выбирают из других соображений. Метод неопределенных коэффициентов (как и метод разностной аппроксимации) применйм к уравнениям с непрерывными и достаточное число раз дифференцируемыми коэффициентами и решениями.
Из-за сравнительной громоздкости он применяется реже двух ранее описанных методов. К р а е в ы е у с л о в и я. Остановимся на записи разностной схемы в нерегулярных узлах (на границе илн вблизи нее). В этих узлах для записи разностных уравнений необходимо привлекать краевые условия. Например, в разностных схемах (16) и (18) для уравнения теплопроводности и, = йи„„ нерегулярными являются граничные узлы и =О, и =уо'. Лля первой краевой И-1 и ы задачи и(0, Т)=111(~), и(а, ()=ро(() зот лппгоксимлция уравнение (18) можно написать при и =0: 1 /г —,, (Уо — Уо) =9 (У-т 3Уо+Уд. Заменим в левом краевом условии (30) производную симметричной разностью: 1 2а (Ут Исключая из последних двух уравнений фиктивную точку, получим разностный аналог краевого условия: 1 Ь Ь(У1 УО) Р1 (~ ) + 2 (Уо Уо)' (33) Заметим, что это уравнение содержит только одно значение с нового слоя у„ т. е.
оно явное. Метод уменьшения невязки менее нагляден, но более универсален. Выразим и(х„т) при помощи формулы Тейлора: и (х„1) = и (х„1)+Й и, (х„()+- Ь'и„„+... На основании краевого условия (30) положим и,(х„1)=р,(1), а из уравнения теплопроводности (15а) найдем и, -ийй. Под- ставляя эти величины в формулу Тейлора, получим и(х1, 1) = и(хо, 1)+Ьрг Е+,— и1+". А и(х) — )(х) =О, хеп б, (34а) Ри(х) — р(х)=0, хенГ. (34б) Введем в области О+Г сетку с шагом и, состоящую из множества внутренних (регулярных) узлов ю„и множества граничных (нерегулярных) узлов у„. Заменим задачу (34) в регулярных узлах разностным аналогом уравнения.(34а): А ьуь (х) — ~рь (х) = О, х ев ы„, (35а) Заменяя здесь и~ — (у,— у,)!т, снова придем к краевому условию (33).
В последнем способе можно учесть большее число членов ряда Тейлора и получить краевые условия не только нормальной, но и повышенной точности. 5. Аппроксимация и ее порядок. Пусть имеется область 6 переменных х = (х„..., х ) с границей Г и поставлена корректная задача для некоторого уравнения с граничными условиями: УРАВНЕНИЯ В ЧАСТНЫХ ПРОИЗВОДНЫХ )ГЛ. 1Х (37) (ь ) 112 Ци (х)Ц =(((р,(х) и„'(х)+р, (х) и'(х)) с)х~ а р„(х) ) О, р, (х) ~ О. (38в) Употребляются и другие нормы. Напомним (см.
главу !), что из локальной близости функций следует их среднеквадратичная близость; поэтому Ц Ц, называют более сильной, чем Ц Цс, Нетрудно проверить, что энергетическая норма (Ззв) сильней Ц Ц . Выбор той или иной нормы в конкретной задаче определяется двумя соображениями. Желательно, чтобы разностное решение у а в нерегулярных узлах — разностным аналогом краевых усло- вий (34б): )с»у» (х) )(» (х) =О х е— = 7» (35б) (индексом й отмечены величины, определенные только на сетке; мы будем опускать его там, где это не вызовет недоразумений).
Близость разностной схемы (35) к исходной задаче (34) будем определять по величине невязки: )р» (х) = (Аи — )) — (А»и — гр„), х ен ю», У» (х) =(Рль — р) — Я»и — )(») х ен у». Определение. Разностная схема (35) аппраксимирует за- дачу (34), если Ц 1Р Це»-ьО, ЦУ Цх„-ьО пРи )т-+О; (38) аппроксимация имеет р-й порядок, если )) РЦ»=О()гл), Цу Ц,„=О()ьл) при )г — О. Обсудим вопрос о выборе норм в этом определении.
Функции и (х), ((х), р(х) определены обычно на отрезке а(х =(т или во всех точках некоторой области пространства большего числа измерений, Для них можно ввести такие нормы, как чебышевская (локальная): Ци(х) Цс= шах !и(х)), (38а) а(к(Ь или гильбертова (среднеквадратичная): гь ч1/2 ))и (х)))с, = р (х) и'(х) с(х~, р (х) ) О (38б) (выражения написаны для одномерного случая). Часто исполь. зуют связанные с оператором А энергетические нормы, напоми- нающие формулы для полной энергии колебательной системы, например: % г1 АППРОКСИМАЦИЯ ) Аг 'цг !!у!1,= шах !у„!, !!у!1г, =( ~', р„у,"-/г„~ 0(л~м го= 1 (39) В выборе разностных аналогов норм существует некоторый произвол. Например, сумму в (39) козино брать по л=б, 1, ..., У вЂ” 1, что соответствует выбору другой квадратурнсй формулы для интеграла (386). Этим пользуются, определяя сеточные нормы так, чгобы облегчить доказательство сходнмости. Как надо понимать Ь -ь ОУ Лля равномерной сетки это не требует пояснений.
На неравномерных сетках рассматривают совокупность шагов Ьа как некоторую сеточную функцию и вводят какую-либо норму шага, например: ум — ~ 'пг Ь=!!Ь„!1,=снах Ь„или Ь=!!Ьз!!л =~ ~„"Ь,',/ л —. о Эту норму считают «величиной шага» в определениях аппроксимации, порядка аппроксимации н т. д. Если невязку оценивают в 11 11„то аппроксимацию называют локальлог". Для уравнений с достаточно гладкими решениями наличие локальной аппроксимации и ее порядок легко проверяются; в таких задачах нередко ограничиваются установлением локальной аппроксимации. Однако наиболее сильные результаты по сходимости разностных схем связаны с использованием более слабых норм для невязки (но сильных норм для решения).
Замечание 1. Факт наличия или отсутствия аппроксимации и порядок аппроксимации зависят не только от операторов А и А„"), но также от классов, к которым принадлежат и(х), )(х), и от выбора норм. Чем сильнее норма или чем шире классы функций, тем ниже, вообще говоря, порядок аппроксимации *) Операторы краевых условий для краткости обычно будем опускать. было близко к точному решению в возможно более сильной норме; например, в задачах на разрушение конструкций малость деформаций в 11 11с, не гарантирует сохранения конструкции, а малость в 11 1!с — гарантирует. С другой стороны, чем слабее 11 11„, тем легче построить сходящуюся в этой норме разностную схему и исследовать ее. Заметим, что функции и(х), )(х), р(х) принадлежат, вообще говоря, разным классам.
Например, если и(х) есть четырежды дифференцируемая функция и А = (г(г)с(хг), то ) (х) является дважды дифференцируемой функцией. Поэтому каждую из этих функций можно оценивать в своей норме: 11 11в, 11.!1), 11 11„. ФуНКцИИ уь, Чга, уа ОПрЕдЕЛЕНЫ ТОЛЬКО На СЕТКЕ, ПОЭтОМу дпя них надо ввести сеточные нормы 11 11 „, 11.1!ть, 11 1!х„. Их вводят так, чтобы при й — 0 они переходили в выбранные 11 11„, 11 1!), 11 11„. За разностные аналоги чебышевской и гильбертовон норм можно принять соответственно зш (ГЛ. (Х УРАВНЕНИЯ В ЧАСТНЫХ ПРОИЗВОДНЫХ (последнее видно по замечанию 2 к п. 3, если оценивать невязку в (( ((с).
Замечание 2. Как правило, решение и(х) исходной задачи (34) неизвестно, так что использовать его для получения невязки затруднительно. В этом случае берут достаточно широкий класс У функций и(х), которому и (х) заведомо принадлежит (обычно это класс функций, непрерывных вместе с достаточным числом своих производных).
Если на всех функциях класса У имеется аппроксимация порядка р: ((Ап(х) — Аьп(х)+(ра(х) — )(х)ЦА=О(йл) при й- О, то, очевидно, аппроксимация на решении и (х) имеет порядок ие ниже р. В подобных случаях аппроксимация на решении и (х) может иметь порядок выше р. В замечании 3 к п. 3 мы видели, что для уравнения из=йи„„явная разностная схема (18) при т = = из((бя) имеет в классе сколь угодно гладких функций п(х) аппроксимацию 0(6я), а на решений — более высокого порядка (четвертого, как нетрудно проверить). Сл у ч а и м н о г и х и е р е м е н ц ы х имеет некоторые особенности. Определение аппроксимации остается в основном прежним; надо только требовать стремления к нулю шагов по всем переменным. Порядок аппроксимации может быть разный по разным переменным.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
