Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

05 семестр / Книги и методические указания / Панаиотти С.С. Лекции по гидромеханике / Панаиотти С.С. Лекции по гидромеханике. Часть 1 - Учебное пособие. - Калуга КФ МГТУ им. Н.Э. Баумана, 2003

.pdf
Источник:
Скачиваний:
258
Добавлен:
04.03.2014
Размер:
1.17 Mб
Скачать
= cos(V, y) ,

При таком движении частицы в своем мгновенном вращении поворачиваются вокруг касательных к линиям тока, как показано на рис. 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 ∂ϕ является функцией координат и времени. Одна-

2t

ко уравнения (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 контур займёт

новое положение

AB′ . Пусть циркуляция вектора скорости в

момент времени 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