новая папка 1 / 292807
.pdfАДИАБАТИЧЕСКОЕ ДВИЖЕНИЕ
СЖИМАЕМОГО ГАЗА
Методические указания к курсовой работе
Министерство образования и науки Российской Федерации
Уральский федеральный университет имени первого Президента России Б. Н. Ельцина
А. С. Носков, Н. В. Служеникина
АДИАБАТИЧЕСКОЕ ДВИЖЕНИЕ СЖИМАЕМОГО ГАЗА
Методические указания к курсовой работе
Рекомендовано методическим советом УрФУ для студентов, обучающихся по направлениям: 270000 – Строительство, 280700 – Техносферная безопасность,
140100 – Теплоэнергетика и теплотехника,
141100 – Энергетическое машиностроение
Екатеринбург Издательство Уральского университета
2014
УДК 533.21:536-34(075.8)
Составители: А. С. Носков, Н. В. Служеникина Научный редактор – канд. техн. наук, доц. А. В. Некрасов
Адиабатическое движение сжимаемого газа : методические указания / сост. А. С. Носков, Н. В. Служеникина. – Екатеринбург : Изд-во Урал. ун-та, 2014. – 16 с.
Методические указания предназначены для выполнения курсовой работы по дисциплинам «Механика жидкости и газа» и «Гидрогазодинамика». Данные методические указания могут быть использованы и при изучении других дисциплин, содержащих раздел «Газодинамика».
В методические указания входят следующие темы:
–«Адиабатическое движение с трением сжимаемого газа в цилиндрической трубе»,
–«Адиабатическое движение сверхзвукового потока газа в цилиндрической трубе со скачком уплотнения».
Для выполнения заданий приведены варианты исходных данных и
методика расчета.
© Уральский федеральный университет, 2014
2
ВВЕДЕНИЕ
Курсовая работа включает в себя задания на расчет движения сжимаемого газа в цилиндрической трубе с трением. Номер задания студент выбирает по указанию преподавателя. Результаты курсовой работы должны быть изложены в виде пояснительной записки. Графики следует выполнять с обязательным указанием масштаба по осям координат. Под рисунком должна быть подпись, поясняющая его содержание. В записке студент должен указать ход расчета и его результаты. Все массовые однотипные вычисления сводятся в таблицу.
КОНТРОЛЬНЫЕ ЗАДАНИЯ
Задание 1 Адиабатическое движение с трением сжимаемого газа
в цилиндрической трубе
(Схема движения газа в трубе приведена на рисунке)
Газ через цилиндрическую трубу истекает в пространство с давлением pн. Движение газа в трубе адиабатическое с трением. Заданы: длина L и диаметр трубы d; материал, из которого изготовлена труба (для определения величины шероховатости);
температура торможения газа T0.
Требуется определить давление торможения p0в на входе в трубу для подачи заданного массового расхода газа G и изменения параметров газового потока по длине трубы (скорости V, температуры T, давления p, плотности ρ).
Варианты численных значений исходных данных приведены в таблице 1. Для новой стальной трубы величину абсолютной шероховатости принять = 0,04 мм, для медной трубы = 0,01 мм.
b |
|
L |
a |
|
|
|
|
O |
d |
|
p |
|
|
|
x |
b |
|
|
a |
Рис. Схема движения газа в трубе:
a-a – сечение на выходе из трубы в пространство с давлением pн, b-b – сечение на входе в трубу
3
Таблица 1 Варианты численных значений исходных данных к заданию 1
№ вар. |
Вид газа |
Труба |
L, м |
d, мм |
pн, кПа |
T0, К |
1 |
воздух |
новая стальная |
1,0 |
10 |
100 |
300 |
2 |
воздух |
новая стальная |
1,5 |
20 |
120 |
300 |
3 |
воздух |
новая стальная |
2,0 |
30 |
100 |
300 |
4 |
воздух |
новая стальная |
5,0 |
40 |
100 |
300 |
5 |
воздух |
новая стальная |
5,0 |
10 |
250 |
300 |
6 |
воздух |
новая стальная |
2,0 |
15 |
140 |
300 |
7 |
воздух |
новая стальная |
1,5 |
5 |
200 |
500 |
8 |
воздух |
новая стальная |
3,0 |
20 |
170 |
400 |
9 |
воздух |
новая стальная |
6,0 |
35 |
300 |
1000 |
10 |
воздух |
новая стальная |
4,0 |
20 |
120 |
600 |
11 |
воздух |
новая стальная |
3,5 |
12 |
150 |
300 |
12 |
воздух |
медная новая |
1,4 |
6 |
110 |
700 |
13 |
воздух |
медная новая |
1,0 |
8 |
190 |
800 |
14 |
водород |
новая стальная |
1,5 |
18 |
120 |
300 |
15 |
водород |
новая стальная |
2,0 |
32 |
100 |
400 |
16 |
водород |
новая стальная |
1,2 |
100 |
100 |
300 |
17 |
водород |
медная новая |
2,0 |
15 |
150 |
300 |
18 |
водород |
медная новая |
1,0 |
7 |
120 |
400 |
19 |
водород |
медная новая |
0,7 |
5 |
200 |
300 |
20 |
метан |
новая стальная |
1,0 |
5 |
200 |
300 |
21 |
метан |
новая стальная |
1,5 |
10 |
140 |
400 |
22 |
метан |
новая стальная |
1,2 |
7 |
250 |
600 |
23 |
метан |
медная новая |
2,0 |
16 |
120 |
300 |
24 |
метан |
медная новая |
3,0 |
18 |
160 |
500 |
25 |
метан |
медная новая |
0,5 |
6 |
110 |
700 |
26 |
воздух |
новая стальная |
3,0 |
10 |
80 |
300 |
27 |
воздух |
новая стальная |
2,4 |
12 |
420 |
400 |
28 |
воздух |
новая стальная |
1,8 |
20 |
90 |
500 |
29 |
воздух |
медная новая |
0,9 |
15 |
180 |
300 |
30 |
воздух |
медная новая |
2,1 |
16 |
330 |
350 |
31 |
воздух |
медная новая |
1,6 |
8 |
270 |
700 |
32 |
кислород |
медная новая |
1,2 |
10 |
80 |
300 |
33 |
кислород |
медная новая |
2,6 |
12 |
310 |
300 |
34 |
кислород |
медная новая |
1,8 |
8 |
150 |
300 |
35 |
кислород |
медная новая |
0,9 |
12 |
250 |
300 |
|
|
4 |
|
|
|
|
Задание 2 Адиабатическое движение сверхзвукового потока газа
в цилиндрической трубе со скачком уплотнения
(Схема движения газа в трубе приведена на рисунке, с. 3)
Поток газа поступает в стальную цилиндрическую трубу со сверхзвуковой скоростью Vb, движется по ней адиабатически при наличии трения и затем истекает во внешнее пространство, давление в котором равно pн. Задан коэффициент скорости на входе в трубу λb, длина L и диаметр трубы d, температура торможения газа T0.
Требуется определить изменение скорости V, давления p, давления торможения p0, плотности ρ и температуры T газового потока по длине трубы.
Варианты численных значений исходных данных приведены в таблице 2. Величину абсолютной шероховатости для новой стальной трубы принять = 0,04 мм.
Таблица 2 Варианты численных значений исходных данных к заданию 2
№ |
λb |
L, см |
d, мм |
Вид газа |
pн, кПа |
T0, K |
1 |
1,5 |
7,2 |
8 |
воздух |
100 |
300 |
2 |
1,4 |
8,8 |
12 |
воздух |
160 |
300 |
3 |
1,6 |
32 |
20 |
воздух |
100 |
500 |
4 |
1,75 |
90 |
30 |
воздух |
200 |
400 |
5 |
1,5 |
22 |
18 |
воздух |
140 |
600 |
6 |
2,0 |
32 |
10 |
воздух |
100 |
400 |
7 |
2,4 |
110 |
16 |
воздух |
250 |
300 |
8 |
2,2 |
147 |
24 |
воздух |
220 |
300 |
9 |
1,9 |
119 |
32 |
воздух |
130 |
600 |
10 |
2,1 |
142 |
28 |
воздух |
270 |
800 |
11 |
2,3 |
75 |
14 |
воздух |
80 |
700 |
12 |
2,0 |
50 |
14 |
водород |
300 |
300 |
13 |
1,8 |
60 |
20 |
водород |
260 |
400 |
14 |
1,6 |
70 |
36 |
водород |
200 |
600 |
15 |
2,1 |
28 |
8 |
водород |
240 |
500 |
16 |
2,3 |
94 |
16 |
водород |
120 |
700 |
17 |
2,4 |
195 |
25 |
водород |
90 |
300 |
18 |
2,2 |
185 |
30 |
водород |
170 |
400 |
|
|
|
|
5 |
|
|
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Окончание табл. 2
№ |
λb |
L, см |
d, мм |
Вид газа |
pн, кПа |
T0, K |
19 |
1,5 |
61 |
40 |
метан |
170 |
300 |
20 |
1,7 |
75 |
30 |
метан |
220 |
500 |
21 |
2,3 |
130 |
20 |
метан |
140 |
300 |
22 |
1,8 |
43 |
16 |
метан |
110 |
500 |
23 |
2,0 |
46 |
12 |
метан |
180 |
400 |
24 |
2,4 |
200 |
26 |
метан |
90 |
250 |
25 |
2,1 |
174 |
33 |
метан |
120 |
400 |
МЕТОДИЧЕСКИЕ УКАЗАНИЯ К ВЫПОЛНЕНИЮ КУРСОВОЙ РАБОТЫ
При выполнении курсовой работы используются стандартные газодинамические функции [1] π(λ), ε(λ), τ(λ), q(λ), y(λ), зависящие от коэффициента скорости λ = V / Vк (их аналитические выражения приведены в Приложении), V – скорость движения газа, Vk – критическая скорость газа:
Vk |
2kRT0 |
, |
(1) |
||
(k 1)M |
|||||
|
|
|
|
||
где T0 – температура торможения газового потока, R = 8314 |
Дж |
||||
моль К |
|
– универсальная газовая постоянная, М – молекулярная масса газа,
k Cp – отношение теплоемкости газа при постоянном давлении (Сp)
Cv
и постоянном объеме (Сv). Для двухатомных газов (воздух, водород, кислород) k = 1,40, для многоатомных (метан) k = 1,33. Численное значение газодинамических функций для данного k приведено в
«Прикладной газовой динамике» .
Расход газа G, текущего по трубе, с помощью газодинамических функций q и y выражается в виде:
G m |
р0 q |
, |
(2) |
|
|||
|
T0 |
|
Абрамович Г. Н. Прикладная газовая динамика. 5-е изд., перераб. и доп. М. : Наука, 1991.
6
G m |
р y |
, |
(3) |
|
|||
|
T0 |
|
где p, p0 – давление и давление торможения газа, ω – площадь поперечного сечения трубы,
|
|
|
|
2 |
|
k 1 |
|
|
|
|
M , |
|
|
|
k 1 |
|
|
|
|
m N |
N |
k |
|
|
|
. |
(4) |
||
|
|
||||||||
|
R |
|
k 1 |
|
|
|
|
Критический расход газа Gk (расход при λ = 1) определяется по формуле
Gk m |
pн |
y 1 , |
(5) |
|
|||
|
T0 |
|
где pH – давление во внешней среде (в которую истекает газ из трубы). В формулах (2), (3), (5) давление измеряется в Па, ω – в м2.
Закон сохранения массы для движения газового потока в трубе с постоянной площадью поперечного сечения позволяет, используя формулы (2) и (3), записывать соотношение между давлением и коэффициентом скорости в любых двух сечениях трубы 1 и 2:
q(λ1) p01 = q(λ2) p02 = y(λ1) p1 = y(λ2) p2. |
(6) |
При адиабатическом течении газа с трением в трубе диаметром
d коэффициент скорости изменяется от значения λ1 |
в первом сечении |
|||||
до λ2 во втором сечении, находящемся от первого на расстоянии l |
||||||
вниз по потоку, по уравнению: |
|
|
|
|
||
φ(λ1) = φ(λ2) + χ , |
(7) |
|||||
φ – газодинамическая функция |
|
|
|
|
||
φ(λ) = λ–2 + 2 lnλ , |
(8) |
|||||
χ – приведенная длина трубы |
|
|
|
|
||
|
2k |
|
|
l |
, |
(9) |
k 1 |
|
|||||
|
|
d |
|
где ζ – коэффициент сопротивления трения.
Коэффициент сопротивления трения вычисляется по формуле Альтшуля:
|
|
|
68 0,25 |
(10) |
ζ 0,11 |
|
. |
||
d |
|
Re |
|
|
|
|
7 |
|
|
Число Рейнольдса вычисляем по формуле |
|
Re Gdμω . |
(11) |
Из формулы (11) следует, что Re изменяется при движении газа по трубе постоянного сечения только за счет зависимости динамической вязкости газа μ от температуры газового потока Т. Эта зависимость достаточно слабая, и ей можно пренебречь. Поэтому значения μ для подсчета значений Re можно выбирать при температуре торможения газового потока. Значения μ определяются по формуле Сатерленда:
μ |
|
T |
3/ 2 |
Tн С |
(12) |
||
|
|
||||||
μ |
|
|
|||||
H |
T |
|
T С . |
||||
|
|
H |
|
|
Значения μH и С для различных газов приведены в таблице 3.
Таблица 3
Значения μH |
и С в формуле (12) |
||
Газ |
|
μH ∙ 105 , Па ∙ с |
C, К |
воздух |
|
1,72 |
122 |
водород |
|
0,83 |
83 |
кислород |
|
1,92 |
138 |
метан |
|
1,04 |
198 |
ЗАДАНИЕ 1
При адиабатическом движении газа в трубе постоянного сечения дозвуковой поток, при наличии трения, ускоряется. Поэтому максимальная скорость течения достигается газовым потоком в выходном сечении. Если G < Gk , то коэффициент скорости на выходе из трубы λa < 1, при G = Gk λa = 1. Если давление торможения на
входе в трубу p0b таково, что G > Gk , то λa = 1.
При этом расход газа увеличивается за счет увеличения плотности газового потока на выходе из трубы, вследствие того, что статическое давление на выходе из трубы больше давления во внешнем пространстве pH. Если G > Gk , то из (3) и (5) следует:
pa pн |
G |
. |
(13) |
G |
|||
|
k |
|
|
8 |
|
|
Характерным для течения газа в трубе при G > Gk является то, что распределение скорости и температуры по длине трубы совпадает с распределением этих параметров по длине трубы при G = Gk.
ЗАДАНИЕ 2
При адиабатическом течении газа с трением в трубе постоянного сечения сверхзвуковой поток тормозится. Если приведенная длина трубы χ больше
кр b 1, |
(14) |
то непрерывное торможение сверхзвукового потока невозможно, т. к. при χ = χкр λa достигнет значения, равного 1. Тогда на некотором расстоянии от входа в трубу xс поток затормозится до λ′ > 1, и в этом сечении образуется скачок уплотнения, после которого поток –
дозвуковой 1 1.
С другой стороны, если приведенная длина трубы χ больше
|
|
|
1 |
|
1, |
(15) |
|
|
|
|
|||
|
|
|||||
|
max |
|
|
|
|
|
|
|
|
b |
|
|
то скачок уплотнения должен произойти во входном участке трубы,
т. к. при χ = χmax λ′ = λb, поэтому сверхзвуковое течение газа в трубе возможно лишь при χ < χmax.
Если χкр < χ < χmax, то положение скачка уплотнения зависит от величины давления торможения на входе в трубу. Коэффициент
скорости перед скачком уплотнения λ′ определяется с помощью нахождения значения вспомогательной функции
|
b a , |
(16) |
||
где |
|
1 |
. |
(17) |
|
|
|||
|
|
|
|
|
Нахождение значения λ по известному значению Φ производится численно.
Минимальное значение λa,min (pob)min определится из условия, что скачок уплотнения находится во входном сечении трубы (χc = 0,
λ′ = λb).
Тогда из (6), (16), (17) cледует, что
9