2-Гидромеханика 1
.pdf61
Физическаяинтерпретация
Длявсехточекданнойлиниитокаилиданнойвихревойлинии (адляпотенциальногоивинтовогопотоковдлявсехточекпотока) сумматрехудельныхэнергий– потенциальнойэнергииположения, давленияикинетической– естьвеличинапостоянная.
Энергетический смысл
z – удельная потенциальная энергия положения,
ρ pg – удельная потенциальная энергия давления, z + ρ pg – удельная потенциальная энергия,
u2
2g – удельнаякинетическаяэнергия.
Частица, имеющаямассуM ивесG = Mg идвижущаясясоско-
ростьюu, обладаеткинетическойэнергией Mu2 2 . Удельнаякинетическаяэнергия– энергия, отнесеннаякединицевеса, равна
Mu2 = Mu2 = u2 .
2G 2Mg 2g
Вуравнении Бернулли каждый член выражает собой часть механическойэнергииединицывесажидкости. Суммавсехчленов выражает суммарную (потенциальную и кинетическую) механическуюэнергию, отнесеннуюкединицевесапротекающейжидкости.
Энергия, отнесенная к единице веса жидкости, исчисляемая относительно произвольно выбранной горизонтальной плоскости, называетсяудельнойэнергией.
УравнениеБернуллипоказывает, чтоудельнаяэнергиявпотоке, присущаякаждойединицевесаневязкойжидкостисустановившимся движением, состоящаяизкинетическойипотенциальной, остается неизменной.
Впотенциальном и винтовом потоках невязкой жидкости удельнаяэнергияраспределенаравномерноповсемупотоку. Впотоке жевихревом(кромевинтового) удельнаяэнергияраспределенапо потоку неравномерно и сохраняет постоянное значение лишь на отдельныхлинияхтоковивихревыхлиниях.
62
Скоростная трубка (трубка Пито)
|
|
|
|
∆h = u |
2 |
Скоростнаятрубкапредназна- |
|
|
|
|
|
|
ченаглавнымобразомдляизмерения |
||
|
|
|
|
|
2g |
||
h |
|
= |
p1 |
u2 |
|
|
скоростногонапора. ТрубкаПито |
1 |
ρg |
+ |
|
|
состоитиздвухизогнутыхподуглом |
||
|
|
2g |
= p2 |
|
|||
|
|
|
|
h2 |
|
трубок, изкоторыхтрубка1 входит |
|
u |
|
|
u1=0 |
ρg |
|
во внутреннюю полость трубки 2 |
|
|
|
|
|
||||
|
|
|
1 |
u2=u 2 |
|
|
(рис. 4.4). |
|
|
|
|
|
Трубка 1 открытым концом |
||
|
|
|
|
|
|
Рис. 4.4 |
устанавливается против потока |
|
таким образом, чтобы скорость потока была направлена по оси отверстия. Вэтомслучаеприобтеканиижидкостьютрубкивточке1 происходитуменьшениескоростидонуля(u1= 0).
Благодаряэтомупротекающаяжидкостьвтрубке1 поднимается навысоту
h1 = ρp1 + u2 . g 2g
В трубке 2 имеется группа отверстий, расположенных в поперечномсечениинаопределенномрасстоянииотпереднегоконца. Жидкостьзаполняеттрубку2 иподнимаетсянавысоту
h2 = ρpg2 =ρ pg1 .
Изрисункавидно, что
h |
− h |
|
p |
|
u2 |
p |
|
u2 |
= ∆h, |
|
= |
1 |
+ |
|
− |
2 |
= |
|
|||
|
|
|
|
|||||||
1 |
2 |
|
ρ g |
|
2g ρ |
g |
|
2g |
|
|
|
|
|
|
|
|
откудаскоростьпотокажидкости
u = 2g∆h.
Длятогочтобыучестьвлияниевязкостижидкостиинарушений впотоке,вызываемыхустановкойтрубки,a такжеитообстоятельство, чтоточки1 и2 находятсянанекоторомконечномрасстояниидругот друга, необходимополученнуюформулудляскоростипредставитьв виде
u =ϕ 2g∆h,
гдекоэффициентϕ =1 1,04 определяетсяопытнымпутем.
63
Глава V
Динамика реальной вязкой жидкости
Преждемырассматривалижидкости, вкоторыхотсутствовала вязкость, т. е. поверхностныесилымоглибытьтольконормальными.
Реальнаяжидкостьвтойилиинойстепениобладаетвнутренним сцеплением и способна оказывать то или иное сопротивление касательнымусилиям.Величинавязкостизависитотродажидкостии оттемпературыжидкости.Силатренияявляетсядобавочнойсилой, влияющейнадвижениереальнойжидкости.
Дляреальнойжидкостидифференциальныеуравнениядвижения получаются более сложными, чем для идеальной. Для несжимаемой жидкостиониносятназваниеуравненийНавье–Стокса
F |
− |
1 |
|
∂ p |
+ ν∆ |
u= |
|
dux |
, |
||||||||
|
|
|
|
||||||||||||||
едx |
|
|
ρ ∂ |
x |
|
x |
|
dt |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|||||||||
F |
− |
|
1 ∂ p |
+ ν∆ |
u= |
|
duy |
|
|
, |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||
едy |
|
|
ρ ∂ |
y |
|
y |
|
dt |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|||||||||
F |
− |
1 |
|
|
|
∂ p |
+ ν∆ |
u= |
duz |
|
, |
||||||
|
|
|
|
|
|
||||||||||||
едz |
|
|
|
ρ ∂ |
z |
|
z |
|
dt |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
гдеν – кинематическаявязкость, ∆ u = |
∂ 2u |
x+ |
∂ |
2u |
x+ |
∂ |
2u |
x |
, ∆ – опе- |
||
|
|
|
|
|
|
|
|||||
x |
∂ x2 |
∂ |
y2 |
∂ |
z2 |
||||||
ратор Лапласа. |
|||||||||||
|
|
|
|
|
|
|
|
|
|
||
ИнтегрированиеуравненийНавье– Стоксапредставляетсобой |
чрезвычайно сложную задачу даже для простейших случаев течений жидкости. Этими вопросами занимаются специалисты по гидромеханике. Мы получим уравнение движения для потока реальной жидкости (уравнение Бернулли) другим путем. Проанализируем полученноеранееуравнениеБернуллидляпотоканевязкойжидкостии внесемвнегопоправки,учитывающиевязкостьжидкости.
64
Уравнение Бернулли для струйки реальной жидкости в установившемся потоке
УравнениеБернуллидляструйкиидеальнойжидкости
z |
+ |
p1 |
+ |
u12 |
= z |
2 |
+ |
p2 |
+ |
u22 |
. |
|
|
|
|
||||||||
1 |
|
ρ g |
|
2g |
|
ρ g |
|
2g |
|||
|
|
|
|
|
|
|
Удельная энергия в струйке реальной жидкости при установившемсядвиженииуменьшаетсяпомерепродвиженияжидкостиот одного сечения струйки до другого, так как часть механической энергиипревращаетсявтепловуюэнергиюблагодарясиламтрения, возникающимвжидкости.
ВсечленыуравненияБернуллиимеютлинейнуюразмерность, поэтому удельную энергию, затрачиваемую на преодоление сопротивлениядвижению, обозначаютсимволомhтр– потерянапора натрение.
УравнениеБернуллидляструйкиреальнойжидкостипринимает
вид
z |
+ |
p1 |
+ |
u12 |
= z |
|
+ |
p2 |
+ |
u22 |
+ h |
, |
ρ g |
2g |
|
|
2g |
||||||||
1 |
|
|
|
2 |
ρ g |
|
тр |
|
гдеhтр – работавсехсилтрения, отнесеннаякединицевесавязкой жидкостиприперемещенииеёотпервогосечениядовторого.
Распределение давления по живому сечению потока при установившемся прямолинейном движении
Прирешениигидравлическихзадачвозникаетнеобхо-димость распространить уравнение Бернулли на поток в целом. Для этого предварительнотребуетсянайтизаконраспределениядавленияпо живомусечениюпотока. Этотвопросможетбытьрешендляживых сеченийпотоков, находящихсявусловияхплав-нойизменяемости. Приплавноизменяющемсядвиженииуголрасхожденияструек(линий токов) крайненезначителен, арадиусихкривизнывесьмавелик.
Рассмотримпрямолинейноедвижениежидкости. Еслинаправитьось0x вдольдвижения, топроекциискоростейбудутравны
ux= u, uy= 0, uz= 0,
ауравнениенеразрывностиупростится
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
65 |
|
|
|
|
|
|
|
|
|
|
|
∂ ux |
= |
∂ u |
= 0 . |
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
∂ x |
|
∂ x |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Дифференциальныеуравнениядвиженияидеальнойжидкости |
|||||||||||||||||||||||||||||||
(уравненияЭйлера) примутвид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
∂ p |
= ρ Fедx |
, |
|
∂ p |
|
= ρ Fедy , |
|
∂ p |
= ρ Fедz , |
||||||||||||||||||||
|
|
|
∂ y |
|
|||||||||||||||||||||||||||
|
|
∂ x |
|
|
∂ z |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂ u |
|
|
|
|||||||||||
так как при установившемся движении |
|
|
= 0 и правая часть |
||||||||||||||||||||||||||||
|
|
∂ t |
|||||||||||||||||||||||||||||
развернутогопервогоуравненияЭйлера |
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
du |
x |
|
= |
∂ u |
x |
+ |
∂ |
u |
x |
u |
x |
+ |
|
∂ |
uy |
u |
y |
+ |
∂ u |
z |
|
u |
z |
= 0. |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
dt |
|
|
|
|
|
|
|
∂ |
y |
|
|
z |
||||||||||||||||||
|
|
∂ t |
|
∂ x |
|
|
∂ |
|
|
|
|
Этиуравненияпоказывают, чтодавлениепоживомусечениюв установившемся, плавно изменяющемся движении идеальной жидкости распределяется по закону гидростатики.
В реальных потоках вязкой жидкости отклонения от этого законанезначительны, иимипренебрегаютисчитают, чтодавление поживомусечениюизменяетсяпогидростатическомузакону, т. е.
z + ρ pg = const .
Уравнение Бернулли для установившегося потока реальной жидкости
ВыведемуравнениеБернуллидляпотокареальнойжидкости (река, канал, водопроводит. д.).
Какбылопоказановыше, вразличныхточкахживогосечения такогопотока
z + ρ pg = const.
u2
Удельнаякинетическаяэнергия 2g будетразличнойдляраз-
ныхточексеченияпотокавязкойжидкости, таккакu = u(x, y, z). Вычислимсреднеезначениеудельнойкинетическойэнергиипо
сечению. Средняяскоростьврассматриваемомсеченииv=Q/S. Местныескоростиuбудутотличатьсяотсреднейнавеличину±∆ u :
u = v ± ∆ u .
66
Подсчитаемрасходпотокареальнойжидкости
|
Q = ∫udS = Sv, |
|
или |
S |
|
Q = ∫(v ± ∆ u)dS= v∫ dS+ ∫±∆ udS= vS, |
||
S |
S |
S |
значит |
|
. |
|
∫ ± ∆ udS = 0 |
|
|
|
S
Кинетическая энергия массы жидкости, протекающей через сечениеS вединицувремени,
∫±∆
S
как
|
|
|
u2 |
|
ρ |
3 |
|
ρ |
|
3 |
K = ∫ρ udS |
|
= |
|
∫ u |
dS = |
|
|
∫ (v ± ∆ u ) dS= |
||
2 |
2 |
|
2 |
|||||||
S |
|
|
|
S |
|
|
S |
|||
= |
ρ |
∫(v3 + 3v∆ u2± 3v∆2 u± ∆ u3 ) dS, |
||||||||
|
2 |
S |
|
|
|
|
|
|
|
|
u3dS – пренебрегаем ввиду его малости, 3v2 ∫±∆ udS= 0 – (так
∫±∆ udS= |
0 ), v3 ∫ dS = v3S . |
|
|
|
|
|
S |
|
|
|
|
|
|
|
|||
S |
S |
|
|
|
|
|
|
|
В результате получим |
|
|
|
|
|
|
||
|
|
3 |
|
|
∫ ∆ |
u2dS |
||
|
|
ρ v S |
|
|
S |
|
|
|
|
K = |
|
1+ 3 |
|
|
. |
||
|
2 |
|
Sv |
2 |
||||
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
Для определения удельной кинетической энергии разделим выражениенавесжидкости, протекающейвединицувремени, ρ gSv:
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
∫ ∆ |
u2dS |
|
|
2 |
|
||
|
|
|
|
|
|
|
|
K |
= |
v |
|
|
1+ 3 |
S |
|
|
|
= |
α v |
|
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
ρ gSv |
|
|
|
Sv |
2 |
2g |
|||||||||
|
|
|
|
|
|
|
|
2g |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∫ ∆ |
u2dS |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
где α = |
|
1+ |
3 |
S |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
Sv |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
67 |
|
Отношениедействительногозначениякинетическойэнергиик |
|||||||||
условновычисленнойвпредположении, чтоскоростивовсехточках |
||||||||||
живогосеченияравнысреднейскорости, |
|
|
||||||||
|
|
|
|
∫u3dS |
|
|
|
|||
|
|
|
|
α = S |
Sv3 |
|
|
|
||
|
|
|
|
|
|
|
|
|||
называетсякоэффициентомкинетическойэнергииα . |
||||||||||
|
Удельнаяэнергияпотока |
p |
|
|
α |
v2 |
|
|||
|
|
|
|
|
|
|
||||
|
|
|
H = z + ρ g + |
2g . |
||||||
|
В турбулентных потоках, встречающихся на практике, при- |
|||||||||
нимаютα =1 (дляламинарного– α |
= 2). |
|
|
|||||||
|
Уравнение Бернулли примет вид для потока реальной |
|||||||||
жидкости |
p |
α 1v12 |
|
|
p |
|
α |
v2 |
||
|
z1 + |
= z2 + |
+ |
|||||||
|
1 + |
2g |
|
2 |
|
2 2 + hтр. |
||||
|
|
ρ g |
|
ρ |
g |
|
|
2g |
||
|
|
|
|
h1-2 |
|
|
|
|
|
Такимобразом, мыуста- |
α |
v 2 |
|
|
2 |
|
новили, что уравнение Бер- |
||||
|
|
α v |
|
|||||||
1 |
1 |
|
|
2 |
|
нуллидляцелогопотокавяз- |
||||
2 g |
|
|
2 |
|
|
|||||
|
|
2 g |
|
|
||||||
|
|
|
|
|
|
|||||
|
p |
|
|
p |
|
|
кой жидкости по своему по- |
|||
|
1 |
|
|
|
|
строению аналогично урав- |
||||
ρ g |
|
|
2 |
|
|
|||||
|
|
ρ g |
|
|
||||||
|
1 |
|
|
|
|
|
нениюБернуллидляэлемен- |
|||
|
|
|
|
2 |
|
|
тарной струйки. Мы как бы |
|||
|
z1 |
|
|
|
|
увеличили элементарную |
||||
|
|
|
z2 |
|
|
|||||
|
|
|
|
|
|
струйку до размеров целого |
||||
|
l |
|
|
|
|
|
||||
|
|
|
|
|
|
потока. Новым элементом |
||||
|
|
|
|
|
|
|
||||
|
Рис. 5.1 |
|
|
|
|
|
здесьявляютсякоэффициенты |
|||
|
|
|
|
|
|
|
α |
1 |
и α 2, учитывающие влия- |
|
ниенеравномерностираспределенияскоростейпоживомусечению |
||||||||||
потоканаегокинетическуюэнергию. |
|
|
|
|
||||||
|
Отношениепотерьнапораhтр |
кдлинеl, накоторойэтапотеря |
||||||||
происходит, называют гидравлическим уклоном |
|
|
|
|
z |
+ |
p |
+ |
α 1v12 |
− |
|
z |
|
+ |
p |
+ |
α |
v |
2 |
|
|
||
|
|
|
|
1 |
|
|
|
|
2 |
|
2 2 |
|
|
|
||||||||
|
|
|
|
|
2 |
|
|
|
|
|
||||||||||||
|
h |
|
|
1 |
|
ρ g |
|
2g |
|
|
|
ρ g |
|
|
2g |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
I = |
тр |
= |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
l |
|
|
|
|
|
|
|
l |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
68
Глава VI
РЕЖИМЫ ДВИЖЕНИЯ ЖИДКОСТИ
Ламинарное и турбулентное движение жидкости
Основной задачей гидродинамики является исследование закономерностей, характеризующихпотоквцелом. Однакопроблема гидравлических сопротивлений и многих других явлений требует ясногопредставленияохарактередвиженияотдельныхчастиц.
Еслижидкостьнеимеетвсвоемсоставечастицтвердоготелаи вязкость ее подчиняется закону Ньютона, для данной жидкости возможенлиболаминарныйрежим, либотурбулентный.
Существованиедвухпоследнихрежимовдвиженияжидкости было доказано в 1883 г. Рейнольдсом путем постановки опыта на специальнойустановке, схемакоторойпоказананарисунке.
1
3
2 v
От водопровода
Слив
Рис. 6.1
Рейнольдс наблюдал движение воды в стеклянных трубках различныхдиаметров. Онвпускализсосуда1 спомощьюкапилляра 2 поосистекляннойтрубки3 струйкукраскиисоздавалразличные скорости течения воды v. Он заметил, что при малых скоростях подкрашеннаяструйкавыделяласьввидетонкойнити, несмешиваясь с остальной массой жидкости. Тем самым Рейнольдсом было установленоналичиеструйчатого, илислоистого, движенияпотока,
названногоимламинарнымпотоком.
69
Придальнейшемповышениисреднейскоростипотокапосле определенногоее ее значения параллельность струек нарушалась, происходилоперемешиваниеокрашенныхчастицсостальноймассой жидкостиипотокповсемусечениюприобреталодинаковыйцвет, чтоуказывалонанепрерывноеобразованиевихрейвпотоке. Такой потокРейнольдсомбылназвантурбулентным.
Движение жидкости, при котором отсутствует пульсация скорости, приводящаякперемешиваниючастиц, называетсяламинарным. Движениежидкостиспульсациейскоростей, приводящейк перемешиваниючастицпотока, называетсятурбулентным.
В процессе опытов Рейнольдсом было замечено, что после нарушенияламинарностипотоканесразуустанавливаетсявполне устойчивыйтурбулентныйрежим, асуществуетобластьперехода. Эта областьимбылаограниченавведениемпонятиянижнейиверхней критических скоростей. Скорость, при которой происходило нарушениеструйчатостипотока, Рейнольдсназвалнижнейкритическойскоростьюvнк, аскорость, послекоторойпотокстановился устойчивым турбулентным, – верхней критической vвк. Им также установлено, чтоверхняякритическаяскоростьпревосходитзначения нижнейкритическойскоростив6 раз, т. е. vвк= 6 vнк.
ДляопределениярежимадвиженияРейнольдсомнаосновании анализаопытныхданныхиисходяизтеоретическихсоображенийдля потокапроизвольнойформыпоперечногосечениябылаустановлена зависимостьввидеформулы
Re = vlµρ ,
гдеl – характерныйлинейныйразмерпоперечногосеченияпотока. Зависимостьдлякруглойтрубысзаполненнымсечениемимеет
вид |
|
|
|||
Re = |
vdρ |
|
, |
||
|
|
|
|||
µ |
|
|
|||
или |
|
|
|||
Re = |
vd |
, |
|||
ν |
|||||
|
|
|
|
гдеRe – безразмернаявеличина– числоРейнольдса, v – средняяскоростьтеченияжидкости,
d – диаметртрубы,
ν – кинематическаявязкость.
Есливкачествелинейногоразмерапринятьсоответствующую гидравлическую характеристику, то формула перепишется в следующемвиде:
70
Re |
|
= |
vd |
, |
Re |
|
= |
vR |
, |
Re |
|
= |
vh |
, |
d |
ν |
r |
ν |
h |
|
|||||||||
|
|
|
|
|
|
|
ν |
|
гдеd – диаметртрубы, R – гидравлическийрадиус, h – глубинаканала. ВтеориикораблявводитсявеличинаL – длинасудна.
Рейнольдсом было предложено два значения Re, соответствующихнижнейиверхнейкритическойскорости– 2000 и12000.
Впоследствии, в 1925 г., на основании многочисленных экспериментальныхработШиллерпредложилвместодвухпределовдляRe однозначение, аименно– 2320.
Это значение и принимается за Reкр
Reкр = |
vd |
= 2320, |
Rer (кр) = |
vR |
= 580. |
|
ν |
ν |
|||||
|
|
|
|
ВламинарномпотокеприRe <2320 возмущения, возникающие от посторонних причин, затухают. В потоке при Re >2320 и в турбулентномнезатухают, аещебольшеразвиваются.
Всоответствии с этим величина критической скорости при достаточной длине трубы с гладкими стенками может быть с достаточной точностью для практики определена по формуле
vкр = 2320d ν .
Анализ экспериментальных работ по определению потери напора в прямых трубопроводах, проведенных различными исследователями, притом с разнообразными жидкостями, показывает зависящее от режима движения различное влияние скорости на величину гидравлических потерь. Так, например, потери напора при ламинарном режиме
hтр= b1v
ипритурбулентномрежиме
hтр= b2vn,
гдеn – показательстепени; онколеблетсявпределахот1,75 до2. Следовательно, нет плавного перехода от одного режима к
другому. Действительно, приламинарномрежимеправноединице, а вслучаетурбулентногорежимапоказательстепенипрезковозрастает, хотябыдосвоегоминимальногозначения– до1,75.