Математика / Курсовая работа: Дослідження збіжності рішень для диференціальних рівнянь у частинних похідних, отриманих методом сіток

Курсовая работа: Дослідження збіжності рішень для диференціальних рівнянь у частинних похідних, отриманих методом сіток

Метод сіток

Ідея методу сіток відома давно ще за часів Ейлера. Однак, практичне використання цього методу наштовхувалося на серйозні труднощі, тому що одержання з його допомогою досить точного розв’язку крайової задачі звичайно приводило до колосальних систем алгебраїчних рівнянь, на розв’язок яких при ручному розрахунку були потрібні роки. Положення різко змінилося з появою електронних обчислювальних машин. Метод сіток допускає зручну реалізацію на ЕОМ, тому що застосування його зазвичай зводиться до масової повторюваності однорідних циклів. В даний час метод сіток є одним з найбільш ефективних методів розв’язку лінійних диференціальних рівнянь. Метод сіток (інакше метод скінчених різниць) для наближеного розв’язку крайових задач двовимірних диференціальних рівнянь полягає в наступному:

У плоскій області G, у якій розшукується розв’язок, будується сіткова область Gh, що складається з однакових осередків і наближає дану область G;

Задане диференціальне рівняння заміняється у вузлах побудованої сітки відповідним скінчено-різницевим рівнянням;

На підставі граничних умов установлюються значення шуканого розв’язку в граничних вузлах області Gh

Розв’язавши отриману систему скінчено-різницевих рівнянь, ми знайдемо значення шуканої функції у вузлах сітки, тобто будемо мати чисельний розв’язок нашої задачі. Вибір сіткової області здійснюється в залежності від конкретної задачі, але у всіх випадках контур Гh сіткової області Gh варто вибирати так, щоб він якнайкраще апроксимував контур Г заданої області G. Сітка будується таким чином, щоб вузли (xi,yi) сітки Sh або належали області G, або відстояли від її границі Г на відстань меншому, ніж h. Точки (вузли) сітки Sh називаються сусідніми, якщо вони розміщені один від одного в напрямку осі Ох або осі Оу на відстань, що дорівнює кроку сітки h. Вузол Ah сітки Sh називається внутрішнім, якщо він належить області G, а всі чотири сусідніх з ним вузла – множині Sh; інакше він називається граничним. Граничний вузол сітки Sh називається вузлом першого роду, якщо він має сусідній внутрішній вузол цієї сітки, інакше граничний вузол називається вузлом другого роду. Внутрішні вузли і граничні вузли першого роду сітки Sh називаються розрахунковими точками. Граничні вузли другого роду не входять в обчислення і можуть бути вилучені із сітки.

На перший погляд процедура застосування методу сіток, що складається з трьох етапів, може здатися простою і легко реалізованою. Однак насправді це не так. Через велику розмаїтість типів і розмірів сіток, видів рівнянь у часткових похідних, граничних і початкових умов, можливих кінцево-різницевих апроксимацій цих рівнянь і методів їхнього розв’язку, чисельне розв’язку рівнянь у часткових похідних вимагає модифікацій алгоритму при розгляді кожного конкретного приклада.

Стійкість скінчено-різницевої схеми для розв’язку рівнянь параболічного типу (рівняння теплопровідності)

Як приклад рівняння параболічного типу розглянемо рівняння теплопровідності:

де u=u(x,t) – температура, t – час, - довжина стрижня.

Для простоти покладемо, а=1.

Початкові і крайові умови:

и..

При використанні скінчено-різницевої схеми для розв’язку крайової задачі виникає питання про стійкість такої схеми. Під цим розуміють наступне: скінчено-різницева схема називається стійкою, якщо малі похибки в процесі розв’язку загасають або у всякому разі залишаються малими при необмеженому збільшенні номера поточного шару.

Для рівняння

(3)