
- •Глава 1. Предмет и методы механики сплош- ной среды
- •1.1 Предмет механики сплошной среды 3
- •Глава 2. Кинематика деформируемой сре- ды
- •Глава 6. Математическая модель гидроди- намики
- •Глава 7. Магнитогидродинамическая модель течений проводящего газа и моде- ли солнечного ветра
- •Глава 8. Модель стационарных течений иде- альной несжимаемой жидкости
- •2 Поток не со-
2 Поток не со-
осей,
здает ни поперечной, ни продольной силы, т.е. тело не испытывает сопротивления.
Таково поперечное обтекание кругового цилиндра идеальной несжи- маемой жидкостью при отсутствии циркуляции. Вязкость вносит прин- ципиальные изменения в эту картину. Наличие трения в погранич- ном слое приводит к появлению продольной силы трения. Кроме то-
го, в областях DB и D0B, где давление вдоль контура растет (см.
рис.
8.8.1), происходит
отрыв
пограничного
слоя
(рис.
9.8.1(слева)),
который
из-за
неустойчивости
течения
превращается
в
дорожку
Кар-
мана
(рис.
9.8.1(справа)).
D
y
y
A
B
A
B
x
x
D'
Á Â
Рис. 9.8.1. Безциркуляционное обтекание кругового цилиндра при наличии вязкости
Cx
100
40
20
10
4
2
|
|
|||||||
|
||||||||
|
||||||||
|
||||||||
|
||||||||
|
||||||||
|
||||||||
|
||||||||
|
||||||||
|
||||||||
|
|
|
|
|
|
|
|
|
1.0
0.4
0.2
0.1 -1 Re
10
100
10
102
10 3
104
10 5
106
Рис. 10.8.1. Сопротивление цилиндра как функция числа Рейнольдса
Давление в зоне отрыва падает, что создает основную часть про- дольной силы силы лобового сопротивления X . На рис. 10.8.1 по- казана зависимость коэффициента сопротивления цилиндра от чис- ла Рейнольдса. При достижении критического числа Рейнольдса за
счет турбулизации течения в пограничном слое профиль скорости становится становится более полным, что приводит к смещению точ- ки отрыва по направлению к B и следовательно к уменьшению силы сопротивления.
Рассмотрим теперь, как изменится обтекание, если на симметрич- ное течение наложить циркуляцию против часовой стрелки. Величи- на циркуляции пока остается неопределенной. В этом случае ком- плексный потенциал, по (8.25), имеет вид (V = 0)
µ
w = U z +
a2 ¶ iΓ
z − 2π
· ln z, (8.28)
где U скорость на бесконечности, направленная вдоль оси Ox, и Γ
циркуляция скорости при обходе цилиндра в положительном на- правлении.
Чтобы получить картину течения, определим положение крити- ческих точек точек торможения. Они определяются из квадратного уравнения
v¯кр =
µ dw ¶
U a2 iΓ 1
2
= − − ·
z
= 0,
и следовательно
dz кр
iΓ
zкр = ±
2
z
кр
s µ Γ
−
2π кр
¶2
+ a2. (8.29)
4πU
4πU
При |Γ| > 4πU a обе точки лежат на мнимой оси: одна
|z(1)
Ї Ї sµ ¶2
Ї Ї
(1)
Γ
Γ
кр | = Ї Ї +− a2 (a < |zкр | < ∞)
Ї 4πU Ї
вне цилиндра, другая
4πU
|z(2)
Ї Ї sµ ¶2
Γ
Ї Ї − a
(0 < |z(2)| < a)
Γ
кр | = Ї Ї2 кр
Ї 4πU Ї −
4πU
внутри цилиндра (см. C и C 0 на рис. 11.8.1).
|
Линии тока при этом имеют вид, изображенный на том же рисун- ке. Если уменьшать |Γ , то точка C опускается, а C 0 поднимается,
| | U
так что при Γ = 4π a они сливаются на поверхности цилиндра,z(1)
(2)
(1)
(2)
кр = zкр = ia. Наконец, при |Γ| < 4πU a, имеем |zкр | = |zкр | = a
обе точки расположены на поверхности цилиндра. Картина линий тока в этом случае показана на рис. 12.8.1.
y
C
C’
A
B
x
Рис. 11.8.1. Циркуляционное обтекание кругового цилиндра
y
A B x
Рис. 12.8.1. Циркуляционное обтекание кругового цилиндра
Вычислим продольную X и поперечную Y силы, отнесенные к единице длины цилиндра. Воспользуемся интегралом Бернулли, учи- тывая, что давление торможения p0 = const. Получим (см. рис. 4.8.2)
I
X = −
C
I
p cos(n, x) ds =
d
2π
ρa Z
2
0
2π
ρa Z
2
|v|
2
cos θ dθ,
(8.30)
так как
Y = −
C
I
p cos(n, y) ds =
d 2
0
0
I|v|
sin θ dθ,
n, x) ds =
C C
p sin(n, x) ds = 0.
d
Здесь ~n внешняя нормаль, ds элемент дуги контура, так что
n, x) = cos θ, cos(d) = sin θ, ds = a · dθ. Согласно (8.27) и (8.19),
cos( d
n, y
на контуре
так что
ρa
2
|v|
2π
Z
|
= |v¯ 2
µ Γ ¶2
= 2U sin θ − 2πa ,
2π
U Γ Z
X = 4U 22
0
sin2 θ cos θ dθ 2
πa
0
sin θ cos θ dθ+
Γ2
+ 4π2a2
2π
2π Z
−
cos θ dθ
0
= 0,
2π
(8.31)
ρa Z
Y = 4U 22
0
sin3 θ dθ − 2
U Γ Z
πa
0
sin2 θ, dθ+
Γ2
+ 4π2a2
2π Z
sin θ dθ
0
= −ρvΓ.
Таким образом при поперечном обтекании цилиндра с циркуля- цией (при V = 0)
X = 0, Y = −ρU Γ (8.32) Чтобы определить направление поперечной силы Y, нужно век-
тор скорости на бесконечности повернуть на 90◦ навстречу цирку-
ляции, т.е. возникающая при циркуляционном обтекании цилиндра сила сопротивления перпендикулярна вектору скорости на бесконеч- ности и повернута в направлении, обратном направлению циркуля- ции.