FDTD電磁場検出器

1. はじめに

Digital voltmeter measuring an electric potential field, illustrating how a detector samples a field at a given location — 電位場を測定する電圧計と同様に、検出器もまたFDTDシミュレーションにおいて電磁場をサンプリングすることで同様の役割を果たす。

FDTDシミュレーションにおいて、検出器は計算領域内に配置される仮想的な測定プローブとして機能します。これらは光学的シミュレーションおよび無線周波数（RF）シミュレーションの両方に用いられます。波長や用途は大きく異なる場合もありますが、電場および磁場の測定原理自体は同一です。??にある電圧計が適切な例となります。電圧計が電気的電位を測定するのに対し、FDTD検出器は電磁波がシミュレーション領域内の特定領域を通過する際にその電磁場をサンプリングするのです。

検出器は、FDTDメッシュによって定められた点において電磁場の値をサンプリングする。検出器を用いれば時間経過に伴う信号の変化を追跡することも可能であり、三次元シミュレーションにおいては表面全体にわたる電磁エネルギーの分布を測定することもできる。OghmaNanoにおいては、検出器の表面は3Dビュー上で格子状の平面として表示され、これは??に示されている。解析の目的に応じて、測定された電磁場の値から透過率や反射率、あるいは物体によって散乱されるパワーを算出することも可能である。

FDTDシミュレーションの実行中、検出器は時間依存の電磁場データを取得可能である。またOghmaNanoでは複素周波数領域の電磁場データを直接積分することも可能であり、全時間履歴を保存した上でFFT処理を行うことなく所定の周波数におけるスペクトル測定が実施できる。このような機能により、同一の検出器フレームワークが光学素子、アンテナ、およびRF散乱問題にも適用可能となる。

検出器の配置や設定に関する一般的な情報については、光学検出器の項をご参照ください。本ページでは、FDTD法による測定内容、およびそれらがどのようにしてスペクトルデータやパワー値として得られるかについて詳述します。

2. 検出器の設定

光学用リボンから??に示す検出器エディタを開くか、3Dビュー上で検出器面を右クリックして??に示すように「オブジェクト編集」を選択します。シミュレーションでは複数の検出器を使用することも可能で、例えばデバイスの手前と奥にそれぞれ検出器を配置したり、RF散乱体に向けた検出器を用いることもできます。

検出器エディタの画面は??に示されている。その「Configure」タブでは検出器の形状や位置、さらには正規化スペクトルといった各種オプションを設定できる。参照用または入力用の測定値に対して出力値を正規化することは、光源のスペクトルが平坦でない場合でも透過率の解析に非常に有効である。検出器表面を通過するパワーを測定する際には表面の向きも重要となる。得られる結果の符号はその法線方向を基準に定義される。

??内の「周波数メッシュ」タブでは、複素電磁場情報が積分される周波数を指定する。光学応用においてはこれらの周波数は通常波長として表され、RFシミュレーションではMHzまたはGHz単位で指定される。この周波数サンプリングは検出器の設定項目であり、幾何形状や電磁場を解像するための空間的なFDTDメッシュと混同してはならない。

Right-click menu on a detector in the 3D view showing the Edit object option — 検出器には、右クリックして「オブジェクトの編集」を選択することでもアクセス可能です。

Optical ribbon in OghmaNano showing the Optical Detectors button used to open the detector editor — 光学用リボンから光学検出器エディタを開く。

Optical detectors editor showing detector list, geometry configuration, and normalisation options between detectors — 検出器の設定パネルには、形状や位置、正規化の各種オプションが表示される。

Frequency mesh tab defining wavelength range and sampling points for FDTD spectral accumulation — 「周波数メッシュ」タブは、FDTD積分器用の波長範囲を定義する際に使用されます。

3. FDTDの周波数領域積分器

FDTD法では電磁場が連続するタイムステップごとに更新されていきますが、多くの場合、デバイスの特性は周波数の関数として記述する方が適切です。すべてのタイムステップを保存して後から処理する代わりに、検出器を用いればシミュレーション実行中に所定の周波数における応答を算出することも可能です。これは複素フーリエ積分器を利用することで実現されます。

時刻\(t_n=n\Delta t\)においてサンプリングされた電磁場成分に対し、角周波数\(\omega=2\pi f\)における累積値は以下の通りである

\[\widetilde{E}(\omega)=\sum_{n=0}^{N-1} E(t_n)\,e^{-i\omega t_n}\,\Delta t\]

各時間ステップにおいて、選択された周波数に対して複素値の寄与が加算される。その周波数で振動する成分は位相的に整合して蓄積される一方、他の周波数における寄与は十分に長い観測時間の間に互いに打ち消し合う傾向がある。チルダ記号は複素値を表す。 複素フーリエ積分まだ正規化されていない正弦波状の電場振幅である。同様の計算を電場の3つの成分すべてに対して行うことができる。 \((E_x,E_y,E_z)\) そして磁場 \((H_x,H_y,H_z)\).

この積分器は両方の値を保持する 振幅と位相この違いは重要である。透過スペクトルは電場強度に基づいて求められるが、方向性を持つ電磁パワーを計算するには電場と磁場間の位相関係が必要となる。十分な時間にわたって蓄積された定常な連続波応答においては \(T\)複素のピーク電場振幅はおおよそ \(2\widetilde{E}/T\)同様の式が他の場合にも適用される \(\widetilde{H}\)有限な観測時間や過渡現象がこの近似に影響を与える可能性があるため、励起条件に応じて積分区間を適切に選定する必要がある。

各要求周波数ごとに個別の積分器が用意されており、大規模な時間領域データを保存することなくスペクトル情報を得ることが可能である。周波数点の設定は「周波数メッシュ」タブにて行う（??）。この数値計算手法の詳細については、FDTD理論を参照されたい。

4. 検出器の出力結果

シミュレーション実行後、検出器の結果はOghmaNanoの出力タブから直接開くことが可能である。利用可能なファイルは有効化された検出器測定項目に応じて異なり、時間領域測定データ、周波数領域スペクトル、さらに較正済みの散乱シミュレーションにおいては散乱パワーデータも含まれる。

Output tab in OghmaNano showing a detector represented as a CCD-style icon — 検出器の出力結果は、出力タブ内でCCD風のアイコンとして表示される。

Expanded detector output folder showing files such as spectra and time-domain data — シミュレーション終了後に検出器によって生成されるファイル。

シミュレーション実行後、検出器の出力結果は出力タブ内に小型のCCD風アイコンとして表示されます（??）。このアイコンをダブルクリックすると、??に示す出力ファイルが表示されます。

Raw electric field spectrum versus wavelength obtained from FDTD frequency-domain accumulator — 検出器から得られた電磁場の原始スペクトル。

Normalised electric field spectrum showing transmission calculated by dividing output detector by input detector — 正規化処理済みのスペクトルであり、透過率の情報を示している。

Time-domain power signal recorded by an FDTD detector showing pulse propagation through the simulation — 検出器によって記録された時間領域におけるパワー値。

ファイル lam_E.csv および lam_E_norm.csv は、周波数領域におけるFDTD積分器から生成されるものである。前者には電場に関連する未処理のスペクトルが記録され、後者には入力測定値を基準とした正規化済みスペクトルが記録されている。これらは透過率や波長選択性を持つデバイスの解析に有用である。対応する例は ?? および ?? に示されている。

ファイル power.csv には時間領域における検出器信号が記録されており、その様子は ?? に示されている。時間領域および周波数領域の結果はいずれも同一の電磁場シミュレーションを表しているが、周波数領域のスペクトルはこのプロットされたパワー波形にFFT処理を施すことで得られるわけではない。実際には、第3節で説明するように、シミュレーション中に電場および磁場から直接積分処理が行われてスペクトルが生成されるのである。

RF散乱シミュレーションにおいては、scattered_power.csvという出力結果も得られる場合がある。これは単に物体の有無によって測定された総パワーの差異を示すものではない。この値は、メインシミュレーションと基準となる較正シミュレーションとを比較することで得られる複素散乱電場および磁場を用いて算出される。その測定方法について以下に説明する。

5. 散乱電磁パワーの測定方法

電磁波がゴルフボールや金属物体、その他の三次元構造物に遭遇すると、生成される電磁場は元の入射波と物体によって散乱された波の合成波となる。散乱現象を測定するため、OghmaNanoでは物体が存在しない状態での較正シミュレーションを実施する。この際、光源、計算領域、検出器の位置、周波数サンプリング条件、積分時間はすべて同一に設定される。その後、物体が存在する状態で再度シミュレーションを行う。この比較処理により、検出器が記録した入射波成分から散乱電磁場成分が分離される。

各検出器位置および周波数において、較正済みの電磁場成分から対象物が存在する場合に記録された電磁場成分を差し引く。

\[\widetilde{\mathbf E}_{\mathrm{sc}}=\widetilde{\mathbf E}_{\mathrm{object}}-\widetilde{\mathbf E}_{\mathrm{cal}}\qquad \widetilde{\mathbf H}_{\mathrm{sc}}=\widetilde{\mathbf H}_{\mathrm{object}}-\widetilde{\mathbf H}_{\mathrm{cal}}\]

ここで、\(\widetilde{\mathbf E}\)および\(\widetilde{\mathbf H}\)は第3節で説明した複素フーリエ累積値である。電場の大きさではなく複素電場そのものを減算することで位相情報が保持され、これが散乱電場を分離するために不可欠となる。計算には6つの直交座標系における電場成分すべてが用いられる。

電磁エネルギーの流れはポインティングベクトルによって記述される。複素正弦波状の電場振幅の場合、その時間平均値は以下の通りである。

\[\langle\mathbf S_{\mathrm{sc}}\rangle=\frac{1}{2}\operatorname{Re}\!\left(\mathbf E_{\mathrm{sc}}\times\mathbf H_{\mathrm{sc}}^{*}\right)\]

外積はエネルギーの流れの方向を表し、複素共役である\((*)\)は電場および磁場の相対位相を表す。局所的なポインティングベクトルの単位はW/m²である。定常な正弦波状の電磁場においては、フーリエ積分器を直接用いて検出器表面にわたって積分することで

\[\boxed{P_{\mathrm{sc}}(f)\simeq\frac{2}{T^2}\sum_{k} \operatorname{Re}\!\left(\widetilde{\mathbf E}_{\mathrm{sc},k}\times \widetilde{\mathbf H}_{\mathrm{sc},k}^{*}\right) \cdot\hat{\mathbf n}\,\Delta A_k}\]

この式において、\(T\)は電磁場の積分時間を表し、\(k\)は検出器のサンプリング点を指定する。また、\(\Delta A_k\)はその点が示す面積、\(\hat{\mathbf n}\)は検出器表面の法線ベクトルである。\(2/T^2\)という係数は、積分された電磁場を正弦波の振幅に変換した後、それらの時間平均値を乗算することで得られる。検出器表面全体にわたってポインティングベクトルの法線成分を加算することで、散乱パワースペクトルが求められ、その値はワット単位でscattered_power.csvとして保存される。

💡 符号の解釈について：散乱パワーとは検出器を通過する有符号のフラックスであり、電磁場強度の非負の尺度ではない。負の値は、実質的に散乱パワーが検出器の法線方向とは逆方向に流れていることを意味し、検出器が物体に向いている場合には光源方向への後方散乱を示唆する。また報告される値は、全パワーから空気のみの場合のパワーを差し引いた値とも異なる。なぜなら両パワーには入射波と散乱波の干渉項が含まれており、電磁場を先に差し引いた場合にはそのような干渉項は存在しないからである。

この較正手法は、周囲のシミュレーション条件を変更せずに、異なる形状からの散乱現象の周波数依存性を比較する際に特に有用である。具体的な適用例としては、10 GHzにおけるゴルフボールからの電波散乱やストーンヘンジからの電波散乱が挙げられる。信頼性の高い定量的結果を得るためには、十分な空間・時間分解能、一貫した較正設定、そして適切な検出器配置が不可欠である。