05 семестр / Книги и методические указания / Панаиотти С.С. Лекции по гидромеханике / Панаиотти С.С. Лекции по гидромеханике. Часть 1 - Учебное пособие. - Калуга КФ МГТУ им. Н.Э. Баумана, 2003
.pdf
При таком движении частицы в своем мгновенном вращении поворачиваются вокруг касательных к линиям тока, как показано на рис. 1.23. Направляющие косинусы векторов щ иV
cos (ω, x ) = |
ω |
x |
|
, cos (V, x ) = |
V |
x |
|
|
||||||
|
|
|
|
|
|
; |
||||||||
|
|
|
|
|
|
|
|
|||||||
|
|
ω |
|
V |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ω |
x |
|
|
Vy |
|
|||||||
cos (ω, y ) = |
|
|
, cos (V, y ) = |
|
|
|
|
; |
||||||
|
|
|
|
|
|
|
|
|||||||
|
|
ω |
|
|
V |
|
||||||||
cos (ω, z ) = |
|
ω |
x |
, cos (V, z ) = |
V |
|
|
|
|
|||||
|
|
|
|
z |
. |
|||||||||
|
|
|
|
|
|
|||||||||
|
|
ω |
V |
|
|
|
|
|||||||
Так как cos(ω, x) = cos(V, x), cos(ω, y)
то
|
|
|
ω |
x |
= |
ωy |
= |
ω |
z |
. |
|
|
|
|
|
|
|
|
|||
|
|
|
Vx |
Vy |
|
Vz |
||||
Переписывая (1.58) в виде |
|
|
|
|
|
|
|
|||
|
V 2 |
|
|
|
|
|
|
|
|
|
d U − P − |
|
= |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
cos(ω, z) = cos(V, z),
(1.60)
=−2 (Vy ωz − Vz ωy )dx + (Vz ωx −Vxωz )dy + (Vxωy − Vy ωx )dz
иимея в виду (1.60), получаем
U − P − |
V |
2 |
= C (для всех точек потока). |
(1.61) |
|
|
|
||||
2 |
|||||
|
|
|
|||
На рис. 1.24 показана вихревая пелена за крылом. Аналогичная пелена появляется за лопастью осевой гидромашины.
Рис. 1.24. Вихревая пелена за крылом по данным [10]
41
Вихревые нити «подвешены» к задней кромке крыла и при движении крыла сбегают с него, как нить с веретена. Движение частиц жидкости, составляющих пелену — винтовое. Поэтому
для любых точек пелены справедливо соотношение (1.61). Интеграл Лагранжа. Теперь интегрируем уравнения (1.47)
для неустановившегося потенциального движения жидкости, при котором ωx = ωy = ωz = 0 . В этом случае уравнения (1.47)
принимают вид:
∂ |
V |
2 |
∂V |
|
∂ |
|
|
|
V 2 ∂Vy |
|
|||||||
|
U − P − |
|
|
|
= |
x |
, |
|
|
U − P − |
|
= |
|
, |
|||
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
∂t |
|
|
|
|
|
2 |
|
∂t |
|
||
∂x |
|
2 |
|
∂y |
|
|
|
|
|||||||||
|
|
|
|
∂ |
|
|
V 2 |
|
∂V |
|
|
|
|
||||
|
|
|
|
|
|
U − P − |
|
|
= |
z |
. |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
2 |
|
∂t |
|
|
|
|
||
|
|
|
∂z |
|
|
|
|
|
|
|
|||||||
Частные производные от потенциала скорости ϕ по координа-
там равны: |
∂ϕ |
= V , |
∂ϕ |
= V |
|
, |
∂ϕ |
= V , |
|
|
y |
|
|||||
|
∂y |
x |
∂y |
|
∂z |
z |
||
|
|
|
|
|
||||
Поскольку значение частной производной не зависит от порядка дифференцирования, то
∂V |
∂ |
∂ϕ |
|
∂ |
∂ϕ |
|
∂Vy |
|
∂ |
|
∂ϕ |
|
∂V |
z |
|
∂ |
∂ϕ |
||||||||
x |
= |
|
|
|
|
= |
|
|
|
|
, |
|
= |
|
|
|
|
, |
|
= |
|
|
|
. |
|
|
|
|
|
|
∂t |
|
|
∂t |
|
|
|
||||||||||||||
∂t |
∂t |
∂x |
|
∂x |
∂t |
|
|
∂y |
∂t |
|
|
|
∂z |
∂t |
|||||||||||
Подставляя эти производные в предыдущие уравнения, найдем:
∂ |
V |
2 |
|
|
∂ϕ |
∂ |
V |
2 |
|
∂ϕ |
||||
|
U − P − |
|
|
− |
|
= 0 , |
|
U − P − |
|
|
− |
|
= 0 , |
|
|
|
|
|
|
|
|
|
|||||||
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
∂x |
|
∂t |
∂y |
2 ∂t |
||||||||||
∂ |
V |
2 |
|
|
∂ϕ |
|
|
|
|
|
|
|
||
|
U − P − |
|
|
− |
|
|
= 0 . |
|
|
|
|
|
(1.62) |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
∂z |
|
|
∂t |
|
|
|
|
|
|
|
||||
В общем случае неустановившегося движения четырехчлен
U − P − V 2 − ∂ϕ является функцией координат и времени. Одна-
2∂t
ко уравнения (1.62) показывают, что в данном случае он не зависит от координат и поэтому в данный момент времени
U − P − |
V 2 |
− |
∂ϕ |
= C (t) (для всех точек потока). (1.63) |
|
2 |
∂t |
||||
|
|
|
42
Здесь C(t) — функция интегрирования, которая зависит от времени t и не зависит от координат x, y, z. В данный фиксированный момент времени эта функция постоянна.
Выясним физический смысл производной ∂ϕ в уравне-
∂t
нии (1.63). Выберем в потоке жидкости две произвольные точки 1 и 2 и соединим их дугой длиной S (рис. 1.25). Дифференциал потенциала скорости равен дифференциалу циркуляции
d ϕ = |
∂ϕ |
dx + |
∂ϕ |
dy + |
∂ϕ |
dz = V dx + V |
dy + V dz = d = V dS . |
||
|
|
|
|||||||
|
∂x |
|
∂y |
∂z |
x |
y |
z |
s |
|
|
|
|
|
|
|
||||
Рис. 1.25. К вычислению производной потенциала скорости по времени
Интегрируя это уравнение в данный момент времени от 1 до 2, получим:
2
ϕ2 − ϕ1 = ∫Vs dS .
1
Производная по времени от этого выражения
∂ϕ∂t 2
|
∂ϕ |
|
∂ |
2 |
|
|
− |
|
|
= |
|
∫ |
Vs dS . |
|
∂t |
|||||
|
∂t |
|
|
|||
|
1 |
|
|
1 |
|
|
Так как S и t являются независимыми переменными, то выполним вначале дифференцирование под знаком интеграла, а затем интегрирование.
43
|
∂ϕ |
|
|
∂ϕ |
|
2 ∂V |
||||
|
|
|
|
− |
|
|
= |
∫ |
s |
dS . |
|
|
|
|
|||||||
|
∂t |
2 |
|
∂t |
|
∂t |
||||
|
|
|
|
1 |
|
1 |
|
|
||
Производная ∂Vs — ускорение в данной точке или проекция
∂t
инерционной силы, отнесённой к единице массы жидкости. По-
этому ∂Vs dS является элементарной работой инерционной си-
∂t
2 ∂V
лы на бесконечно малом пути dS , а ∫ s dS — работа при пе-
1 ∂t
ремещении вдоль кривой |
S от точки 1 до точки 2 . Следова- |
|||||||||
|
|
∂ϕ |
|
∂ϕ |
|
|
∂ϕ |
|||
тельно, |
|
|
|
− |
|
|
— |
это работа, а |
|
— энергия единицы |
|
|
|
||||||||
|
|
∂t 2 |
|
∂t 1 |
|
|
∂t |
|||
массы жидкости в данной точке потока.
Уравнение (1.63) можно истолковать следующим образом. При неустановившемся потенциальном течении полная энергия единицы массы жидкости в любых двух точках потока в данный момент времени одинакова.
В частном случае движения несжимаемой жидкости в поле
силы тяжести U = −gz + const , а P = p + const . Поэтому уравне-
с
ния (1.56), (1.59) и (1.61) примут вид: |
|
|
|
|
|
|
|
|||||||
|
p |
V 2 |
|
|
|
p |
2 |
|
V 2 |
|
||||
gz + |
1 |
+ |
|
1 |
= gz |
|
+ |
|
+ |
|
2 |
. |
(1.64) |
|
|
|
|
|
|
|
|||||||||
1 |
с |
2 |
|
2 |
|
с |
2 |
|
|
|||||
Причем, в случаях потенциального и винтового движений точки 1 и 2 — любые две точки потока. В случае вихревого движения точки 1 и 2 выбираются на данной линии тока. При тех же условиях интеграл Лагранжа ( 1.63) запишется так
|
|
p V 2 |
|
∂ϕ |
|
p |
2 |
V 2 |
|
∂ϕ |
|
|||||
gz1 |
+ |
1 |
+ |
1 |
+ |
|
|
= gz2 + |
|
+ |
2 |
+ |
|
|
. (1.65) |
|
с |
2 |
|
с |
2 |
|
|||||||||||
|
|
|
|
∂t |
|
|
|
∂t |
2 |
|||||||
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
||
В этом выражении точки 1 и 2 — любые две точки в области неустановившегося потенциального потока.
44
1.11.5. ВИХРЕВОЕ ДВИЖЕНИЕ НЕВЯЗКОЙ ЖИДКОСТИ
Теорема Томсона. Циркуляция скорости по замкнутому жидкому контуру в (1) идеальной, (2) баротропной жидкости при наличии массовых сил, (3) обладающих однозначным потен- циалом во всё время движения жидкости остаётся неизменной.
Проведём в движущейся жидкости в момент времени t контур AB (рис. 1.26). Этот контур считается жидким. Каждая частица M этого контура движется со скоростью V и за промежу-
ток времени |
t пройдёт расстояние dS = V t , переместившись |
|
в положение |
M ′ . |
В момент времени t′ = t + t контур займёт |
новое положение |
A′B′ . Пусть циркуляция вектора скорости в |
|
момент времени t |
равна , а в момент времени t′ она равна ′ . |
|
Тогда производная по времени от циркуляции по жидкому контуру
d |
= lim |
′ − |
. |
(1.66) |
|
|
|||
dt t →0 t |
|
|||
Необходимо доказать, что d = 0 . dt
Рис. 1.26. К доказательству теоремы Томсона о постоянстве циркуляции по замкнутому жидкому контуру
Положение произвольной частицы жидкости M на дуге AB в момент времени t можно задать её расстоянием S до точки A (рис. 1.26). В таком случае декартовы координаты частицы
45
x = x(S, t), y = y(S, t), z = z(S, t) .
Пусть в данный момент времени t проекции элемента дуги дS на оси координат равны дx, дy, дz, а проекции перемещения dS частицы M за промежуток времени t на оси координат есть dx, dy, dz . Циркуляция в момент времени t
B
= ∫(Vxдx + Vyдy + Vzдz ),
A
а её производная по времени
d |
= |
d |
B (V |
дx + V |
дy + V дz ) = B |
d |
(V |
дx + V |
дy + V дz ), |
||
|
|
|
|||||||||
dt |
dt |
∫ |
x |
y |
z |
∫ dt |
x |
y |
z |
||
|
|
|
|
|
|
||||||
|
|
|
A |
|
|
|
A |
|
|
|
|
вследствие независимости операций интегрирования по длине дуги и дифференцирования по времени.
Выполняя дифференцирование по времени под знаком интеграла, получим:
d |
= |
B |
|
|
d (дx ) |
|||||||
|
|
|
|
|
Vx |
|
|
|
||||
|
|
|
|
|
|
|
|
|||||
dt |
|
|
|
|
∫A |
|
|
dt |
||||
|
B |
dV |
x |
|
|
|
dVy |
|||||
+ |
∫ |
|
|
|
|
|
дx + |
|
|
|||
|
dt |
dt |
||||||||||
|
|
|
|
|
||||||||
|
A |
|
|
|
|
|
|
|
|
|
|
|
+ Vy |
d (дy ) |
+ Vz |
d (дz ) |
+ |
||||
|
|
|
|
|||||
dt |
dt |
|||||||
|
|
|
|
(1.67) |
||||
|
|
|
|
. |
|
|
||
дy + |
dVz |
дz |
|
|
|
|||
|
|
|
|
|||||
|
|
dt |
|
|
|
|
||
|
|
|
|
|
|
|||
Так как скорость |
d (дx ) |
|
удлинения отрезка |
|
дS вдоль оси Оx |
||||||||||||||||||||
|
dt |
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
равна разнице скоростей на его концах дVx , то |
|
|
|
|
|||||||||||||||||||||
|
|
|
d (дx ) |
= дVx , |
|
d (дy ) |
= дVy , |
|
d (дz ) |
= дVz , поэтому |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
dt |
|
|
|
|
|
dt |
|
|
dt |
|
|
|
|
|
|
||||||
|
|
d (дx ) |
+ V |
|
d (дy ) |
+ V |
|
d (дz ) |
= V |
|
|
|
+ V |
|
|
+ V дV = |
V 2 |
||||||||
V |
x |
|
|
y |
|
|
|
|
|
|
дV |
x |
дV |
y |
д |
|
. |
||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
dt |
|
dt |
|
z |
|
dt |
x |
|
|
|
y |
|
z z |
|
2 |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Рассмотрим второе подынтегральное выражение. В соответствии с уравнениями Эйлера (1.38) — (1.40)
dV |
1 ∂p |
|
dVy |
|
1 ∂p |
|
dV |
|
|
1 ∂p |
|
|||||||||
x |
дx = Xдx − |
|
|
|
|
дx, |
|
дy = Yдy − |
|
|
|
|
дx, |
|
z |
дz = Zдz − |
|
|
|
дz. |
|
с ∂x |
|
с ∂y |
|
|
с ∂z |
||||||||||||||
dt |
|
dt |
|
|
dt |
|
|
|
||||||||||||
46
Так как в соответствии с условиями теоремы внешние массовые силы обладают потенциалом U , а жидкость баротропна, то
|
|
|
|
|
|
|
|
1 |
|
∂p |
|
|
|
|
∂p |
|
|
∂p |
|
дp |
|
Xдx + Yдy + Zдz = дU и |
|
|
|
|
дx + |
|
дy + |
|
дz = |
|
. |
||||||||||
|
∂x |
∂y |
∂z |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
с |
|
|
|
|
|
|
|
с |
|||||
Следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
d |
|
B |
|
|
|
|
дp |
|
|
|
V 2 |
|
|
|
|
|||
|
|
|
|
|
∫ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
= |
|
дU |
− |
с |
|
− |
д |
2 |
|
. |
|
|
|
|
|||
|
|
|
|
A |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Если контур замкнут, то точки |
A и |
B совпадают и при одно- |
|||||||||||||||||||
значном потенциале U |
|
|
∫ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
d |
= 0 или = |
(V |
dx + V |
|
dy + V dz ) = const , |
(1.68) |
||||||||||||||
|
|
|
|||||||||||||||||||
|
|
y |
|||||||||||||||||||
|
dt |
|
|
|
x |
|
|
|
|
|
|
|
z |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
что и доказывает теорему Томсона.
Следствием этой теоремы являются ещё две теоремы.
Вторая теорема Гельмгольца о сохранении вихревых тру-
бок. При тех же предположениях, что и в теореме Томсона,
вихревая трубка во всё время движения состоит из одних и тех же частиц жидкости.
Возьмём на поверхности S вихревой трубки бесконечно малый контур L (рис. 1.27а). По формуле (1.33)
Рис. 1.27. Во все время движения вихревая трубка состоит из одних и тех же частиц жидкости, а ее интенсивность постоянная
47
d L = 2ωndS .
Так как вектор щ касается поверхности S , то во всякой точке этой поверхности нормальная составляющая вектора угловой скорости вращения частицы ωn = 0 и циркуляция d L = 0 . В
следующий момент времени t′ жидкий контур L займёт положение L′ . Бесконечно малый контур L′ будет лежать на поверхности S′ , образованной теми же частицами жидкости, которые раньше составляли поверхность S . По теореме Томсона
d L′ = d L = 0.
Но тогда из формулы (1.33) следует, что
ω′ dS′ = 0
n
Бесконечно малый контур L′ можно взять в любом месте поверхности S′ . В таком случае во всякой точке поверхности S′
ω′ = 0 .
n
Это означает, что поверхность S′ есть поверхность вихревой трубки. Каждая индивидуальная вихревая трубка перемещается в пространстве вместе с частицами, её составляющими. Таким же свойством сохраняемости обладают и вихревые линии. Вихревые нити с «нанизанными» на них частицами сохраняются во всё время движения, как показано на рис. 1.27б.
Третья теорема Гельмгольца о сохранении интенсивности вихревых трубок. При тех же предположениях, что и в тео-
реме Томсона, интенсивность любой вихревой трубки во всё время движения остаётся постоянной.
Интенсивность I вихревой трубки в момент времени t равна циркуляции скорости по контуру l, то есть I = l , а в момент
времени t′ интенсивность I ′ = ′ . Но так как по теореме Томсо- |
|
l |
|
на l = l′ , то |
|
I = I ′ = const. |
(1.69) |
На рис. 1.28 показана фотография пары развивающихся крупных атмосферных вихрей — циклонов над Индийским океаном.
48
Рис. 1.28. Пара атмосферных вихрей [12]
Интенсивность обоих вихрей одинакова. Циркуляция скорости по контуру, охватывающему вихри в начальный момент их развития, равна нулю. Так как по теореме Томсона циркуляция должна сохраняться равной нулю и во все последующие моменты времени, то вихри будут вращаться в противоположные стороны, что и подтверждается фотографией. Теорема Томсона позволяет объяснить возникновение циркуляции вокруг крыла и пары вихрей — при разгоне и остановке крыла (рис. 1.29).
Рис. 1.29. Пара вихрей при разгоне и остановке крыла по Прандтлю [7]
Однако, если какое-либо из трёх условий теоремы Томсона нарушается, то циркуляция скорости по замкнутому жидкому контуру может изменяться и в жидкости образуются или исче-
49
зают вихри. Так, при движении маловязкой жидкости вблизи твёрдых поверхностей её нельзя считать идеальной. В этих местах возможно образование вихрей. На рис. 1.30 показан вихрь, развивающийся за острым ребром в потоке жидкости.
Рис. 1.30. Вихрь за острым ребром по Прандтлю. Течение слева направо [7]
На рис. 1.31 схематически показаны два вихря, возникшие в подводящей камере осевого насоса.
Рис. 1.31. Вихри в подводящей камере осевого насоса
50
