Обтекание колеблющегося крыла потоком идеальной несжимаемой жидкости
5
Сетка координатных линий свободной вихревой поверхности в
текущий расчетный момент времени ,
q
t
соответствующая введенной
на несущей поверхности сетке и разбивающая правую половину сво-
бодной вихревой поверхности
1
1
2
0
( ( )
0, 0
)
q
b t
b
b d
на элемен-
ты (рис. 2)
1
1 1 2 2 2
1
1
,
b b b b b b
,
(10)
определяется следующими формулами:
1 1
0
2
2
2
2
2
2
1
2
1
1
( )
,
1 ,
0, 1, ..., ,
1, 2, ...;
( ), 0
,
,
1, ...,
1, 0
1;
( ),
,
,
2 ,
2,
3,...,
1,
,
1, ..., 1;
1
0.
j
i
b b t
ct
c
t
q q
b g a
b n
j j
n
a
b g a n b d d m n m n i
n n
m n i m m
a
(11)
При этом
1
2
1
2
1
2
1
2
( )
, ( )
,
a
a
g a
m n g a
а шаг по времени
t
считается постоянной величиной.
Рис. 2. Расчетная схема
Сетка координатных линий свободной вихревой поверхности соот-
ветствует сетке координатных линий несущей поверхности в том смыс-
ле, что координатные линии
2
const
b
свободной вихревой поверхно-
