Математическая энциклопедия - параболического типа уравнение

численные методы решенияметоды решения уравнений параболич. типа на основе вычислительных алгоритмов. Для решения П. т. у. часто применяются приближенные численные методы, рассчитанные на использование быстродействующих ЭВМ. Наиболее универсальным является метод сеток (конечно разностный метод).

Ниже рассмотрен метод сеток на примере уравнения теплопроводности

(1) с краевыми условиями 1-го рода

Вводится равномерная сетка узлов (xi, tn).

и обозначения

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

(2)

где s числовой параметр и сеточная аппроксимация функции f( х, t), напр. ).

Система уравнений (2) решается по слоям, т. е. для каждого n=0, 1, ... по известным значениям ,

находятся новые значения .

Если s=0 (явная схема), то выражаются явным образом через . Если же (неявные схемы), то относительно , i=l, 2, . . ., N-1, возникает система уравнений, имеющая трехдиагональную матрицу. Эта система уравнений решается прогонки методом. Недостатком явной схемы является сильное ограничение на шаг т, возникающее из условия устойчивости, а именно . Неявные схемы при абсолютно устойчивы, т. е. устойчивы при любых шагах hи т. Известны и другие разностные схемы для уравнения (1) (см. [1], [2]).

Если схема (2) устойчива и аппроксимирует f(x, t), то при h,решение разностной задачи сходится к решению u( х i, tn).исходной задачи (см. [1]). Порядок точности зависит от параметра а. Так, симметричная схема имеет второй порядок точности по t и по h, т. е. для каждого пвыполняется оценка

с константой М, не зависящей от hи t. При s=0,5-h²/(12t) и специальном выборе схема (2) имеет точность О(t²+h⁴). Для остальных s точность O(t+h²).

При решении П. т. у. на больших отрезках временя существенное значение имеет асимптотич. устойчивость разностной схемы. Решение уравнения (1) с f(x, t)=0, и 1 (г)=ц 2 (г)=0 ведет себя при как е ^-l1t, l1=p²/l². Этим свойством обладает не всякая разностная схема, аппроксимирующая уравнение (1). Напр., симметричная схема (s=0,5) асимптотически устойчива при условии . Построены схемы 2-го порядка точности, асимптотически устойчивые при любых t, h, однако они не входят в семейство (2) (см. [3]).

Краевые условия 3-го рода

аппроксимируются с порядком 0(h²).разностными уравнениями

Для уравнения теплопроводности в цилиндрич. и сферич. координатах

вводятся сетки (ri,tn), где tn=nt,n=0,1, . . .,

и рассматриваются разностные уравнения (см. [4])

где

Для решения П. т. у. с переменными коэффициентами применяются однородные консервативные разностные схемы, выражающие на сетке законы сохранения, присущие исходному дифференциальному уравнению (см. [1], [5]). При построении разностных схем для П. т. у. с переменными коэффициентами применяются методы: баланса, вариационный, конечных элементов. Напр., для уравнения

используются разностные схемы с весами

где

Доказана сходимость и получены оценки погрешности однородных консервативных разностных схем для П. т. у. как в случае непрерывных, так и в случае разрывных коэффициентов (см. [1]).

Для решения систем уравнений параболич. типа, содержащих одну пространственную переменную, применяются те же разностные схемы, что и для одного уравнения. Вектор у ^п+1 решения на новом слое находится матричной факторизации методом.

При решении квазилинейного П. т. у. используются неявные абсолютно устойчивые разностные схемы. Решение yⁿ⁺¹ разностного уравнения находится итерационным методом с использованием прогонки. Напр., для уравнения

(3) применяют чисто неявную разностную схему

к-рую решают с помощью итерационного метода (s номер итерации):

Решение на каждой итерации находится методом прогонки. При решении нелинейного П. т. у. нашли также применение схемы, основанные на идее Рунге Кутта метода. Так, для уравнения (3) используется двухэтапный метод