Совершенствование вычислительных алгоритмов для решения уравнений Навье-Стокса на структурированных сетках - page 9

1. Задание начальных значений
u
(0)
,
v
(0)
и
p
(0)
.
2. Задание
n
p
= 0
.
3. Получение начального приближения к решению уравнений
Навье–Стокса из решения системы
 
M
h
(
u
) =
p
(0)
ij
p
(0)
i
1
j
h
x
,
M
h
(
v
) =
p
(0)
ij
p
(0)
ij
1
h
y
.
(20)
4. Задание
n
p
=
n
p
+ 1
.
5. Коррекция давления
p
(
n
p
)
ij
.
6. Решение уравнений движения (18) и (19)
 
M
h
(
u
) =
p
(
n
p
)
ij
p
(
n
p
)
i
1
j
h
x
,
M
h
(
v
) =
p
(
n
p
)
ij
p
(
n
p
)
ij
1
h
y
.
7. Проверка сходимости, если необходимо — возврат к п. 4.
Далее совокупность действий п. 4–7 будет называться
итерацией
по давлению
, а
n
p
— соответствующим счетчиком итераций.
Модифицированный сегрегированный алгоритм (МСА), с включе-
нием вспомогательной задачи, представим в виде следующей после-
довательности действий.
1. Задание начальных значений
u
(0)
,
v
(0)
,
(
p
xy
)
(0)
,
(
p
x
)
(0)
и
(
p
y
)
(0)
.
2. Задание
n
p
= 0
.
3. Получение начального приближения (
P
-итерации)
 
 
M
h
(
u
) =
p
x
i
p
x
i
1
h
x
(
p
xy
)
(0)
ij
(
p
xy
)
(0)
i
1
j
h
x
,
1
Z
0
u
(
x, y
)
dy
=
. . .
;
 
M
h
(
v
) =
p
y
j
p
y
j
1
h
y
(
p
xy
)
(0)
ij
(
p
xy
)
(0)
ij
1
h
y
,
1
Z
0
v
(
x, y
)
dx
=
. . .
(21)
4. Задание
n
p
=
n
p
+ 1
;
5. Коррекция давления
(
p
xy
)
(
n
p
)
ij
.
86
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2008. № 2
1,2,3,4,5,6,7,8 10,11,12,13,14,15,16,17
Powered by FlippingBook