
Gidrogazovaya_dinamika / Учебники по гидрогазовой динамике / ghukov-a
.pdf
После интегрирования по напору от Н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).
Определим время, в течение которого уровни в резервуарах выравниваются.
paт |
|
|
|
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