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

книги из ГПНТБ / Криогенные поршневые детандеры

..pdf
Скачиваний:
45
Добавлен:
24.10.2023
Размер:
10.95 Mб
Скачать

Глава III. ТЕОРЕТИЧЕСКИЕ ОСНОВЫ ИССЛЕДОВАНИЯ РАБОЧИХ ПРОЦЕССОВ КРИОГЕННЫХ ПОРШНЕВЫХ ДЕТАНДЕРОВ С РЕАЛЬНЫМ РАБОЧИМ ТЕЛОМ

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

Это — процессы

определяющие. При их анализе для криоген­

ных детандеров

необходимо найти такие соотношения между

параметрами, Которые в доступной форме учитывали бы реаль­ ность рабочего тела.

Физическая обстановка в цилиндре расширительной машины существенно усложнена процессами массо- и теплообмена. Ана­ лиз этих явлений с учетом специфических особенностей низко­ температурных машин должен давать количественную и каче­ ственную их оценку и устанавливать степень их влияния на к. п. д.

1.Изменение температуры реального газа

вравновесных адиабатном и политропном

процессах

Теоретически процессы внутреннего расширения и обратного сжатия могут рассматриваться равновесными (обратимыми). В зависимости от условий они могут быть адиабатными или политропными. Уравнения равновесного адиабатного (изоэнтропного) процесса в переменных Т v и Т р имеют вид:

Исходя из

обобщенного

уравнения

состояния

реальных

газов:

,

. В

 

 

 

 

 

 

pv

С ,

D

Е

Н

,

(47)

—— = z = 1

-\------

-------------------------- v*

v*

vs

I»8

RT

 

v

 

 

где z — коэффициент сжимаемости;

В,

С, D,

Е, Н — вириаль-

ные коэффициенты, являющиеся функциями температуры.

43

Определим

производные

в уравнение (46), найдем

\ дТ

 

Л * - ‘).Гг + г ^

 

[

V дТ ]v

у х

I

T { k - 1) z + т

( f ) . = “• =

pkx

 

здесь

 

и /"-^-Л и, подставив их

/ о \ дТ ) Р

(48)

дг

T (k — l) z f OiJ

~дТ

pkx

 

щ -

l f -)

- ^ = ! L

( - £

) .

(49)

 

 

\ dp h

pkXx

\ dTv i

ijp

 

R

R

выражены

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

Отношения —

и

 

CV

Майера:

+(-Sr)J+'-f(-H(-£).}-**- (50)

так:

k — 1

R

k — \

k x

CV

X

Величина as полностью определена значениями p и T, так как все величины, входящие в уравнение (48), являются функ­ циями давления и температуры. При увеличении температуры и уменьшении давления <xs, как правило, увеличивается.

Проинтегрировать уравнения (48) можно только в том слу­ чае, если известны функциональные связи z = z(p, Т) и k =

=k(p, Т). Например, при постоянных z и k, х = г. Тогда

к— 1

(51)

т » - т к = ти

При постоянном коэффициенте сжимаемости г, не равном единице, газ не является в обычном смысле идеальным, но величина дифференциального эффекта дросселирования обра­ щается в нуль [см. уравнение (49)]. В действительности значе­ ния z и k в различных областях состояний несомненно различны.

44

Вместе с тем для большинства газов, используемых в криоген­

ной технике (кроме водорода), изменение температуры

в изо-

энтропных

процессах, с

точностью

вполне достаточной

для

практических

расчетов,

может быть вычислено по форму­

лам (51) при постоянном значении k.

Весьма существенно,

что

уравнения

(51)

действительны для этих реальных газов

в ши­

рокой области состояний, в том числе и вблизи пограничной кривой. Эта замечательная особенность проиллюстрирована табл. 6, в которой для разных газов приведены значения вели­

чины р, представляющей отношение действительной

разности

температур (Гн— TK)R в изоэнтропном

процессе, определенной

по тепловым диаграммам,

к

разности температур (Гн—

Т'к)ур.(51), вычисленной по уравнению (51):

 

 

 

ДТрасч

'

н

'к;ур.(51)

 

 

Учитывая точность, с которой можно определить

изменение

температур в процессах s =

const

 

по

тепловым

диаграммам,

констатируем, что данные табл. 6

убедительно

подтверждают

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

уравнениям

(51) для

реальных газов в широкой области состояний.

 

 

Однако исключением является водород. Если для гелия, неона, азота, воздуха вероятная погрешность определения тем­ ператур по уравнению (51) не превышает 2—3%, то для водо­ рода при постоянном значении k = 1,407, она может достигать 20—30% (в области температур ниже 130° К).

С физической стороны эта особенность водорода объясняется следующим. Известно, что все виды движений молекул (посту­ пательное, вращательное и колебательное) вносят определен­ ный вклад в величину внутренней энергии, а следовательно, и в теплоемкость. На каждую степень свободы при поступатель­ ном движении центра массы и на каждую «вращательную» степень свободы приходится кинетическая энергия, равная llzkT. Полная энергия, приходящаяся на одну колебательную степень свободы, равна kT. Однако в соответствии с квантовой теорией молекула может находиться лишь в стационарных со­ стояниях, соответствующих определенным значениям энергии. Молекула может воспринимать энергию только определенными порциями (квантами) и при этом она переходит из одного ста­ ционарного состояния в другое. Для вращательного и колеба­ тельного движений дискретные уровни энергий вполне опреде­ ленны. Поэтому могут быть определены так называемые харак­ теристические температуры Гвр и Ткол (температуры «вырожде­ ния»), ниже которых вклад в теплоемкость вращательного или

колебательного движений не существенен (табл. 7).

 

Для

водорода характеристическая

температура Твр весьма

высока. Поэтому при температурах,

близких к 85° К и

более

низких,

теплоемкость cv водорода

уже практически

равна

45

Т а б л и ц а 6

гн к

300

250

200

180

150

120

гн к

300

250

200

180

150

120

Значения отношений разностей действительных и расчетных температур (р), вычисленные по уравнению (52)

 

Воздух. рк =

1

ат (абс.),

k = const = 1,4

 

 

 

 

 

 

Рн , ат (абс.)

 

 

 

80

60

50

 

30

15

10

5

3

 

 

 

 

 

 

 

 

1,007

1,009

1,008

1,001

0,996

1,005

1,008

1,02

_

_

._

 

1,007

1,001

1,003

0,998

1,02

_

_

_

 

1,01

1,01

1,01

1,02

_

_

_

 

1,013

1,013

1,018

1,03

_

_

_

 

1,02

1,03

 

1,03

 

 

 

 

 

 

 

 

 

Азот. рк =

1

ат (абс.),

k = const =

1,4

 

 

 

 

 

 

Рн , ат (абс.)

 

 

 

80

60

50

 

30

15

10

5

2,5

1,008

1,008

1,008

1,001

0,995

1,002

1,005

0,998

___

___

 

1,007

1,002

1,0

0,998

1,003

___

 

1,011

1,0

1,01

1,01

_

____

 

1,015

1,01

1,02

1,01

_

___

 

1,01

1,01

 

1,02

Неон. рк = 1 ат (абс.), k = const = 1,68

рн, ат (абс.)

гн к

200

80

50

30

15

10

6

4

 

300

0,99

0,995

0,99

0,99

1,0

0,99

0,99

0,99

250

0,99

0,99

0,99

0,99

0,99

0,99

1,0

0,99

200

0,99

0,99

1,0

1,0

1,0

1,0

1,0

180

0,99

1,0

0,997

1,0

0,99

1,0

1.0

150

 

0,99

0,99

0,99

0,99

1,0

1.0

------

100

___

___

___

0,98

0,99

0,99

1,0

77

___

___

___

0,98

0,99

0,99

47

0,99

 

 

Гелий.

рк = 1

ат (абс.),

k = const = 1,66

 

 

г„К

 

 

 

р н ,

ат (абс.)

 

 

 

200

100

60

40

30

20

10

5

 

300

1,005

1,007

1,007

1,01

1,01

1,01

1,01

1,001

250

1,006

1,01

1,008

1,01

1,01

1,0

1,0

1,01

200

1,008

1,01

1,005

1,01

1,01

1,01

1,0

1,005

150

1,008

1,01

1,01

1,01

1,0

1,0

1,01

0,99

100

1,007

1,01

1,02

1,02

1,02

1,02

0,99

1,01

77

1,003

1,01

1,01

1,02

1,02

1,02

1,02

1,02

50

1,01

1,007

1,008

1,008

1,01

1,02

1,03

1,02

20

0,98

0,99

0,99

0,99

0,99

0,99

 

 

 

 

 

 

 

 

10

— •

1,007

46

Т а б л и ц а 7

Су

Характеристические температуры двухатомных газов

R

Газ

Г - К

г

°к

 

 

 

 

вр

 

 

 

 

 

н2

85,4

6100

 

 

 

о2

2,86

3340

 

 

 

n 2

2,07

2230

Рис. 7. Молярная

теплоемкость

 

 

 

 

некоторых газов

 

теплоемкости cv

одноатомного газа

(см.

рис. 7). Для дру­

гих многоатомных газов

это «вырождение»

происходит при

существенно более низкой температуре.

 

 

Чтобы избежать больших ошибок при вычислении изме­

нения температур водорода по уравнению

(51),

необходимо

значения показателя Пуассона k выбирать как средние по про­ цессу s = const. В этом случае достигается необходимая точ­ ность.

В табл. 8 приведены значения &ср и р для нормального водо­ рода иН2, равновесный состав которого: 25% />Н2 и 75% оН2. В детандерах установок для ожижения водорода может расши­ ряться уже практически полностью конвертированный водород: >95% рН2. Теплоемкости ср и cv пароводорода, ортоводорода и нормального водорода различны (табл. 9). Вероятные значения

показателя Пуассона для пароводорода (й)Рн2

в

реальном

газовом состоянии

можно приближенно найти

по

величине

этого показателя

для нормального водорода (6)„н2

и отно-

(^о)рН шению 7т~— - при определенной температуре (см. табл. 9):

'лН,

( * о ) 'pH, ( £ ) РН , ~ ( £ ) л Н г лН,

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

жет быть выражено через величину теплоемкости политропного процесса с:

dq = cdT.

(53)

47

Т а б л и ц а 8

Значения р, вычисленные по уравнению (52) для нормально-

ft + А

го водорода; р к = 1 ат (абс.); ft = (ftCp b =Const = ~

— ~

гн к

300

250

200

180

150

90

33

 

 

 

 

 

, а т (абс.)

 

 

 

 

 

100

 

80

|

60

40

|

20

|

10

2

ftcp = 1,527

1,520

 

1,501

1,496

 

1,468

1,452

1,412

II О

00 00

0,989

 

0,994

0,980

 

1,0

0,972

0,985

1,552

 

1,547

 

1,54

1,524

 

1,497

1,475

1,427

0,988

 

0,995

 

0,993

0,985

 

1.0

«

1,0

0,980

1,59

 

1,59

 

1,57

1,559

 

1,535

1,52

1,457

0,995

 

0,993

 

1,0

1,0

 

0,997

«

1,0

« 1,0

1,61

 

1,607

 

1,593

1,579

 

1,546

1,548

1,472

0,995

 

0,995

 

0,998

0,999

 

0,990

0,99

« 0,999

1,66

 

1,65

 

1,636

1,623

 

1,592

1,573

1,494

0,991

 

0,990

 

0,993

0,998

 

0,995

1,01

« 1,0

 

 

 

1,757

 

1,742

1,680

1,62

 

 

 

1,02

 

1,02

1,101

0,995

 

 

 

 

 

 

 

 

 

 

 

 

 

1,766

1,72

1,65

 

 

 

 

 

 

 

1,03

1,02

1,01

 

 

 

 

 

 

 

1,822

1,75

1,65

 

 

 

 

 

 

 

 

1,03

0,96

ftcp = 1 > 7 8

II CCL

о

СО

Т а б л и ц а

9

Теплоемкости и показатель Пуассона для водорода

 

в идеальном газовом состоянии

 

 

 

СР_

 

СУ

 

^0

 

 

 

R

 

R

 

<*о)рН,

Т К

 

 

 

 

 

п Н 2

рН г

п Н 2

pHг

лН,

pH,

(*0>лН,

 

 

300

3,4515

3,585

2,4513

2,585

1,409

1,387

0,986

200

3,2667

3,895

2,2657

2,895

1,441

1,347

0,935

150

3,0403

3,9130

3,0413

2,913

1,490

1,343

0,902

100

2,7073

3,2370

1,7093

2,237

1,585

1,447

0,913

90

2,6420

3,028

1,6430

2,028

1,610

1,492

0,927

77

2,5704

2,7812

1,5703

1,7812

1,638

1,561

0,952

70

2,5441

2,6763

1,5441

1,6763

1,650

1,597

0,968

60

2,5178

2,5713

1,5178

1,5713

1,658

1,638

0,989

50

2,5

1,5

1,5

1,5

1,667

1,667

1,0

40

2,5

1,5

1,5

1,5

1,667

1,667

1,0

30

2,5

1,5

1,5

1,5

1,667

1,667

1,0

48

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

Уравнение равновесного политропного процесса в перемен­ ных Т и р имеет вид

cdT = cpd T - T ( ^ d p ;

(54)

для реального газа, выполнив те же операции, что и при выводе уравнения (48), получим

& L - ^ - r 3 r ^ [ ‘+T&).]-

(55)

 

 

pkx

(-t)

 

 

 

 

 

 

 

Обозначив

 

 

 

 

 

 

 

1—

 

 

 

m = k

 

с р -

 

(56)

 

1—k ■

 

 

 

 

 

наидем

 

 

 

 

 

 

 

Т(т— 1) г + Т

( П \

(57)

V др Л

 

рт%

 

 

 

 

Интегрирование этого уравнения при постоянных т и z дает

т — 1

т

ТК

- (

\

— ) “ ;

 

 

Рк /

тп— 1

m

тн- т к = т ■.'-(to

Для количества тепла q получим

<7=

т~ к ср(ТК

Т„) = JUzzl. су(ТкТи).

 

k(m— 1)

ТП'— 1

(58)

(59)

При сделанных допущениях показатель политропы m полу­ чает физический смысл коэффициента, который определяет, ка­ кая часть введенного тепла перешла в увеличение внутренней

энергии ^Ди = —— qj и какая часть превращена в работу

1—k

газа

4 Заказ 1397

49

В действительности распределение тепла между Ли и Д/ на протяжении процесса не остается постоянным и поэтому опи­ сание реальных процессов уравнением политропы представ­ ляется прямой аппроксимацией. Простая математическая форма уравнения политропы (pvm = const), позволяющая в первом приближении описать сложные процессы, является решающим его преимуществом. Действительно, многолетний опыт приме­ нения этого уравнения показал, что кроме тепловых взаимодей­ ствий показателем т могут учитываться специфические свой­ ства реального газа, переменное его количество, присутствие таких трудно контролируемых факторов, как трение и натека­ ние газа, присутствие окислительных реакций и других необра­ тимых процессов. В таких случаях показатель т является про­ сто эмпирическим коэффициентом и никакого физического смысла не имеет, а само понятие политропного процесса не имеет ничего общего с равновесной политропой. При анализе реальных процессов надо иметь в виду, что если простота мате­ матической формы политропы является преимуществом, то неопределенность показателя т как функции многочисленных параметров — объективным недостатком.

2. Изменение температуры реального газа

внеравновесном адиабатном и политропном процессах

Процессы впуска и выхлопа, протекающие почти мгновенно при открытии клапанов являются неравновесными (необрати­ мыми). Наиболее вероятной теоретической моделью этих про­ цессов является модель неравновесного адиабатного процесса. Тепловые взаимодействия газа со стенками исключены по усло­ вию, и неравновесность обусловлена тем, что при изменении объема силы давления газа на контрольной поверхности систе­ мы не уравновешены силами противодавления.

Для анализа процесса удобно воспользоваться схемой рис. 8. Первоначальные параметры газа в цилиндре Гн и рн. Цилиндр герметически закрыт задвижкой с газовым поршнем. После того как задвижка будет освобождена, она начнет двигаться без трения в выпускной трубе. Давление газа на задвижку бу­ дет постепенно падать. В то же время на газовый поршень будет действовать постоянная сила противодавления, так как по условию газ вытекает в область постоянного давления рк. Когда в момент времени Тк давление в цилиндре достигнет величины рк, задвижка с газовым поршнем остановится. Спустя некоторое время после окончания процесса во всех частях системы установится равновесная температура Тк, которая мо­ жет быть определена. Действительно, для этого типично нерав-

50

Рис. 8. Мгновенное и равновесное распределение температур в процессе выхлопа

новесного процесса, исключая из рассмотрения явления, происходящие в газовом поршне, имеем (для 1 кг газа)

для реального газа

“ к— ын~ Р к К — v„y,

(60)

 

 

 

 

 

 

du = cvdT

dv\

(61)

поэтому при постоянных г и cv

 

 

u„— uK= cv(Tll— TK) = p llRz ( —s ~

(62)

учитывая (50), найдем

 

V

Рк

 

 

£—1

 

 

т

_Т

= Т

 

 

1н

* К

1

н

 

 

 

 

 

 

 

(63)

Т к

Ри

4*

51

Понижение температуры газа при одинаковых Тв, рн и рк в неравновесном адиабатном процессе меньше, чем в процессе

изоэнтропном. Следует обратить внимание

на

то,

что

при

Др = (рв — pK)-*-dp будут соблюдены

условия

квазиравнове­

сия. В этом случае дифференциальные

эффекты

процесса

вы­

хлопа и изоэнтропного процесса будут неразличимы:

 

 

Тяр ^ я 1) h

+ 4 a

r H )

J -

( б 4 )

Существенно следующее: при принятых выше допущениях все выражения для рассматриваемого процесса [уравнения (63)] можно легко получить из уравнения (64), если положить

а выхл = (

)

= c o n s t = ( а 5) н =

- ^ - . - ^

- .

( 6 5 )

\

dp

Jвыхл

Рн

*

 

Таким образом определяется математическая форма урав­ нения процесса выхлопа. Учитывая связь авыхл и as (65) и из­ ложенное в предыдущем параграфе, можно полагать, что точ­ ность вычисления температур реальных газов (кроме водорода) в процессе выхлопа по уравнениям (63) при постоянных z и 4 будет вполне удовлетворительной. Для водорода, как и при вы­ числении температур в равновесном адиабатном процессе, следует пользоваться средними значениями показателя Пуас­

сона. Исходя из уравнения

(65), нетрудно построить линии

«выхл = (a«)H= const для

процесса

выхлопа

в координатах

s — Т при заданных рв и Тв

(рис. 9).

Из этого

рисунка ясно,

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

Если процесс выхлопа рассматривать во времени, то неизбежно приходится учитывать «температурное расслоение» или возникновение градиента температуры в потоке выходяще­ го из цилиндра газа (см. рис. 8). В течение всего процесса в оди­ наковые моменты времени ti, п , тк, температуры в выпускной

трубе Г1в; Г2в; Т кв будут выше, чем температура газа

в ци­

линдре Г 1ц; Г2ц; Т кц. Температура газа, находящегося

в ци­

линдре в момент окончания процесса выхлопа 7’кц может быть

в первом приближении определена по уравнению адиабаты Пуассона (51), т. е. как для равновесного адиабатного процес­ са. Это верно потому, что на контрольной поверхности подси­ стемы, содержащей массу газа, остающуюся в цилиндре к мо­ менту окончания выпуска (контур а на рис. 8), можно считать выполняющимися условия квазиравновесия. Средняя темпера­ тура вышедшего газа к моменту окончания выхлопа (Тв)ср может быть просто вычислена из условия, что после полного смешения всех порций газа, оставшихся и вышедших из ци-

52