Метод скінченних різниць

Методи скінченних різниць базуються на дискретизації функції на сітці

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

$\frac{\partial u}{\partial t}=Au+f,(1)$

$(x,t)\in \mathrm{\Omega },$,

$lu=g,(x,t)\in \partial D\times T;u(x,0)=u_0(x)$

де $\mathrm{\Omega }=D\times T$ — циліндрова область інтеграції $t\in T=[0,t{}_0]$ — межа області $\mathrm{\Omega }$, $D$ — її основа, $u$— шукана вектор-функція $f$ і $g$— задані вектор-функції, $x$ — просторовий векторний аргумент, $A$ і $l$ — оператори (не обов'язково обмежені), то найпростіша схема інтеграції початкового рівняння має вигляд:

$\frac{u^{n+1}-u^n}{\tau }={\mathrm{\Lambda }}_1{u}^{n+1}+{\mathrm{\Lambda }}_0{u}^n+F(2)$,

$(x,t)\in {\mathrm{\Omega }}_h$

Тут $u^n$ — сіткова функція, що є розв'язком різницевого рівняння, ${\mathrm{\Lambda }}_1,{\mathrm{\Lambda }}_0,\lambda$— різницеві оператори, залежні від параметрів $\tau$, $h$ сітки, $t=n\tau$, ${\mathrm{\Omega }}_h$ — сіткова область, що апроксимує деяким чином область $\mathrm{\Omega }$, $\partial D\times T$— її границя $F$і $G$— сіткові функції, що апроксимують функції $f$і $g$ відповідно. Окремим випадком схеми (2) є схема з вагами, коли ${\mathrm{\Lambda }}_1=\alpha \mathrm{\Lambda }$ ${\mathrm{\Lambda }}_0=(1-\alpha )\mathrm{\Lambda }$, $\alpha$ — ваговий коефіцієнт.

Схема (2) називається двошаровою, оскільки вона зв'язує між собою значення $u^n$ і $u^{n+1}$ різницевого розв'язку на двох тимчасових шарах $n\tau$; можливі також і багатошарові схеми. Якщо оператор $E-\tau {\mathrm{\Lambda }}_1$, де $E$ — одиничний оператор, оборотний, то схема (2) може бути представлена у вигляді, що розв'язується $u^{n+1}=\sigma u^n+\mathrm{\Phi }(3)$ де оператор $\sigma$ називаєстья оператором кроку різницевої схеми і враховує крайові умови, а $\mathrm{\Phi }$ — функція, залежна від $F$і $G$. Кажуть, що оператор $\mathrm{\Lambda }\left(\tau \right)$, залежний від параметра $\tau$, апроксимує (приблизно) оператор $A$, якщо $\Vert [\mathrm{\Lambda }\left(\tau \right)-A]u\Vert =\epsilon n\left(\tau \right)\to 0$ при $\tau \to 0$. Тут $u\in U$ — деяке еталонне сімейство функцій, на якому перевіряється апроксимація (наприклад, сімейство достатньо гладких функцій).

Схема (2) називається коректною, або стійкою, якщо ${\Vert \sigma \Vert }_B=1+O(\tau )$, де ${\Vert \sigma \Vert }_B$ означає норму оператора $\sigma$ в деякому банаховому просторі. $B$ (див. простір абстрактний у функциональному аналізі), яке може залежати від $h$.

Схема (2) апроксимує рівняння (1), якщо ${\mathrm{\Lambda }}_1+{\mathrm{\Lambda }}_0\propto l$ і $\lambda \propto l$. Для лінійних систем рівнянь встановлені теореми збіжності, що стверджують, що збіжність різницевого рішення до розв'язку початкового рівняння виходить з апроксимації і коректності (стійкості) різницевої схеми.

Якщо властивості апроксимації, стійкості і збіжності мають місце лише при деякому співвідношенні між параметрами сітки $\tau$, $h$, де $h=h(\tau )$, то їх називають умовними. Якщо ж ці властивості справедливі при будь-якому співвідношенні між $\tau$ і $h$, то їх називають абсолютними.

Схема (2) називається явною, якщо ${\mathrm{\Lambda }}_1\equiv 0$, і неявною, якщо ${\mathrm{\Lambda }}_1\ne 0$.

Схеми, що абсолютно збіжні, існують лише в класі неявних схем. Як правило, при відповідному виборі параметрів схеми (наприклад вагових коефіцієнтів) неявні схеми є абсолютно стійкими, вони допускають скільки завгодно великий крок $\tau$. Але перетворення оператора $E-\tau {\mathrm{\Lambda }}_1$ ускладнює алгоритм. У разі одновимірних задач неявні схеми реалізують методом факторизації; вони є достатньо економічними. Для багатовимірних задач неявні економічні схеми одержують за допомогою дробових кроків методу, який зводить багатовимірні задачі до послідовності одномірних або простіших задач.

Для розв'язання стаціонарних задач застосовують метод стаціонування (встановлення), в якому стаціонарний розв'язок розглядається як межа нестаціонарного розв'язку із стаціонарним (або що встановлюються) краєвими умовами. Відповідно до цього стаціонарну задачу вирішують ітераційним методом, аналогічним різницевому методу інтеграції (2). На відміну від нестаціонарного випадку, оператор $\sigma$ для ітераційного процесу повинен бути сильно стійкий, тобто повинен задовольняти умові ${\Vert \sigma \Vert }_B=1+\epsilon (h)$, $\epsilon (h)>0$. При розв'язку нелінійних задач, особливо в механіці суцільного середовища, застосовують комбінації схем інтеграції з ітераційними методами (т.з. ітерації по нелінійності).

Дробний порядок

Нехай $f(x)$ є скалярною функцією із $\mathrm{D}\mathrm{o}\mathrm{m}(f)\subset R,$ а $f_k=f(x_k),$ де $x_k=a+kh$ (де $h=\mathrm{c}\mathrm{o}\mathrm{n}\mathrm{s}\mathrm{t}>0,k=0,\pm 1,...$) - точки із $\mathrm{D}\mathrm{o}\mathrm{m}(f).$

Перша різниця у точці $x_k$ є прирощенням

$\mathrm{\Delta }{f}_k=f_{k+1}-f_k.$

Кінцева (права) різниця порядку $m$ функції у точці $x_k$ визначається рекурентною формулою

${\mathrm{\Delta }}^m{f}_k={\mathrm{\Delta }}^{m-1}{f}_{k+1}-{\mathrm{\Delta }}^{m-1}{f}_k=\sum _{j=0}^m(-1{)}^j\left(\begin{array}{c}m\\ j\end{array}\right)f_{k+m-j},m=1,2,...;$

${\mathrm{\Delta }}^1{f}_k=\mathrm{\Delta }{f}_k,{\mathrm{\Delta }}^0{f}_k=f_k.$

Скінченні різниці ${\mathrm{\Delta }}^m{f}_k,m=1,2,...$ є випередженням. Запізненням ж є

$\mathrm{\nabla }{f}_k=f_k-f_{k-1}=\mathrm{\Delta }{f}_{k-1},{\mathrm{\nabla }}^m{f}_k={\mathrm{\Delta }}^{m}f-k-m,m=1,2,...$

Запізнювані різниці називаються також лівими різницями. Центральною різницею є

${\mathrm{\Delta }}^m{f}_{k-\frac{m}{2}}.$

Кінцевою (правою) різницею порядку $a\in ℂ$ у точці $x_k$є

${\mathrm{\Delta }}^{a}f-k=\sum _{j=0}^m(-1{)}^j\left(\begin{array}{c}a\\ j\end{array}\right)f_{j+m-j},m=1,2,...$

Відношення $\frac{{\mathrm{\Delta }}^a{f}_k}{h^a}(h^a=\exp (a\ln k)),k=0,\pm 1,\pm 2,...$ є різницевою похідною (дробного) порядку $a$ у точці $x_k$.

Нехай $x\in \mathrm{D}\mathrm{o}\mathrm{m}(f)$ - довільне число, більше від $a$,

$x_k=a+kh,h={\delta }_n=(x-a)/n,n=1,2,...$

При виборі точок $x_k$ ($k=0,\pm 1,\pm 2,...$) й числа $h$, які є вузлами й кроком, мають місце вирази

${\mathrm{\nabla }}^m{f}_n={\mathrm{\nabla }}^m{f}_{n-m}=\sum _{j=0}^m(-1{)}^j\left(\begin{array}{c}m\\ j\end{array}\right)f_{n-j}=\sum _{j=0}^m(-1{)}^j\left(\begin{array}{c}m\\ j\end{array}\right)f(x-j{\delta }_n),m=1,2,...$

Звідси

${\mathrm{\nabla }}^{m}f(x)=\sum _{j=0}^m(-1{)}^j\left(\begin{array}{c}m\\ j\end{array}\right)(x-j{\delta }_n),m=,1,2,...$

Відношення

${\mathrm{\Delta }}_n^{a}f(x)=\frac{\underset{n}{\overset{a}{\mathrm{\nabla }}}f(x)}{{\delta }_x^a},n=1,2,...$

є різницевою похідною функції у точці $x.$

Див. також

• Метод скінченних елементів
• Метод рухливих клітинних автоматів

Література

• Шаблон:ЕК, Н. Н. Япепко.
• Годунов С. К., Рябенький В. С. Введение в теорию разностных схем. М., 1962 [библиогр. с. 272—274];
• Яненко Н. Н. Методы дробных шагов решения многомерных задач математической физики. Новосибирск, 1967 [библиогр. с. 189—193];
• Самарский А. А. Введение в теорию разностных схем. М., 1971 [библиогр. с. 538—550];
• Рихтмайер Р.  Д., Мортон К. Разностные методы решения краевых задач. Пер. с англ. М., 1972 [библиогр. с. 381—413].