ВУЗ:
Составители:
Рубрика:
29
Алгебраическое уравнение (63) совместно с разностным уравнением не-
разрывности решается методом прогонки. При этом решение уравнения (63)
представляется в форме
.
,, jjijji
BuAu (64)
Прогоночные коэффициенты
jj
BA , определяются соотношением:
.;
1
1
1 jjj
jjj
j
jjj
j
j
bAa
Bad
B
bAa
c
A
(65)
Расчет структуры пограничного слоя осуществляется последовательно, на-
чиная с сечения i = 1, при этом все параметры потока в предыдущем сечении i =
0 известны из граничных условий для х = 0. В каждом расчетном сечении сна-
чала определяются прогоночные коэффициенты
jj
BA , во всех узлах, начиная с
j = 1. Эту операцию называют прямой прогонкой. На поверхности стенки при j
= 0 значения прогоночных коэффициентов находятся из граничных условий для
y = 0. Так, при расчете пограничного слоя с граничными условиями для скоро-
сти u получаем .0;0
00
BA Далее во всех точках j, начиная с точки с наи-
большим номером, обратной прогонкой с использованием уравнений (65) нахо-
дят значения скорости u и других параметров.
По результатам расчета структуры пограничного слоя в интересующих се-
чениях потока определяются коэффициент трения
f
c
;
0
2
y
f
y
u
u
c
(66)
толщина вытеснения
*
, толщина потери импульса
**
пограничного слоя и
соответствующие числа Рейнольдса Re
*
, Re
**
;Re;1
*
*
0
*
u
dy
u
u
(67)
;Re;1
**
**
0
**
u
dy
u
u
u
u
(68)
безразмерные скорости
φ и координаты η
.;
*
*
yv
v
u
(69)
Более подробно численный метод интегрирования дифференциальных
уравнений пограничного слоя изложен в работе [2].
Программа расчета
Рассмотренный численный метод расчета пограничного слоя реализован
в виде кода «Программа для расчета эффективности пленочного охлаждения
поверхности, обтекаемой высокоскоростным потоком с воздействиями» на ал-
горитмическом языке Microsoft Visual C++. Программа зарегистрирована в рее-
стре программ для ЭВМ (Свидетельство о государственной регистрации про-
грамм для ЭВМ №2011613213). Окно компилятора Microsoft Visual C++ пред-
ставлено на рис 5.