
- •1.1 Идеальные жидкость и газ.
- •1.1.1 Напряжения в идеальной жидкости (газе).
- •1.1.2 Уравнение движения Эйлера для идеальной жидкости (газа).
- •1.1.4 Полная система уравнений движения баротропной идеальной жидкости (газа).
- •1.1.5 Полная система уравнений движения идеальной несжимаемой жидкости.
- •1.1.6 Начальные и граничные условия.
- •1.2 Первые интегралы уравнений движения идеальной жидкости.
- •1.2.1 Интеграл Бернулли при установившемся течении баротропной жидкости.
- •1.3 Потенциальное течение идеальной жидкости (газа).
- •1.3.1 Теорема Лагранжа о сохраняемости потенциального течения во времени.
- •1.3.2 Общие нелинейные уравнения потенциального течения идеальной баротропной жидкости (газа).
- •1.3.3 Возмущенное движение сжимаемой жидкости (газа).
- •1.3.4 Потенциальное течение однородной несжимаемой жидкости.
- •1.3.5 Установившееся плоское движение сжимаемой жидкости (газа).
- •1.3.6 Неустановившееся одномерное течение.
- •1.4 Примеры решения простейших задач гидро и газо динамики с применением интеграла Бернулли.
- •1.4.1 Задача об истечении жидкости из бочки.
- •1.4.2 Задача о переливе жидкости через плотину.
- •1.4.3 О сути подъёмной силы при обтекании различных профилей.
- •1.4.4 Течение в трубе переменного сечения. Эффект пульверизатора.
- •1.4.5 Кавитация.
- •1.5 Задача о движении сферы в бесконечной идеальной несжимаемой жидкости.
- •1.5.1 Задача об обтекании сферы идеальной несжимаемой жидкостью.
- •1.5.2 Парадокс Даламбера.
- •1.6 Контрольные вопросы и задачи по теме "Идеальная жидкость (газ)".

или даже отрицательным. В этих местах в жидкости образуются пузырьки, заполненные парами жидкости или же газами, растворенными в ней. Такое явление называется кавитацией. Подробнее с явлением кавитации можно познакомиться в книге [2]
1.5Задача о движении сферы в бесконечной идеальной несжимаемой жидкости.
Пусть абсолютно твердая сфера радиуса a движется поступательно со скоростью V (t) относительно неподвижной системы декартовых координат в неограниченном объёме идеальной несжимаемой жидкости. Массовые силы положим равными нулю, т.е.
~
F = 0.
Обозначим неподвижные координаты через x1; x2; x3. Свяжем с телом декартову систему координат x1; x2; x3, которая движется поступательно вместе со сферой в направлении оси x1. Оси подвижных координат направим по осям неподвижной системы координат, а центр расположим в центре сферы.
p
Обозначим через r = xixi = x21 + x22 + x23 расстояние от начала координат до некоторой точки подвижного пространства, а через # угол между лучем Or и осью Ox1, причем cos # = x1=r (рис. 5).
Рис. 5: Движение сферы в идеальной жидкости
При движении сферы частицы жидкости отодвигаются сферой в направлении нормали к поверхности тела со скоростью V cos #. Этим действием сфера приводит в движение всю жидкость. Пусть ~v(x) = [v1(x); v2(x); v3(x)] вектор скорости жидкости относительно подвижных координат. На границе сферы проекция вектора скорости ~v жидкости на нормаль к сфере будет равна скорости точки сферы в направлении нормали
|
|
|
vn = v~n r=a = vr r=a = V cos # |
(1.66) |
Граничное условие (1.66) называется условием непроницаемости (поток жидкости внутрь тела равен нулю). Условие непроницаемости диктует характер взаимодействия тела с жидкостью. При 0 6 # < =2 жидкость продвигается вперёд и в стороны. На экваторе, т.е. при # = =2 радиальная компонента вектора скорости vr равна нулю, тангенциальная компонента вектора скорости v# = V (t). При =2 < # 6 жидкость увлекается за сферой и, тем самым, обеспечивается безотрывность обтекания.
Другое условие на движение жидкости состоит в том, что на бесконечности жидкость покоится, т.е.
lim vr ! 0 (1.67)
r!1
14

В нашем случае массовые силы отсутствуют, несжимаемая однородная жидкость является баротропной жидкостью с функцией давления P = p=%0 + const:, движение жидкости начинается из состояния покоя, т.е. при t = 0 вихри в жидкости отсутствуют. Следовательно, все условия теоремы Лагранжа о сохраняемости потенциальности течения выполнены.
Течение жидкости вызванное движением в ней жесткой сферы будет потенциальным, т.е. ~v = grad ', при этом
vr = v~n = grad ' ~n = |
@' |
ni = |
@' |
; |
v# = v~ = |
1 @' |
; |
(1.68) |
|||
|
|
|
|
|
|
||||||
|
@r |
|
r @# |
||||||||
|
@xi |
|
|
|
|
|
Где ~n и ~ единичные векторы нормали и касательной к поверхности сферы. Потенциальное течение несжимаемой однородной жидкости описывается уравнением Лапласа (1.54)
' = 0 |
(1.69) |
и интегралом Коши-Лагранжа в подвижных координатах (1.42), в котором следует положить U = 0, P = p=%0 + const:, ~vпер = [V (t); 0; 0], а функцию f(t) приписать к потенциалу. В результате интеграл Коши-Лагранжа примет вид:
@' |
|
@' |
|
1 |
|
p |
|
(1.70) |
|||
|
|
|
|
V (t) + |
|
|
jgrad 'j2 |
+ |
|
+ const: = 0 |
|
@t |
@x1 |
2 |
%0 |
Константу выбираем из условия на бесконечности для давления p(x; t) ! p1. По-
r!1
сле этого из интеграла Коши-Лагранжа найдём распределение давления
@' |
@' |
1 |
(1.71) |
||||
p(x; t) = p1 %0 |
|
+ %0 |
|
V (t) %0 |
|
jgrad 'j2 |
|
@t |
@x1 |
2 |
Искомый гармонический потенциал должен удовлетворять условию непроницаемо-
сти (1.66) и условию на бесконечности, т.е. |
|
|
|
|
|
|
|
|||||
@' |
r=a = |
V |
cos |
# ; |
lim ' |
! |
0 ; |
lim |
@' |
! |
0 |
(1.72) |
@r |
|
|
r!1 |
|
r!1 |
@r |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Задача (1.69), (1.72) является классической задачей Неймана для внешности шара. Её решением, убывающем на бесконечности, является гармоническая функция
' = A@x1 |
r |
= A r3 = A |
r2 |
(1.73) |
|
@ |
1 |
|
x1 |
cos # |
|
Константу A найдём из условия непроницаемости
|
@r |
r=a |
= @r A |
r2 |
r=a |
|||||
@' |
|
@ |
|
|
cos # |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
Отсюда |
|
|
|
|
|
|
|
V a3 |
|
|
|
|
|
|
|
|
|
A = |
|
||
|
|
|
|
|
|
|
|
|||
и потенциал скорости принимает вид: |
2 |
|
||||||||
|
|
|
||||||||
|
|
|
'(x; t) = |
a3V (t) cos # |
||||||
|
|
|
|
2r2 |
= |
2A cos # |
= V cos # |
(1.74) |
|
a3 |
||||
|
|
|
= |
a3V (t)x1 |
(1.75) |
2r3 |
15

Далее найдем компоненты вектора скорости |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
a3V (t) |
|
x1 |
|
|
a3V (t) |
1 + |
3x2 |
|
|
a3V (t) |
|
|
|||||||||||||||
v1 |
= ';1 |
= |
|
|
|
|
;1 |
= |
|
|
|
|
1 |
|
= |
|
|
|
|
1 + 3 cos2 |
# |
|||||||||
2 |
|
r3 |
|
2r3 |
|
r2 |
|
2r3 |
||||||||||||||||||||||
|
|
|
a3V (t)x |
1 |
1 |
|
|
|
|
3a3V (t) x1x2 |
|
|
|
3a3V (t) |
|
|
||||||||||||||
v2 |
= ';2 |
= |
|
|
|
|
|
;2 |
= |
|
|
|
|
|
= |
|
|
|
|
|
sin # cos # |
|
||||||||
2 |
|
|
|
r3 |
2 |
|
|
r5 |
|
|
|
2r3 |
|
|
||||||||||||||||
|
|
|
a3V (t)x |
1 |
1 |
|
|
|
|
3a3V (t) x1x3 |
|
|
|
3a3V (t) x3 |
|
|
||||||||||||||
v3 |
= ';3 |
= |
|
|
|
|
|
;3 |
= |
|
|
|
|
|
= |
|
|
|
|
|
|
cos # |
|
|||||||
2 |
|
|
|
r3 |
2 |
|
|
r5 |
|
|
|
2 |
|
|
r4 |
|
(1.76)
(1.77)
(1.78)
Найдём радиальную и тангенциальную составляющие вектора скорости в плоскости x1x2
@' |
a3V (t) cos # |
|
|
@' |
|
a3V (t) sin # |
|
) j~vj = |
a3V (t) |
|||||||||
vr(r; #) = |
|
= |
|
|
|
; |
v#(r; #) = |
|
= |
|
|
|
|
|
|
|||
@r |
|
r3 |
|
r@# |
|
r3 |
|
r3 |
||||||||||
На поверхности сферы, т.е. при r = a |
|
|
|
|
|
|
|
|
(1.79) |
|||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
j~v(a)j = q |
|
= V (t) |
(1.80) |
||||||||||
vr(a; #) = V (t) cos # ; |
|
v#(a; #) = V (t) sin # |
) |
|
vr2 + v#2 |
Отсюда видно, что линии тока начинаются и заканчиваются на поверхности сферы, поскольку
vr(a; #) = V (t) cos # ; |
v#(a; #) = V (t) sin # |
||||
vr(a; #) = V (t) cos # ; |
v#(a; #) = V (t) sin # |
||||
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 6: Скорость жидкости в симметричных точках и линии тока (помечены красным цветом)
Давление на поверхности сферы определяется по формуле (1.71), которая для случая постоянной скорости движения сферы принимает вид:
@' |
|
% v2 a |
) |
|
% V 2 |
3 cos2 # 2 |
|
|
p(a; #) = p1 + %0 |
|
(a)V |
0 ( |
= p1 + |
0 |
(1.81) |
||
@x1 |
2 |
|
2 |
16

1.5.1Задача об обтекании сферы идеальной несжимаемой жидкостью.
Пусть теперь жидкость обрекает неподвижную сферу. Систему координат оставим прежней Ox1x2x3 с центром O в центре сферы. Жидкость течет из бесконечности вдоль оси x1 против её направления, т.е. скорость жидкости на бесконечности равна V . Такое течение, также как и в предыдущем случае является потенциальным течением.
Потенциал течения обозначим через . Постановка задачи об обтекании неподвижной сферы жидкостью состоит из уравнения Лапласа
|
|
= 0 ; |
|
|
|
(1.82) |
|
условия непроницаемости |
|
|
|
|
|
|
|
wn = w~n r=a = wr |
r=a = |
@ |
r=a |
= 0 |
(1.83) |
||
@r |
|||||||
|
|
|
|
|
|
|
|
и условия на бесконечности |
|
|
|
~ |
|
|
(1.84) |
|
rlim grad |
|
|
|
|
||
|
|
! V k1 |
|
||||
|
!1 |
|
|
|
|
|
|
Рассмотрим случай, когда сфера отсутствует, а течение есть. Тогда такое течение описывается потенциалом
'0 = V x1 |
(1.85) |
Очевидно, что уравнение Лапласа и условие на бесконечности для этого потенциала удовлетворяются.
Если теперь на течение, вызванное движением сферы вместе с осями координат наложить движение жидкости противоположное движению сферы, то сфера остановится2, а среда будет обтекать сферу. Потенциал такого суммарного течения является гармонической функцией и он равен
|
|
= ' + '0 = |
a3V (t)x |
1 |
V (t)x1 = V r + |
|
a3 |
|
cos # |
|
(1.86) |
||||||||||||||||
|
|
2r3 |
|
2r2 |
|
||||||||||||||||||||||
Вычислим компоненты вектора скорости |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
@ |
= V 1 |
|
a3 |
cos # ; |
|
|
|
|
|
@' |
|
|
|
1 + |
a3 |
sin # |
(1.87) |
||||||||||
wr(r; #) = |
|
|
v#(r; #) = |
|
= V |
|
|||||||||||||||||||||
@r |
r3 |
r@# |
2r3 |
||||||||||||||||||||||||
Проверяем условие непроницаемости |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
@ |
|
|
|
|
|
|
|
a3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
@r |
|
r=a = V cos # 1 r3 r=a = 0 |
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Условие на бесконечности |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
a3V x1 |
|
|
~ |
|
|
|
|
~ |
|
|
|
|
||||
|
|
r!1 |
|
|
= r!1 |
|
|
|
2r3 |
|
|
V k |
1 = |
|
V k |
1 |
|
|
|
||||||||
|
|
lim grad |
|
|
lim grad |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Следовательно, все условия задачи (1.82)–(1.84) выполняются, а потенциал |
, пред- |
ставленный формулой (1.86) является решением задачи (1.82)–(1.84) об обтекании жесткой недеформируемой сферы идеальной несжимаемой жидкостью.
2по крайней мере очевидно, что скорость лобовой точки сферы станет равной нулю
17