Пространственная локализация тепловых возмущений в нелинейном процессе теплопроводности - page 5

y
(
k
)
n
y
n
τ
=
a
2
2
y
(
k
)
n
1
2
y
(
k
)
n
+
y
(
k
)
n
+1
h
2
+
+
a
2
2
y
n
1
2
y
n
+
y
n
+1
h
2
f
(
k
1)
n
, k
= 1
,
2
, . . .
;
y
(0)
n
=
y
n
.
Опишем итерационную процедуру нахождения решения на вре-
менн´ом слое
t
m
+1
.
Шаг 1
(
k
= 1
). За нулевое приближение
f
(0)
n
принимаем значение
f
n
с предыдущего временн´ого слоя
t
m
и определяем первое прибли-
жение
y
(1)
n
решения
y
n
.
Шаг 2
(
k
= 2
). Используя значение
y
(1)
n
, вычисляем
f
(1)
n
=
f y
(1)
n
,
а затем — второе приближение
y
(2)
n
решения
y
n
и т.д.
Шаг
k
. Величины
y
(
k
)
n
,
n
= 1
, N
1
, могут быть найдены методом
прогонки из решения трехдиагональной системы линейных алгебраи-
ческих уравнений
2
a
2
+
h
2
τ
y
(
k
)
1
+
a
2
y
(
k
)
2
=
= 2
a
2
h
2
τ
y
1
a
2
(
y
0
+
y
2
) + 2
h
2
f
(
k
1)
1
a
2
y
(
k
)
0
;
a
2
y
(
k
)
n
1
2
a
2
+
h
2
τ
y
(
k
)
n
+
a
2
y
(
k
)
n
+1
= 2
a
2
h
2
τ
y
n
a
2
(
y
n
1
+
y
n
+1
) + 2
h
2
f
(
k
1)
n
, n
= 2
, N
2;
a
2
y
(
k
)
N
2
2
a
2
+
h
2
τ
y
(
k
)
N
1
=
= 2
a
2
h
2
τ
y
N
1
a
2
(
y
N
2
+
y
N
) + 2
h
2
f
(
k
1)
N
1
a
2
y
(
k
)
N
.
(9)
Элементы матрицы системы (9) удовлетворяют условиям диаго-
нального преобладания, которые являются достаточными условиями
того, что в формулах для определения прогоночных коэффициентов
ни один из знаменателей не обращается в нуль, а обратная прогонка
устойчива по входным данным [5].
Для выбора значений шагов
h
и
τ
может быть применена методика,
изложенная в работе [6].
Пример численного расчета.
Приведем результаты численного
расчета для задачи вида
∂u
∂t
=
2
u
∂x
2
48
u
1
/
2
,
0
< x <
1
, t >
0;
u
(
x,
0) = 0
,
0
< x <
1;
u
(0
, t
) = 1
, u
(1
, t
) = 0
, t >
0
.
(10)
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2013. № 4
31
1,2,3,4 6,7
Powered by FlippingBook