Кейс: тепловой двигатель работает в режиме конечного времени между двумя резервуарами с конечной теплоёмкостью — как оптимизировать цикл для максимальной мощности и какой компромисс с КПД возникает
Кратко — модель, оптимизация, результат и компромисс. 1) Модель (минимально необходимая постановка) - Два резервуара: начальные температуры Th0>Tc0T_{h0}>T_{c0}Th0>Tc0 и удельные теплоёмкости Ch, CcC_h,\;C_cCh,Cc. - Тепловой двигатель эндореверсивен (внутри обратим), контакты описываются законом теплопередачи (линейно): Q˙h=Kh(Th(t)−Thw(t)),Q˙c=Kc(Tcw(t)−Tc(t))\dot Q_h=K_h\big(T_h(t)-T_{hw}(t)\big),\qquad \dot Q_c=K_c\big(T_{cw}(t)-T_c(t)\big)Q˙h=Kh(Th(t)−Thw(t)),Q˙c=Kc(Tcw(t)−Tc(t)). - Внутренняя обратимость даёт связь потоков: Q˙hThw=Q˙cTcw\dfrac{\dot Q_h}{T_{hw}}=\dfrac{\dot Q_c}{T_{cw}}ThwQ˙h=TcwQ˙c. - Температуры резервуаров меняются по теплоёмности: T˙h=−Q˙hCh,T˙c=+Q˙cCc\dot T_h=-\dfrac{\dot Q_h}{C_h},\qquad \dot T_c=+\dfrac{\dot Q_c}{C_c}T˙h=−ChQ˙h,T˙c=+CcQ˙c. - Мощность в момент времени: P(t)=Q˙h(1−TcwThw)P(t)=\dot Q_h\Big(1-\dfrac{T_{cw}}{T_{hw}}\Big)P(t)=Q˙h(1−ThwTcw). 2) Как оптимизировать для максимальной мощности - Общая идея: в каждом момент оптимально выбирать управляющие переменные Thw(t),Tcw(t)T_{hw}(t),T_{cw}(t)Thw(t),Tcw(t) так, чтобы максимизировать мгновенную мощность P(t)P(t)P(t) при данных Th(t),Tc(t)T_h(t),T_c(t)Th(t),Tc(t). Затем интегрировать во времени, пока резервуары не исчерпаются. Формально это задача управления, решаемая методом Понтрягина; для линейной теплопередачи решение даёт простую условную стратегию: поддерживать внутренние температуры (или их отношение) постоянными на интервале и выбирать их из условия максимума P при данных Th,TcT_h,T_cTh,Tc. - Частный, хорошо известный случай — симметричные проводимости Kh=KcK_h=K_cKh=Kc и фиксированные резервуарные температуры: при стационарном режиме оптимум по Thw,TcwT_{hw},T_{cw}Thw,Tcw даёт Curzon–Ahlborn эффективность ηCA=1−TcTh.
\eta_{CA}=1-\sqrt{\dfrac{T_c}{T_h}}. ηCA=1−ThTc.
То есть внутренние температуры выбираются так, что эффективность при максимальной мощности равна ηCA\eta_{CA}ηCA. - Для конечных теплоёмкостей резервуаров оптимальное управление фактически сводится к последовательному применению локально оптимальных внутренних температур (решение «оптимум по моменту времени»), затем интегрированию мощности по времени с учётом убывающих Th(t),Tc(t)T_h(t),T_c(t)Th(t),Tc(t). 3) Формальные условия стационарного оптимума (резюме) - Из вариационного условия (оптимизация по Thw,TcwT_{hw},T_{cw}Thw,Tcw) получаем систему алгебраических уравнений: {∂P∂Thw=0,∂P∂Tcw=0,где P=Kh(Th−Thw)(1−TcwThw) и Kh(Th−Thw)Thw=Kc(Tcw−Tc)Tcw.
\begin{cases} \dfrac{\partial P}{\partial T_{hw}}=0,\\[4pt] \dfrac{\partial P}{\partial T_{cw}}=0, \end{cases} \quad\text{где }P=K_h(T_h-T_{hw})\Big(1-\dfrac{T_{cw}}{T_{hw}}\Big) \text{ и } \dfrac{K_h(T_h-T_{hw})}{T_{hw}}=\dfrac{K_c(T_{cw}-T_c)}{T_{cw}}. ⎩⎨⎧∂Thw∂P=0,∂Tcw∂P=0,гдеP=Kh(Th−Thw)(1−ThwTcw)иThwKh(Th−Thw)=TcwKc(Tcw−Tc).
Решение даёт оптимальные внутренние температуры как функции Th,Tc,Kh,KcT_h,T_c,K_h,K_cTh,Tc,Kh,Kc. 4) Компромисс мощность ⇄ КПД - Общий принцип: максимум мощности достигается при существенной генерируемой необратимости (большие перепады температур между резервуаром и рабочим телом), поэтому КПД при максимальной мощности меньше карно: ηMP<ηC=1−TcTh.
\eta_{MP}<\eta_C=1-\dfrac{T_c}{T_h}. ηMP<ηC=1−ThTc.
- В предельном симметричном случае получают Curzon–Ahlborn: ηMP=ηCA=1−Tc/Th\eta_{MP}=\eta_{CA}=1-\sqrt{T_c/T_h}ηMP=ηCA=1−Tc/Th. Для малых ηC\eta_CηC это даёт универсальную первую аппроксимацию ηMP≈ηC2+O(ηC2).
\eta_{MP}\approx\dfrac{\eta_C}{2}+O(\eta_C^2). ηMP≈2ηC+O(ηC2).
- Для конечных теплоёмкостей точное значение ηMP\eta_{MP}ηMP зависит от соотношения Ch/CcC_h/C_cCh/Cc и от Kh/KcK_h/K_cKh/Kc: асимметрия смещает оптимум, и эффективность при максимуме может быть выше или ниже ηCA\eta_{CA}ηCA. При стремлении к обратимому (медленному) режиму мощность → 0, КПД → ηC\eta_CηC. 5) Практические рекомендации - Если цель — максимальная мощность за конечное время до выравнивания резервуаров — в каждый момент используйте локально оптимальные Thw,TcwT_{hw},T_{cw}Thw,Tcw (решение выше) и распределяйте контактные времена согласно проводимостям и теплоёмкостям. - Если важен КПД — уменьшайте теплоперепад (работайте медленнее), приближаясь к обратимому циклу; при этом мощность падает примерно обратно пропорционально времени цикла. - Для количественных оптимальных значений (Thw(t),Tcw(t)T_{hw}(t),T_{cw}(t)Thw(t),Tcw(t), времена контактов) решите систему ОДУ и алгебраических условий или примените численную оптимизацию (Pontryagin), задав конкретные Ch,Cc,Kh,Kc,Th0,Tc0C_h,C_c,K_h,K_c,T_{h0},T_{c0}Ch,Cc,Kh,Kc,Th0,Tc0. Коротко: оптимизация сводится к локальному максимуму мгновенной мощности при каждой паре температур резервуаров (эндореверсивная модель) — в симметричном пределе это даёт Curzon–Ahlborn эффективность ηCA=1−Tc/Th\eta_{CA}=1-\sqrt{T_c/T_h}ηCA=1−Tc/Th. Компромисс: больше мощности — меньше КПД; КПД при максимальной мощности всегда меньше Карно и приближается к нему только при бесконечно медленной (нулевой мощности) работе.
1) Модель (минимально необходимая постановка)
- Два резервуара: начальные температуры Th0>Tc0T_{h0}>T_{c0}Th0 >Tc0 и удельные теплоёмкости Ch, CcC_h,\;C_cCh ,Cc .
- Тепловой двигатель эндореверсивен (внутри обратим), контакты описываются законом теплопередачи (линейно):
Q˙h=Kh(Th(t)−Thw(t)),Q˙c=Kc(Tcw(t)−Tc(t))\dot Q_h=K_h\big(T_h(t)-T_{hw}(t)\big),\qquad \dot Q_c=K_c\big(T_{cw}(t)-T_c(t)\big)Q˙ h =Kh (Th (t)−Thw (t)),Q˙ c =Kc (Tcw (t)−Tc (t)).
- Внутренняя обратимость даёт связь потоков: Q˙hThw=Q˙cTcw\dfrac{\dot Q_h}{T_{hw}}=\dfrac{\dot Q_c}{T_{cw}}Thw Q˙ h =Tcw Q˙ c .
- Температуры резервуаров меняются по теплоёмности:
T˙h=−Q˙hCh,T˙c=+Q˙cCc\dot T_h=-\dfrac{\dot Q_h}{C_h},\qquad \dot T_c=+\dfrac{\dot Q_c}{C_c}T˙h =−Ch Q˙ h ,T˙c =+Cc Q˙ c .
- Мощность в момент времени: P(t)=Q˙h(1−TcwThw)P(t)=\dot Q_h\Big(1-\dfrac{T_{cw}}{T_{hw}}\Big)P(t)=Q˙ h (1−Thw Tcw ).
2) Как оптимизировать для максимальной мощности
- Общая идея: в каждом момент оптимально выбирать управляющие переменные Thw(t),Tcw(t)T_{hw}(t),T_{cw}(t)Thw (t),Tcw (t) так, чтобы максимизировать мгновенную мощность P(t)P(t)P(t) при данных Th(t),Tc(t)T_h(t),T_c(t)Th (t),Tc (t). Затем интегрировать во времени, пока резервуары не исчерпаются. Формально это задача управления, решаемая методом Понтрягина; для линейной теплопередачи решение даёт простую условную стратегию: поддерживать внутренние температуры (или их отношение) постоянными на интервале и выбирать их из условия максимума P при данных Th,TcT_h,T_cTh ,Tc .
- Частный, хорошо известный случай — симметричные проводимости Kh=KcK_h=K_cKh =Kc и фиксированные резервуарные температуры: при стационарном режиме оптимум по Thw,TcwT_{hw},T_{cw}Thw ,Tcw даёт Curzon–Ahlborn эффективность
ηCA=1−TcTh. \eta_{CA}=1-\sqrt{\dfrac{T_c}{T_h}}.
ηCA =1−Th Tc . То есть внутренние температуры выбираются так, что эффективность при максимальной мощности равна ηCA\eta_{CA}ηCA .
- Для конечных теплоёмкостей резервуаров оптимальное управление фактически сводится к последовательному применению локально оптимальных внутренних температур (решение «оптимум по моменту времени»), затем интегрированию мощности по времени с учётом убывающих Th(t),Tc(t)T_h(t),T_c(t)Th (t),Tc (t).
3) Формальные условия стационарного оптимума (резюме)
- Из вариационного условия (оптимизация по Thw,TcwT_{hw},T_{cw}Thw ,Tcw ) получаем систему алгебраических уравнений:
{∂P∂Thw=0,∂P∂Tcw=0,где P=Kh(Th−Thw)(1−TcwThw) и Kh(Th−Thw)Thw=Kc(Tcw−Tc)Tcw. \begin{cases}
\dfrac{\partial P}{\partial T_{hw}}=0,\\[4pt]
\dfrac{\partial P}{\partial T_{cw}}=0,
\end{cases}
\quad\text{где }P=K_h(T_h-T_{hw})\Big(1-\dfrac{T_{cw}}{T_{hw}}\Big)
\text{ и } \dfrac{K_h(T_h-T_{hw})}{T_{hw}}=\dfrac{K_c(T_{cw}-T_c)}{T_{cw}}.
⎩⎨⎧ ∂Thw ∂P =0,∂Tcw ∂P =0, где P=Kh (Th −Thw )(1−Thw Tcw ) и Thw Kh (Th −Thw ) =Tcw Kc (Tcw −Tc ) . Решение даёт оптимальные внутренние температуры как функции Th,Tc,Kh,KcT_h,T_c,K_h,K_cTh ,Tc ,Kh ,Kc .
4) Компромисс мощность ⇄ КПД
- Общий принцип: максимум мощности достигается при существенной генерируемой необратимости (большие перепады температур между резервуаром и рабочим телом), поэтому КПД при максимальной мощности меньше карно:
ηMP<ηC=1−TcTh. \eta_{MP}<\eta_C=1-\dfrac{T_c}{T_h}.
ηMP <ηC =1−Th Tc . - В предельном симметричном случае получают Curzon–Ahlborn:
ηMP=ηCA=1−Tc/Th\eta_{MP}=\eta_{CA}=1-\sqrt{T_c/T_h}ηMP =ηCA =1−Tc /Th . Для малых ηC\eta_CηC это даёт универсальную первую аппроксимацию
ηMP≈ηC2+O(ηC2). \eta_{MP}\approx\dfrac{\eta_C}{2}+O(\eta_C^2).
ηMP ≈2ηC +O(ηC2 ). - Для конечных теплоёмкостей точное значение ηMP\eta_{MP}ηMP зависит от соотношения Ch/CcC_h/C_cCh /Cc и от Kh/KcK_h/K_cKh /Kc : асимметрия смещает оптимум, и эффективность при максимуме может быть выше или ниже ηCA\eta_{CA}ηCA . При стремлении к обратимому (медленному) режиму мощность → 0, КПД → ηC\eta_CηC .
5) Практические рекомендации
- Если цель — максимальная мощность за конечное время до выравнивания резервуаров — в каждый момент используйте локально оптимальные Thw,TcwT_{hw},T_{cw}Thw ,Tcw (решение выше) и распределяйте контактные времена согласно проводимостям и теплоёмкостям.
- Если важен КПД — уменьшайте теплоперепад (работайте медленнее), приближаясь к обратимому циклу; при этом мощность падает примерно обратно пропорционально времени цикла.
- Для количественных оптимальных значений (Thw(t),Tcw(t)T_{hw}(t),T_{cw}(t)Thw (t),Tcw (t), времена контактов) решите систему ОДУ и алгебраических условий или примените численную оптимизацию (Pontryagin), задав конкретные Ch,Cc,Kh,Kc,Th0,Tc0C_h,C_c,K_h,K_c,T_{h0},T_{c0}Ch ,Cc ,Kh ,Kc ,Th0 ,Tc0 .
Коротко: оптимизация сводится к локальному максимуму мгновенной мощности при каждой паре температур резервуаров (эндореверсивная модель) — в симметричном пределе это даёт Curzon–Ahlborn эффективность ηCA=1−Tc/Th\eta_{CA}=1-\sqrt{T_c/T_h}ηCA =1−Tc /Th . Компромисс: больше мощности — меньше КПД; КПД при максимальной мощности всегда меньше Карно и приближается к нему только при бесконечно медленной (нулевой мощности) работе.