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

книги / Эксергетические расчеты технических систем

..pdf
Скачиваний:
11
Добавлен:
20.11.2023
Размер:
19.71 Mб
Скачать

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

Целевое назначение рассматриваемой установки — отвод из охлаждаемого объекта некоторого количества теплоты Сохл. Исходя из понятия о приведенной эксергетической холодопроизводительности [33] и определения эксергии холо­ да (теплоты при температурах ниже окружающей среды), целевым назначением установки можно также считать сообщение охлаждаемому объекту приведенной эксергетической производительности едохл, определяемой по зависимости

^?охл = Ф°хл ^ Тохл

l) = ФохлТе.

(9.29)

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

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

для зоны III

е3 =

Е 3(eQoxn,

А Т в);

z3i =

2 31 (е<гохл* ©в,

а Т в),

(9.30а)

 

 

 

 

г33 =

^32 (^Qcxn’ ®в>

А Т в);

^32 =

Е 32(^Сохл*

А^в)’

 

для зоны I I

 

 

 

 

 

 

 

 

ег =

Е 2(ез,

A T s);

^21 — ^21 (^3» ®и» A^s)»

(9.306)

е 22 =

Е 22 (е3, ®и> A T J,

^22 =

^22 (*3, ®И» Д

 

 

 

для зоны /

 

 

 

 

 

 

 

 

еи =

^11 (е2> ®к» &Twy,

Zn =

Z 11

®к» Д7*w)9

(9.30в)

е1з =

Егз (е2> ®к»

г12 =

^ 1 2 (е2>

Д7*ш);

^12 =

^12 (e2t ®К» ^ 7 w)>

zi3 =

Z u (е2, 0 К, ДТ’ц,).

 

В выражениях (9.30) строчные буквы ei, zu v±2 обозначают численное зна­ чение величины, а прописные Еь Zit V12 — функциональную зависимость.

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

П = [Цэл (£ц + £13 + ^22 + езт) + V l 2 + Z11 Z12 +

 

+ 213 + z21 + Z22+ ^3i “Ь z32] тг-

(9.31)

Расход электроэнергии на привод двигателей компрессора, насосов охлажда­ ющей воды и хладоносителя, вентилятора, расход охлаждающей воды, а также стоимость компрессора, конденсатора, насосов охлаждающей воды и хладоноси-

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

дП

= 0;

дП

= 0;

дП

= 0;

дП

= о;

а п

= 0.

двк

дАТи

дви

dATs

dATB

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

Использование термоэкономики позволяет представить сложную энергети­ ческую систему как последовательное соединение элементарных процессов (зон). Применяя методику построения характеристики объекта по характерис­ тикам элементарных процессов при их последовательном соединении [61], вы­ ражения (9.30) можно представить в виде

 

 

 

 

е3 = Е3(Qoxnf Е0хл, E0tC, 0 В, АТа)\

 

 

 

 

 

 

&32 =

E32(QQM,

Еохл,

Го.с, ©в»АГв);

 

 

 

 

 

 

^31 =

^31 (Сохл1Т 0хл,

Т 0'С,

0 В, А Т в)\

 

 

 

 

 

 

Z3 2 =

Z 3 2 №охл,

Тохл,

7\).с»

0 В,

А Т в );

 

 

 

^2 “

Е 3 [^*з (Сохл,

Т охл, Т 0 'С1

0 В, А Ев),

0 И, АТ^],

(9.32)

 

е 22

 

=

Е22 [Е3(Сохл,

Е охлу ТОЛ9 0 в*А Г в),

0и, A T J;

 

 

Z21

Z2i [Е3(Сохл, Т0хл, Е0шС, 0 В, ДГВ),

0 И, ATS]',

 

 

е22 “

^22 [Е3 ((дохл, Тохлу Т0ьСу 0 В, АТъ),

0 И, AEJ,

 

*11 = Еп 2 [Е3(Сохл , Еохл, Ео.С, ©

в ,

АГв),

0 И( ATJ,

0• к> ATwy,

е13 = Егз

2 [Е3(Сохл1

 

Ео.С, ©

в ,

АГв),

0 И,

ATS],

©к. A7U;

 

 

 

, ТбхЛ,

 

 

 

 

 

 

»!2 :=

IV [Е2 [Е3охл, Еохл , Ео.С, ©

в ,

АГв),

0„,. ATJ,

©К. А Тш);

= z n 2 [Е3(Сохл, Т охл, Ео.С, ©

в .

АТв), 0 Н, AT’sl. ©к. АТш);

*12 =

Z12[Е2 [Е3охл, ЕОХЛ, Т о.С, ©в, АГВ), 0„,

ATt], ©к, ATwy,

г13 =—^132 [£'з(Фохл» ЕОХЛ, •7Vc © в ,

АГв),

0„,

A7\],

© К .

ьтш).

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

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

е3= Е.6 (eQoxn, 0„, ДТв)\ е2= Е2 (е3, 0 И, ATs).

Для нахождения условного экстремума используется метод неопределен­ ных множителей Лагранжа. Для этого составляется выражение функции Лаг­ ранжа, которая представляет собой сумму исходной функции приведенных затрат (9.31) и произведения уравнений связи, записанных в неявной форме, на некоторые, неизвестные вначале, множители Лагранжа А*.

Такое выражение для рассматриваемой задачи представляется в виде

L = {цЭл [<£*11 (^2»©к»

w) + El3(e2l ©к, ATW)

4~ Б22(е31©и» ^£s) 4~

+

£32 (е<?охл, ©в, А£в)] +

u.wV12{e2i ©к. ATw) +

2 u (£?a, ©к,

+

 

4- 2 12(^2, 0 К, ДTw) + Z13(e2, 0 К, ДTw) + Z21(e3, 0 И, ДТ5) +

 

+

Z22(е3, 0 И, ATS) +

Z31 (eQoxJ1, ©в, ДГВ) + Z32(eQoxjj, 0 В, Д Г В) +

 

+ Х2[£2 (в3, 0 И, ДTs) -

е2] + А3 [£3 (е0охл,

0 В, ДГВ) - е3]} тг.

(9.33)

Выражение для приведенных затрат (9.31) и выражение (9.33) с учетом зависи­ мостей (9.30) полностью тождественны.

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

связи) и приравнять их к нулю.

Частные производные по et позволяют вычислить значения неопределенных

множителей Лагранжа

 

 

 

 

 

~дё— = Тр де2~ ^ эл (2Г11

13) + ЦиУ12 4" 2ц + Z12+

Z13

А2в2] =

0.

Так как тр ф 0, то нулю должно быть равно значение производной, откуда

А2 =

[цэл (£ц + £ 13) 4"

12 4" 2 U + Z12 +

2 13];

 

^ —= тг

— [цЭл£г2 4“ 22I + Z22+

А2£ 2 А3в3] =

0 (тР Ф- 0),

(9.34а)

откуда

 

 

 

 

 

 

 

Я8 =

(Цэл^и + 221 + Z*2+

 

 

(9.346)

Частные производные по оптимизирующим переменным после аналогичных преобразований позволяют получить следующие выражения:

-gg^- [Цэл (£*11 + £ is) +

12 + Z 11 +

Z 12 +

Zisi =

df T~ [цЭл (Еи + Е13) +

nwv 12 + Zn +

z 12+

Z13] = 0;

' 0 0 " (Д эл ^2 2 ” t“ Z 21 “ Ь ^ 2 2 "Ь ^ 2^ 2) =

(Дэл^22 + ^21 + ^22 + ^2^ 2) = O',

(дэл£ 32 + Zsi + z 32+ К Ез) = о;

(Дэл£32 + 2,1 + Z32 + X3E3) = 0.

(9.35)

Система уравнений (9.35) устанавливает связь между диссипацией энергии и затратами в каждой зоне модели при определенных значениях экономических показателей цэл, Ц®, А,2, А,3. Система уравнений (9.34а) и (9.346) позволяет оп­ ределить действительное значение экономических показателей Х2 и ^з-

Если выражения (9.34) записать в форме конечных приращений

л*

А [иЭл (£ ц + EU ) + uwVi2 + 2 ц + 2 i2 + 2 13]

 

а*

- ^ (Цэл^22 + 2Zi + 222 -f А*Е2)

 

3 =

----------------- -------------------- ’

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

С помощью выражений (9.34а) и (9.346) в общем случае определяют ско­

рость изменения затрат при изменении количества эксергии. Однако поскольку зависимости £ u , £ 13, V12, Zn, Z12, Z13— линейные функции e2i а выражения

E 22t z21, Z22l E2 — линейные функции e3, то графические изображения произ­ водной и функции будут совпадать, что обеспечивает равенство А2 = К- и А3 =*

= А3. Следовательно, в рассматриваемом случае неопределенные множители Лагранжа Я* выражают стоимость единицы потока эксергии, выходящего из каждой зоны термоэкономической модели, найденную с учетом всех отчислений от стоимости оборудования и потерь эксергии в рассматриваемой зоне.

Аналогично определяется стоимость единицы приведенной эксергетической производительности. Очевидно, что

Ьо=

(Цэл^зг Z31-J- Z32-f- X3E3).

 

^охл

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

Соотношения, положенные в основу математического описания

термоэкономической модели

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

исходные данные:

холодопроизводительность установки — Q0x*, кДж/с; температура охлаждаемого объекта — Тохл, К; температура окружающей среды — То.с, К;

температура поступающей охлаждающей воды — TWx, К;

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

потери давления в аппаратах и трубопроводах — Apif Па; вид и свойства рабочего тела; вид и свойства используемого хладоносителя;

тариф на электроэнергию — ЦэЛ, руб./(кВт ч); тариф на воду — nw>руб./м3.

Кроме того, в соответствии со сделанными допущениями нужно задать вели­ чины АТПи Д71ж.

Эксергетическая производительность определяется зависимостью (9.29). Температура воздуха, поступающего в воздухоохладитель, может быть принята равной температуре охлаждаемого объекта ТВх = Т0Хл• Исходя из этого находится средняя температура воздуха, циркулирующего через воздухо­

охладитель:

7*

-{- Т

а'р

 

Г Г ___

В 1

1 в 2 ___ 'П

В

 

Потери эксергии при необратимом переходе теплоты Q0XJl с уровня,

соответству­

ющего температуре охлаждаемого объекта

Т0хЛ, на уровень, соответствующий

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

Thi рассчи­

тываются по зависимости

 

 

 

 

Следовательно, эксергия, которой обладает поток воздуха*, используемый для отвода теплоты от охлаждаемого объекта (без учета теплоты, вносимой вен­ тилятором воздухоохладителя), определяется выражением

е в = е<?охл + d 0 = Qoxn( —

\ т

где eqохл — приведенная «нетто» холодопроизводительность установки. Исполь­ зуя эту величину, можно получить выражение для вычисления эксергии, под-

* Здесь и далее рассматривается только термомеханическая эксергия потока.

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

 

 

^ОХЛ________ Яд_____

(9.37>

 

£ .4 2 -----

Те

(СвРвА^вЛв

ЯВ)Т]В

 

 

 

где Нв — напор,

развиваемый

вентилятором

воздухоохладителя, кПа; св

удельная массовая

теплоемкость воздуха, кДж/(кг К); рв — плотность воз­

духа, кг/м3; т]в — КПД вентилятора; т|в.э — КПД электродвигателя вентиля­ тора. Выражение (9.37) справедливо, если вся мощность, подводимая к электро­ двигателю вентилятора, в конечном счете передается воздуху. Это требует некоторого увеличения расхода воздуха для сохранения неизменным принято­ го охлаждения в воздухоохладителе АТв.

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

Выражение (9.37) является развернутой формой представления зависи­ мости для эксергии e32t записанной в системе уравнений (9.30) и в выражении лагранжиана (9.33) в общем виде.

Теплота Q B.O, которая отводится от воздуха в воздухоохладителе, с учетом мощности, потребляемой электродвигателем вентилятора (если этот электро­ двигатель расположен в пределах охлаждаемого объекта),

B O

Qoxn

"f" £;32-

(9.38)

Q .

 

Если электродвигатель вентилятора вынесен за пределы охлаждаемого объек­ та, выражение (9.38) принимает вид

 

Q B.O (Зохл "f" ^згЛв.э»

(9.38а)

Количество эксергии,

которой обладает поток воздуха, с учетом теплоты*

вносимой вентилятором воздухоохладителя, находится по выражению

 

т

П (

 

^°'С

 

£в.о —■Чв.0 I

Т

дт*

 

 

I

--------

 

 

\

охл

2

 

Потери эксергии при необратимом переходе теплоты QB<0 с уровня, соот­ ветствующего средней температуре воздуха Тв, на уровень, соответствующий средней температуре промежуточного хладоносителя TS = TB— 0 В, запи­ сываются как

Эксергия е3, которая подводится потоком промежуточного хладоносителя к зоне III термоэкономической модели (см. рис. 9.8), представляет собой сум­

му эксергии, которой обладает поток воздуха, и потерь эксергии при передаче теплоты

е3 = ев.» ■+■^в.о = Q B .O (

ту;

1V

(9.39)

\Г 0ХЛ----- 2 ^ - 9 *

I

 

Выражение (9.39) с учетом уравнений (9.29), (9.37), (9.38) и (9.38а) можно •представить в виде зависимости от приведенной|холодопроизводительности «нет­ то». При размещении электродвигателя вентилятора воздухоохладителя в ох­ лаждаемом помещении оно будет следующим:

С9рхл

|~ СвРвДТа'ПвПв.з + (1 — Уз) Нв

(9.40а)

Дтв - в .

L

свр>вАТВТ1В

Нв

 

 

 

л при размещении электродвигателя вентилятора воздухоохладителя

вне ох­

лаждаемого помещения таким:

 

 

 

 

ез

 

___^ВрвДТУПв

(9.406)

 

(^ВрвДТвЧв

Яв)

Выражения (9.40а) и (9.406)

представляют собой развернутую форму за­

висимости для эксергии е3, приведенной в системе уравнений (9.30) и в выраже­ нии лагранжиана (9.33) в общем виде.

Используя выражение (9.39), можно вычислить эксергию, подводимую к электродвигателю насоса промежуточного хладоносителя от внешнего источни­

ка, т. е. потребляемую им мощность

 

 

 

 

^22

 

Я,

 

 

5р5ДГ5т), — H s) т]и

(9.41)

 

 

 

 

где

Hs — напор, развиваемый насосом

хладоносителя, кПа;

с, — удельная

массовая теплоемкость хладоносителя,

кДж/(кг

К); ps— плотность хладоно­

сителя, кг/м3; т), — КПД насоса хладоносителя;

т)5Э— КПД электродвигате­

ля насоса.

 

 

 

су

Выражение (9.41) также справедливо, если мощность, подводимая к насо­

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

передаваемую

промежуточному хладоносителю. Это требует для сохранения постоянства из­ менения температуры хладоносителя А Г некоторого увеличения его расхода.

Теплота, эквивалентная потерям мощности в насосе, подводится к проме­ жуточному хладоносителю, вызывая его подогрев и приводя к некоторому пе­ рераспределению температур. В связи с незначительным влиянием на последую­ щий расчет это перераспределение температур во внимание не принимается. Зависимость (9.41) является развернутой формой записи выражения для эксергии е22, приведенного в системе уравнений (9.30) и в выражении лагран­ жиана (9.33) в общем виде.

Теплота, которая отводится от промежуточного хладоносителя в испари­ теле, определяется «брутто» холодопроизводительностью установки

Qop = ^ в*° "Ь e22r\ss>.

(9.42)

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

Приведенная «брутто» холодопроизводительность будет

J>р _

п°

 

 

 

Со =

,бр

Атл

 

 

У о

 

 

Потери эксергии при необратимом переходе теплоты Qop с уровня, соот­

ветствующего средней температуре

хладоносителя

Ts, на

уровень, соответ­

ствующий температуре кипения рабочего тела Т0 =

Ts — 0

И, находятся по за­

висимости

 

 

 

 

du = Qop

ев —еи

Эксергия е2» которая (с учетом потерь) подводится рабочим телом к зоне II гермоэкономической модели (см. рис. 9.8), представляет собой сумму

е2

бр

 

(9.43)

= *0бр + dH= Qo

АГН

 

 

 

Для дальнейшего использования эксергия е2 должна быть выражена как функция потока эксергии е3, выходящего из зоны II:

е2

 

— 1

с5р5АТ$т)5___

АТ п

(cspsA7’st|s H s )

 

'

■—

 

(9.43а)

Зависимость (9.43) представляет собой развернутую форму выражения для эк­ сергии е2>входящего в систему уравнений (9.30) и в выражение лагранжиана (9.33).

Определенные по зависимостям (9.42) и (9.43а) «брутто» холодопроизводи­

тельность установки QoP и эксергия е2, которая подводится к зоне //, позво­ ляют перейти к величинам, характеризующим процессы, осуществляемые ра­

бочим телом.

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

<7о = h\l) h3

(здесь обозначения точек соответствуют рис. 9.9). Эта методика применима и для математического описания термоэкономической модели. Однако для полу-

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

 

Яо ==Г0

СХ(Т’з

Г0),

Рис. 9.9. Цикл одноступенчатой холодиль­

где сх — средняя

изобарная теплоем­

ной машины на Ts-диаграмме

кость жидкого рабочего тела в интерва-

 

ле температур Т3 — Т0, кДж/(кг К); г0 — теплота парообразования, кДж/кг. Теплота парообразования и теплоемкость жидкости с достаточно высокой точностью могут быть найдены по эмпирическим зависимостям [15], справедли­ вым только для выбранного рабочего тела. Такие зависимости позволяют по­ лучить простые выражения, дающие допустимую для инженерных расчетов

точность.

Массовый расход рабочего тела Gai обеспечивающий получение «брутто» холодопроизводительности, определяется выражением

G„ =

Qop

Чо

и с учетом уравнения (9.43) может быть представлен в виде функции потока эксергии

 

т-------- ;------------- .

(9-44)

 

J [л> — сх (Т* Т3)]

 

где Гд = Тк — АТЖ, г0 = г0 (Г0), сх = сх

т°. j ; Тк — температура

конден-

сации; Т0 — температура кипения.

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

— ^2а -- hi.

Работа действительного процесса находится из соотношения

г д е г)н — индикаторный КПД компрессора.

268

Соседние файлы в папке книги