= У AaU, Aa = ka-, fea = const, (67)

X = {Xx, X, Xp}.

Аппроксимируем это уравнение симметричной неявной схемой, которую назовем исходной:

\{У-У)=\ 2 Ла(Р+г/), (68)

а= 1

где Аа - разностные операторы, аппроксимирующие А с погрешностью 0{Ы); обычно для них используют формулы (57), соответствующие г = 2. Благодаря симметричной форме исходная схема имеет погрешность

\ а /

Однако эта схема неэкономична, потому что не найдено хорошего алгоритма вычисления у.

Наряду с исходной схемой построим локально-одномерную схему. Введем промежуточные слои и на каждом слое в правой части (68) вместо возьмем рАц; в левой части поставим

шаг т/р. Обозначим решение на промежуточных шагах через Ша (а=1, 2.....р). Тогда функции Wa будут удовлетворять разно-

каждом слое составим схему типа (59), неявную по одному направлению и явную по остальным. Во-первых, такая схема несимметрична и имеет аппроксимацию лишь 0{т). Во-вторых, она оказывается условно устойчивой при kx<.h и, тем самым, неэкономичной.

Экономичные многомерные разностные схемы можно строить локально-одномерным методом, также используя промежуточные слои. Эти схемы имеют лишь суммарную аппроксимацию. На промежуточных слоях они вообще не аппроксимируют исходное дифференциальное уравнение; но погрешности аппроксимации промежуточных слоев при суммировании гасят друг друга так, что на целом слое аппроксимация есть. При этом разностное решение следует сравнивать с точным только на целых слоях, не придавая промежуточным слоям самостоятельного смысла.

Рассмотрим многомерное параболическое уравнение (66); для простоты ограничимся случаем анизотропной теплопроводности с постоянными коэффициентами:

Операторы попарно перестановочны; операторы Ад получаются тоже попарно перестановочными. Последовательно применяя (70) и используя перестановочноегь операторов, нетрудно установить следующее равенство:

П (£ - { тЛ„) = {П (£ +1 тЛ„) wi.

\a=l

а= 1

Раскроем произведения операторов и положим Wi = у, Wp= у. Пренебрегая членами высокого порядка по т, получим запись схемы (69):

la а, р j

la а. Р j

стным уравнениям и начальным условиям следующего вида:

\{wa-Wa)==Ya{wa + Wa), а = 1,2,...,р; (69а) Wi = y, W2 - W1, и0з=щ,..., Wp=Wpi, у = Wp. (696)

Поскольку Лц - одномерные операторы, то каждая Wa является решением одномерной разностной схемы; поэтому схему (69) называют локально-одномерной. Исследуем ее.

Устойчивость. Каждое уравнение (69а) является одномерной неявной симметричной схемой типа схемы (6) при 0 = 1/2-Последняя схема безусловно устойчива, так что ошибка начальных данных не возрастает ни на одном промежуточном слое. Следовательно, схема (69) также безусловно устойчива и позволяет вести расчет с шагом xrh.

Вычисление разностного решения несложно. Каждое уравнение (69а) решается одномерной прогонкой. По тем же причинам, что и в случае схемы (6), прогонка устойчива, а разностное решение у существует и единственно.

Для нахождения решения на новом целом слое надо выполнить прогонки по всем р направлениям. Это требует Юр действий на каждую точку сетки независимо от величин шагов ha-Таким образом, локально-одномерная схема экономична.

Аппроксимацию исследуем, сравнивая схему (69) с исходной. Для этого перепишем (69) в следующем виде:

-l-xAa)wa=(E+хАа]ш>а, = Ша-1- (70)

а= 1

Следовательно, погрешность аппроксимации симметричной локально-одномерной схемы (69) на целых слоях есть

\ а=1

Заметим, что для получения погрешности аппроксимации О (т*) в граничных условиях надо к естественным граничным условиям добавлять поправки типа (646).

Сходимость схемы (69), как следует из сказанного выше, является безусловной с погрешностью O/t-j-aV

Замечание. В некоторых случаях расщепление многомерной задачи на последовательность одномерных бывает точным. Например, многомерный перенос по характеристике точно эквивалентен последовательности,одномерных переносов по проекциям этой характеристики на координатные плоскости.

Остановимся на некоторых усложнениях задачи (67), Переменные коэффициенты ka{x, t) приводят к тому, что операторы Лц становятся неперестановочными и Лц -тоже. В этом случае погрешность аппроксимации схемы (69) возрастает

до 0/т+2а]- Поэтому ДЛЯ уравнения

ди dt

= 2 A,u + f{x, t), А,и = -[К{х, t)\ (72)

нередко ограничиваются чисто неявной локально-одномерной схемой

(ша-Ша) = ЛаШа + фа, tOa = tt-a-i, 2 4»=/ (73а)

01=1

с естественными граничными условиями

[ша-1г(х. t]y = Q\ (736)

(-1/)=2«( + )- 2 ЛаЛр(-у) + 0(т=). (71)

а= I а, р

На решениях с непрерывными пятыми производными двойная сумма в (71) есть 0(т:), поэтому (71) отличается от исходной схемы только членами О(т). Но погрешность аппроксимации исходной схемы равна