Численное моделирование дозвукового отрывного обтекания осесимметричных тел
7
ним вихревых многоугольников. Скорость в контрольных точках на
поверхности тела находим как предельное значение градиента потен-
циала двойного слоя
1
1
( )
Г ( )
( ).
2
N
k k
k
V C V
w C V C
Давление
p
определялось из интеграла Бернулли:
2 2
1
.
2
p p
V V
(12)
В соответствии с (12) безразмерные коэффициенты давления
2
2(
) / (
)
p
c
p p
V
рассчитываются по формуле
2
1 ( / ) .
p
c
V V
Результаты численных расчетов.
По предложенной методике
произведено моделирование обтекания осесимметричных тел раз-
личной формы. На данном этапе исследований было признано целе-
сообразным осуществлять построение эквивалентного тела наиболее
простой конфигурации. Поэтому поправка на толщину вытеснения
пограничного слоя не производилась и передний участок поверхно-
сти эквивалентного тела совмещался с поверхностью обтекаемого
тела. За линией отрыва потока, которая считалась совпадающей с ли-
нией донного среза, хвостовой участок поверхности эквивалентного
тела, моделирующий влияние спутного следа, представлял собой
продолжение поверхности обтекаемого тела. Длина указанного
участка поверхности эквивалентного тела обозначалась как
сл
.
l
Предполагалось, что нормаль к поверхности эквивалентного тела
при переходе через линию отрыва потока меняется неперывным об-
разом. Это условие равносильно требованию об отсутствии излома
поверхности эквивалентного тела в окрестности донного среза.
Дополнительно предполагалось, что если спереди к линии отрыва
потока примыкает сужающийся или цилиндрический участок экви-
валентного тела, то и хвостовой участок поверхности эквивалентного
тела должен быть сужающимся или по крайней мере нерасширяю-
щимся.
Угол атаки принимался равным нулю, что дало возможность счи-
тать эквивалентное тело осесимметричным.
