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

Методическое пособие 749

.pdf
Скачиваний:
6
Добавлен:
30.04.2022
Размер:
7.44 Mб
Скачать

Выражая в уравнении энергии (4.6) отношение p через уравнение состо-

яния (4.1), получим

v2

 

kRT

 

const .

(4.8)

2

k 1

 

 

 

Последнее равенство показывает, что при отсутствии теплообмена с внешней средой увеличение скорости вдоль струйки приводит к падению температуры газа и наоборот.

Используя формулу для скорости звука, уравнение энергии (4.6) можно представить в виде

v2

 

a2

 

const ,

(4.9)

2

k 1

 

 

 

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

Выражение (4.9) позволяет выяснить смысл постоянной в правой части уравнения энергии. Действительно, в покоящемся газе v 0 и скорость звука

достигает здесь своей наибольшей величины ao . Следовательно, const

a2

и

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 1

 

уравнение энергии может быть представлено в виде

 

 

 

 

v2

 

a2

 

a2

 

 

 

 

 

 

 

0

.

(4.10)

 

 

k 1

 

2

 

 

k 1

 

 

 

Наконец, если использовать понятие энтальпии, или теплосодержания, газа i , рассматриваемое в термодинамике

i cpT

k

RT ,

 

 

k 1

 

 

то уравнение энергии (4.8) приобретает вид

v2

i const .

(4.11)

2

 

 

Таким образом, потенциальная энергия газа выражается в различных формах уравнения энергии (4.6)-(4.11) с помощью различных взаимосвязанных параметров – давления, температуры, скорости звука, энтальпии. Ниже показа-

120

ны примеры применения уравнения энергии в различных формах записи для решения задач одномерного течения.

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

Выясним зависимость скорости течения от площади поперечного сечения потока. Для газа уравнение неразрывности (2.7), или уравнение постоянства массового расхода при установившемся течении имеет вид

vf const .

(4.12)

Логарифмируя и дифференцируя это равенство, получим

d

 

dv

 

df

0 ,

(4.13)

 

v

f

 

 

 

 

откуда

 

 

 

 

 

 

d

 

 

 

dv

 

 

 

 

 

 

 

 

dF

 

 

 

 

 

 

 

 

1

 

 

 

 

 

.

(4.14)

F

v

 

 

dv

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

Из уравнения энергии в дифференциальной форме (4.5) имеем

dvv v12 dp ,

что даёт после подстановки в уравнение (4.14)

dF

 

dv

 

d

 

 

v2

1 .

F

v

dp

 

 

 

Поскольку согласно ранее приведенной формуле

dp

a2 , имеем

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dF

 

dv

 

2

 

 

dv

M 2

1 .

 

 

 

 

 

v

 

1

 

 

(4.15)

 

F

v

 

2

 

 

 

 

a

 

 

 

v

 

 

 

Из уравнения (4.15) следует, что изменение скорости dv при изменении сечения dF происходит по-разному для дозвукового и сверхзвукового течения. В дозвуковом потоке ( v a , M 1, рис. 4.1, а) знаки dv и dF в (4.15) противоположны: уменьшение сечения в конфузорном канале приводит к возрастанию

121

скорости, а в диффузорном канале скорость по потоку уменьшается. При сверхзвуковом течении ( v a , M 1, рис. 4.1, б) в конфузорном канале поток замедляется, в диффузорном – ускоряется.

а

б

Рис. 4.1. Изменение скорости при разных числах Маха

Чтобы пояснить полученные результаты, которые для сверхзвукового течения кажутся неожиданными, сопоставим уравнения (4.15) и (4.14). Имеем

 

 

d

 

 

 

 

 

 

 

M 2

 

 

.

(4.16)

 

 

 

 

dv

 

v

Поскольку левая часть равенства (4.16) всегда положительна, ясно, что знаки d и dv всегда противоположны: рост скорости приводит к уменьшению

плотности. Но при дозвуковом течении ( M 2 1) скорость изменяется более быстро, чем плотность

dv d . v

При сверхзвуковом течении, наоборот, более быстро уменьшается плот-

ность

dv d . v

Для получения сверхзвуковых скоростей газа в технике используется сопло Лаваля, схема которого приведена на рис. 4.2, принцип действия которого ясен из приведённых рассуждений. В дозвуковом потоке, поступающем в суживающуюся часть сопла Лаваля, скорость увеличивается. Если в наименьшем сечении сопла не достигается скорость, равная скорости звука, то в расширяющейся части происходит её уменьшение; скорость по длине сопла изменяется по кривой 1 (рис. 4.2). Если перепад давления достаточно велик, чтобы в наименьшем сечении скорость течения сравнялась со скоростью звука, то при дальнейшем расширении поток переходит в сверхзвуковой и скорость его изменяется по кривой 2.

122

Рис. 4.2. Изменение скорости по длине сопла

Сопло Лаваля имеет широкое применение, являясь составной частью реактивных двигателей, сопловых аппаратов некоторых турбин (в которых рабочие лопатки обтекаются сверхзвуковым потоком), сверхзвуковых аэродинамических труб и т.д. Более полная теория сопла учитывает влияние трения на стенках и волновых явлений на выходе потока.

4.4. Параметры изоэнтропического торможения газа

При торможении газа его кинетическая энергия переходит в потенциальную, при этом давление, плотность и температура возрастают.

В случае полного торможения (остановки) потока, например в точке раздвоения струйки на передней поверхности обтекаемого тела, p, ,T достигают

максимальных для данного потока величин – параметров торможения p0 , 0 ,T0 .

Определим эти величины для адиабатного изоэнтропического процесса торможения, при котором давление и плотность газа связаны соотношением (4.2).

Применим уравнение (4.9) к сечениям струйки «на бесконечности», т.е. там, где на поток не сказывается искажающее влияние обтекаемого тела, и в точке торможения

v2

a2

 

 

a2

 

 

 

 

 

0

.

 

k

 

 

2

 

1

k 1

Как и следовало ожидать, мы получили частный случай уравнения (4.10).

 

 

a2

 

 

 

Разделив последнее выражение на

 

 

 

 

, получаем

k

 

1

 

 

 

 

 

 

 

k 1 v2

1

a2

 

 

 

 

 

0

.

 

 

 

a2

 

2

 

 

 

a2

 

 

 

 

 

 

 

123

Принимая во внимание, что

v

 

M , где М– число М для невозмущён-

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

 

T

 

 

ного потока, и что согласно формуле

0

 

0

, имеем

 

a2

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

T0

 

1

k 1

M 2 .

(4.17)

 

 

 

 

 

T

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

Повышение температуры газа у поверхности тела, обтекаемого при больших числах M , называется аэродинамическим нагревом. Отметим, что термометр, помещённый в поток газа, показывает температуру, очень близкую к температуре торможения.

Используя зависимости (4.3), связывающие температуру адиабатного процесса с давлением и плотностью, получим

 

 

 

k 1

 

k

 

 

 

p

 

 

2

k 1

 

0

0

1

 

 

 

M

;

 

p

 

 

 

 

 

2

 

 

 

Зависимость температуры, давления M представлена графически на рис. 4.3.

 

 

k 1

 

k

 

 

 

 

2

k 1

 

 

1

 

 

 

M

.

(4.18)

2

 

 

 

 

 

 

 

и плотности торможения от числа

Рис. 4.3. Зависимость теплофизических параметров от числа Маха

Расчёт по формулам (4.17) и (4.18) показывает, что при M 0,2 (для воздуха при 15 ºС это соответствует скорости 68 м/с) сжимаемость газа приводит к поправкам в плотности торможения на 2 %, в давлении и температуре – порядка 1 %. Ввиду малости этих поправок ими пренебрегают, считая газ при малых скоростях несжимаемой жидкостью. В задачах, не требующих высокой точно-

124

сти решения, можно считать газ несжимаемым и при больших числах M (порядка 0,3).

Выражения (4.17) и (4.18) являются ещё одной формой записи уравнения

(4.6).

4.5. Истечение газа из резервуара. Максимальная и критическая скорости

Исследуем истечение газа из резервуара, где он находился под давлением p0 , в среду с давлением p .

Применяя к сечениям струйки газа в резервуаре, где скорость близка к нулю, и в сжатом сечении уравнение энергии в форме (4.6), имеем

 

v2

 

 

 

k

 

 

p

 

 

k

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

,

 

 

2

k

1

k 1 0

 

 

 

 

 

 

 

 

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

2k

 

 

p0

 

 

 

 

p

 

 

0

 

 

 

 

 

1

 

 

 

.

 

k

 

0

 

p0

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

Выражая отношение

0

через отношение

p

с помощью уравнения

 

p0

 

 

 

адиабаты (4.2) и используя уравнение состояния, получим формулу Сен-Венана и Ванцеля для скорости адиабатного истечения газа

 

 

 

 

 

 

 

k 1

 

 

 

2k

 

k

 

 

 

 

 

 

 

v

RT0

1

 

p

 

 

 

.

(4.19)

 

 

 

 

 

 

k 1

 

 

p0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если газ вытекает в пустоту

p 0 ,

то достигается максимальная ско-

рость истечения

v

 

2kRT0

 

a

 

2

 

.

(4.20)

 

 

max

k 1

0

 

k 1

 

 

 

 

 

 

 

В частности, если в пустоту вытекает воздух при температуре 15 ºС, то vmax 760 м/с . При постепенном уменьшении давления в среде, в которую вы-

текает газ, начиная от p p0 , согласно формуле (4.20) растёт скорость истече-

ния (рис. 4.4).

125

Рис. 4.4. Изменение линии скорости от внешнего давления

Возрастание скорости в соответствии с уравнением энергии в форме (4.9) приводит к уменьшению местной скорости звука a . Наконец, при достаточно малом давлении среды

 

 

 

 

 

k

 

 

 

 

2

 

 

 

 

 

 

 

k 1

 

 

p pкр

p0

 

 

 

,

(4.21)

 

 

 

k 1

 

 

 

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

тических значений, определяемых формулами

 

 

 

 

1

 

 

 

 

 

 

 

 

2

 

 

k 1

 

; Tкр T0

2

 

 

 

кр

0

 

 

 

 

 

.

(4.22)

 

 

k 1

 

k 1

 

 

 

 

Величину критической скорости легко определить из уравнения энергии в форме (4.9) или (4.10), если принять v a aкр , то

a

 

2k

 

p0

 

a

 

 

2

 

.

 

(4.23)

 

 

 

 

 

кр

k 1

0

 

 

k 1

 

 

 

 

 

0

 

 

 

 

В частности, для воздуха, имеющего температуру 15

ºС,

aкр 0,91 м/с ,

a0 310 м/с . При дальнейшем уменьшении противодавления

p

скорость исте-

чения остаётся неизменной и равной aкр .

 

 

 

 

 

 

 

Постоянство скорости (и расхода) при

p pкр можно объяснить следую-

щим образом. Представим себе (рис. 4.5), что газ вытекает из резервуара 1 в вакуумную камеру 2 через трубу, давление в которой регулируется краном K . При p pкр скорость v aкр и при открытии крана волны разрежения от него,

распространяясь навстречу струе, соответственно увеличивают скорость исте-

126

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

не меняет скорость течения v aкр .

Рис. 4.5. Схема течения газа

Величина критической скорости aкр остаётся постоянной вдоль струйки.

Поэтому удобно измерять скорость течения в долях этой величины; так вводится безразмерная скорость газа, иногда называемая коэффициентом скорости

 

v

.

(4.24)

 

 

aкр

 

Здесь в знаменателе безразмерная скорость – величина постоянная вдоль струйки, тогда как в выражении (4.4) знаменатель является переменной величиной. При v aкр имеем M 1. Связь между этими величинами очевидна из

отношения

M

 

aкр

 

a

 

2

 

 

 

 

0

 

 

 

.

 

a

a

 

k 1

 

 

 

 

Используя уравнение энергии в форме (4.10), получим

M

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 1

 

k 1

 

2

 

 

 

 

 

2

 

2

 

 

 

 

Газодинамические функции для воздуха k 1, 4

 

 

 

k 1

 

 

 

 

 

 

 

 

k 1

 

 

 

 

k

 

1

 

2

 

 

 

 

2

 

k 1

 

 

 

 

 

 

 

, 1

 

 

 

 

 

 

 

,

k 1

 

 

k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 1

 

 

 

 

 

 

1

 

k 1

1

 

 

 

 

 

 

 

 

 

2

 

k 1

 

k 1

 

.

1

 

 

 

 

 

 

 

 

 

, q

 

 

 

 

 

 

k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

127

Переходя в формулах для параметров торможения газа (4.17) и (4.18) от числа M к , получим соотношения

 

 

 

k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 1

 

 

 

k

 

 

T

1

 

2

;

 

p

 

 

 

 

 

2

 

k 1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

;

T0

k 1

 

 

 

p0

 

k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.25)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

0

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

Формулы (4.25) дают изменение параметров газа вдоль струйки в зависимости от скорости. Они носят название газодинамических функций. Их численные значения для различных k и (или M ) сведены в таблицы и графики газодинамических функций. В частности, величины газодинамических функций для воздуха k 1, 4 даны в таблице. В таблицах газодинамических функций

даются также значения приведённого расхода q (отношение удельного рас-

хода массы газа в произвольном сечении струйки к расходу в критическом сечении. Использование таблиц газодинамических функций облегчает расчёты.

 

 

 

 

 

Таблица

 

 

 

 

q

M

 

1

2

3

4

5

6

 

0,00

1,0000

1,0000

1,0000

0,0000

0,0000

 

0,05

0,9996

0,9986

0,9990

0,0788

0,0457

 

0,10

0,9983

0,9942

0,9959

0,1571

0,0914

 

0,15

0,9963

0,9870

0,9907

0,2344

0,1372

 

0,20

0,9933

0,9768

0,9834

0,3102

0,1830

 

0,25

0,9896

0,9640

0,9742

0,3842

0,2290

 

0.30

0,9850

0,9485

0,9630

0,4557

0,2760

 

0,35

0,9796

0,9303

0,9497

0,5243

0,3228

 

0,40

0,9733

0,9097

0,9346

0,5897

0,3701

 

0,45

0,9663

0,8868

0,9178

0,6515

0,4179

 

0,50

0,9583

0,8616

0,8991

0,7091

0,4663

 

0,55

0,9496

0,8344

0,8787

0,7623

0,5152

 

0,60

0,9400

0,8053

0,8567

0,8109

0,5649

 

0,65

0,9296

0,7745

0,8332

0,8543

0,6154

 

0,70

0,9183

0,7422

0,8082

0,8924

0,6668

 

0,75

0,9063

0,7086

0,7819

0,9250

0,7192

 

0,80

0,8933

0,6738

0,7543

0,9518

0,7727

 

0,85

0,8796

0,6382

0,7256

0,9729

0,8274

 

0,90

0,8650

0,6019

0,6959

0,9879

0,8833

 

0,95

0,8496

0,5653

0,6653

0,9970

0,9409

 

128

Окончание таблицы

1

2

3

4

5

6

1,00

0,8333

0,5283

0,6340

1,0000

1,0000

1,05

0,8163

0,4913

0,6019

0,9969

1,0609

1,10

0,7983

0,4546

0,5694

0,9880

1,1239

1,15

0,7796

0,4184

0,5366

0,9735

1,1890

1,20

0,7600

0,3827

0,5035

0,9531

1,2566

1,25

0,7396

0,3479

0,4704

0,9275

1,3268

1,30

0,7183

0,3142

0,4374

0,8969

1,4002

1,35

0,6962

0,2816

0,4045

0,8614

1,4769

1,40

0,6733

0,2505

0,3720

0,8216

0,5575

1,45

0,6496

0,2209

003401

0,7778

1 ,6423

1,50

0,2250

0,1930

0,3088

0,3707

1,3721

1,55

0,5996

0,1669

0,2784

0,6807

1 ,8273

1,60

0,5733

0,1427

0,2489

0,6282

1,9290

1,65

0,5463

0,1205

0,2205

0,5740

2,0380

1,70

0,5183

0,1003

0,1934

0,5187

2,1555

1,75

0,4896

0,0821

0,1677

0,4630

2,2831

1,80

0,4600

0,0660

0,1435

0,4075

2,4227

1,85

0,4296

0,0520

0,1210

0,3530

2,5766

1,90

0,3983

0,0399

0, 1002

0,3002

2,7481

1,95

0,3662

0,0297

0,0812

0,2497

2,9414

2,00

0,3333

0,0214

0,0642

0,2024

3,1622

2,05

0,2996

0,0147

0,0491

0,1588

3,4190

2,10

0,2650

0,0096

0,0361

0,1198

3,7240

2,15

0,2296

0,0058

0,0253

0,0857

4,0961

2,20

0,1933

0,0032

0,0164

0,0570

4,5674

2,25

0,1563

0,00151

0,00966

0,0343

5,1958

2,35

0,1183

0,00057

0,00482

0,0175

6,1033

2,40

0,0796

0,00014

0,00170

0,0063

7,6053

2,449

0,0400

0,128 10 4

0,00032

0,0012

10,957

 

0

0

0

0

4.6. Одномерное течение газа с трением

4.6.1.Изотермическое течение в трубах

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

Будем учитывать потери напора на трение вдоль трубы по формуле Дарси, тогда потери на участке трубы длиной dx составят

129