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

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

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

Рис. 9. Кривые сиыхл = а .н = const в координатах Т — s для воздуха при Т„ = 267 К и различных рв

линдра, их температура должна быть равной Тн. Так, при сде­ ланных допущениях

(Тв)

______Рн

Рк______

( 66)

 

 

Ря~рк(^)

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

I + (*-!) —

___________ Ри_

(67)

k

 

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

53

Щ*/молъ

16000

1S500

15000

m oo

m oo

15500

13000

12500

12000

11500

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

Рк

т

k T RX.TН

( 68)

=

Рн

1К--

1 + W,

 

 

 

 

во

Рассматривая процесс

впуска

времени,

необходимо

также,

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

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

(68).

В диаграмме Т — s условные

линии для неравновесных

адиа-

 

 

 

 

 

 

 

1 Неравномерное поле температур в

Рис.

10. Изменение

температуры

объеме,

в который впускается газ, мож­

но в определенной конструкции

стабили­

воздуха при неравновесном

зировать

во

времени

и,

использовав

по­

 

 

сжатии:

 

 

 

 

 

 

вышение

температуры

сжимаемого газа

/ — сжимаемый

воздух

изолирован

для передачи энергии газа в окружаю­

при сжатии

от

входящего

газа [Урав­

нение

(67)1;

2

воздух

сжимаемый

щую

среду,

получить

таким

образом

и входящий

постоянно

смешиваются:

охлаждающее устройство. Это было по­

Гвх “

Гн -

203° К - Р„х - Р к - 5 ат

казано

Гиффордом

и

Лонгфортом

в

[уравнение (68)]

 

 

 

 

статье об изобретенной ими пульсирую­

 

 

 

 

 

 

 

щей

трубе.

 

 

 

 

 

54

батных процессов расширения и сжатия, разумеется, не совпада­ ют (рис. 10). Хотя линия 2 на рис. 10 расположена левее изоэнтропы, суммарная энтропия для взаимодействующих частей воздуха в результате осуществления этого процесса увеличивает­ ся. Входящий воздух переходит в состояние К из состояния О, а сжимаемый — из состояния Я.

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

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

не­

равновесной политропы: например, положив, что q с'ДГ

при

постоянных 2 и т', найдем для процесса расширения:

 

Ь _

т'

 

Гк

1)

(69)

Рн

 

m'—k

CP( T K — T U); пг

 

Я = k(m'~ 1)

 

 

В формулах (69) величина т' имеет смысл эмпирического коэффициента. Для процессов сжатия подобные выражения приобретают форму, аналогичную уравнениям (67) и (68),

сзаменой в них k на т'.

3.Процессы в адиабатной термомеханической системе с переменной массой

Количество рабочего тела в пространстве под поршнем детан­ дера за время цикла меняется. Даже в процессах расширения {23) и обратного сжатия (56) при закрытых клапанах количе­ ство газа в рабочем объеме не остается постоянным, так как существуют неизбежные утечки и натекание газа и специфи­ ческие «мертвые» объемы между цилиндром и поршнем,

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

вмертвом пространстве.

Любая термодинамическая система с переменным количе­ ством рабочего тела может быть проанализирована и описана с привлечением методов классической термодинамики. Однако способы решения задач для систем с переменной массой могут

быть разными. Иногда для решения

задачи удобно

выбрать

подсистемы или дополнительные системы с постоянной

массой

и рассмотреть их взаимодействие [5].

В других случаях более

простым оказывается прямой путь исследования системы с пе­ ременной массой [48].

55

Рис. 11. Внешнеадиабатная система с переменной массой

На рис. 11 изображе­ на внешнеадиабатная си­ стема, контрольная по­ верхность которой совпа­ дает с границами рабо­ чего объема над порш­

нем.

В сечении аа

в си­

стему

втекает

поток dG

с температурой

Та,

При

движении поршня параметры системы: р, Т, G, V меняются. Для составления уравнения, связывающего изменения этих парамет­ ров, выполним следующие операции.

Дифференцируя

уравнение

состояния и деля его на G.,

найдем

dp

dV

dT

dz

 

dG

(70)

G

~ p

V

~T~

~ 7~ '

 

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

idG = dU + pdV.

(71)

В более общем случае, когда в систему входят или из си­ стемы уходят несколько потоков, левая часть уравнения (71) принимает вид: hidG + dQ (если система неадиабатная). Рас­ сматривая рабочее тело при постоянных z, cv и ср, найдем

dU = cyd(GT) = cvTdG + cvGdT\ idG = cpTadG.

Уравнение (71) представим теперь так:

(kTa- T ) - ^ - = dT + ( k - l ) T ^ - .

(73)

После подстановки в это выражение уравнения (70) найдем искомое соотношение, связывающее текущие параметры газа в системе:

( Т а - Т ^ - Т е —

dp

(74)

Р

56

В зависимости от условия задачи, температура Та может быть переменной или постоянной величиной. Полученное выра­ жение является общим. Из него можно получить известные частные решения. Например, если температура газа Та в сече­ нии аа остается все время равной текущей температуре газа Т, что может быть в двух случаях — когда газ в систему не посту­ пает или когда газ уходит из системы, то уравнение (74) пре­ образуется в известное уравнение равновесной адиабаты Пуассона:

-^ - = —

.

(75)

р

k — l

Т

v

Если температура Та не равна

температуре газа

в системе,

то происходит необратимое смешение. Уравнение (74) остается

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

если теплота смешения равна нулю.

Например,

для процесса впуска газа в какой-либо объем dV =

= 0 и Та =

Гвх = const. Тогда из уравнения (74) найдем

 

 

L

dp

■ ^ + ТВХ^ ~ = о.

(76)

 

k

Р

р

Т

 

или, приводя к виду, удобному для интегрирования, имеем диф­ ференциальное уравнение процесса впуска:

 

-Р -

(

\тТ

■jd T .

 

(77)

 

 

 

 

Интегрируя выражение

(77) в пределах от ри до рк,

получим

 

Т =-

kTВХ^н Рк

 

(78)

 

 

Ри

 

 

1К

 

kTBX4- Тя

Рк.

 

 

 

 

 

 

 

 

 

 

Ри

 

 

 

 

 

 

 

 

 

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

(68).

се

В другом типичном процессе с переменной массой — процес­

наполнения при постоянном давлении dp = 0

и Гл = Ги =

=

const. В этом случае

из уравнения (74) найдем дифферен­

циальное уравнение процесса наполнения при р =

const.

 

 

dV

Твх

dT

 

 

 

V

 

Т и к - Т

Т

 

(79)

 

 

dT

 

 

dV

dT

 

 

 

V

 

т

т - т в

 

 

 

Интегрируя выражение (79)

в пределах от Ун До Ук, полу­

чим

 

 

 

 

 

 

Т

= Т

 

 

(80)

 

1 К

 

1 вх

 

 

 

Vk- V h + v

Т .

57

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

4. Волновые явления в поршневом детандере

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

бых возмущений, т. е. скорость звука по Лапласу а — V kRT увеличивается

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

ную волну сжатия — так называемую

ударную волну. При расширении газа

во время выхлопа температура падает

и волны разрежения в цилиндре не

могут образовать ударных волн. В ударной волне энергия рассеивается в ре­ зультате столкновения молекул. Кинетическая энергия молекул, как частиц потока, переходит в энергию теплового (хаотического) движения. Таким об­ разом, процесс сжатия газа при прохождении ударной волны также является неравновесным адиабатным процессом, сопровождающимся увеличением энтропии. Изменение параметров газа определяется неравновесной или удар­ ной адиабатой Гюгонио:

 

1+

к'— 1

рн

 

------ • —

7 Н

__

к

+

1

р к

 

 

 

 

 

(81)

7 к

|

к

 

1

Р к

 

 

к

+

1

р н

Возникнув при впуске, ударная волна многократно отражается от порш­ ня и днища цилиндра и постепенно затухает во время процесса наполнения

12. Возникновение ударных волн при впуске не является неизбежным. Кроме

Рвх

отношения давлений------- , решающую роль играют такие конструктивные

Ре

факторы, как профиль впускных каналов и динамический режим открытия клапана. Отвлекаясь от этих конструктивных факторов, и считая, что проис­ ходит мгновенное открытие клапана, можно найти число М (отношение мест­

ной скорости газа к местной скорости звука) при входе газа в цилиндр [49].

58

В первом приближении, при равенстве коэффициентов Пуассона для газа при

параметрах до клапана и в цилиндре (kBZ ~ ke) , связь М и

Рвх

------- имеет вид

 

Р«

- 2 - ю - ± = ±

Рвх

Л + 1

k + \

Pt

 

(82)

1

 

 

 

Изменения температуры

и давления

при прохождении ударной волны

в зависимости от числа М и начальных параметров определяются следующим

образом:

 

2k

 

 

 

 

Рх—Рв — Ри fe+1 (М2— 1);

 

(83)

 

2(k -

1 ) (Af2— 1)'(1 + Ш 2).

 

 

 

(k + l)2Af2

 

 

Скорость распространения ударной волны Су и скорость спутного, т. е.

следующего за ударной волной потока Сс, находятся так:

 

 

Су —Мвц,

 

 

 

 

А+ 1 (*-£г) ви

(84)

 

Сс = — - —

 

 

 

Влияние

Рвх на относительную

амплитуду

колебаний давления _Др_

 

Ре

 

 

Рн

и температуры

ДГ

пространстве

над поршнем при прохож­

газа в рабочем

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

рнх

Pt

Т а б л и ц а 10

Относительные амплитуды колебаний давления

итемпературы газа при прохождении ударной волны

врабочем пространстве над поршнем

м

дГ

Ар

дг

 

 

Д 7 -

АР

АТ

г н

•° н

Г н

^ВХ

м

Т н

Рн

Тн

 

 

 

У р а в н е н и я

 

Pt

 

У р а в н е н и я

 

( 8 2 )

( 8 3 )

( 8 3 )

( 6 7 )

 

( 8 2 )

( 8 3 )

( 8 3 )

( 6 7 )

 

= 1 , 6 6

 

 

 

 

* = 1 ,4

 

 

2

1.12

0,117

0,317

0,398

2

1,15

0,099

0,375

0,287

5

1,33

0,32

0,96

1,59

5

1,4

0,182

1,12

1,145

10

1,54

0,583

1,712

3,58

10

1 . 6

0,388

1,82

2,57

20

1,825

0,853

2,910

7,55

20

1,825

0,551

2,72

5,43

50

2,38

1,60

5,83

19,47

50

2,14

0,805

4,165

14

59

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

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

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

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

Ъ

60

Рис. 13. Периодическое изменение тем­ ператур в системе газ — стенка

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

 

дТ(г; т)

 

д2Т (п т)

_!_

дТ{г; т)

 

 

 

 

дх

- [ ■

 

дг2

 

+

г

 

дг

 

 

 

(85)

 

_ х dT(R0; х)

[Tr_ T(Ro. т)].

 

 

 

 

 

 

 

 

 

 

 

дг

 

 

 

 

 

 

 

 

 

 

 

 

 

dT(R; т)

--0

(изоляция);

 

 

 

 

 

 

 

дг

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т (г;

0) = Г ср,

 

 

 

 

 

 

 

где Тг—ТСр + ДТр cos сот

 

— температура

газа; г,

R„

и R — соответст­

 

2

 

 

 

 

 

 

 

 

 

 

 

 

венно текущий, внутренний и внешний радиусы цилиндра;

а — коэффициент

теплообмена на внутренней

поверхности

цилиндра;

Я,

а — соответственно

коэффициенты теплопроводности

и

температуропроводности;

т — время;

ш —

круговая

часгота; Гср — средняя

температура

газа;

ДГг — полная

амплитуда

изменения температуры газа.

 

 

 

 

 

 

изменения

температур

Конечно, аппроксимация действительного закона

косинусоидой не для всех циклов одинаково удовлетворительна.

Она

тем

надежнее,

чем меньше с0 и больше Ь0\

при больших с0

и малых

Ь0 действи­

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

61

Решая задачу для случая установившихся температурных волн, получим

 

 

 

Т(Г,

т) - Г ,ср = A cos(PdFo— Af) =

 

 

 

 

 

 

 

ДГг/2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

 

 

 

 

 

Ni ~ N~i M

 

 

( 86)

 

=

(NiN_i)

2 cos

PdFo— arctg

 

 

 

Ni+N -i) Г

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K,(fe УЛр5 )/,(р

/ Г р З) + / , ( й 1^rPd)/C0(p УГрд)

 

 

 

[/<,(* VTPd)lt{VTPd) + /,(fty7pd)K o( УТр З)]—

 

 

 

 

Bi

 

 

 

 

 

 

 

 

 

кгра)/с,(/гра)],

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

причем p = —

; k-

------; Bi = —

 

 

 

 

 

 

 

to

 

 

Rg — критерий Био; Pd = —

RZ— критерий

 

Ro

 

Ro

 

 

a,

 

 

 

 

 

 

 

a

 

 

 

Предводителева;

Fo = ——

— число

Фурье;

/v (z)

и

K v (z )— функции

Бес-

селя

 

 

Rl

 

 

 

первого

и второго

рода

соответствующего

от чисто мнимого аргумента

порядка v; AL; определяется аналогично

N и только

V <

заменяется везде

на

V I-

температуру

на

зеркале цилиндра,

приняв

р = 1,

тогда

ам­

 

Определим

плитудная и фазовая характеристики выразятся соответственно так:

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

Л* =

1-

Bi*~

+

<№>■

]

1.

 

 

 

«(87)

 

 

 

Af* =

arctg Г—

 

2

 

 

 

 

 

 

 

 

 

 

 

L 2Bi*— e

J

 

 

 

 

 

 

 

 

Модифицированное число Био для

стационарно-периодического состояния

 

 

 

 

 

Bi* =

а

 

 

 

 

 

 

 

 

(88)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V ХсрцШ ’

 

 

 

 

 

 

 

где

с, рц — соответственно удельная теплоемкость

и

плотность

материала

цилиндра.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Коэффициенты 6,

е и <р выражаются уравнениями

 

 

 

 

 

 

 

g2 + е2

 

 

 

ed + gb

<р = 2

eb— gd

 

 

 

 

 

d2 + b2

е _

 

d2 + 62

 

d2+ b2

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b = her'(ftz)bei(z) +

hei'(Az)ber(z) — ber'(fez) hei (z) — bei'(fez) her(z);

 

 

d =» her'(£z)ber(z)— hei'(fe)bei(z) — ber'(fez)her(z) + bei'(fez)hei(z);

 

 

e = her'(Az)ber'(z) — hei' (Az)bei'(z)«—ber'(ftz)her'(z) + bei'(fez)hei'(z);

 

 

g = her'(fez)bei'(z) + her' (fez)ber'(z)— ber'(ftz)hei'(z) —

 

 

 

 

 

 

 

— bei'(fcz)her'(z);

 

 

 

 

 

 

 

здесь [her(z), her'(z),

ber(z)

и т. д.— функции Кельвина]

 

 

 

 

 

 

 

 

 

 

z =

/ P d .

 

 

 

 

 

 

 

 

62