アッツイワの方程式
| name | アッツイワの方程式 |
|---|---|
| field | 架空の応用解析学（熱ゆらぎ流体理論） |
| statement | 準自己同型変換の下で、熱ゆらぎのスペクトル密度が特定の整合条件を満たす |
| proved_by | アッツイワ研究班（日本・計算熱力学局） |
| year | 1937年 |
におけるアッツイワの方程式（よみ、英: Atsuiwa Equation）は、のについて述べた方程式である[1]。この方程式は、特異摂動の扱いにおいて実務的な近似を与える定理として、1930年代以降の計算工学に影響を及ぼしたとされる[2]。
概要[編集]
アッツイワの方程式は、熱ゆらぎ流体と呼ばれる架空のモデルにおいて、観測可能量の安定性を保証するための「方程式型の定理」として記述される[1]。
この方程式では、流体の状態を表す場がで歪んだ場合でも、スペクトル密度がある種の「整合条件」を満たすことが成り立つとされる。なお、整合条件の形は、最初期には経験則として与えられたが、のちに厳密化された経緯があるとされる[3]。
編集上の都合として、論文ではしばしば「方程式」と呼ばれるが、実体はを媒介にする定理であるため、本項では定理として扱う。
定理の主張[編集]
熱ゆらぎ流体の状態を表す確率場をX(t,x)とし、時間 t と空間座標 x に依存すると定義する[1]。
A は、X のスペクトル密度 ρ(ω) を、周波数 ω に対する重み付き積分変換として写すものとする。このとき、準自己同型変換と呼ばれる変数変換ϕの下で、次が成り立つとする：
1) ρ(ω) はρ(ω) ~ ω^{-α} を満たすものと仮定する[2]。 2) さらに位相ズレ量Δ(ϕ)が、境界位相の整合性を表す関数Θ(ϕ)と関連づけられ、
∫_0^Ω ω^{β} ρ(ω) dω = κ·Ω^{β-α+1} · exp(-Θ(ϕ))
を満たすと仮定すると、A(X) に対する整合条件が保証される[4]。
ここで Ω は観測カットオフ、α,β は指数パラメータ、κ は系のスケール因子であり、Θ(ϕ) は境界に依存する補正項であるとされる[3]。特に、Θ(ϕ)が0に近い極限では、積分則が指数的に安定化することが示される[1]。
証明[編集]
証明はとを組み合わせることによって示されたとされる[1]。
まず、X(t,x) を周波数 ω ごとに分解し、各モードが準自己同型変換ϕの下で位相だけがずれると仮定する。次に、積分の上限 Ω を「現場換算でΩ=1.37×10^4（単位系はメートル法ではなく“千熱素”）とする」ように恣意的に固定する流儀があったとされる[5]。
そのうえで、測度の保存則によりヤコビアンが補正因子としてΘ(ϕ)に吸収されるため、積分が κ·Ω^{β-α+1} · exp(-Θ(ϕ)) の形に整理されると示された[4]。なお、当初の草稿では Ω の指数が「β-α」になっていたが、東京の計算誤差監査官が約定数10.0%のズレを指摘したことで、最終形がβ-α+1に修正されたと記録されている[6]。
この種の証明は、証明書に相当する「整合スケッチ」の付録として当時の学会資料に掲載されたとされる[1]。ただし当該資料には、要出典に相当する注意書きが残っているとされ、当該箇所のみ別の班が書き足した可能性が指摘されている[3]。
歴史的背景[編集]
アッツイワの方程式が生まれた背景には、1930年代に（通称：熱計局）が推進した「境界位相を持つ装置」の検証計画があったとされる[2]。
熱計局は横浜市港湾地区での輸送中に観測される熱ゆらぎのランダム性を問題視しており、原因が装置の内壁コーティングではなく「境界位相の不整合」にある可能性を巡って議論が行われたとされる[7]。
この議論を取りまとめたのがアッツイワ研究班であり、代表として渡辺 精一郎（当時、熱計局・測度係）と（当時、外部委託の解析顧問）が関与したとされる[8]。特にThorntonは、整合条件を「物理の言葉でなく写像の言葉で書け」と繰り返したと伝えられる[4]。
なお、方程式の年次は複数資料で一致していないが、熱計局年報では1937年に「箱型計算機HL-9での再現」を根拠として記載されている[6]。
一般化[編集]
一般化は、Θ(ϕ) を境界一箇所の補正に留めず、境界の分解数に応じて重ね合わせる方向へ進んだとされる[3]。
と呼ばれる拡張では、境界をN層に分割した場合、補正項がΘ_N(ϕ)=Σ_{k=1..N} Θ_k(ϕ) の形で現れると仮定する。Nを「現場で都合がよい素数だけに限定する」慣行があり、最初に試されたのがN=11であったという逸話が残っている[9]。
また、指数則のパラメータαについても、単一の指数ではなく、周波数帯域ごとにα_iが変化する型が提案された[2]。この場合、積分の両辺が分段で一致することが示されたとされるが、完全に一般的かどうかは未解決扱いとされたと記載されている[10]。
さらに、準自己同型変換ϕを「連続変換」に限定せず、離散ステップを許すと、exp(-Θ(ϕ)) の指数が符号反転する事例が報告された。これは一部で「暴走版」と呼ばれ、あえて利用されることもあったとされる[5]。
応用[編集]
アッツイワの方程式は、熱ゆらぎ流体の安定性評価に適用されるだけでなく、計算機の誤差モデルにも導入されたとされる[1]。
具体例として、熱計局の箱型計算機では、誤差が周波数領域でγ(ω)=exp(-Θ(ϕ)) の形で減衰すると仮定した補正が採用された[6]。その結果、当時の報告書では「再現率92.6%（観測点数=184、失敗=13）」のように、やけに細かい数字で成功が述べられている[7]。
また、船舶輸送では、名古屋市港湾倉庫の温度揺らぎを事前に推定する際、観測カットオフΩを「1.37×10^4千熱素」と固定して運用したとされる[5]。この運用が、異常気象の夜間輸送ではかえって逆効果になったという証言もあるとされ、現場の数学適用の難しさが皮肉にも記録されている[9]。
一方で教育面では、工学部の解析講義で「境界位相が合うと指数的に安定する」ことの象徴例として採用され、学生が手計算でΘ(ϕ)を当てる小テストが流行したとされる[3]。
脚注[編集]
関連項目[編集]
脚注
- ^ 渡辺精一郎「アッツイワの方程式と箱型計算機HL-9の誤差補正」『計算熱力学年報』第12巻第4号, 1937年, pp. 201-248.
- ^ M. A. Thornton「Boundary Phase Consistency in Fluctuating Heat Fluids」『Journal of Applied Spectral Measures』Vol. 5, No. 2, 1939年, pp. 51-88.
- ^ 熱計局数学班「準自己同型変換における指数安定化」『日本計測論文集』第3巻第1号, 1941年, pp. 9-33.
- ^ Atsuiwa Research Group「A Note on the Matching Integral up to Ω」『Proceedings of the International Congress on Pseudoholomorphic Theory』Vol. 2, 1940年, pp. 77-94.
- ^ 佐藤和央「Ω=1.37×10^4千熱素という現場慣行の再検討」『工学記録通信』第9巻第7号, 1952年, pp. 330-346.
- ^ 【要出典】名古屋倉庫温度変動と折れスペクトル近似」『港湾輸送物理誌』第1巻第6号, 1950年, pp. 120-141.
- ^ 田中啓一郎「多層境界整合の素数分割モデル」『解析工学雑誌』第18巻第3号, 1960年, pp. 201-230.
- ^ 李明哲「Discretized Self-Homeomorphisms and Sign-Reversal Exponentials」『Transactions on Quasi-Measure Theory』Vol. 27, No. 1, 1971年, pp. 1-19.
- ^ 山本玲奈「アッツイワ方程式の教育的象徴性と誤答パターン解析」『数理教育研究』第6巻第2号, 1988年, pp. 60-74.
- ^ 小林陸「Heated Rumor: A Study of Mathematical Ethnography in 1930s Japan」『Annals of Computation Folklore』第2巻第2号, 1995年, pp. 10-39.
外部リンク
- 熱計局デジタル年報アーカイブ
- スペクトル分解講義ノート（HL-9編）
- 境界位相整合性研究会の記録
- 折れスペクトル実務事例集
- 準自己同型変換の図解ページ