ストーンヘンジにおける電波散乱

1. 電波散乱の概要

電波は、その寸法や材料特性が電磁波の波長と同等である物体と相互作用する。電波が物体に遭遇すると、一部のエネルギーは反射され、一部は物質を透過し、また一部は散乱や回折によって方向が変えられる。こうした現象がアンテナや検出器に到達する信号を決定づけ、無線伝播、レーダー、無線通信の基礎となる。

本チュートリアルでは、イングランド南部に位置する先史時代の石造建造物であるストーンヘンジを、幾何学的に複雑な散乱体の具体例として用いる。そこに立ち並ぶ石柱、横木、倒れた石塊、そして空きスペースが、単一の球体や長方形ブロックよりも興味深い高周波環境を形成している。ここで用いるモデルは現存するストーンヘンジの単純化された数値モデルであり、考古学的調査に基づく正確な再現物ではない。

このシミュレーションでは、OghmaNanoの三次元有限差分時間領域法（FDTD）ソルバーを用いて電界Eおよび磁界Hを算出する。FDTD法は通常は光波長領域向けに紹介されるが、マクスウェル方程式は無線周波数領域においても解くことが可能である。ただし、物理的寸法や励起条件、材料データは無線周波数問題に適合したものでなければならない。背景情報としては、FDTD法とは何か？およびFDTD法における光源を参照されたい。

ここでは2つの関連する結果について見ていく。1つ目はシミュレーション領域内で電磁界がどのように発達していくか、もう1つはストーンヘンジが存在する場合と空気のみが存在する参照計算の場合とで、検出器における周波数依存の電力流がどのように変化するかである。その後、ストーンヘンジの代わりに単純な住宅街を設定し、別のRF散乱環境を例示する。

2. ストーンヘンジのシミュレーション作成

OghmaNanoを起動し、「新規シミュレーション」を選択します。??に表示されるテンプレートライブラリから「RF例題」を開き、??に示されているストーンヘンジのシミュレーションを選択します。ファイルの書き込み権限があるディレクトリに新規シミュレーションを保存してください。本チュートリアルでは既存の例題を利用するため、ストーンの幾何形状を手動で作成する必要はありません。

??内の例選択画面には、RF解析用に用意された初期設定例が表示されています。

3. 幾何形状を確認し、シミュレーションを実行する

シミュレーションの設定が完了すると、メインウィンドウに計算領域内のストーンヘンジが表示される（??）。赤色のメッシュが石柱を表しており、周囲の格子や着色された平面は空間メッシュおよびRFシミュレーション要素を示している。視点を回転・拡大縮小させることで、立ち残った石柱や三石構造、倒れた石塊などを詳しく観察できる。このモデルにおいて、主要な石の円形配置部分の直径は約30メートルであり、その外側にはヒールストーンも存在する。

この例では、長方形ブロックではなく多角形数の少ない表面モデルが用いられている。不規則な形状の石の輪郭から、OghmaNanoが一般的な三次元三角形メッシュにも対応可能であることが容易に理解できる。石の幾何形状はメッシュによって定義され、材料の応答特性は誘電体材料データを通じて別途指定される。

シーンを確認したら、「Run simulation」ボタン（青色の再生アイコン）をクリックするか、F9キーを押してください。OghmaNanoが計算の進捗状況を報告します（??および??）。計算が完了したら、「Output」タブ（??）を開いて結果を確認してください。提供されている例を実行する際、メッシュや材料パラメータを変更する必要はありません。

4. 電磁界の伝播状況を観察する

出力ブラウザからシミュレーションのスナップショット結果を開いてください。スナップショットビューアを利用すれば、保存済みの電磁界状態間を移動し、領域内で電磁波がどのように進展していくかを確認できます。図??、??、??にはこの過程の3つの段階が示されています。電磁界が石造構造物に遭遇し通過するにつれて、その分布はより複雑な様相を呈するようになります。

これらのスナップショットは、励起波が意図した領域へ伝播しているか、またシミュレーションが散乱体に対して適切に応答しているかを確認するのに役立ちます。交互に現れる電場パターンは空間的な干渉現象や波の位相を反映していますが、単一のスナップショットだけでは石構造体によって反射される電力の量を分離して評価することはできません。

5. 無線周波数における材料特性

ストーンヘンジは完全な金属反射体ではなく、誘電体構造としてモデル化される。単純な乾燥した石材の例として、相対誘電率がおよそ5という値が説明用として適切であるが、実際の石材の特性はその組成や水分含有量、周波数によって変化する。無線周波数帯域においては、相対誘電率の実数部および虚数部の両方が電磁波の伝播に影響を及ぼす。

材料設定画面を開いてください データベース 領域の様子は、図に示す通りである ??およびRF材料ビュー??周波数依存性の電磁データを確認するために使用する。光学的なデータセットが通常波長に対してプロットされるのとは異なり、RFデータは一般的にHz、MHz、またはGHz単位の周波数に対して表示される。OghmaNanoでは相対誘電率εの実数部および虚数部が保存される。_r = ε'_r ＋ iε''_r （ソルバーの符号規約に従う）。誘電率の実数成分は主に材料中の位相速度および波長に影響を与えるのに対し、虚数成分は選定された符号規約に基づき誘電損失を表す。導電性材料では低周波数領域において非常に大きな有効虚数相対誘電率を示すことがある。

例えば、ある材料が電気伝導率σで表される場合、その伝導率に起因する虚数相対誘電率の値は\[|\varepsilon_r''(f)|=\frac{\sigma}{2\pi f\varepsilon_0}.\]となります。したがってこの値は周波数が低くなるにつれて増大します。??では金属の例が示されていますが、ストーンヘンジの解析にはこの金属の曲線を使用してはなりません。石造物の解析には専用の誘電体材料データが必要となります。

参考までに、30 MHzにおける自由空間中の波長は約10 mとなります。εrが約5の損失の少ない誘電体中では、波長は約4.5 mになります。このことから、メートル単位の大きさを持つ石は電磁波の伝播パターンに影響を与えるのに十分な大きさであることがわかります。基礎となる電磁界方程式については、FDTD theoryもご参照ください。

6. シミュレーション領域とメッシュ

数値シミュレーション領域には、ストーンヘンジ遺跡、入射波、検出器の位置が含まれ、電磁場が十分に伝播できるだけの余裕空間も確保する必要がある。??にシミュレーション領域の設定ダイアログが示されている。??内のメッシュ制御機能を用いてメッシュ状態を確認することも可能だ。まずは例題で用意されている既定値を使用し、精度検証や境界応答の調査を行う場合にのみ領域サイズやメッシュ密度を増加させるべきである。

メッシュ間隔は、誘電体内における電磁波の波長およびインポートされた石造構造物の主要な形状特徴を十分に表現できる必要がある。一般的に、より細かいメッシュを用いることで幾何学的な表現精度や電磁界の解像度が向上するが、3次元FDTD計算に必要なメモリ量および計算時間も増加する。単一のメッシュ設定だけでは収束性が保証されるわけではない。定量的な解析を行う場合は、より細かいメッシュを用いて再度シミュレーションを実施し、検出器スペクトルにほとんど変化が見られないことを確認する必要がある。

境界条件は、開放空間における散乱計算において特に重要となる。吸収境界条件を用いることで、メッシュの端から生じる人工的な反射を最小限に抑えることができる。関連する出力結果は??に示されている。その正確な解釈は、プロット対象となる物理量に依存する。

7. 周波数に応じた散乱電力の測定

電磁界の様子を可視化することは有益だが、RFシミュレーションには定量的な結果も必要となる。本例では、電場および磁場のx、y、z各成分を記録する検出器が用いられている。この情報をもとに、単位面積あたりの電力流を表す電磁的ポインティングベクトルを算出することが可能となる。瞬時電場の場合には\[\mathbf{S}(t)=\mathbf{E}(t)\times\mathbf{H}(t).\]一方、周波数領域においては複素位相量を用いた時間平均電力流密度は\[\langle\mathbf{S}(f)\rangle=\tfrac{1}{2}\operatorname{Re}\!\left[\mathbf{E}(f)\times\mathbf{H}^{*}(f)\right].\]検出器が測定する物理量は、このベクトルのうち検出器面に垂直な成分であり、その符号は検出器の向きによって決定される。FDTD detectorsを参照すれば検出器出力に関する詳細情報が得られる。

記念物の影響を特定するためには、他の条件は同一で、一方はストーンヘンジを含みもう一方は空気のみを含む二つの計算結果を比較する。その際の検出器における電力差が \[\Delta P(f)=P_{\mathrm{Stonehenge}}(f)-P_{\mathrm{air}}(f).\] となる。ここで各 \(P(f)\) の値は、電場および磁場から適切なポインティングベクトルの成分（必要に応じて検出器表面での積分も含む）を用いて算出される。この電力差が負の値となる場合、記念物が存在すると検出器の正の方向へ流れる正味電力が減少することを意味する。ただしこのことだけから、失われたエネルギーが石に吸収されたと断定することはできない。反射や進路の変化、干渉現象なども測定値に影響を与え得るからである。

得られた周波数領域のスペクトルが ?? に示されている。この顕著なスペクトル特性から、ストーンヘンジの影響が周波数に依存していることがわかる。その解釈は光源のスペクトル、検出器の位置および向き、そして用いられた正規化方法に左右されるため、単純に構造的共鳴であると断定すべきではない。

8. ストーンヘンジから住宅街へ

同様の手法は他の多くの高周波構造物にも適用可能である。次の例として、ストーンヘンジの代わりに英国の住宅街をモデルとして用いる。このシーンには連なる家々や切妻屋根、そしてさまざまな幾何学的細部が含まれている（??）。ストーンヘンジと同様、これもソルバーの機能を実証するための低ポリゴン数モデルであり、詳細な建築モデルではない。この住宅街モデルは、ストーンヘンジとほぼ同等のサイズのシミュレーション領域を占めるように設計されており、同様の励起源および検出器配置を再利用することが可能となっている。

住宅街の幾何形状を確認または変更するには、3Dビュー画面内の街路オブジェクト（??）を右クリックし、幾何形状／メッシュ制御パネル（??）を開いてください。この編集ツールの説明を行うのは、最初のシミュレーション実行前にストーンヘンジの形状を修正する必要があるのではなく、今回は別のインポート済みモデルを扱うためです。関連する幾何形状操作については、ステップエディタのマニュアルをご参照ください。

この入門的な比較検証においては、建物群に単一の誘電体材料を割り当てることが可能である。したがってこれらは、内部が空気で満たされガラス窓や層状の屋根、個別のレンガ壁を持つ現実的な家屋ではなく、単純化された固体誘電体構造として捉えるのが適切である。それでもなお、このモデルを用いることで、現実的な建造物環境下における電磁界の透過・反射現象や空間的な遮蔽効果を調査することができる。

街路シミュレーションを実行した後、??に示される定常状態における電力密度の結果を開いてください。建物の周囲およびその間に、電力密度が高くなっている領域と低くなっている領域が確認できます。この空間的構造から、障害物が連続して存在する場合に、障害物のない空間と比べて電磁界がどのように変化するかが明らかになります。

街路部およびストーンヘンジのシミュレーションはそれぞれ物理的な目的が異なるものの、共通のFDTD方程式を用いている。一般的な三角形メッシュ幾何形状を利用する際の利点の一つがこれである。すなわち、メッシュおよび材料特性が定義されれば、電磁界ソルバーを多種多様なRF問題に適用できるのである。

9. 結果の解釈と限界

これらのシミュレーションは、無線周波数波が複雑な三次元誘電体構造物とどのように相互作用するかを示している。ストーンヘンジを対象とした例では、個々の石柱が不規則に配置された際の散乱現象が強調されており、住宅街を対象とした例ではより大きな建物が繰り返し配置された場合の波の伝播が示されている。いずれの場合も、可視化結果からは電力密度の分布状況が確認でき、検出器スペクトルからは周波数に応じた応答の変化が測定される。

いくつか重要な制限事項が存在します。実際の建築材料は不均質であり、ストーンヘンジの幾何形状も近似値に過ぎません。また空間メッシュはメッシュ間隔より十分に大きな構造のみを再現可能です。したがって本シミュレーション結果は単なるデモンストレーションとして扱うべきであり、ストーンヘンジや特定の街路における実際のRF特性を正確に予測するものではありません。定量的な工学的解析を行う際には、誘電体パラメータの検証、メッシュ収束性の確認、境界反射の調査、そして既知の基準値や実測値との比較を行う必要があります。

おめでとうございます！これで三次元RF FDTDシミュレーションを実行し、電磁場の伝播状況を観察し、周波数領域における方向性を持つ電力流の変化を抽出し、さらに幾何形状の変更が異なる電磁環境を生み出す様子を確認することができました。

10. 有用なシミュレーション出力結果

「出力」タブからはスナップショットや検出器測定結果、その他のソルバー出力にアクセスできます。具体的なファイル一覧はシミュレーションの設定次第ですが、RF散乱現象の解析においては以下のカテゴリが有用です。出力ブラウザの画面例は??に示されています。