Оцените пределы эффективности холодильной машины и теплообменника, опираясь на второй закон термодинамики и цикл Карно; какие реальные факторы (трение, теплопередача, утечки) приводят к отклонениям от идеальной эффективности и как их количественно учесть в расчёте?
Теоретический предел (цикл Карно) - Для холодильной машины (охлаждение при температуре TcT_cTc, сброс тепла при ThT_hTh, в К): COPCarnot,ref=TcTh−Tc\displaystyle \mathrm{COP}_{\text{Carnot,ref}}=\frac{T_c}{T_h-T_c}COPCarnot,ref=Th−TcTc. - Для теплового насоса (полезный тёплый поток при ThT_hTh): COPCarnot,hp=ThTh−Tc\displaystyle \mathrm{COP}_{\text{Carnot,hp}}=\frac{T_h}{T_h-T_c}COPCarnot,hp=Th−TcTh. (Все температуры в К; это следствие второго закона: при обратимом цикле энтропия не создаётся.) Учет реальных необратимостей через энтропию (общая формула) - Второй закон для рабочего цикла в стационарном режиме: S˙gen=Q˙HTh−Q˙LTc\displaystyle \dot{S}_{\text{gen}}=\frac{\dot{Q}_H}{T_h}-\frac{\dot{Q}_L}{T_c}S˙gen=ThQ˙H−TcQ˙L. - Энергетический баланс: W˙=Q˙H−Q˙L\displaystyle \dot{W}=\dot{Q}_H-\dot{Q}_LW˙=Q˙H−Q˙L. Отсюда выражение для реального COP холодильника с учётом ненулевой генерации энтропии: W˙=Q˙L(ThTc−1)+ThS˙gen\displaystyle \dot{W}=\dot{Q}_L\Big(\frac{T_h}{T_c}-1\Big)+T_h\dot{S}_{\text{gen}}W˙=Q˙L(TcTh−1)+ThS˙gen, COPreal=Q˙LW˙=1ThTc−1+ThS˙genQ˙L=TcTh−Tc+TcThS˙genQ˙L\displaystyle \mathrm{COP}_{\text{real}}=\frac{\dot{Q}_L}{\dot{W}}=\frac{1}{\frac{T_h}{T_c}-1+\dfrac{T_h\dot{S}_{\text{gen}}}{\dot{Q}_L}}=\frac{T_c}{T_h-T_c+\dfrac{T_cT_h\dot{S}_{\text{gen}}}{\dot{Q}_L}}COPreal=W˙Q˙L=TcTh−1+Q˙LThS˙gen1=Th−Tc+Q˙LTcThS˙genTc. Таким образом любое добавочное S˙gen>0\dot{S}_{\text{gen}}>0S˙gen>0 снижает COP относительно Карно. Главные реальные факторы и как их количественно включить 1) Нисходящий КПД компрессора / трение (механические необратимости) - Модель: изоэнтропная эффективность компрессора ηs\eta_sηs. Для массового расхода m˙\dot{m}m˙: W˙comp,act=m˙(h2s−h1)ηs\displaystyle \dot{W}_{\text{comp,act}}=\frac{\dot{m}(h_{2s}-h_1)}{\eta_s}W˙comp,act=ηsm˙(h2s−h1), где h2sh_{2s}h2s — энтальпия после идеального изоэнтропного сжатия. Разность W˙comp,act−m˙(h2s−h1)\dot{W}_{\text{comp,act}}-\dot{m}(h_{2s}-h_1)W˙comp,act−m˙(h2s−h1) входит в E˙d\dot{E}_{\text{d}}E˙d и через неё в S˙gen\dot{S}_{\text{gen}}S˙gen. 2) Теплопередача через конечный тепловой поток (конечный ΔT → тепловые необратимости) - Для теплообменника с коэффициентом теплопередачи UAUAUA: Q˙=UAΔTlm\displaystyle \dot{Q}=UA\Delta T_{\text{lm}}Q˙=UAΔTlm, где ΔTlm=ΔT1−ΔT2ln(ΔT1/ΔT2)\displaystyle \Delta T_{\text{lm}}=\frac{\Delta T_1-\Delta T_2}{\ln(\Delta T_1/\Delta T_2)}ΔTlm=ln(ΔT1/ΔT2)ΔT1−ΔT2. - Энтропия, создаваемая при передаче Q˙\dot{Q}Q˙ между температурами Th∗T_h^\astTh∗ и Tc∗T_c^\astTc∗: S˙gen,ht=Q˙(1Tc∗−1Th∗)\displaystyle \dot{S}_{\text{gen,ht}}=\dot{Q}\Big(\frac{1}{T_c^\ast}-\frac{1}{T_h^\ast}\Big)S˙gen,ht=Q˙(Tc∗1−Th∗1). Для малых ΔT\Delta TΔT можно приближённо: S˙gen,ht ≈Q˙ΔTTm2\displaystyle \dot{S}_{\text{gen,ht}\,}\approx\dot{Q}\frac{\Delta T}{T_m^2}S˙gen,ht≈Q˙Tm2ΔT. 3) Давление и потери в трубопроводах (гидравлические потери) - Потери давления Δp\Delta pΔp требуют доп. насосной/компрессорной работы W˙fric≈V˙Δp\dot{W}_{\text{fric}}\approx\dot{V}\Delta pW˙fric≈V˙Δp или включаются в снижение массового расхода; эта работа идёт в тепло (увеличивает S˙gen\dot{S}_{\text{gen}}S˙gen). 4) Утечки и смешение (массовые потери и утечка хладагента) - Модель: доля утечки α\alphaα → эффективный холодопроизводительность Q˙Leff=(1−α)Q˙L\dot{Q}_L^{\text{eff}}=(1-\alpha)\dot{Q}_LQ˙Leff=(1−α)Q˙L; при утечках часто требуется доп. подкачка/перекачка, увеличивается W˙\dot{W}W˙. 5) Недоминирование теплообменника (конечный NTU) — метод NTU/ε - Определение: NTU=UACmin\displaystyle NTU=\frac{UA}{C_{\min}}NTU=CminUA, Ci=m˙icp,iC_i=\dot{m}_i c_{p,i}Ci=m˙icp,i, C=CminCmaxC=\frac{C_{\min}}{C_{\max}}C=CmaxCmin. - Эффективность для counterflow: ϵ=1−exp[−NTU(1−C)]1−Cexp[−NTU(1−C)]\displaystyle \epsilon=\frac{1-\exp[-NTU(1-C)]}{1-C\exp[-NTU(1-C)]}ϵ=1−Cexp[−NTU(1−C)]1−exp[−NTU(1−C)]. Для C=1C=1C=1: ϵ=NTU1+NTU\displaystyle \epsilon=\frac{NTU}{1+NTU}ϵ=1+NTUNTU. - Тогда Q˙=ϵ Cmin(Th,in−Tc,in)\dot{Q}=\epsilon\,C_{\min}(T_{h,in}-T_{c,in})Q˙=ϵCmin(Th,in−Tc,in). Низкая ϵ\epsilonϵ уменьшает Q˙L\dot{Q}_LQ˙L и повышает относительную долю необратимостей. 6) Фазовые явления, недопущение идеального теплообмена (градусные потери при кипении/конденсации) - Модельируются через локальные тепловые сопротивления и температурные перепады на интерфейсах → увеличивают ΔTlm\Delta T_{\text{lm}}ΔTlm и S˙gen\dot{S}_{\text{gen}}S˙gen. Количественный алгоритм расчёта реальной эффективности 1. Для каждого компонента (компрессор, конденсатор, испаритель, трубопроводы) составить энергетический баланс и второе начало → найти S˙gen,i\dot{S}_{\text{gen},i}S˙gen,i. Примеры формул: - Компрессор: S˙gen,comp=W˙act−W˙revT0\dot{S}_{\text{gen,comp}}=\frac{\dot{W}_{\text{act}}-\dot{W}_{\text{rev}}}{T_{0}}S˙gen,comp=T0W˙act−W˙rev (приблизительно) или через реальные и идеальные энтальпии. - Теплообменник (с заданными температурами потоков): S˙gen,hex=m˙hcp,hlnTh,outTh,in+m˙ccp,clnTc,outTc,in\displaystyle \dot{S}_{\text{gen,hex}}=\dot{m}_h c_{p,h}\ln\frac{T_{h,out}}{T_{h,in}}+\dot{m}_c c_{p,c}\ln\frac{T_{c,out}}{T_{c,in}}S˙gen,hex=m˙hcp,hlnTh,inTh,out+m˙ccp,clnTc,inTc,out. 2. Суммировать: S˙gen=∑iS˙gen,i\displaystyle \dot{S}_{\text{gen}}=\sum_i\dot{S}_{\text{gen},i}S˙gen=i∑S˙gen,i. 3. Подставить в общую формулу COP (или в формулу COP теплового насоса) через найденную S˙gen\dot{S}_{\text{gen}}S˙gen: COPreal=TcTh−Tc+TcThS˙genQ˙L\displaystyle \mathrm{COP}_{\text{real}}=\frac{T_c}{T_h-T_c+\dfrac{T_cT_h\dot{S}_{\text{gen}}}{\dot{Q}_L}}COPreal=Th−Tc+Q˙LTcThS˙genTc. Альтернативный практический подход: разбить необратимости на эксергетические потери и учесть их как дополнительные «встроенные» работы/потери: - Компрессор: учесть ηs\eta_sηs → повысит W˙\dot{W}W˙. - Теплообменники: через UAUAUA и NTU → уменьшат Q˙L\dot{Q}_LQ˙L. - Давление/утечки: прямое добавление потерь к работе/уменьшение Q˙L\dot{Q}_LQ˙L. Краткое резюме - Карно даёт верхнюю границу: COPCarnot,ref=TcTh−Tc\mathrm{COP}_{\text{Carnot,ref}}=\dfrac{T_c}{T_h-T_c}COPCarnot,ref=Th−TcTc. - Любая генерация энтропии (S˙gen>0\dot{S}_{\text{gen}}>0S˙gen>0) снижает реальный COP согласно формуле через S˙gen\dot{S}_{\text{gen}}S˙gen. - Реальные источники необратимостей: конечная теплопередача (UA, NTU, ε), эффективность компрессора/трение, гидравлические потери, утечки, фазовые и конструктивные ограничения — все их можно включить количественно через модели (η_s, UA→NTU→ε, дополнительные работы/потери и расчёт S˙gen\dot{S}_{\text{gen}}S˙gen).
- Для холодильной машины (охлаждение при температуре TcT_cTc , сброс тепла при ThT_hTh , в К):
COPCarnot,ref=TcTh−Tc\displaystyle \mathrm{COP}_{\text{Carnot,ref}}=\frac{T_c}{T_h-T_c}COPCarnot,ref =Th −Tc Tc .
- Для теплового насоса (полезный тёплый поток при ThT_hTh ):
COPCarnot,hp=ThTh−Tc\displaystyle \mathrm{COP}_{\text{Carnot,hp}}=\frac{T_h}{T_h-T_c}COPCarnot,hp =Th −Tc Th .
(Все температуры в К; это следствие второго закона: при обратимом цикле энтропия не создаётся.)
Учет реальных необратимостей через энтропию (общая формула)
- Второй закон для рабочего цикла в стационарном режиме:
S˙gen=Q˙HTh−Q˙LTc\displaystyle \dot{S}_{\text{gen}}=\frac{\dot{Q}_H}{T_h}-\frac{\dot{Q}_L}{T_c}S˙gen =Th Q˙ H −Tc Q˙ L .
- Энергетический баланс: W˙=Q˙H−Q˙L\displaystyle \dot{W}=\dot{Q}_H-\dot{Q}_LW˙=Q˙ H −Q˙ L .
Отсюда выражение для реального COP холодильника с учётом ненулевой генерации энтропии:
W˙=Q˙L(ThTc−1)+ThS˙gen\displaystyle \dot{W}=\dot{Q}_L\Big(\frac{T_h}{T_c}-1\Big)+T_h\dot{S}_{\text{gen}}W˙=Q˙ L (Tc Th −1)+Th S˙gen ,
COPreal=Q˙LW˙=1ThTc−1+ThS˙genQ˙L=TcTh−Tc+TcThS˙genQ˙L\displaystyle \mathrm{COP}_{\text{real}}=\frac{\dot{Q}_L}{\dot{W}}=\frac{1}{\frac{T_h}{T_c}-1+\dfrac{T_h\dot{S}_{\text{gen}}}{\dot{Q}_L}}=\frac{T_c}{T_h-T_c+\dfrac{T_cT_h\dot{S}_{\text{gen}}}{\dot{Q}_L}}COPreal =W˙Q˙ L =Tc Th −1+Q˙ L Th S˙gen 1 =Th −Tc +Q˙ L Tc Th S˙gen Tc .
Таким образом любое добавочное S˙gen>0\dot{S}_{\text{gen}}>0S˙gen >0 снижает COP относительно Карно.
Главные реальные факторы и как их количественно включить
1) Нисходящий КПД компрессора / трение (механические необратимости)
- Модель: изоэнтропная эффективность компрессора ηs\eta_sηs . Для массового расхода m˙\dot{m}m˙:
W˙comp,act=m˙(h2s−h1)ηs\displaystyle \dot{W}_{\text{comp,act}}=\frac{\dot{m}(h_{2s}-h_1)}{\eta_s}W˙comp,act =ηs m˙(h2s −h1 ) , где h2sh_{2s}h2s — энтальпия после идеального изоэнтропного сжатия. Разность W˙comp,act−m˙(h2s−h1)\dot{W}_{\text{comp,act}}-\dot{m}(h_{2s}-h_1)W˙comp,act −m˙(h2s −h1 ) входит в E˙d\dot{E}_{\text{d}}E˙d и через неё в S˙gen\dot{S}_{\text{gen}}S˙gen .
2) Теплопередача через конечный тепловой поток (конечный ΔT → тепловые необратимости)
- Для теплообменника с коэффициентом теплопередачи UAUAUA: Q˙=UAΔTlm\displaystyle \dot{Q}=UA\Delta T_{\text{lm}}Q˙ =UAΔTlm , где ΔTlm=ΔT1−ΔT2ln(ΔT1/ΔT2)\displaystyle \Delta T_{\text{lm}}=\frac{\Delta T_1-\Delta T_2}{\ln(\Delta T_1/\Delta T_2)}ΔTlm =ln(ΔT1 /ΔT2 )ΔT1 −ΔT2 .
- Энтропия, создаваемая при передаче Q˙\dot{Q}Q˙ между температурами Th∗T_h^\astTh∗ и Tc∗T_c^\astTc∗ :
S˙gen,ht=Q˙(1Tc∗−1Th∗)\displaystyle \dot{S}_{\text{gen,ht}}=\dot{Q}\Big(\frac{1}{T_c^\ast}-\frac{1}{T_h^\ast}\Big)S˙gen,ht =Q˙ (Tc∗ 1 −Th∗ 1 ). Для малых ΔT\Delta TΔT можно приближённо: S˙gen,ht ≈Q˙ΔTTm2\displaystyle \dot{S}_{\text{gen,ht}\,}\approx\dot{Q}\frac{\Delta T}{T_m^2}S˙gen,ht ≈Q˙ Tm2 ΔT .
3) Давление и потери в трубопроводах (гидравлические потери)
- Потери давления Δp\Delta pΔp требуют доп. насосной/компрессорной работы W˙fric≈V˙Δp\dot{W}_{\text{fric}}\approx\dot{V}\Delta pW˙fric ≈V˙Δp или включаются в снижение массового расхода; эта работа идёт в тепло (увеличивает S˙gen\dot{S}_{\text{gen}}S˙gen ).
4) Утечки и смешение (массовые потери и утечка хладагента)
- Модель: доля утечки α\alphaα → эффективный холодопроизводительность Q˙Leff=(1−α)Q˙L\dot{Q}_L^{\text{eff}}=(1-\alpha)\dot{Q}_LQ˙ Leff =(1−α)Q˙ L ; при утечках часто требуется доп. подкачка/перекачка, увеличивается W˙\dot{W}W˙.
5) Недоминирование теплообменника (конечный NTU) — метод NTU/ε
- Определение: NTU=UACmin\displaystyle NTU=\frac{UA}{C_{\min}}NTU=Cmin UA , Ci=m˙icp,iC_i=\dot{m}_i c_{p,i}Ci =m˙i cp,i , C=CminCmaxC=\frac{C_{\min}}{C_{\max}}C=Cmax Cmin .
- Эффективность для counterflow:
ϵ=1−exp[−NTU(1−C)]1−Cexp[−NTU(1−C)]\displaystyle \epsilon=\frac{1-\exp[-NTU(1-C)]}{1-C\exp[-NTU(1-C)]}ϵ=1−Cexp[−NTU(1−C)]1−exp[−NTU(1−C)] .
Для C=1C=1C=1: ϵ=NTU1+NTU\displaystyle \epsilon=\frac{NTU}{1+NTU}ϵ=1+NTUNTU .
- Тогда Q˙=ϵ Cmin(Th,in−Tc,in)\dot{Q}=\epsilon\,C_{\min}(T_{h,in}-T_{c,in})Q˙ =ϵCmin (Th,in −Tc,in ). Низкая ϵ\epsilonϵ уменьшает Q˙L\dot{Q}_LQ˙ L и повышает относительную долю необратимостей.
6) Фазовые явления, недопущение идеального теплообмена (градусные потери при кипении/конденсации)
- Модельируются через локальные тепловые сопротивления и температурные перепады на интерфейсах → увеличивают ΔTlm\Delta T_{\text{lm}}ΔTlm и S˙gen\dot{S}_{\text{gen}}S˙gen .
Количественный алгоритм расчёта реальной эффективности
1. Для каждого компонента (компрессор, конденсатор, испаритель, трубопроводы) составить энергетический баланс и второе начало → найти S˙gen,i\dot{S}_{\text{gen},i}S˙gen,i . Примеры формул:
- Компрессор: S˙gen,comp=W˙act−W˙revT0\dot{S}_{\text{gen,comp}}=\frac{\dot{W}_{\text{act}}-\dot{W}_{\text{rev}}}{T_{0}}S˙gen,comp =T0 W˙act −W˙rev (приблизительно) или через реальные и идеальные энтальпии.
- Теплообменник (с заданными температурами потоков):
S˙gen,hex=m˙hcp,hlnTh,outTh,in+m˙ccp,clnTc,outTc,in\displaystyle \dot{S}_{\text{gen,hex}}=\dot{m}_h c_{p,h}\ln\frac{T_{h,out}}{T_{h,in}}+\dot{m}_c c_{p,c}\ln\frac{T_{c,out}}{T_{c,in}}S˙gen,hex =m˙h cp,h lnTh,in Th,out +m˙c cp,c lnTc,in Tc,out .
2. Суммировать: S˙gen=∑iS˙gen,i\displaystyle \dot{S}_{\text{gen}}=\sum_i\dot{S}_{\text{gen},i}S˙gen =i∑ S˙gen,i .
3. Подставить в общую формулу COP (или в формулу COP теплового насоса) через найденную S˙gen\dot{S}_{\text{gen}}S˙gen :
COPreal=TcTh−Tc+TcThS˙genQ˙L\displaystyle \mathrm{COP}_{\text{real}}=\frac{T_c}{T_h-T_c+\dfrac{T_cT_h\dot{S}_{\text{gen}}}{\dot{Q}_L}}COPreal =Th −Tc +Q˙ L Tc Th S˙gen Tc .
Альтернативный практический подход: разбить необратимости на эксергетические потери и учесть их как дополнительные «встроенные» работы/потери:
- Компрессор: учесть ηs\eta_sηs → повысит W˙\dot{W}W˙.
- Теплообменники: через UAUAUA и NTU → уменьшат Q˙L\dot{Q}_LQ˙ L .
- Давление/утечки: прямое добавление потерь к работе/уменьшение Q˙L\dot{Q}_LQ˙ L .
Краткое резюме
- Карно даёт верхнюю границу: COPCarnot,ref=TcTh−Tc\mathrm{COP}_{\text{Carnot,ref}}=\dfrac{T_c}{T_h-T_c}COPCarnot,ref =Th −Tc Tc .
- Любая генерация энтропии (S˙gen>0\dot{S}_{\text{gen}}>0S˙gen >0) снижает реальный COP согласно формуле через S˙gen\dot{S}_{\text{gen}}S˙gen .
- Реальные источники необратимостей: конечная теплопередача (UA, NTU, ε), эффективность компрессора/трение, гидравлические потери, утечки, фазовые и конструктивные ограничения — все их можно включить количественно через модели (η_s, UA→NTU→ε, дополнительные работы/потери и расчёт S˙gen\dot{S}_{\text{gen}}S˙gen ).