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

книги из ГПНТБ / Васильев В.К. Термодинамические основы исследовательского проектирования судовых энергетических установок

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

(для обозначения критической скорости и других параметров потока при этой скорости примем подстрочный индекс «к»);

максимальная скорость потока смакс, определяемая форму­ лами (296) или (297);

скорость звука в среде, определяемой параметрами затормо­

женного потока с подстрочным индексом «О».

Соответственно получаются четыре скоростных характеристики

потока:

формулами

 

 

1)

число М, выражаемое

 

 

 

 

с

_

с

(312)

 

 

y j T T T W

 

2)

скоростной коэффициент %

 

 

 

 

 

 

с

 

(313)

 

 

 

 

 

3)

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

 

 

 

 

 

С

 

 

 

(314)

 

смакс

 

 

 

 

 

 

 

 

4)

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

 

 

 

 

 

с

______с

 

 

(315)

 

 

 

 

 

Из последних формул видно, что лишь одна характеристика — число М — определяется текущими параметрами потока. Остальные характеристики относят местную скорость к характерным скоростям, которые выражаются параметрами торможения. Если поток тепло- и энергоизолирован, то указанные характерные скорости сохраняют постоянное значение вдоль всего потока.

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

Используя формулу (213), связывающую местное значение энталь­ пии i с местным значением произведения pv, получаем

k — 1

 

с2

М2

2

k

 

I

226

Прибавляя к обеим частям этого равенства единицу, будем иметь

1+

. |

^

 

 

k 1

2

('о

'

2

i

~ i ~~ Т

Как видно, число М2 связано

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

торможения Т 0.

Эта связь дает возможность выразить и остальные параметры тормо­ жения в любой точке процесса расширения через число М:

k — 1 М2;

2

k

(316)

 

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

1

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

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

Здесь не будем выводить этой зависимости, а возьмем ее из [16];

(№ - l ) - f =

4 - - ± i L , -

± d L ,

+ Д р { - % - ) / L r . (317)

Для идеального газа

^ •

Подставив эту величину

в уравнение (317),

получим

 

 

(М2— 1)-^- = dF_

 

(318)

 

F

 

 

15*

227

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

или в дифференциальной форме

dc____ dT

М2

с

Т

После несложных преобразований из уравнения (318) получим

(319)

Это уравнение определяет влияние внешних воздействий на число М. Если воздействия заданы, то можно найти число М, а затем по уравнениям, определяющим параметры потока, получить и сами эти параметры. Однако проинтегрировать уравнение (319) можно только в том случае, если известны связи внешних воздействий с параметрами потока. Если эти связи неизвестны или их вообще не существует (воздействия произвольны), то исчезает определенность

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

иуравнение (319) становится пригодным только для получения ка­ чественных зависимостей.

§ 32. ИНТЕГРИРОВАНИЕ УРАВНЕНИЙ ОДНОМЕРНОГО ПОТОКА

Изолированное геометрическое воздействие* Изолированное гео­ метрическое воздействие может иметь место в случае течения рабо­ чего агента через прямоосный канал переменного сечения при тепло­ вой и энергетической изоляции потока от внешней среды. При этом предполагается отсутствие внешнего и внутреннего трения. Таким условиям удовлетворяют идеализированные процессы течения через сопловые каналы турбинных проточных частей и по различным участкам прямоосных труб и аппаратов. Используя уравнение (319), следует принять в нем LT = 0 и Lr — 0. Тогда уравнение получит вид

М2 _

1

dM2

k — \ М2

(320)

М2

— s —

М а

 

228

Переменные здесь разделены, и уравнение легко интегрируется.

В интегральной форме оно будет выглядеть так:

 

 

к—1

Ь_ 1

 

 

1+ i _ L Ms

 

(ЛЫг

k+1

 

(321)

)

к _ I

\ МF

 

причем здесь интегрирование выполнено в пределах от начальных параметров процесса расширения с подстрочным индексом «1» до конечных параметров без индекса (текущих).

Уравнение (321) показывает, что в процессе расширения число М однозначно связано с проточной частью канала. Г1рименив закон

обращения воздействий

[16]

к уравнению (320),

получим при М = 1

(т. е.

при достижении критиче­

 

 

 

 

 

 

ской скорости) dF=0.

Это зна­

 

 

 

 

 

 

чит, что площадь поперечного

 

 

 

 

 

 

сечения канала должна принять

 

 

 

 

 

 

минимальное значение,

а удель­

 

 

 

 

 

 

ный массовый расход через эту

 

 

 

 

 

 

минимальную площадь

должен

 

 

 

 

 

 

стать

максимальным.

 

 

 

 

 

 

 

 

 

 

В докритической зоне про­

 

 

 

 

 

 

цесса течения (М < 1 )

знак dF

Рис. 35. Форма прямоосного канала для

должен быть отрицательным, а

 

расширяющегося потока.

в сверхкритической (М >

1)—

 

 

 

ходу процесса рас­

положительным. Следовательно, канал по

ширения в докритической

зоне

должен быть сходящимся

(с умень­

шающимися площадями поперечных сечений),

в

сверхкритической

же — расходящимся (рис.

35).

 

 

 

 

нетрудно

перейти

От числа М, определяемого формулой (321),

к определению

значений

текущих

параметров

 

процесса

расшире­

ния в рассматриваемом случае.

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

 

7’0 =

T1( l + - Ь ± М ? ) = T ( l +

^

z

l M2) ,

 

откуда

 

 

 

ь _1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1+ —— — М 2

Т 1

 

 

 

 

 

 

 

^

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ь _1

 

Заменив в

уравнении

(321)

отношение

1 + ~ ~ 2 ~ М2

отноше-

-------г— ,-----

нием температур T J T ,

получим

k-i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т

 

/

m1F1 \ k

 

 

 

(322)

 

 

 

т х

-

\

MF

)

 

 

 

 

 

 

 

 

 

 

229

Используя зависимости изоэнтропийного процесса расширения (217), найдем и другие относительные параметры потока:

 

 

 

 

k

 

 

 

 

 

2k

 

 

 

р

(

Т \ k-l

~

(

MjFj

\ k+i .

 

 

Pi

 

П )

V

м F

)

2

 

 

 

 

 

 

I

 

 

М1F1

\

'

 

 

р

 

Т \ k-l

/

k+i

(3-23)

 

Pi =

(

Тг )

=

\

МF

)

 

 

 

 

 

 

1

 

 

 

 

 

k-i

 

 

 

/'

т

\ 2

М

/

MxFi

\

k+i

М

 

 

. Тг )

Мх -

V

MF

)

 

 

Mj

 

Положив в уравнениях (322)

и (323)

М = 1

и определив значе­

ние FMHHиз уравнения (321) при М =

 

1, получим по уравнениям (322)

и (323)

значения критических

параметров

процесса расширения:

 

 

 

 

 

 

к ~

х ж\\

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

Рк

 

 

 

k

1

w 2 \

k— I

!

 

 

Pi

 

 

 

T T T Ml)

 

(324)

 

 

 

 

 

 

 

 

 

I

 

 

 

 

 

 

 

 

 

 

 

 

 

Рк

 

 

 

k

1

, X2 \ ~ ^

 

 

 

Pi

 

 

 

 

 

 

 

 

I

 

 

 

 

 

 

 

 

 

 

 

 

 

_£к

 

 

 

 

 

k + 1 «

T

 

 

Cl

 

 

 

 

 

0

4

 

Если можно пренебречь начальной скоростью потока,

считая

Мх = 0,

то уравнения

(324)

примут вид:

 

 

 

 

(325)

2k

PtPi = VkpA = aK.

сК V - k + 1

Нетрудно заметить, что в обоих случаях критическая скорость представляет собой скорость звука, рассчитанную по критическим параметрам.

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

230

невозможно,

расчет

критических параметров следует

выполнять

по формулам

(324).

 

 

Комбинированное

воздействие — геометрическое и

воздействие

трения. Такое комбинированное воздействие соответствует работе прямоосного соплового канала с учетом сопротивления процессу

течения. Уравнение

(319)

при этом перепишется так:

 

(Мг- 1

) ^

= 2 А ^ _ 2 А А dLr.

(326)

Здесь, как и в предыдущем случае, легко отделяются переменные. Выразив элементарную работу трения через коэффициент потери энергии элементарного действительного (адиабатного с трением) процесса расширения:

 

dLr = ~ l d i = -

l sdis,

(327)

в результате интегрирования

уравнения (326) получим

 

 

м-1

 

2k

 

(4

г

 

<з28>

В последних

выражениях

введено

равенство:

 

 

1 +

М2 =

- у - .

(329)

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

комбинированное

воздействие, удобно преобразовать

уравнение

(289) и при dh = 0 записать такое уравнение энергетического баланса:

dQa — d(^i Н— -f-dLT.

В это уравнение не вошла работа сил трения, так как кинетиче­ ская энергия потока dLr, затраченная на преодоление сил трения, перешла в эквивалентное количество теплоты dQr, которое осталось в потоке и добавилось к внешнему теплообмену dQa.

Таким образом, если надо явно ввести в уравнение энергетиче­ ского баланса воздействие трения, то к правой части уравнения следует прибавить величину dLn а к левой части — величину dQr. Очевидно, подобная операция не нарушит энергетического баланса, поскольку dLr = dQr. Однако она окажет влияние на параметры потока. Величины г и с в уравнении будут не такими, как при отсут­ ствии трения. В данном случае поток предполагается теплоизолиро­ ванным от внешней среды (dQa = 0) и принимается отсутствие меха­ нического воздействия (dLr = 0). С учетом сказанного будем иметь

dQr = di0 + dLr.

(330)

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

231

тенциальной энергии при постоянстве их суммы. Это было ясно и без уравнения (330), которое, однако, позволяет обратить внима­ ние на самый характер теплообмена. Это не внешний, а внутренний теплообмен, поскольку внешняя теплоизоляция потока остается. Кроме того, уравнение (330) показывает, что при внешней энерго­ изоляции потока его кинетическая энергия снижается из-за расхо­ дования некоторого ее количества (dLr) на преодоление сопротивле­ ний движению внутри потока. Таким образом, можно прийти к за­ ключению, что при постоянстве i 0 из-за трения снижается скорость с

 

 

потока

и

увеличивается

его текущая эн­

 

 

тальпия.

 

 

 

 

 

 

 

Как видно, воздействие трения свя­

 

 

зано с

необратимостью

действительного

 

-dp

процесса

расширения в рассматриваемых

1

/

условиях.

Так

как все формулы

класси­

 

УС

ческой термодинамики выводятся при ус­

/

ловии обратимости (равновесности) про­

 

 

цессов,

то выявленная

здесь

необрати­

 

 

мость лишает нас возможности

использо­

p-ip

 

вать закономерности классической термо­

 

 

динамики, что, несомненно,

затруднит

 

 

исследование

рассматриваемого

случая.

 

 

В § 29 было разъяснено,

что причина нео­

 

 

братимости кроется в необратимости про­

 

 

цесса механического трения (затрата ме­

Рис. 36. Элементарный про-

ханической энергии на преодоление сопро­

тивления движению зависит не от

самого

цесс

расширения

потока

и закономерностей его

движения,

 

 

аот конструктивного оформления каналов

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

можно отнести к любым необратимым (неравновесным) процессам

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

иот их закономерностей.

При анализе процесса расширения с учетом внешних воздей­ ствий (геометрическом и трения) обратим внимание на то, что урав­ нение (330) написано для элементарных процессов изменения состоя­ ния рабочего агента. На рис. 36 в диаграмме i—s изображен эле­ ментарный процесс расширения АВС, протекающий от давления р до давления рdp. Предположим, что процесс — теплоизолирован­ ный, но с трением. Обозначим через dis изменение энтальпии, соот­ ветствующее изоэнтропийному процессу при ds = 0. Тогда в соот­ ветствии с уравнением (206) получим dis = vdp. Но при отсутстствии трения величина vdp измеряет отданную во внешнюю среду работу, рассчитанную на единицу массы рабочего агента, включая

232

работу перемещения d (pv) и увеличение кинетической энергии потока. Эту работу мы обозначали dLT. Если процесс расширения происходит с трением, то действительная отданная работа будет меньше располагаемой работы:

dLT= — r\nvdp = r\ndis.

(331)

Из этих соотношений можно установить, что величина т]п является отношением действительно отданной на сторону работы dLT к изме­ нению энтальпии dis, соответствующему при изоэнтропийном про­ цессе изменению давления на ве­

личину dp. Такое отношение на­ зывается политропным коэффи­ циентом полезного действия про­ цесса расширения:

 

ч . —

тгг-

<332)

 

 

Знак минус в формулах (331)

 

 

и (332) обусловлен тем, что в про­

 

 

цессе

изоэнтропийного

расшире­

 

 

ния изменение энтальпии dis всег­

 

 

да отрицательно, а удельная рабо­

 

 

та и к. п. д. всегда положительны.

 

 

В

действительности

адиабат­

 

 

ный процесс с трением может быть

 

 

осуществлен только в том случае,

 

 

если рабочий агент имеет физи­

 

 

ческие

свойства

реальных газов

 

 

и паров. Поэтому к реальным про­

Рис. 37. Одностадийный

процесс рао

цессам нельзя применять формулы,

ширения в диаграмме i—s.

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

 

в этом слу­

идеального газа. В частности уравнение состояния (204);

чае следует учесть коэффициент сжимаемости а — а (Т,

р) и взять

уравнение состояния (222).

Предположим, что имеется процесс расширения реального рабочего агента от давления р г до конечного давления р (рис. 37). Как уже отмечалось, следует условно принять, что теплота трения Qr сооб­ щается рабочему агенту непрерывно и обратимо в процессе его рас­ ширения. Обозначим указанный теплообмен Qr0бРАбстрагируясь от внутреннего источника этой теплоты и полагая, что она сообщается

как бы извне, получим процесс с приращением энтропии ds = .

Линия такого процесса АС будет отклоняться от изоэнтропы АВ. Конечную точку С процесса расширения на изобаре р можно полу­ чить двумя способами. Действительный процесс расширения пойдет от точки А к точке С по линии АС с постоянным увеличением энтро­ пии. Но можно достичь точки С последовательно двумя процессами, из которых первым будет изоэнтропийный процесс АВ, а вторым —

£33

изобарный процесс ВС, в течение которого рабочему агенту сооб­ щается теплота Qr0бр и энтропия его увеличивается на

J d s = j Щ&В-.

вв

Замена действительного процесса АС совокупностью последо­ вательно идущих процессов АВ и ВС удобна и часто используется при расчетах. Однако при этом не учитывается часть возвращенной

теплоты потери Qr,

измеряемая площадью треугольника АВС на

 

 

рис. 38,

где

процесс

дан

в

 

диаграмме

 

 

Тs.

Эта теплота

сообщается

рабочему

 

 

агенту

еще в процессе

расширения и ис­

 

 

пользуется для совершения добавочной по­

 

 

лезной

работы.

 

 

 

 

 

 

 

 

 

При последовательно идущих процес­

 

 

сах АВ и ВС располагаемая

работа

про­

 

 

цесса

расширения

сообщается

рабочему

 

 

агенту

в

изоэнтропийном

процессе

АВ

 

 

в виде разности энтальпий Ais. Изобар­

 

 

ный нагрев

теплотой трения

происходит

 

 

в изобарном

процессе

ВС,

за

счет

чего

 

 

увеличивается энтальпия рабочего агента

 

 

и его энтропия.

Разность

энтальпий

то­

 

 

чек А и С начала и конца

процесса рас­

 

 

ширения будет, очевидно,

на

величину

О

S

Qr0бр

меньше

разности

энтальпий

A/s

в начале и конце изоэнтропийного про­

Рис. 38. Одностадийный про­

цесса АВ

(см. рис.

37).

Если

обозначить

цесс расширения в диаграм­

ме Т —s.

 

разность энтальпий

начала и

конца

дей­

(адиабата с трением)

 

ствительного

процесса

расширения

АС

через Ai (рис. 37),

то отношение

 

 

 

 

 

Ж

-

ч.

 

 

 

 

 

 

 

<333>

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

В конечных процессах расширения (при конечной разности давле­ ний процесса расширения) эти два к. п. д. имеют различное значе­ ние. В адиабатном процессе с трением АС (см. рис. 38) нагрев рабо­

чего агента теплотой трения начинается с точки А.

Процесс идет

по наклонной линии АС. Как известно,

площадь,

ограниченная

линией процесса и изоэнтропами его начала и конца,

в диаграмме

T —s измеряет теплообмен в ходе процесса.

Если расширение идет

от давления р г до р (рис. 38), то в пределах этого процесса рабочий

234

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

-

Д£

_

 

Л/

(334)

1111 =

Ats +

Д ’

115

= ~ К 1 Г '

 

Рассмотрим значение политропного к. п. д. более подробно. Обращаясь к рис. 37, напишем расчетную формулу для величины Ai:

р

р 1

 

Ai = — j

r\nvdp = т)п

J

vdp.

(335)

P i

р

 

 

 

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

нее значение т]п, получая его по формуле (334).

Нетрудно показать, как выведена формула (334). Вычислим интеграл в формуле (335), для чего рассмотрим круговой процесс, изображенный на рис. 38. В него входят изоэнтропийное расшире­ ние 12, изобарный теплоотвод 23, изоэнтропийное сжатие 3—4 ■и изобарный нагрев 41. Полученная в этом цикле работа в диа­ грамме Т—s изобразится заштрихованной площадкой 1—2341. С другой стороны, та же работа равна разности между работой рас­ ширения vdp и работой сжатия v'dp, т. е. площадь 1234— 1 —

= vdp v'dp,. или

 

vdp = площадь 123—4— 1 + v'dp.

(336)

Отсюда можно перейти к интегральным зависимостям. Сумма площадей 1234— 1 составляет площадь треугольника АВС,

которая была нами обозначена через А. Интеграл J v’dp есть не

что иное, как Ats; это видно из формулы (206) при ds = 0. Тогда интегральная форма зависимости (336) будет

Pi

 

J vdp = A{'s +А

(337)

р

и, наконец, из уравнения (335) получим

А/ = г)п (Ais + А).

(338)

235

Соседние файлы в папке книги из ГПНТБ