통계물리 1 보충 자료: 카르노 엔진 너머
Chambadal–Novikov–Curzon–Ahlborn 엔진
유한한 시간 안에 일률을 내야 하는 카르노 엔진입니다. 아래에서 작동 온도를 움직이며 효율과 일률 사이의 절충 관계를 직접 확인한 뒤, 유도 과정을 단계별로 따라가 보세요.
모형 세우기
준정적으로 작동하는 카르노 엔진은 \(\eta_C = 1 - T_c/T_h\)에 도달하지만, 이는 각 등온 과정에 무한한 시간이 걸리기 때문이며 그 결과 일률은 0입니다. 유한한 속도로 열을 주고받으려면 작동 물질이 고온 열원보다 차갑고(\(T_{hw} < T_h\)) 저온 열원보다 뜨거워야(\(T_{cw} > T_c\)) 합니다. Schroeder를 따라, 열 흐름은 양쪽에서 같은 열전도 계수 \(K\)를 갖는 선형 열전도 법칙을 따르고 두 등온 과정은 같은 시간 \(\Delta t\) 동안 일어난다고 가정합니다.
\[ \frac{Q_h}{\Delta t} = K\,(T_h - T_{hw}), \qquad \frac{Q_c}{\Delta t} = K\,(T_{cw} - T_c). \]엔진 내부는 모두 가역적이고, 비가역성은 오직 열을 전달하는 두 벽에서만 생깁니다. 그래서 이 모형을 내부가역(endoreversible) 엔진이라고 부릅니다.
이름에 관하여. \(\eta = 1 - \sqrt{T_c/T_h}\)라는 결과는 1975년 논문의 저자 이름을 따서 흔히 Curzon–Ahlborn 효율이라 불리며, 이 논문은 Schroeder 교과서 연습문제 4.6의 바탕이기도 합니다. 그러나 같은 식이 이미 1957년 Pierre Chambadal과 Igor Novikov의 원자력 발전소 효율 연구에 등장했으므로, 이 페이지에서는 엔진을 Chambadal–Novikov–Curzon–Ahlborn 엔진이라 부릅니다. 기호 \(\eta_{CA}\)는 인용 논문들과 맞추기 위해 그대로 사용합니다.
P–V 도표로 본 실제 작동 순환
벽 안쪽에서 작동 물질은 보통의 카르노 순환을 합니다. 다만 \(T_h\)와 \(T_c\)가 아니라 \(T_{hw}\)와 \(T_{cw}\) 사이에서 작동합니다. 아래 도표는 단원자 이상기체(\(\gamma = 5/3\))에 대해 이 순환을 그린 것으로, \(Nk = 1\)이고 고온 등온 과정이 시작될 때의 부피를 \(V_1 = 1\)로 두었습니다. 일률 \(P\)와 구별하기 위해 압력은 \(p\)로 씁니다.
순환은 네 과정으로 이루어집니다. \(T_{hw}\)에서의 등온 팽창(1→2), 단열 팽창(2→3), \(T_{cw}\)에서의 등온 압축(3→4), 단열 압축(4→1)입니다. 각 등온 과정은 같은 시간 \(\Delta t\) 동안 진행되므로 흡수하는 열은 \(Q_h = K\Delta t\,(T_h - T_{hw}) = NkT_{hw}\ln(V_2/V_1)\)이고, 이것이 팽창비를 결정합니다. 단열 과정은 \(V_3/V_2 = V_4/V_1 = (T_{hw}/T_{cw})^{1/(\gamma-1)}\)을 정합니다.
\(\Delta t\)가 고정되어 있으므로 둘러싼 넓이는 일률에 비례합니다(\(P = W/2\Delta t\)). 따라서 순환의 넓이는 정확히 최대 일률 지점에서 가장 크고, 두 극한에서는 모두 사라집니다. 단락 극한에서는 두 등온선이 합쳐져 순환이 하나의 등온선 위로 납작해집니다. 열은 많이 지나가지만 넓이가 없습니다. 준정적 극한에서는 고정된 \(\Delta t\)에서 \(Q_h \to 0\)이므로 등온 과정이 점으로 줄어들어 순환이 단열선 위로 무너집니다. 교과서의 카르노 순환이 이 극한에서도 유한한 넓이를 갖는 것은 \(\Delta t \to \infty\)를 허용하기 때문입니다.
최대 일률에서의 효율 유도
엔트로피 균형
한 순환이 끝나면 작동 물질은 처음 상태로 돌아오므로 엔트로피 변화가 0입니다. 내부가 가역적이므로 \(T_{hw}\)에서 흡수한 엔트로피는 \(T_{cw}\)에서 방출한 엔트로피와 같아야 합니다.
\[ \frac{Q_h}{T_{hw}} = \frac{Q_c}{T_{cw}}. \]열전도 법칙을 대입하면 (\(K\Delta t\)는 약분됩니다)
\[ \frac{T_h - T_{hw}}{T_{hw}} = \frac{T_{cw} - T_c}{T_{cw}} \quad\Longleftrightarrow\quad \boxed{\;\frac{T_h}{T_{hw}} + \frac{T_c}{T_{cw}} = 2\;} \]이를 저온 쪽 작동 온도에 대해 풀면 아래에서 쓸 형태를 얻습니다.
\[ T_{cw} = \frac{T_c\,T_{hw}}{2T_{hw} - T_h}. \]엔진이 일을 하려면 \(T_{cw} < T_{hw}\)여야 하고, 이는 \(\tfrac12(T_h + T_c) < T_{hw} \le T_h\)에 해당합니다. 아래 끝에서는 \(T_{cw} = T_{hw}\)(단락), 위 끝에서는 \(T_{hw} = T_h\), \(T_{cw} = T_c\)(준정적 카르노)가 됩니다.
일률
단열 과정에 걸리는 시간을 무시하면 한 순환은 \(2\Delta t\) 동안 진행됩니다. 열역학 제1법칙에 의해 \(W = Q_h - Q_c\)이므로
\[ P = \frac{Q_h - Q_c}{2\Delta t} = \frac{K}{2}\Big[(T_h - T_{hw}) - (T_{cw} - T_c)\Big]. \]1단계의 제2법칙 관계식으로 \(T_{cw}\)를 소거하면
일률 최대화
\(T_h\)와 \(T_c\)를 고정하고 미분합니다. \(\dfrac{d}{dT_{hw}}\dfrac{T_{hw}}{2T_{hw}-T_h} = \dfrac{-T_h}{(2T_{hw}-T_h)^2}\)를 이용하면
\[ \frac{dP}{dT_{hw}} = \frac{K}{2}\left[-1 + \frac{T_h T_c}{(2T_{hw} - T_h)^2}\right] = 0 \quad\Longrightarrow\quad 4T_{hw}^2 - 4T_h T_{hw} + T_h^2 - T_h T_c = 0. \]이 2차방정식의 해는 \(T_{hw} = \tfrac12\big(T_h \pm \sqrt{T_h T_c}\big)\)입니다. 음의 근은 \(\tfrac12(T_h+T_c)\)보다 작아서 \(T_{cw}\)를 음수로 만들므로 물리적이지 않습니다. \(d^2P/dT_{hw}^2 = -2K T_hT_c/(2T_{hw}-T_h)^3 < 0\)이므로 양의 근이 최댓값입니다.
각 작동 온도는 해당 열원의 온도와 기하평균 \(\sqrt{T_hT_c}\)의 정확히 중간에 놓입니다. 다시 대입하면 최대 일률은
\[ P_{\max} = \frac{K}{4}\left(\sqrt{T_h} - \sqrt{T_c}\right)^2 \]이며, 이는 Curzon–Ahlborn 논문의 식 (19)에서 \(\alpha = \beta = K\)로 둔 것과 같습니다.
최대 일률에서의 효율
효율은 제1법칙과 1단계의 엔트로피 균형으로부터 바로 나옵니다.
\[ \eta = 1 - \frac{Q_c}{Q_h} = 1 - \frac{T_{cw}^*}{T_{hw}^*} = 1 - \frac{\sqrt{T_c}\,(\sqrt{T_c} + \sqrt{T_h})}{\sqrt{T_h}\,(\sqrt{T_h} + \sqrt{T_c})} \quad\Longrightarrow\quad \boxed{\;\eta_{CA} = 1 - \sqrt{T_c/T_h}\;} \]석탄 화력 증기 터빈의 경우 \(T_h = 600\,^\circ\mathrm{C} = 873.15\ \mathrm K\), \(T_c = 25\,^\circ\mathrm{C} = 298.15\ \mathrm K\)이므로
\[ \eta_{CA} = 1 - \sqrt{\tfrac{298.15}{873.15}} \approx 0.416, \qquad \eta_C = 1 - \tfrac{298.15}{873.15} \approx 0.659. \]최대 일률에서 작동 온도는 \(T_{hw}^* \approx 691.7\ \mathrm K\)(418.5 °C), \(T_{cw}^* \approx 404.2\ \mathrm K\)(131.0 °C)입니다. 약 42%인 Curzon–Ahlborn 값은 66%인 카르노 한계보다 실제 석탄 발전소의 효율(약 40%)에 훨씬 가깝습니다. 실제 발전소는 준정적 극한에 다가가는 것이 아니라 일률을 내도록 설계되고 운전되기 때문입니다.
엔트로피는 어디서 생기는가
작동 물질의 엔트로피는 매 순환마다 원래 값으로 돌아오지만 열원은 그렇지 않습니다. 고온 열원이 잃는 \(Q_h/T_h\)는 엔진이 흡수하는 \(Q_h/T_{hw}\)보다 작고, 저온 열원이 얻는 \(Q_c/T_c\)는 엔진이 방출하는 \(Q_c/T_{cw}\)보다 큽니다. 따라서 한 순환 동안 생성되는 전체 엔트로피는
\[ \Delta S_{\rm tot} = \frac{Q_c}{T_c} - \frac{Q_h}{T_h} > -\frac{Q_h}{T_{hw}} + \frac{Q_c}{T_{cw}} = 0. \]이를 제1법칙과 결합하면 간결한 에너지 수지(Gouy–Stodola 관계)를 얻으며, 데모의 색 막대가 바로 이것을 보여줍니다.
\[ W = \eta_C\,Q_h - T_c\,\Delta S_{\rm tot}. \]벽에서 생성된 엔트로피는 그만큼 \(T_c\)를 곱한 양의 일을 잃게 만듭니다. 슬라이더의 양 끝 두 극한에서 이 절충 관계가 분명해집니다.
준정적 극한(\(T_{hw}\to T_h,\ T_{cw}\to T_c\)). 온도 차이가 사라져 엔트로피가 생성되지 않고 \(\eta \to \eta_C\)가 됩니다. 하지만 유한한 열을 옮기려면 열전도 법칙에 의해 \(\Delta t \to \infty\)가 되어야 하므로 일률은 0으로 갑니다.
단락 극한(\(T_{hw}\to T_{cw}\)). 열은 최대한 빨리 흐르지만, 내부의 카르노 엔진이 같은 두 온도 사이에서 작동하게 되어 일을 전혀 하지 못하고 열을 고온에서 저온으로 그냥 통과시킵니다. 전지를 도선으로 단락시킨 것과 비슷합니다. 일률과 효율이 모두 0이 됩니다.
최대 일률 지점은 이 두 개의 0 사이의 절충점이며, 그 효율은 신기하게도 열원의 온도에만 의존합니다.
실제 발전소와의 비교
Curzon과 Ahlborn은 논문의 Table I에서 자신들의 공식을 실제 가동 중인 발전소와 비교했습니다. 아래 값은 표에 적힌 온도로 다시 계산한 것이며, 맨 위에는 Schroeder가 예로 든 전형적인 석탄 화력 터빈을 추가했습니다.
| 발전원 | \(T_{c}\) (°C) | \(T_{h}\) (°C) | \(\eta_{C}\) | \(\eta_{CA}\) | 실측 \(\eta\) |
|---|
실측 효율은 Curzon과 Ahlborn(1975)의 Table I에서 가져왔고, 전형적인 석탄 발전소의 40%는 Schroeder가 인용한 값입니다. 원 논문의 표와 약간 다른 것은 섭씨–켈빈 변환의 반올림 차이 때문입니다.
단순화 가정이 결과를 바꿀까?
Schroeder 교과서의 각주에 따르면, 최종 효율을 얻는 데 두 열전도 계수가 같다는 가정도, 두 등온 과정의 시간이 같다는 가정도 필요하지 않습니다. Curzon과 Ahlborn은 고온 쪽 \(\alpha\)와 저온 쪽 \(\beta\)를 다르게 두고 등온 과정의 시간 \(t_1, t_2\)도 자유롭게 두었습니다. 두 작동 온도에 대해 모두 최대화하면
\[ T_{hw}^* = C\sqrt{T_h},\quad T_{cw}^* = C\sqrt{T_c},\quad C = \frac{\sqrt{\alpha T_h} + \sqrt{\beta T_c}}{\sqrt{\alpha} + \sqrt{\beta}},\quad \frac{t_1}{t_2} = \sqrt{\frac{\beta}{\alpha}} \]이므로 다시 \(T_{cw}^*/T_{hw}^* = \sqrt{T_c/T_h}\)가 됩니다. \(\alpha = \beta\)이면 최적 시간이 저절로 같아지므로, 이 경우 Schroeder의 \(\Delta t\)가 같다는 가정은 아무런 손해도 없습니다. 아래 위젯은 열전도 계수가 다른데도 시간을 같게 강제하면 어떻게 되는지도 확인합니다. 최적점이 이동하고 일률은 줄어들지만, 최대 일률에서의 효율은 여전히 \(\eta_{CA}\)입니다.
덧붙임: 양자점 열기관
Curzon–Ahlborn의 결과는 최대 일률에서의 효율에 관한 긴 연구 흐름의 출발점이 되었습니다. Van den Broeck은 평형 근처(작은 \(\eta_C\))에서 선형 비가역 열역학을 이용해, 열 흐름이 일을 만드는 흐름에 정확히 비례하는 "강결합(tight coupling)" 엔진이라면 \(\eta_{\rm op} \approx \tfrac12\eta_C\)여야 함을 보였습니다. \(\eta_{CA} = \tfrac12\eta_C + \tfrac18\eta_C^2 + \tfrac1{16}\eta_C^3 + \cdots\)이므로 Curzon–Ahlborn 엔진도 이를 만족합니다.
이상훈(강의 담당 교수), 엄재곤, 박형규의 2018년 연구는 이 "\(\tfrac12\)-보편성"이 얼마나 견고한지를, 상상할 수 있는 가장 단순한 엔진 중 하나로 살펴보았습니다. 에너지 \(E_{\rm QD}\)(\(\mu_1\) 기준)의 단일 양자점 준위가 고온 도선(\(T_1\))과, 화학 퍼텐셜이 더 높은(\(\Delta\mu = \mu_2 - \mu_1\)) 저온 도선(\(T_2\))을 잇는 구조입니다. 전자 하나가 통과할 때마다 고온 도선에서 열 \(E_{\rm QD}\)를 가져와 일 \(\Delta\mu\)를 하므로
\[ J = \tfrac12\left[\frac{1}{e^{E_{\rm QD}/T_1}+1} - \frac{1}{e^{(E_{\rm QD}-\Delta\mu)/T_2}+1}\right],\qquad \dot W = J\,\Delta\mu,\qquad \eta = \frac{\Delta\mu}{E_{\rm QD}} \]입니다(\(k_B = 1\)). 이 엔진은 구조적으로 강결합 조건을 만족합니다. 내부가역 엔진에서 \(T_{hw}\)가 하던 역할을 이제 두 개의 손잡이가 맡습니다. 게이트 전압이 \(E_{\rm QD}\)를, 소스–드레인 전압이 \(\Delta\mu\)를 정합니다. 어떤 손잡이를 돌릴 수 있느냐가 결과를 크게 바꾼다는 것이 핵심입니다.
두 손잡이 모두 자유(전역 최적화). \(\eta_{\rm op} = \tfrac12\eta_C + \tfrac18\eta_C^2 + 0.0775\,\eta_C^3 + \cdots\)로, 2차까지는 Curzon–Ahlborn과 같고 3차에서야 CA 계수 \(1/16\)과 달라집니다.
게이트 고정, \(\Delta\mu\) 조절. \(\tfrac12\)은 유지되지만 2차 계수가 \(E_{\rm QD}\)에 의존합니다: \(\eta_{\rm op} = \tfrac12\eta_C + \tfrac{E_{\rm QD}}{16T_2}\tanh\!\big(\tfrac{E_{\rm QD}}{2T_2}\big)\eta_C^2 + \cdots\).
바이어스 고정, 게이트 조절. 여기서는 보편성이 깨집니다: \(\eta_{\rm op} \approx \eta_C - (T_2/\Delta\mu)\,\eta_C^2\)로, 선형 계수가 \(\tfrac12\)이 아니라 1입니다. 전개의 기준점 \(E_{\rm QD}\to\infty\)는 흐름이 0이지만 평형이 아닌 점이므로 선형 응답 이론이 적용되지 않습니다. 역학적 흐름은 선형이 아니라 본질적 특이성(essential singularity)을 가지며 사라집니다. 흐름이 \(J_2 \propto (X_2 + \xi' X_t)X_2^{\,n}\)처럼 행동한다면 \(\eta_{\rm op} = \frac{n+1}{n+2}\eta_C\)를 얻는데, 보통의 선형 경우 \(n=0\)은 \(\tfrac12\)을, 본질적 특이성은 \(n\to\infty\)처럼 작용하여 1을 줍니다. 게이트 전압을 조절한 Josefsson 등(2018)의 양자점 실험에서도 \(\eta_{CA}\)보다 훨씬 높은 최대 일률 효율이 보고되어 이 그림과 부합합니다.
아래 곡선은 \(\dot W\)를 수치적으로 최대화하여 실시간으로 계산한 것입니다(\(T_2 = 1\)).
정리하면, "최대 일률에서의 효율"은 엔진마다 하나로 정해진 숫자가 아닙니다. 수송 법칙(선형 열전도는 Curzon–Ahlborn을 줍니다), 어떤 변수를 최적화하는지, 그리고 전개의 기준이 되는 상태가 진짜 평형인지에 따라 달라집니다. Curzon–Ahlborn 엔진에서는 시간을 같게 강제해도 \(\eta_{CA}\)가 그대로이지만, 양자점에서는 바이어스 대신 게이트를 고르는 것만으로 선형 차수의 답이 바뀝니다.