エクセルギー
エクセルギー (exergy) とは、『系が外界とのみ熱・仕事を交換しながら、外界と平衡するまで状態変化するとき、系から理論上取り出せる最大の仕事量』のことである。availability、available energy、有効エネルギー などとよばれることもある。
熱力学第二法則によると、熱を仕事に変換するには常にカルノー効率による制約が伴うので、熱に関連したエネルギーを扱うときには、そのうちの有効に仕事に変換できる部分とできない部分とを区別して扱うことが必要である。前者がエクセルギーとよばれる。
熱力学第一法則(エネルギー保存則)は常に成り立つので、外界を含めた拡大系では、エネルギーの量は一定不変である。「省エネルギー」、「エネルギーを節約する」と言うときの「エネルギー」は、「エネルギーの価値」、「エネルギーの質」の意味であり、これがエクセルギーに対応していると考えられる。省エネルギーを念頭において機器の開発や改良を行う際には、エクセルギーを用いた評価法は極めて有力な手段となる。
「熱から最大いくらの仕事が取り出せるか」との問題に最初に答えようとしたのは、サディ・カルノー((N. L. Sadi Carnot))であり(1824)、その四半世紀後に、トムソン(W. Thomson, Baron Kelvin)、クラウジウス(R. J. E. Clausius)らが熱力学を確立する基礎となった。このカルノーの問題意識はエクセルギーに対する問いそのものであり、熱力学第二法則がその答えとなっている。
availability という用語は、1868年に初めてテイト(P. G. Tait) の書籍で「熱の仕事として利用できる度合い」の意味で用いられた。マクスウェル(J. C. Maxwell)やトムソンはこれを available energy とよんでいた。1873年にギブズ(J. W. Gibbs)が available energy の解析的表示を提示し、自由エネルギーの概念を発表した。1889年にグイ(Louis Georges Gouy)が ギブズと同様な表示式を提示し、1898年にストドラ(Aurel Stodola)がエントロピー発生とエクセルギー損失の関係(グイ・ストドラの定理)を発見した。アメリカでは、1930年代にキーナン(J. H. Keenan)が availability と命名して、熱力学的概念の確立と普及に大きく寄与した。
エクセルギー概念が普及するにつれ、日用語と区別できる名称の必要性が認識され、1956 年にユーゴスラビア(現スロヴェニア)の工学者 ゾラン・ラント(Zoran Rant) がエネルギー(en + ergon; 内への仕事)にならって、ラテン語由来の接頭辞 ex (外へ) とギリシア語 ergon (仕事) の合成語として、エクセルギー(exergy)の術語を提案した。現在ではほぼエクセルギーが一般的に使用されている。
熱は物体から他の物体へ移動するエネルギーのひとつの形態であり、物体の温度に密接に関連したエネルギーである。温度 T の物体からとりだした Q の熱から得られる仕事の最大値は、外界(温度 T0)との間でカルノーサイクルを動かしたとして、次式となる。
W = \left( 1 - \frac{T_0}{T} \right) Q  \qquad \qquad \qquad \qquad (1)
外界よりも温度の低い物体は「冷熱」を持っているということもある。この場合、外界を高温熱源、物体を低温熱源とするカルノーサイクルを用いれば、最大
W = \left( 1 - \frac{T}{T_0} \right) Q_0
の仕事が得られる。ただし、Q0 はカルノー機関が温度 T0 の外界から取り出す熱であり、温度 T の物体に吸収される熱量 -Qとの間に
\frac{Q_0}{-Q} = \frac{T_0}{T}
の関係がある。これを用いると、
W = \left( 1 - \frac{T}{T_0} \right) \frac{T_0}{T} (-Q) = \left( 1 - \frac{T_0}{T} \right) Q
となり、前記の Q > 0 の場合の式 (1) と同じ形となる。冷熱の場合は、Q の値は負であり、カルノー効率 (1 - T0/T) も負であり、結果として得られる仕事量は正の値となる。
結局、冷熱の場合も含めて、温度 T の物体が放出する熱流 Q (物体から出る方向を正とする)のエクセルギーは
E = \left( 1 - \frac{T_0}{T} \right) Q       \qquad \qquad \qquad \qquad (2)
で表される。
ある機器から放熱がある場合、放熱の間の物体温度は未知の場合も多いが、物体の温度が一様に低下する(熱容量が一定である)のであれば、放熱前後の温度 T1 と T2 の対数平均値
T_m = \frac{T_1 - T_2}{\ln (T_1/T_2)}
を用いて、
E = \left( 1 - \frac{T_0}{T_m} \right) Q     \qquad \qquad \qquad \qquad (3)
とすることができる。
ある状態の物質の持つエネルギーは物質の出入のない系(閉じた系)では内部エネルギー、そうでない系(開いた系)ではエンタルピーであるが、この中には、仕事に変換できる部分とできない部分が含まれている。物質が外界と平衡するまで状態変化を行うとき、そこから取り出せる仕事の最大値が物質の持つエクセルギーである。
対象とする物質(系)が任意の状態変化をするとき、系は外界に対して仕事をし、同時に熱を放出する。外界に対して行った仕事は上記の仕事の一部となるのに加えて、外部へ放出した熱も、その後仕事に変換できる可能性がある。
そこで、右図のように系と並行して動作し、系が放出する熱を受け取って仕事を取り出す補助熱機関を考える。系自身が行う仕事 Wi を内部仕事とよび、補助熱機関が行う仕事を外部仕事とよぶ。系から取り出せる仕事は、内部仕事と外部仕事を合わせた有効仕事 Wg = Wi + We であると考えることができる。ある状態の物質の持つエクセルギーは、外界と平衡するまで状態変化を行うときに取り出せる有効仕事の最大値である。
このような有効仕事に対して、次の定理が成立する。
『与えられた二つの状態間で変化が行われる場合、変化が可逆の場合に有効仕事が最大(符号も含めた代数的意味で最大)となる。可逆であれば、有効仕事は変化の経路に依存せず、両端の状態のみにより定まる。』
(証明)
状態1から2へ至る二つの状態変化 R と I を考える。状態変化 R は外部に対して有効仕事 WgR を行い、外界へ -QR の熱を放出し,状態変化 I は外部に対して有効仕事 WgI を行い、外界へ -QI の熱を放出するとする。
状態変化 R が可逆変化であるとき、R の向きを逆にして、1から I に沿って2へ行き、R に沿って 1 へ戻るサイクルを考えると、有効仕事と熱の出入りは 右図 のようになる。
サイクル後は系および補助熱機関は元の状態に戻るので、熱力学第一法則(エネルギー保存則)より
W_{gI} + (-Q_I) = W_{gR} + (-Q_R)
(-Q_R) - (-Q_I) = W_{gI} - W_{gR}
が成立する。
もし、WgI > WgR であれば、外界から (-QR) - (-QI) の熱を受け取り、それをすべて仕事 WgI - WgR に変換する第二種永久機関となる。したがって、熱力学第二法則に矛盾しないためにはWgI ≤ WgR でなければならない。つまり、可逆変化 R の有効仕事は任意変化 I の有効仕事より必ず大きくなる(等号は下記のように、双方共に可逆変化の場合に対応)。
I も可逆変化である場合には、上の R と I を入れ替えた結果も成立するので、WgI = WgR でなければならず、可逆変化の有効仕事は、経路のいかんにかかわらず等しくなる。(証明終わり)
以上より、ある状態の物質のもつエクセルギーは、現在の状態から外界と平衡するまで任意の可逆経路に沿って変化したときに行う有効仕事として求まる。
閉じた系の微小変化では熱力学第一法則より次式が成立する。
dQ = dU + dW_i
ここで、系は外界で囲まれているとすると、系が行う内部仕事 dWi には、外界の圧力 p0 に抗して行う排除仕事 p0 dV も含まれているが、これは取り出すことはできないので、後で取り除かなければならない。
この間、系が放出する熱は -dQ であり、これより補助熱機関がdWe = (1 - T0/T) (-dQ) の外部仕事を行うので、この間の有効仕事は次式となる。
dW_g = dW_i + dW_e = dQ - dU + \left(1 - \frac{T_0}{T}\right) (-dQ)
= - dU + \frac{T_0}{T} dQ
ここで、dQ/T = dS は系のエントロピー増加量 dS であるので、
dW_g = - dU + T_0 dS
これを現在の状態 (U, S) から外界と平衡する状態 (U0, S0) まで積分して
W_g = - (U_0 - U) + T_0 (S_0 - S) = (U - U_0) - T_0 (S - S_0)
閉じた系のエクセルギーは、この有効仕事から排除仕事 p0 (V0 - V) を除いた次式となる。
E = W_g - p_0 (V_0 - V) = (U - U_0) - T_0 (S - S_0) - p_0 (V - V_0) \qquad (4)
ヘルムホルツの自由エネルギー F およびギブズの自由エネルギー G は、次式で定義される状態量である。
\begin{align}
F &= U - T S \\
G &= H - T S
\end{align}
\begin{align}
dF &= dU - T dS \\
dG &= dH - T dS
\end{align}
と表される。
一方、等温変化での閉じた系のエクセルギー変化量は、系と外界が熱平衡である(T = T0)として、
dE = dU - T dS + p_0 dV
となる。
系の体積が一定に保たれているときは、 dV = 0 であるので、
dE = dU - T dS   = dF
となる。また、系の圧力が一定に保たれているときは、 p0 = p と置き換えて、U + pV = H を用いれば、
dE = dU - T dS + p dV = dH - T dS = dG
となる。つまり、ヘルムホルツの自由エネルギーの変化量は、等温等積条件におけるエクセルギーの変化量であり、ギブズの自由エネルギーの変化量は、等温等圧条件におけるエクセルギーの変化量である。
自由エネルギーは、温度一定の条件で生じる変化の方向を示すのに用いられる。体積一定の系ではヘルムホルツの自由エネルギーが減少する方向に変化が生じ、圧力一定の系ではギブズの自由エネルギーが減少する方向に変化が生じる。これは、一般の変化では、エクセルギーが減少する方向の変化だけが生じることと合致している。
温度一定に保たれた系で、膨張を伴わない仕事 dWe(例えば、電池等の電気的な仕事、液柱を持ち上げる等の仕事など)を生じるとすると、その非膨張仕事の最大値は、体積一定の系では dWe ≤ -dF となり、圧力一定の系では dWe ≤ -dG となる。一般の変化で系から取り出せる 仕事の最大値 は、系のエクセルギーの減少量 dW ≤ -dE であることと合致している。
化学反応を伴う多くの系では、反応に伴うエネルギー変化(生成エンタルピーの差)は圧力、温度などの物理的変化に伴うエネルギー変化に比べて数桁大きく、後者を無視しても問題ない場合が多い。後者が無視できない場合は、自由エネルギーではなくエクセルギーを用いる必要がある。
つまり、等温等積または等温等圧の条件での自由エネルギーの概念を、一般の変化に拡張したものがエクセルギーであると考えることができる。ただし、自由エネルギーは物質の状態だけで決まる状態量であるのに対して、エクセルギーは、外界条件を固定しなければ値が確定しない点に注意を要する。
実用的な多くの系は、物質の出入を伴う開いた系となる。開いた系の微小変化では、熱力学第一法則 dQ = dH + dWi より、内部仕事は、
dW_i = dQ - dH
と表される。また、補助熱機関が行う外部仕事は、閉じた系と同一であるので、有効仕事は次式となる。
\begin{align}
dW_g &= dW_i + dW_e = dQ - dH + \left(1 - \frac{T_0}{T}\right) (-dQ) \\
&= - dH + \frac{T_0}{T} dQ = -dH + T_0 dS
\end{align}
ただし、dQ/T = dS は可逆変化における系のエントロピー増加量である。
これを現在の状態 (H, S) から外界と平衡する状態 (H0, S0) まで積分すれば、開いた系で物質がもつエクセルギーは次式となる。
E = W_g = H - H_0 - T_0 (S - S_0)                       \qquad \qquad (5)
物質がエネルギー変換に関与する実用的なプロセスのほとんどは流動系であるため、以下では開いた系に限って説明する。また、単位量(質量、標準体積またはモル)あたりのエクセルギー(比エクセルギー)を対象とする。
e = E/G = h - h_0 - T_0 (s - s_0)
ただし、 h は比エンタルピー、s は比エントロピーである。
多くの固体や液体では圧力の影響は無視することができ、狭い温度範囲であれば比熱は一定と見なすことができる。この場合は、dh = c dT,   ds = (c/T) dT = c d(ln T) より、次式より比エクセルギーを求めることができる。
e = c \left[ T - T_0 - T_0 \ln \left(\frac{T}{T_0}\right) \right]
= (h - h_0) \left[ 1 - \frac{T_0}{T - T_0} \ln \left(\frac{T}{T_0}\right) \right]
\qquad (6)
エネルギー(エンタルピー)のうちエクセルギーの占める比率を有効度という。
上記の比熱一定の物質の有効度 λ は
\lambda = \frac{e}{h-h_0}
= 1 - \frac{T_0}{T - T_0} \ln \left(\frac{T}{T_0}\right)
となる。
多くの場合、気体を理想気体(完全ガス)で近似することができる。このとき、比熱 cp, cv は一般に温度の関数になるが、狭い温度範囲であれば、これも一定と見なせる。
圧力が変わらなければ、式(6)に比熱として定圧比熱 cp を用いればよい。
温度に加えて圧力も変化する場合は、h = cp T, s = cp ln T - (cp - cv ) ln p = cp { ln T - [(κ-1)/κ] ln p } となるので、次式から比エクセルギーを求めることができる。
\begin{align}
e &= c_p (T - T_0)
- c_p T_0
\left[ \ln \left(\frac{T}{T_0}\right)
- \frac{\kappa - 1}{\kappa} \ln \left(\frac{p}{p_0}\right) \right] \\
&= (h - h_0)\left\{
1 - \frac{T_0}{T - T_0}
\left[ \ln \left(\frac{T}{T_0}\right)
- \frac{\kappa - 1}{\kappa} \ln \left(\frac{p}{p_0}\right) \right]
\right\}                                                       \qquad (7)
\end{align}
\lambda =
1 - \frac{T_0}{T - T_0}
\left[ \ln \left(\frac{T}{T_0}\right)
- \frac{\kappa - 1}{\kappa} \ln \left(\frac{p}{p_0}\right) \right]
となる。ただし、κ = cp/cvは定圧比熱と定積比熱の比であり、2原子分子気体では約 1.4 となる。
成分 1 , 2 をそれぞれ n1、n2モル含む混合気体(系)を取り上げる。同成分の外界でのモル比を n10 : n20 とする。ただし、n1+n2 = n10+n20 = n とし、系の温度、圧力は外界の温度 T0 、圧力 p0に等しいものとする。
このとき系のエクセルギーは、特定の成分気体だけが通過できる半透膜を用いた思考実験により、次式のように求まる。
E = R T_0 \left(n_1 \ln\frac{n_1}{n_{10} + n_2 \ln\frac{n_2}{n_{20} \right)
\qquad \qquad (8)
絶対湿度 x (kg/kgDA) の湿り空気は、乾き空気(水分を除いた空気) 1 kg の中に x kg の水蒸気を含んでいる。水蒸気(分子量 18.015)と乾き空気(分子量 28.966)のモル比は
n_1 : n_2 = n_1 : (1 - n_1) = \frac{x}{18.015} : \frac{1}{28.966} = 1.608 x : 1
であるので、絶対湿度 x の湿り空気 1kmol 中のそれぞれの kmol 数は
n_1 = \frac{1.608 x}{1 + 1.608 x} , \qquad \qquad n_2 = \frac{1}{1 + 1.608 x}
となる。
一方、外界の絶対湿度を x0 とすると、その 1kmol中の水と乾き空気の kmol 数は
n_{10} = \frac{1.608 x_0}{1 + 1.608 x_0} , \qquad \qquad n_{20} = \frac{1}{1 + 1.608 x_0}
これらを式 (8) に用いて、含まれる乾き空気の質量 28.966 n2 で割ると、温度 T0、圧力 p0 で乾き空気 1kg を含む湿り空気のエクセルギー(比エクセルギー) e は次式となる。
e = R_a T_0 \left[1.608 x \ln \frac{x}{x_0}
+ \left(1 + 1.608 x \right) \ln \frac{1 + 1.608 x_0}{1 + 1.608 x} \right]
\qquad \qquad (9)
ただし、Ra = 0.2870 kJ/(kg K) は,乾き空気 1 ㎏ のガス定数である。
任意の温度、圧力の湿り空気の比エクセルギーは、式(7)と(9)の和となる。
燃料についても、化学的エネルギーを含めたエンタルピー、エントロピーを知ることができれば、式 (5) 等を用いてエクセルギーを求めることができる。しかし、石油、石炭等の化石燃料では化学的組成を含めて正確に知ることは難しく、実用上はいくつかの近似式が用いられている。その例を下表に示す。
外界を含めた全体系のエントロピーの増加量を ΔS* とすると、エクセルギー損失 LW (Loss of Work の略)は次式で表される。
LW = T_0 \Delta S^*
これは グイ=ストドラ(Gouy=Stodola) の定理と呼ばれている。
(証明)
補助熱機関を加えた大きい系が、図 2 のように状態 1 (H1, S1) から状態 2 (H2, S2) へ非可逆変化し、その間に WgI = WiI + WeI の有効仕事をし、外界へ -QI の熱を放出したとする。この間の全体系のエントロピー増加量は次式となる。
\Delta S^* = S_2 - S_1 + \frac{-Q_I}{T_0}
もし状態 1 から状態 2 へ可逆変化をする場合は、有効仕事 WgR = WiR + WeR をし、外界へ -QR の放出をするとすると、この場合の全体系のエントロピー増加量は 0 であるので、次式が成立する。
\Delta S^*_R = S_2 - S_1 + \frac{-Q_R}{T_0} = 0
したがって、
\Delta S^* = - \frac{-Q_R}{T_0} + \frac{-Q_I}{T_0} = \frac{(-Q_I) - (-Q_R)}{T_0}
と表すことができる。
一方、熱力学第一法則(エネルギー保存則)より
W_{gI} + (-Q_I) = W_{gR} + (-Q_R)
(-Q_I) - (-Q_R) = W_{gR} - W_{gI}
が成立する。したがって、
\Delta S^* = \frac{W_{gR} - W_{gI}{T_0}
となり、これより有効仕事の減少量つまりエクセルギー損失は
LW = W_{gR} - W_{gI} = T_0 \Delta S^*
となる(証明終わり)。
物質の出入りを伴う開いた系では、熱力学第一法則(エネルギー保存則)により、出入りするエンタルピー、熱、仕事(工業仕事)の間に保存則が成り立つ。したがって、ある機器に流入するエネルギーと流出するエネルギーをリストアップし、両者の差を求めれば、それは、放熱や漏洩などにより(機器外へ逃げて)失われたエネルギーとなる。このような作業は通常、熱管理、熱勘定、熱収支、等とよばれている。ある機器のエネルギー(エンタルピー)の流れは右図のようになる。
一方、エクセルギーは、機器内で非可逆変化が生じれば、グイ=ストドラの定理にしたがって必ず減少する。出入りする物質、熱、仕事のエクセルギー差を求めれば、それがこの機器におけるエクセルギーの損失 LW である。
ある機器のエクセルギーの流れは右図のようになる。エクセルギー評価は熱力学第二法則に基づく評価ということもできる。
各機器のエクセルギー評価を行うには、以下のような手順に沿えばよい。
外界と平衡する状態を選定する。これには、次節を参考にできる。
系(機器)に流入する物質、熱、仕事を列挙し、流入エネルギー(エンタルピー)を求める。物質のエンタルピーの原点は通常は便宜的(例えば三重点の水等)に定められているので、外界と平衡する状態のエンタルピーとの差を用いるのがよい。
系から流出するエネルギー(エンタルピー)についても同様に行う。
両者の差を求めて、未知のエネルギー損失(放熱、漏洩等)を推定する。
推定した放熱、漏洩等も含めて、流入・流出エネルギーをエクセルギーに換算する。
流入エクセルギーと流出エクセルギーの差を求めて、機器内のエクセルギー損失を推定する。
エンタルピーを用いたエネルギー収支では、その機器内でどれだけの(エクセルギーの)損失があったかは分からない。この種の損失は、熱と動力の変換を伴う設備(熱機関、冷凍機、各種化学プラント等)で、設備を構成する全機器のエネルギー収支を集計した結果、外界へ放出する熱量の増加となって現れる。しかし、どの機器の何を改善すれば損失が減るかは、勘と経験をもとにして試行錯誤を行わなければ知ることができない。
これに対して、エクセルギーを用いた評価では、それぞれの機器ごとに原因別に損失の内訳が明らかになるので、これを用いて見通しよく改善策を立案することができる。ここにエクセルギー評価の最大のメリットがある。
エクセルギーを求めるには基準状態となる「外界と平衡する状態」を確定する必要がある。熱力学で平衡という場合、通常、(1)熱的平衡、(2)力学的平衡、(3)化学的平衡の 3 者を満たすことを意味している。しかし、このうち、関与するプロセスによってはいくつかの条件を満たさなくてもよい場合があり、また、実用上考慮しなくてよい場合もある。以下に実用上妥当な外界条件の決め方を例示する。
(A) [熱的平衡] 温度が外界温度に一致する。 復水式蒸気タービン設備や空調機の室外機などのように、系と外界とが伝熱管等を介して間接的に(熱的にのみ)接触する場合、圧力は外界と同じになる必要はなく、また化学的な平衡も必要ない。
(B) [熱的平衡+力学的平衡] 温度および圧力が外界に一致する。内燃機関、ガスタービン、ボイラなどの燃焼排ガスは外気と直接接触混合するので、温度に加えて圧力も外界と一致する必要がある。
(C) [熱的平衡+力学的平衡+化学的平衡] 温度、圧力が外界に一致することに加えて、拡散・反応によって外界を構成する物質と化学的にも平衡する。前記の燃焼排ガスが外界(大気)と平衡するのは、厳密には排ガス成分の拡散や何らかの化学反応を経た後となるが、現実には、拡散過程やこれらの化学反応から有用な仕事を得る手段がないため、化学的平衡を無視しても問題ない。実用上、このタイプの平衡を考えなければならない場合は多くないと思われる。
右図のようなボイラとタービンを用いた自家発電(熱電併給)設備を例に、エクセルギー評価の方法について解説する(書籍を参考にしているが、大幅に変更簡略化している)。
ボイラで燃料を燃やして過熱蒸気を発生する。燃焼プロセスと熱交換プロセスを検討するために、ボイラを燃焼器と熱交換器に分けて扱う。ボイラで発生した過熱蒸気は長い管路を通って、蒸気タービンへ導かれる。タービンでは動力を取り出して発電すると共に、タービン途中から一部の蒸気を抽気し、またタービン背気(約 5 気圧)を用いて、工場内の種々の加熱用熱源として利用する(熱電併給、コジェネレーション)。
また、比較のために、外界温度の復水器を用いて全て動力として取り出す場合(発電のみ)も並行して計算し、検討の材料とする。
このような自家発電設備は、実際はかなり複雑な構成となるが、ここではエクセルギー評価の概要と利点を紹介することを目的としているので,空気予熱器や脱気器などの副次的な機器はないものとし、できる限り簡略化した構成とする。対象とした設備の各部における蒸気等の条件を下表に示す。タービン出口の条件は、タービン内部効率が 80 % となるように定めた。
また、用いた外界条件を下表に示す。
右図に示すように、C 重油を燃料とする燃焼器(ボイラ)で、燃料の化学的エネルギーを熱エネルギーに変換するプロセスを取り上げる。燃料の C 重油の高位発熱量は Hh = 44200 kJ/kg, 低位発熱量は Hl = 41400 kJ/kg であり、燃料 1 kg/s あたり、20℃ の燃焼用空気(比熱 1.30 kJ/(N m3 K) ) を12.2 N m3/s 供給して燃焼させ、2560 ℃ の燃焼ガス 12.2 N m3/s が発生したとする。燃焼ガスの比熱も空気と同じく 1.30 kJ/(N m3 K) として扱う。
流入エンタルピーは、燃料の化学的エネルギー(低位発熱量) 1×41400=41400 kJ/s と、空気のエンタルピー 12.2×1.30×(293.15 - 293.15) = 0 kJ/s を合わせて、41400 kJ/s である。
流出エンタルピーは、燃焼ガスのエンタルピー 12.2×1.30×(2560-20)=40284 kJ/s であり、差 41400 - 40284 = 1116 kJ/s が、放熱によるエネルギー損失となる(放熱自身は測定困難であるので、エンタルピー収支から推定することになる)。
エンタルピー収支の結果を下表のエンタルピー欄に示す。
ラントの近似式(表 1 )を用いて、燃料のエクセルギーは e = 0.975 × 44290 = 43095 kJ/kg となる。流入エクセルギーは、燃料の 1 × 43095 = 43095 kJ/s だけである。
燃焼ガスが持ち出すエクセルギーは、式 (7) より
E = 12.2 \times 1.30 \times
\left[ 2560 - 20 - (273.15+20) \times \ln\left(\frac{273.15+2560}{273.15+20}\right) \right]