
книги из ГПНТБ / Криогенные поршневые детандеры
..pdfРис. 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