シェルヴィツキの3次元波動式
| 分類 | 波動理論・計算物理（架空） |
|---|---|
| 提唱者 | シェルヴィツキ（Shervitsky） |
| 主な対象 | 音響・弾性・水面揺らぎ（応用） |
| 特徴 | 3次元の位相遅延を明示する補正項 |
| 成立年（通説） | 1937年 |
| 発表媒体（通説） | 『工学音響年報』 |
| 影響分野 | 非破壊検査、都市防災の避難放送（架空） |
| 注目ポイント | 「一次格子でのエネルギー保存」規則 |
シェルヴィツキの3次元波動式（しぇるびつきのさんじげんはどうしき）は、3次元空間における波の伝播を記述するとされる理論式である。流体工学と音響計算の双方に影響を与えたとされ、大学の実験講義でも「簡易版」として触れられることがある[1]。
概要[編集]
シェルヴィツキの3次元波動式は、3次元空間での波の伝播を記述する枠組みとして知られている。形式的には偏微分方程式の一種であるとされるが、実際には「数値格子に優しい補正」を含むため、解析よりも計算機実装で参照されることが多いとされる[1]。
この式は、位相のずれを「遅延ポテンシャル」と呼ばれる量で補う点に特色があるとされる。とくに初期の通信音響研究では、同一の音源でも地形によって波面が歪む現象を説明する材料として用いられたと説明される[2]。
なお、式そのものよりも、導出手順で導入される規則（一次格子エネルギー保存規則、以下「一次保存規則」）が、後年の技術標準にまで波及したとする見方がある。一方で、一次保存規則が実験的妥当性を欠くのではないかという指摘も早い段階から存在したとされる[3]。
仕組み[編集]
3次元波動式は、通常の波動方程式に対し、空間座標の回転に対しては見かけ上不変であるよう調整された項を加える、と説明されることが多い。シェルヴィツキは、位相遅延が「速度場のねじれ」に相当するという直観を採用し、遅延ポテンシャルを速度の回転成分から計算できるとしたとされる[4]。
具体例として、港湾の音響モデルに適用した際、波の屈折を再現するために「格子幅の1.0007倍を基準にする」規約が用いられたと記録されている。これは、当時の計算機が倍精度を採用する前夜であり、小数点以下の丸め誤差がしばしば観測誤差に見えてしまったことへの対症療法だった、とする説明がある[5]。
また、この式では境界条件の選び方が強調される。シェルヴィツキは、境界上でのエネルギー密度の偏りが一定以下になるよう、格子点の外挿係数を「0.5から0.5027へ」と微調整する手順を提示したとされる。ただし、その数値の由来については、後年に「当時の珈琲の抽出比率から逆算した」とする逸話が混じっているとされ、研究者の間で軽い茶化しの対象になったと報告されている[6]。
歴史[編集]
発端：停泊船の“聞こえ方”事件[編集]
シェルヴィツキの3次元波動式が成立した背景として、の旧港で起きたとされる停泊船の通信齟齬が挙げられることがある。1930年代、港湾無線の当直が「同じ周波数なのに命令が遅れて聞こえる」と繰り返し訴え、局地的な遅延現象が問題化したとされる[7]。
調査チームは傘下の臨時委員会を編成し、音響解析を担当する若手としてシェルヴィツキが招かれた。彼は、波面の遅れを2次元断面だけで説明しようとする従来モデルでは不足すると主張し、3次元補正が必要だと説いたとされる[8]。
委員会の記録によれば、1934年の夏に近郊で行われた夜間試験では、送信から到達までの遅延が最大で約12.3ミリ秒（測定装置の時計誤差±0.2ミリ秒を含む）まで変動したと報告された。シェルヴィツキは、この変動を「位相遅延ポテンシャルとして再整理すれば、格子化しても保存量が崩れない」と述べ、一次保存規則の原型を提出したとされる[9]。
数値標準化：一次保存規則の“通行証”[編集]
1937年、シェルヴィツキは式とともに一次保存規則を提案し、同時に「格子幅は全区間で同じであるべき」という設計指針も添えたとされる。この指針が、のちに音響計算の品質検査に使われることになる。研究会では、検査合格の目安として「エネルギー残差がL2ノルムで10^-5未満」とする基準が掲げられた[10]。
ところが、基準を厳密に運用すると計算時間が爆発するため、現場の計算技師の間では「10^-4でも“概ね”合格」とする黙認運用が広がったとされる。さらに、の商業計算会社がこの運用を商材化し、「一次保存規則の通行証」を持つモデルのみを広告表示できるようにしたことが、社会への浸透を加速させたという[11]。
その結果、1940年代には学術から実務へと波及し、非破壊検査装置の校正に3次元波動式が採用されたと説明される。一方で、厳密な意味での保存が成立しない条件が混在していたのではないか、という疑義も同時期に指摘されたとされる。とくに、格子外挿が絡む境界では、理論上の前提が崩れる可能性があるとされ、後の批判につながった[12]。
冷戦期の応用：避難放送の“聞き分け”[編集]
冷戦期、都市防災の文脈で「避難放送の聞こえ遅れ」を補正する研究が行われたとされる。この分野では、単なる音の強さだけでなく、到達タイミングのズレが行動に影響すると考えられた。そこで、シェルヴィツキの3次元波動式を用いた“タイミング補償”が導入されたとされる[13]。
架空の逸話として、のある区役所が導入検討会で「補償アルゴリズムの結果、サイレン到達は平均で1.6秒早まるが、標準偏差は0.09秒で広がる」と試算したとされる。しかし、会議議事録の添付資料には“どこから出たのか不明な係数”が数行だけ手書きで追記されていたという。この手書き係数が、一次保存規則の省略近似に由来する可能性があると後年推測された[14]。
さらに、東西いずれの陣営でも、式の一部が“音響迷彩の計算”へ応用されたと語られることがある。ただし、これは公式資料では否定されており、内部文書の存在だけが噂として残ったともされる。結果として、理論は表向きの防災に位置づけられつつ、影の用途として語られることがあった[15]。
批判と論争[編集]
一次保存規則が、実験室条件と現場条件のギャップに耐えられないのではないか、という批判がある。特に境界条件が複雑な場合、遅延ポテンシャルが実測の位相ずれと一致しないことがあるとされる。反論としては、測定系の帯域制限によって“位相が見えていないだけ”である、と主張されることがあった[16]。
また、式の導出に含まれる係数の由来が曖昧である点が問題視された。ある研究者は、格子幅を「1.0007倍」とする規約が、理論的に必要というより当時の装置較正に合わせた名残だと指摘したとされる[5]。その指摘を受けて、シェルヴィツキの弟子筋とされる人物が「係数は自然の方から合図してくる」と詩的に語った、と伝わるが、同時にこの発言は“出典のない文学”として処理されたと報告されている[17]。
一方で、論争の多くは“数学”ではなく“運用”をめぐるものだった。計算機の丸め誤差が観測と区別できない時代には、厳密な理論よりも再現性が重視される傾向があったとされる。だからこそ3次元波動式は広がったが、後に計算機が高精度化した際には、運用差による再解釈が続出したという[12]。
脚注[編集]
関連項目[編集]
脚注
- ^ Elżbieta K. Wysocka『波面遅延の工学的再構成：1930年代港湾試験の再読』Springer, 1972.
- ^ A. L. Morozov「The “One-Grid Conservation” Heuristic and Its Aftermath」『Journal of Applied Acoustics』Vol. 41, No. 3, 1989, pp. 201-219.
- ^ シェルヴィツキ『三次元波動式の導出と遅延ポテンシャル』工学音響年報, 1937.
- ^ 渡辺精一郎『格子化誤差と保存則の関係：簡易版三次元波動論』養賢堂, 1956.
- ^ Marta R. Nowak「Phase Delay Potentials in Urban Siren Models」『Proceedings of the International Society of Computational Physics』Vol. 12, No. 2, 1994, pp. 77-104.
- ^ 佐久間啓介『防災サイレンの時間設計：聞こえ遅れ補償の数理』日本建築学会, 1961.
- ^ John P. Ketteridge『Boundary Coefficients and Hidden Assumptions』Oxford University Press, 2003.
- ^ Klaus E. Brunner「丸め誤差が“理論”に見える瞬間」『数値計算通信』第5巻第1号, 1978, pp. 33-51.
- ^ Viktor P. Shervitsky『Three-Dimensional Wave Computation for Everyone』Cambridge University Press, 1979.
- ^ Robert L. Thorne『Acoustics Without Grids: A Counterfactual History』Academic Press, 2011.
外部リンク
- シェルヴィツキ式保存則アーカイブ
- 港湾音響試験データベース（暫定）
- 一次保存規則検定室
- 遅延ポテンシャル図書館
- 避難放送タイミング補償研究会