В.Н. Тимофеев
8
( ) ( );
k
k
a w C n C
( ).
b V n C
При проведении вычислений методом дискретных вихрей для
замкнутых поверхностей матрица, составленная из коэффициентов
vk
a
, получилась вырожденной, поэтому систему (12) решали с ис-
пользованием регуляризирующей переменной [6].
После нахождения неизвестных циркуляций
k
Г
скорость потока
в точках, находящихся вне поверхности тела, определяли из соотно-
шения (11). В контрольных точках, лежащих на поверхности тела,
учитывалось, что скорость как градиент потенциала двойного слоя
испытывает разрыв:
1
2
1
2
( )
( )
( )
( )
( ).
v
v
v
v
v
g C
g C
V C
C
C
(13)
Здесь частные производные
1
( )
v
g C
и
2
( )
g C
следует вычислять
по направлениям двух взаимно ортогональных ортов
1
( )
C
и
2
( )
C
,
лежащих в касательной плоскости, проходящей через контрольную
точку
.
v
C
На каждой панели
v
поверхностная плотность
( )
g C
была постоянной и равной
Г
, поэтому в формуле (13) частные про-
изводные определяли численным дифференцированием по значениям
циркуляций
v
-го и соседних с ним вихревых многоугольников. После
нахождения векторов
( )
v
V C
скорость потока в контрольных точках
определяли как предельное значение градиента потенциала двойного
слоя:
1
1
( )
( )
( ),
2
N
k k
k
V C V Г w C V C
а статическое давление
p
находили из интеграла Бернулли:
2 2
1 (
),
2
p p
V V
где
V
модуль вектора скорости в рассматриваемой точке течения;
V
модуль вектора скорости набегающего потока.
Безразмерные коэффициенты давления
2
2(
) /
p
c
p p
V
в
соответствии с интегралом Бернулли рассчитывали по формуле
2
1 ( / )
p
c
V V
. (14)
