Методы моделирования теплоэнергетических процессов. Ртищева А.С. - 19 стр.

UptoLike

Составители: 

18
Коэффициент турбулентного переноса
æ определяется зависимостью,
предложенной Ковальноговым Н.Н., учитывающей влияние фактора
динамической нестационарности, продольного градиента давления и кривизны
обтекаемой поверхности
0
0
21,41
4,9
1
0,4
=
=
+
=
y
x
0
x
0
xx
y
x
y
w
u
x
p
r
w
ww
p
y
w
ρ
ρ
τ
æ
.
(1.25)
Здесь индекс характеризует параметры течения в анализируемом
сечении за пределами пограничного слоя.
Анализ различных разностных схем для решения системы
дифференциальных уравнений пограничного слоя (1.17) – (1.22) показывает,
что наиболее удобными являются неявные шеститочечные схемы. Для
составления такой схемы на координатной плоскости x, y выбираются основная
и две вспомогательные сетки.
Координаты узлов
основной сетки определяются соотношениями:
xix
i
=
;
yjy
i
=
.
(1.26)
Координаты узлов вспомогательных сеток определяются соотношениями:
xix
i
= ; yjy
j
+=
+
2
1
2
1
;
xix
i
+=
+
2
1
2
1
; yjy
i
=
.
(1.27)
В выражениях (1.26) – (1.27) x, yшаги сетки вдоль координатных осей
x и y соответственно.
Тогда конечно-разностная аппроксимация дифференциального уравнения
энергии (1.18) будет иметь следующий вид, где значения
ji,
T
,
1n
ji,
x
w
и т.д. в
момент времени
τ
, а
ji,
T ,
1n
ji,
x
w
и т.д. в момент времени
τ
+
τ
: