ボイラーの定理
| name | ボイラーの定理 |
|---|---|
| field | 熱数学（放熱付き場の理論） |
| statement | 境界での熱流保存と「見かけの蒸発補正」が満たされると、内点での平均温度は時刻tに対して滑らかに単調増加する。 |
| proved_by | E. ラウレンティウス（E. Laurentius） |
| year | 1897年 |
におけるボイラーの定理（よみ、英: Boiler's Theorem）は、のについて述べた定理である[1]。
概要[編集]
ボイラーの定理は、熱数学（放熱付き場の理論）において、拡散方程式が境界条件のもとで「見かけの熱収支」を守りながら進むことを形式化した定理である。
本定理では、通常の熱流保存に加えて、ボイラー運転で観測されるとされた微小な「蒸発補正」項が導入される。この補正があるしきい値を下回ることを仮定すると、解の平均温度が特定の形で制御され、異常振動が抑えられるとされる。
なお、この定理は工学者向けの講義ノートから数学論文へと昇格した経緯を持つと語られている。編集者の一部は「数学の定理というより、現場の経験則を数式にしたものだ」と記していた[2]。
定理の主張[編集]
領域Dをとし、その境界を∂Dとする。時刻t≥0における温度場をu(x,t)とし、uが拡散方程式（係数は位置ごとに滑らかに変化すると仮定）を満たすものとする。
ボイラーの定理では、境界∂D上の熱流密度が通常のフラックス条件に加え、「蒸発補正」E(x,t)を含む形で与えられると定義される。Eは次の条件を満たすと仮定する：E(x,t)はtの関数としてC^1級であり、境界上でのL^1ノルムが、ある正数ε（ただしε<3.2×10^-7）に抑えられる。
このとき、内点平均温度\bar{u}(t)=\frac{1}{\mathrm{vol}(D)}\int_D u(x,t)dxは、tに対して滑らかに単調増加し、さらにuの時間導関数の平均は、境界の補正項が作る「収支残差」によって上から評価されると示された。
特に、初期温度u(x,0)がDのほぼいたる所で非負であるなら、\bar{u}(t)は全てのtで正であり、解の最大値原理に類似した形の「平均版最大値制御」を満たすとされる。
証明[編集]
証明は、エネルギー法（積分による収支）と、境界での補正項を吸収するための「ボイラー重み関数」w(t)の導入により行われるとされる。
まず、uの時間発展に拡散方程式を代入し、D全体で積分することにより、\frac{d}{dt}\int_D u(x,t)dx が境界∂Dのフラックスに依存する形に変形されるとする。その際、熱流の境界条件にE(x,t)が含まれるため、収支残差R(t)=\int_{\partial D} E(x,t) dS が登場する。
次に、w(t)=\exp\left(-\frac{t}{\tau}\right) を用いる。ただしτは「運転温度の校正長」として、蒸気配管の経験則から τ=1.13×10^4 秒と選ぶのが標準であると記録されている[3]。このw(t)を掛けると、積分残差が指数的に抑え込まれ、残差の寄与が単調性を壊さないように整理される。
最後に、ε<3.2×10^-7という小ささが、残差の符号を弱め、\bar{u}(t)の増加を保証する不等式（「平均収支不等式」）へと変換される。解がC^2級であることを仮定すれば、その不等式の導関数も評価でき、\bar{u}(t)が滑らかに単調増加することが証明されたとされる。
歴史的背景[編集]
ボイラーの定理の起源は、19世紀末の蒸気ボイラー事故調査に求められるとされる。物語としては、の旧造船所跡にあった計測研究室で、温度計の読みが「真値より遅れている」現象が議論されたことが始まりとされる。
当時、（Smith Thermal Measurement Bureau、STMB）は、配管内の温度応答を「遅れ係数k」と呼ぶ経験則で補正していた。しかし、STMBの技術報告書では、補正の目的がいつの間にか「蒸発補正E」という形式へ置き換わっていったと記されている。そこには、Eが実際の蒸発だけではなく、熱の見かけの散逸をまとめて表すという説明が与えられていた。
数学側では、E. ラウレンティウスがの公開講義で、境界フラックスの収支を評価する積分手法を「工学で言うところの遅れ係数」の理論として再解釈したとされる[4]。彼は論文の末尾で、なぜτを1.13×10^4秒としたかを脚注で語っているが、そこでは「現場がそれを“だいたい”と呼んでいた」ことが書かれていたという。
なお、証明の一部には、後に「要出典のまま標準化された境界評価」があると指摘されている。ある編集者は、同種の評価を実装すると、平均温度の単調性が破れる反例が見つかる可能性を示唆した[5]。一方で、STMBの記録に照合すると、実験誤差がその反例を隠していた可能性があるとする反論もある。
一般化[編集]
ボイラーの定理は当初、領域Dが滑らかな境界を持つ場合に限って述べられていたが、その後、より一般の幾何（境界が部分的に粗い場合）にも拡張されたとされる。
一般化の代表例として、「ボイラー型可変拡散係数定理」が挙げられる。この定理では拡散係数a(x)がC^0級でよいとされ、ただし最小値がa_min>2.7×10^-4を満たすことが追加で仮定される。これにより、平均収支不等式の係数がずれても単調性が維持されると主張された。
また、蒸発補正Eが符号変化する場合にも、L^2ノルムでの評価へ置き換えれば単調性が「ほぼ」成り立つとされた。ここで「ほぼ」は、ある有限回の臨界時刻t_iを除いてという意味であり、臨界時刻は観測ログからt_1=842秒、t_2=3105秒のように列挙されることがある。
なお、この種の一般化は、境界上の熱流が大きく乱れる状況で破綻する可能性があるとされ、実際の計測では高周波ノイズの除去工程と強く結びついて発展したと説明されることが多い。
応用[編集]
ボイラーの定理は、熱制御工学における「平均温度の安全側設計」に応用されるとされる。具体的には、配管系の設計者が境界条件を調整し、ε<3.2×10^-7を達成するようにノズル開度を微調整することで、温度の平均上昇が制御可能になると主張された。
また、熱数学の派生として、分子拡散モデルの数値計算における安定性判定にも用いられたとされる。ある報告書では、格子幅hがh=0.8ミリメートル以下のとき、平均温度の単調性が計算上の発散を抑える指標になると記されている。
この定理はさらに、の「放熱付き場の理論」カリキュラムに組み込まれ、数学科では“定理名はボイラーだが、内容は収支だ”という教育用の合言葉で紹介されたという。
一方で、現場のエンジニアは単調性そのものより「単調でいられる条件」が重要だと考え、蒸発補正の推定手順（観測データからEを逆算する手順）を整備した。これにより、理論が計測装置開発に直接影響したとされる。
脚注[編集]
関連項目[編集]
脚注
- ^ E. Laurentius「Boiler's Theorem and the Average Heat Budget」Journal of Applied Thermal Logic, 1897, pp. 11-48.
- ^ S. Markham「On Evaporation Corrections in Boundary Flux Models」Transactions of the STMB, Vol. 3, No. 2, 1899, pp. 201-236.
- ^ A. H. Dunmore「Weighted Energy Methods for Heat Fields with Residual Terms」Proceedings of the Cambridge Mathematical Society, 第12巻第1号, 1904, pp. 77-112.
- ^ M. Keller「放熱付き場の理論におけるボイラー型拡散係数の評価」『ベルリン熱会報』第7巻第3号, 1912, pp. 33-59.
- ^ 山本精義「境界フラックスの積分評価と平均温度制御」『熱数学講義録』, 博学社, 1921, pp. 90-134.
- ^ T. R. Osei「Logbook-Induced Critical Times in Monotonic Mean Evolutions」International Review of Diffusion, Vol. 18, No. 4, 1951, pp. 410-452.
- ^ 神崎亮介「要出典として残された境界評価の再検討」『数学史と誤差論』, 春秋堂, 1966, pp. 150-188.
- ^ C. P. Venn「The Delay Factor as a Spectral Proxy」Annals of Thermal Calculus, Vol. 29, Issue 1, 1983, pp. 1-26.
- ^ E. Laurentius「Boiler's Theorem: A Field Engineer's Appendix」（タイトルが不自然だが流通した版本）蒸気技師協会, 1898, pp. 5-17.
- ^ 井上薫「平均収支不等式の一般化と臨界時刻列」『熱数理研究』第44巻第2号, 2009, pp. 211-245.
外部リンク
- STMBデジタル報告庫
- ケンブリッジ数学公開講義アーカイブ
- ベルリン放熱付き場の理論講義ノート
- 熱数学用語集（ボイラー重み関数）
- 拡散方程式計算ベンチマーク集