книги / Эксергетические расчеты технических систем
..pdfПри выборе в качестве целевой функции других показателей (например, суммарной массы или суммарного объема) величина 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, |
©К. А Тш); |
|||||||||
2ц = z n {Е2 [Е3(Сохл, Т охл, Ео.С, © |
в . |
АТв), 0 Н, AT’sl. ©к. АТш); |
||||||||||||
*12 = |
Z12[Е2 [Е3(Сохл, ЕОХЛ, Т о.С, ©в, АГВ), 0„, |
ATt], ©к, ATwy, |
||||||||||||
г13 =—^13 {Е2 [£'з(Фохл» ЕОХЛ, •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). Эта методика применима и для математического описания термоэкономической модели. Однако для полу-