Скачиваний:
116
Добавлен:
29.03.2015
Размер:
1.01 Mб
Скачать

После интегрирования по напору от Н1 до Н2 получим выражение по времени опорожнения

t = 2Ω(H1 H2 ), с.

μω 2g

При выполнении условия H2 = 0 время полного опорожнения резервуара определится из выражения

t =

2Ω

H1

 

=

2ΩH1

=

2V

, с,

 

 

 

 

 

 

 

 

μω

2g

μω 2gH

 

Q

 

 

 

 

 

 

1

 

 

 

где V – объём резервуара, м3; Q – расход жидкости при начальном на-

поре H1, м3/с.

Таким образом, время полного опорожнения резервуара при переменном напоре в два раза больше времени, которое требуется для вытекания из резервуара жидкости при начальном напоре в количестве, равном первоначальному объёму.

Истечение в сообщающихся резервуарах

Рассмотрим случай выравнивания уровней жидкости в сообщающихся резервуарах I и II (рис. 6.11).

Определим время, в течение которого уровни в резервуарах выравниваются.

p

 

 

 

dz1

 

H

H1 pат H2

 

 

z1

 

 

dz2

 

 

 

z2

0

 

 

0

I

Ω1

Ω2

II

Рис. 6.11. Истечение жидкости под уровень при переменном напоре

51

За время dt в результате перетекания жидкости из резервуара I в резервуар II напор над центром тяжести отверстия со стороны первого резервуара уменьшится на величину dz1, а со стороны второго резервуара, наоборот, увеличится на величину dz2. При этом объём жидкости в первом резервуаре уменьшится на величину

− Ω1dz1 = μω2gH dt , м3,

откуда

dt = −

 

Ω1

 

dz1

, с.

μω

 

 

 

 

2gH

 

 

 

 

 

С другой стороны, объём жидкости во втором резервуаре увеличится на величину Ω1dz2 .

Очевидно

−Ω dz

= Ω

2

dz

2

, м3.

 

 

1

1

 

 

 

 

 

 

 

 

 

 

Откуда

 

 

 

 

 

Ω1

 

 

 

 

 

 

 

 

 

dz

2

= −

dz , м.

 

 

 

 

 

 

 

 

 

 

Ω2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dH = dz1 dz2 , м,

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dH = dz +

Ω1

 

dz =

Ω1 + Ω2 dz , м.

 

1

 

Ω2

1

 

 

 

 

 

 

Ω2

1

 

 

 

 

 

 

 

 

 

 

 

Следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dz1 =

 

 

Ω2

 

 

 

dH , м.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ω1 + Ω2

 

 

 

 

 

 

Окончательно имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt = −

 

Ω1Ω2

 

 

 

dH

, с.

Ω1 + Ω2

 

μω

 

 

 

2gH

Проинтегрировав это выражение в пределах разности напоров H1 и H2, получим выражение для времени выравнивания уровней жидкости в резервуарах

t =

2Ω1Ω2

H1

 

 

, с.

μω

 

(Ω + Ω

)

2g

 

1

2

 

 

52

7. ТЕЧЕНИЕ ЖИДКОСТИ И НЕВЯЗКОГО (ИДЕАЛЬНОГО) ГАЗА

Скорость распространения возмущений. Метод малых возмущений. Гидравлический удар. Прямой скачок уплотнения. Характеристики одномерного течения. Движение идеального газа в канале переменного сечения (сопло Лаваля). Сверхзвуковые течения. Косой скачок уплотнения. Особенности двухкомпонентных и двухфазных течений. Течение жидкости при фазовом равновесии. Тепловой скачок и скачок конденсации.

МЕТОДИЧЕСКИЕ УКАЗАНИЯ

Определение скорости распространения малых возмущений в не-

подвижной газовой среде.

Пусть по какой-либо причине в трубе образовалась волна, движущаяся влево со скоростью u (рис. 7.1 а). Причиной образования такой волны может быть, например, ускоренное или замедленное движение поршня. Чтобы изучить стационарное течение, сообщим системе скорость –u . Тогда волна окажется неподвижной, слева от волны газ будет двигаться со скоростью u (вправо), а справа – со скоростью u+∆u (рис. 7.1 б).

Единственное принимаемое допущение – малое значение ∆u по сравнению с u, т.е.

u << 1. u

Другими словами, рассматриваются волны с конечным, но малым изменением параметров.

Чтобы установить связь между параметрами газа до и после волны, применим уравнение неразрывности и количества движения

 

 

 

ρu = (ρ + ρ)(u + u );

(7.1)

 

 

ρu 2 + p = (ρ + ρ)(u +

u)2 + p +

 

 

 

 

p.

 

 

и

 

 

и, ρ

 

 

 

и +

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ +

ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а)

 

б)

 

 

Рис. 7.1. К расчёту скорости распространения малых возмущений

53

Совместное решение уравнений (7.1) при условии u << 1 позво- u

ляет найти зависимость между изменением давления и плотности в волне (∆ p и ∆ρ) и скоростью распространения волны u

u =

p

, м/с.

(7.2)

ρ

 

 

 

Скоростью звука называют скорость распространения волны очень малой интенсивности, т.е.

a =

dp

 

 

 

 

 

, м/с.

(7.3)

 

 

 

 

 

 

 

 

dρ

S

 

 

 

 

 

 

При малой интенсивности волны потери в ней отсутствуют и поэтому распространение звуковой волны (а практически – и волн конечной, но малой интенсивности) можно рассматривать как изоэнтропическое течение. Соответственно этому индекс S в формуле (7.3)

подчёркивает, что производная dp должна браться при постоянной dρ

энтропии.

Известно, что скорость звука зависит от температуры газа. Эту зависимость легко установить, воспользовавшись уравнением

p= const

ρk

иуравнением состояния. Легко убедиться, что

a = kzRT , м/с.

Таким образом, скорость звука зависит от физических свойств газа и его температуры.

Наряду со скоростью звука а, важную роль играет критическая скорость. Смысл введения критической скорости становится ясным, если воспользуемся уравнением энергии

i2* i1* = Q + H

при условии отсутствия подвода энергии и тепла (Q = H = 0), т.е.

i2* = i1* , Дж/кг,

54

или

 

 

 

 

 

 

 

 

 

 

i +

u 2

= i* = const .

(7.4)

 

 

 

 

2

 

 

 

 

 

 

 

 

Выполним несложные преобразования:

 

 

 

 

i = c T =

 

 

k

 

zRT =

 

a2

 

 

 

 

k 1

k 1

p

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

i* = c T * =

 

 

k

zRT * =

a

2

 

 

 

 

 

 

0

.

 

 

 

 

 

p

 

k 1

 

 

k 1

 

 

 

 

 

где a0 – скорость звука при температуре торможения T * .

Реально с такой скоростью звук будет распространяться в резервуаре, из которого происходит истечение газа (и в котором T = T * ), т.е.

a2

 

u2

 

a2

 

 

 

 

0

 

 

k 1

+

2

=

k 1

.

(7.5)

При движении газа возможен случай, когда скорость и станет равной местной скорости звука а. В этом случае её называют критиче-

ской скоростью и обозначают a*

u = a = a* , м/с.

Из выражения (7.5) следует зависимость

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2kzRT *

 

 

 

a =

2

a

 

=

 

 

, м/с.

(7.6)

k +1

 

k +1

*

 

 

0

 

 

 

 

 

Для определения скорости звука в жидкости воспользуемся зависимостью между коэффициентом объёмного сжатия βс , равного отно-

сительному изменению объёма жидкости, поделенному на вызвавший это изменение перепад давления, т.е.

β

 

= −

1

V , 1/Па.

с

 

 

 

V

p

 

 

 

Но согласно уравнению неразрывности имеем

V

= −

ρ

 

ρ

V

55

и, следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

β

 

=

1

 

 

ρ

, 1/Па.

(6.7)

 

 

 

c

ρ

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

Из выражений (7.2) и (7.7) получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a =

 

 

E

 

, м/с,

(6.8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ

 

1

 

 

 

 

 

 

 

 

 

 

 

где Е =

 

 

– модуль объёмной упругости жидкости, Па.

 

β

с

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выражение (7.8) определяет скорость звука в жидкости. Скорость звука в воде при обычной температуре почти в четыре раза больше, чем в воздухе. Из формулы (7.8) следует, что в абсолютно несжимаемой жидкости возмущения распространяются мгновенно (а = ).

Гидравлическим ударом называется явление повышения давления в трубопроводе при резком нанесении возмущения. Такое возмущение может быть вызвано резким открытием или закрытием задвижки. В неблагоприятном случае гидравлический удар может быть причиной разрыва трубопровода.

Метод расчёта гидравлического удара, впервые предложенный Н.Е. Жуковским, заключается в определении величины повышения давления внутри трубопровода в момент возникновения гидроудара.

Повышение давления р , вызванное гидроударом, можно определить, используя уравнение неразрывности

p = −u2 ρ − 2uρ u , Па.

Принимая u = a, и, учитывая, что согласно уравнению неразрывности

ρ= − u ,

ρa

получим формулу Жуковского

p = −ρa u , Па,

(7.9)

где ρ – плотность среды, кг/м3; а – скорость звука в данной среде при данных условиях, м/с; u – разность скоростей до и после волны гидроудара, м/с.

u = u2 u1 , м.

56

Важность изучения явления гидравлического удара следует потому, что, например, если скорость воды в трубопроводе равна 5 м/с, то при скорости звука 1430 м/с повышение давления согласно форму-

ле (7.9) составит более 7 ×106 Па (около 70 атм.). При диаметре трубы d = 0,1 м гидравлический удар вызывает осевое усилие около 5,5 т.

Выражение (7.8) даёт максимальное значение скорости звука, поскольку оно не учитывает деформацию трубопровода во время гидравлического удара.

Сильные возмущения в газе, вызванные взрывом или сверхзвуковым обтеканием тела, также вызывают ударную волну. Ударная волна при обтекании тел или в канале называется скачком уплотнения.

Расчёт параметров скачка уплотнения производят с помощью уравнений неразрывности и движения, однако в данном случае приращение плотности и скорости нельзя считать малыми:

r1и1 = r2и2 ;

r1и12 + р1 = r2и2 + р2 .

Добавив к этим уравнениям уравнение энергии

i + c T = k p ,

p

k -1 r

 

получим замкнутую систему уравнений, позволяющую найти р2, ρ2, и2:

 

 

и и

2

= а2 ;

 

 

 

 

 

1

 

 

 

*

 

 

 

 

r2и2 = r1и1;

 

 

 

 

 

k +1

 

r2

-1

 

 

 

 

 

 

 

 

р

2

 

 

k -1 r

 

 

 

 

=

 

 

 

 

 

1

 

.

(7.10)

 

 

 

 

k +1

 

 

р

 

 

 

 

r2

 

 

1

 

 

 

k -1

- r

 

 

 

 

 

 

 

 

 

 

1

 

 

Следует также отметить, что зависимость (7.10) получена Ранкиным в 1870 г.

Характеристики одномерного течения

Наряду с двумя характерными скоростями – скоростью звука а и критической скоростью а* (см. стр. 54, 55), существует ещё одна ха-

рактерная скорость – максимальная скорость umax , достигаемая при

57

истечении газа в пустоту. Скорость umax

можно определить из уравне-

ния энергии в форме (7.5), приняв а = 0 (Т = 0)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2kzRT *

 

 

 

и

 

=

2

a

 

 

=

 

, м/с.

(7.11)

max

 

k 1

0

 

 

k 1

 

 

 

 

 

 

 

 

 

 

Сопоставляя (7.11) и (7.6), получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

иmax

=

 

 

 

k + 1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а*

 

 

 

k 1

 

 

 

Соответственно трём характерным скоростям вводятся безразмерные характеристики течения:

число Маха М;

безразмерную скорость Θ ;

безразмерную скорость Ψ ;

M = u ; a

Θ = u ; a*

Ψ = u .

umaх

Любая из этих безразмерных характеристик может рассматриваться как основной критерий подобия течений газа с большими скоростями, так как между ними существует однозначная связь. Например, между М и Θ существует зависимость

Θ2 =

 

(k + 1)M

2

.

(7.12)

 

 

 

2

+ (k 1)M 2

 

Принято, однако, считать основным критерием подобия число Маха М, так как оно непосредственно и простым образом связано с некоторыми важными характеристиками течения.

На практике предпочитают пользоваться безразмерной скоростью Θ , так как для течений без теплообмена и подвода энергии a* = const и, следовательно, величина Θ пропорциональна скорости

течения и.

Диапазон изменения чисел М, Θ , Ψ очевиден:

58

0 < М < ;

0< Θ < k +1 ; k 1

0< Ψ <1.

Связь между безразмерной температурой τ и числом Маха М, безразмерной скоростью Θ и безразмерной скоростью Ψ следующая:

τ =

T

=

 

 

 

1

 

;

(7.13)

 

 

 

 

 

 

T *

1 +

k

1

M 2

 

 

 

 

 

 

 

 

2

 

 

 

 

τ = 1

k 1

Θ2 ;

(7.14)

k + 1

 

 

 

 

 

 

 

 

 

 

τ = 1 − Ψ2 .

 

 

 

(7.15)

Движение идеального газа в канале переменного сечения. Сопло Лаваля

Рассмотрим закономерности стационарного движения идеального газа в канале переменного сечения. Воспользуемся уравнениями неразрывности и Бернулли

ρωи = соnst ;

(7.16)

udu =

dp

.

(7.17)

 

 

ρ

 

Прежде всего установим прямую связь между скоростью и площадью проходного сечения. Для этого из уравнения (7.16) найдём приращения функций

du + dρ + dω = 0 . u ρ ω

Далее выполним следующие преобразования

dρ

=

dρ

 

dp

=

1

 

dp

,

ρ

 

 

a2

 

ρ

 

dp ρ

 

 

или, воспользовавшись уравнением Бернулли, получим

dρ = − u du .

ρ a2

(7.18)

(7.19)

(7.20)

59

Подставив (7.18), (7.19) в уравнение неразрывности, окончательно получим

(M 2 1)

1

 

du

=

1

 

dω

.

(7.21)

 

 

 

 

 

u dx ω dx

 

Уравнение (7.21) показывает, что характер изменения скорости в канале переменного сечения зависит от числа М.

Для анализа характера изменения скорости истечения преобразу-

ем выражение (7.21) следующим образом

 

 

 

 

 

du

=

dωω

.

 

 

 

 

 

 

 

 

 

 

 

 

u

M2 1

 

 

Если течение дозвуковое, то

M 2 1 < 0 и отрицательной произ-

водной

dω

соответствует положительная производная

du

. Другими

ω

u

словами, в случае истечения газа при дозвуковых скоростях, как и в случае несжимаемой жидкости, уменьшение площади проходного сечения вызывает увеличение скорости истечения. Отличия по сравнению с течением несжимаемой жидкости только количественные (рис. 7.2).

При сверхзвуковом течении M 2 1 > 0 и положительной произ-

водной

dω

соответствует положительная производная

du

, т.е. увели-

ω

u

чение площади проходного сечения канала вызывает также увеличение скорости истечения газа или жидкости (рис. 7.3).

М < 1

p

1

q

M

Рис. 7.2. Дозвуковое течение

60