ВУЗ:
Составители:
Рубрика:
При аппроксимации оператора переноса B в уравнении (4) считаем
положительными обе компоненты скорости , . При иных знаках
изменения очевидны.
0≥u 0≥v
Полагаем функцию
ϕ
постоянной в ячейке (n,m) в том смысле, что
именно это значение переносится через правую и верхнюю границы ячей-
ки, то есть
()
2/12/1
2/1
2
1
−+
+
+⋅=
mmm
n
m
n
hhuJ
ϕ
()
2/12/1
2/1
2
1
−+
+
+⋅=
nn
m
n
m
n
hhvJ
ϕ
Соответственно
()()
()
(
)
1
2/12/11
2/12/1
22
−
−+−
−+
−++−+=
m
n
m
nnn
m
n
m
n
mm
nn
hh
v
hh
u
B
ϕϕϕϕϕ
3.3.10. Метод решения дискретного уравнения переноса
Обращение начинаем с левого нижнего угла сетки, так как полагаем
, . Если знаки
u
и
v
другие, то начинаем с других углов.
0>u 0>v
При нахождении уже известны (вычислены на предыдущих ша-
гах) и , а значит, можно определить и , . находим
из уравнения
m
n
ϕ
m
n
1−
ϕ
1−m
n
ϕ
m
n
J
2/1−
2/1−m
n
J
m
n
ϕ
() (
2/1
2/11
2
−
−
++=
⎥
⎦
⎤
⎢
⎣
⎡
++
m
n
m
n
m
n
m
nn
mm
n
JJPvhuhS
ϕ
τ
)
.
Для контроля за точностью решения будем вычислять норму невязки
уравнения (3). Итерационный процесс прекращается, когда норма невязки
уменьшилась в требуемое количество раз, например, в 1000 или 10000.
82
Страницы
- « первая
- ‹ предыдущая
- …
- 80
- 81
- 82
- 83
- 84
- …
- следующая ›
- последняя »
