Елементарні поняття з теорії різницевих схем

Перехід від математичної моделі неперервної диференціальної задачі до відповідної їй дискретної моделі та її чисельної реалізації останньої на ЕОМ супроводиться виникненням цілого ряду похибок:

а) похибки апроксимації диференціального рівняння або системи різницевим;

б) похибки апроксимації крайових умов;

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

г) похибки округлень даної ЕОМ.

 Якщо підставити точний розв'язок U деякої храмової задачі в її дискретний аналог (різницеву схему), то цей розв'язок, взагалі кажучи, не буде задовільняти цим скінченно - різницевим співвідношенням. Виникає так звана нев'язка (похибка апроксимації) як для самого диференціального рівняння (чи рівнянь), так і для додаткових (крайових) умов задачі. Близькість різницевої схеми до вихідної задачі якраз визначається по величи­ні цієї нев'язки. Кажуть, що різницева схема апроксимує крайову задачу (диференціальну), якщо величини цих нев'язок прямують до нуля в деякій вибраній нормі при подрібненні сітки (h®0); апроксимація має к -тий порядок, у випадку рівності нев'язок О (h^k ), де О(h^k) =Аh^k, А - деяка стала. Отже, властивість апроксимації означає близькість різницевого оператора до диференціального. Різницеву схему, яка апроксимує крайову задачу для диференціального рівняння в частинних похідних, називають ще узгодженою.

Однак, з цього ще не випливає близькість розв'язків диференціальної крайової задачі і її скінченне різницевого аналогу (різницевої схеми).

Результатом дискретизації задачі, як відмічалось, є система різни­цевих (алгебраїчних) рівнянь. Різницева схема називається стійкою, якщо розв'язок системи різницевих рівнянь неперервно залежить від вхідних даних і ця залежність рівномірна відносно кроку (кроків) сітки h . Ін­шими словами, схема є стійкою, якщо малим змінам її параметрів відпові­дають малі зміни її розв'язків

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

Різницева схема називається коректною, якщо її розв'язок:
1) існує;                                                                                           

2) єдиний при будь-яких вхідних даних;

3) схема стійка.

Основним питанням теорії різницевих схем є питання про їх збіжність .

Під збіжністю розуміють прямування розв'язку різницевої схеми до точного розв'язку вихідної задачі при подрібненні сітки.

Тобто кажуть, що різницевий розв’язок збігається до точного роз­в'язку диференціальної задачі, якщо їх різниця прямує до нуля в деякій нормі при прямуванні кроків сітки до нуля.

Коротко кажучи, збіжність означав близькість різницевого розв'язку до істинного розв'язку вихідної задачі.

При цьому виникає питання про швидкість збіжності різницевої схеми.

Кажуть, що різницева схема збігається з швидкістю О(h^l ), як­що різниця міх її розв'язком і точним розв'язком вихідної задачі є ве­личина порядку О(h^l ).

Доведено, що для того, щоб різницева схема була збіжною, тобто її розв'язок збігався до істинного, необхідно і достатньо, щоб вона апроксимувала вихідну задачу і була стійкою.

Коротко кажуть: "апроксимація і стійкість забезпечують збіжність".

Слід відмітити, що теорія різницевих схем (як результат чисель­ного розв'язування крайових задач методом скінчених різниць) найбільш повно розвинута в роботах вітчизняних математиків: А.А.Самарського, С.К.Годунова, М.М.Яненко, їх учнів та ін.

С.К.Годуновим та його учнями розвинена теорія консервативних різ­ницевих схем тобто таких схем, для яких виконуються різницеві аналоги інтегральних законів збереження, наслідкам яких є основні рівняння математичної фізики.

В основному склалося три основні способи побудови РС на заданому шаблоні :

1.Метод різницевої апроксимації (той, що ми розлядали до сих пір) ;

2.Метод невизначених коефіцієнтів;

3.Інтегро-інтерпаляційний метод (метод балансу);

 

1. Метод різницевої апроксимації полягає в тому, що кожна похідна, що входить в диференціальне рівняння і крайові умови змінюється певним різницевим виразом (включаючи лише вузли шаблону). Саме так були отримані всі розглянуті нами вище різницеві схеми. Цей метод досить простий і додаткових пояснень не потребує.

Метод різницевої апроксимації дозволяє легко складати РС першого чи другого порядку апроксимації на прямокутній сітці для рівнянь з неперервними і достатньо гладкими коефіцієнтами. Однак цей метод важко чи неможливо застосувати в більш складних випадках , а саме для рівнянь з розривними коефіцієнтами, на прямокутних сітках , для рівнянь високого порядку , на нерівномірних сітках і т. д.

⇐ Предыдущая12345678910Следующая ⇒

Поделиться: