佐久間＝服部方程式

佐久間＝服部方程式（さくま=はっとりほうていしき）は、完全黒体から放射された、もしくは熱放射検出器により受信された熱放射、放射束、放射パワーの量を予測するための数学モデルである。 1982年に佐久間史洋、小野晃、服部晋により最初に提案された。1996年、この方程式のさまざまな形式の有用性が研究された。この研究において、ほとんどの応用に最適なプランキアン形式が示された。この研究は3個以下のフィッティング変数を含む佐久間＝服部方程式の10の異なる形式に対して行われた。2008年、BIPM CCT-WG5は960℃未満の放射温度測定の不確かさバジェットにこの方程式を使用することを推奨した。 佐久間＝服部方程式は物体の温度に基づいて熱放射から電磁波信号を与える。信号は電磁流速であってもよいし、この放射を測定する検出器により生成される信号であってもよい。銀点以下では佐久間＝服部方程式を用いた方法が提案されている。一般形は次のようになる。 S ( T ) = C exp ⁡ ( c 2 λ x T ) − 1 , {\displaystyle S(T)={\frac {C}{\exp \left({\frac {c_{2}}{\lambda _{x}T}}\right)-1}},} ここで プランキアン形式は次の置換により実現される。 λ x = A + B T {\displaystyle \lambda _{x}=A+{\frac {B}{T}}} この置換を行うことで次のようなプランキアン形式の佐久間＝服部方程式が得られる。 プランキアン形式は放射温度測定や赤外線温度測定の不確かさバジェットの計算に使用することが推奨されている。また、銀点未満の放射温度計の較正に使用することも推奨されている。 プランキアン形式はプランクの法則と類似している。 S ( T ) = c 1 λ 5 [ exp ⁡ ( c 2 λ T ) − 1 ] {\displaystyle S(T)={\frac {c_{1}}{\lambda ^{5}\left[\exp \left({\frac {c_{2}}{\lambda T}}\right)-1\right]}}} ただし、佐久間＝服部方程式は低温、広帯域の放射温度測定を考える場合に非常に有用である。広いスペクトル帯域でプランクの法則を使用する場合、次のような積分を考慮する必要がある。 S ( T ) = ∫ λ 1 λ 2 c 1 λ 5 [ exp ⁡ ( c 2 λ T ) − 1 ] d λ {\displaystyle S(T)=\int _{\lambda _{1}}^{\lambda _{2}}{\frac {c_{1}}{\lambda ^{5}\left[\exp \left({\frac {c_{2}}{\lambda T}}\right)-1\right]}}d\lambda } この積分により不完全多重対数関数（英語版）が生成されるが、この関数によりこの方程式は扱いにくくなる。標準的な数値処理では指数の幾何級数の不完全積分を展開する。 ∫ 0 λ 2 c 1 λ 5 [ exp ⁡ ( c 2 λ T ) − 1 ] d λ = c 1 ( T c 2 ) 4 ∫ c 2 / ( λ 2 T ) ∞ x 3 exp ⁡ ( x ) − 1 d x {\displaystyle \int _{0}^{\lambda _{2}}{\frac {c_{1}}{\lambda ^{5}[\exp({\frac {c_{2}}{\lambda T}})-1]}}d\lambda =c_{1}({\frac {T}{c_{2}}})^{4}\int _{c_{2}/(\lambda _{2}T)}^{\infty }{\frac {x^{3}}{\exp(x)-1}}dx} （ λ = c 2 / ( x T ) {\displaystyle \lambda =c_{2}/(xT)} , d λ = − c 2 / ( x 2 T ) d x {\displaystyle d\lambda =-c_{2}/(x^{2}T)dx} の置換をしている。）すると、 J ( c ) ≡ ∫ c ∞ x 3 exp ⁡ x − 1 d x = ∫ c ∞ x 3 exp ⁡ ( − x ) 1 − exp ⁡ ( − x ) d x = ∫ c ∞ ∑ n ≥ 1 x 3 exp ⁡ ( − n x ) d x {\displaystyle J(c)\equiv \int _{c}^{\infty }{\frac {x^{3}}{\exp x-1}}dx=\int _{c}^{\infty }{\frac {x^{3}\exp(-x)}{1-\exp(-x)}}dx=\int _{c}^{\infty }\sum _{n\geq 1}x^{3}\exp(-nx)dx} = ∑ n ≥ 1 exp ⁡ ( − n c ) ( n c ) 3 + 3 ( n c ) 2 + 6 n c + 6 n 4 {\displaystyle =\sum _{n\geq 1}\exp(-nc){\frac {(nc)^{3}+3(nc)^{2}+6nc+6}{n^{4}}}} これにより、和をある桁で切り捨てることで概算値が得られる。 上で示した佐久間＝服部方程式は多くの方程式を検討した中で、放射温度計のスケールの補完に最適なフィット曲線を提供することがわかった。 反復計算をすることなく逆佐久間＝服部関数を使用することができる。これはプランクの法則の積分よりも優れている。 1996年の論文では10の異なる形式が研究された。それらを実際の放射測定データに対するフィット曲線の質の高い順に以下の表に示す。 シュテファン＝ボルツマンの法則 プランクの法則 レイリー・ジーンズの法則 ヴィーンの放射法則 ヴィーンの変位則 キルヒホッフの法則 (放射エネルギー) 放射温度計 パイロメーター 細フィラメント高温測定法 サーモグラフィー 黒体 熱放射 放射輝度 放射率 ASTM Subcommittee E20.02 on Radiation Thermometry ^ 銀点とは銀の融点である962℃(961.961 ± 0.017)℃であり、いくつかの温度スケールで較正点として使われている。 安定しており再現性が高いため、赤外温度計の較正に使用されている。 ^ a b F Sakuma, S Hattori, "Establishing a practical temperature standard by using a narrow-band radiation thermometer with a silicon detector", in Temperature: Its Measurement and Control in Science and Industry, vol. 5, edited by J F Schooley, New York, AIP, 421–427 (1982). ^ a b c Sakuma F, Kobayashi M., "Interpolation equations of scales of radiation thermometers", Proceedings of TEMPMEKO 1996, pp. 305–310 (1996). ^ a b c d J. Fischer, P. Saunders, M. Sadli, M. Battuello, C. W. Park, Y. Zundong, H. Yoon, W. Li, E. van der Ham, F. Sakuma, Y. Yamada, M. Ballico, G. Machin, N. Fox, J. Hollandt, M. Matveyev, P. Bloembergen and S. Ugur, "Uncertainty budgets for calibration of radiation thermometers below the silver point" (pdf), CCT-WG5 on Radiation Thermometry, BIPM, Sèvres, France (2008). ^ “2006 CODATA recommended values”. National Institute of Standards and Technology (NIST) (Dec 2003). Apr 27, 2010閲覧。 ^ a b MSL Technical Guide 22 – Calibration of Low Temperature Infrared Thermometers (pdf), Measurement Standards Laboratory of New Zealand (2008). ^ ASTM Standard E2758-10 – Standard Guide for Selection and Use of Wideband, Low Temperature Infrared Thermometers, ASTM International, West Conshohocken, PA, (2010). ^ J Tapping and V N Ojha (1989). “Measurement of the Silver Point with a Simple, High-Precision Pyrometer”. Metrologia 26 (2): 133–139. Bibcode: 1989Metro..26..133T. doi:10.1088/0026-1394/26/2/008.  ^ “Definition of Silver Point - 962°C, the melting point of silver”. 2010年7月26日閲覧。