FDTD法（有限差分時間領域法） — CAE用語解説
FDTD法とは
先生、電磁波シミュレーションでFDTD法ってよく聞くんですけど、一言で言うとどんな手法ですか？
ざっくり言うと、Maxwell方程式を空間・時間ともに有限差分で直接離散化して、電場と磁場を交互に時間更新していく陽解法だよ。1966年にKane Yeeが提案した方法で、正式名称はFinite-Difference Time-Domain法。行列を解かずに時間を刻むだけだから、アルゴリズムが驚くほどシンプルなんだ。
行列を解かないって、FEMとは全然違うアプローチですね。具体的にはMaxwell方程式のどの形を使うんですか？
時間領域の回転方程式を直接使うよ。ソースがない場合はこんな形だ：
上がファラデーの法則、下がアンペールの法則（導電率 $\sigma$ による損失項つき）。この2本を交互に差分更新するのがFDTD法の本質だよ。例えばレーダーの電磁波がステルス機にぶつかって散乱する様子を、まるで動画を撮るように1フレームずつ追いかけられるんだ。
Yeeセルの構造
「Yeeセル」っていう独特な格子構造が肝だと聞いたんですけど、普通の格子と何が違うんですか？
最大のポイントは電場と磁場を空間的に半格子ずらして配置すること。3次元の $\Delta x \times \Delta y \times \Delta z$ の直方体セルを考えると：
- 電場成分（$E_x, E_y, E_z$）はセルの辺の中点に置く
- 磁場成分（$H_x, H_y, H_z$）はセルの面の中心に置く
こうすると、ある $H$ 成分の周りをぐるっと囲む4つの $E$ 成分がちょうど半格子離れた位置にあるから、Maxwell方程式の回転（curl）を中心差分で自然かつ正確に表現できるんだ。
なるほど、空間だけじゃなくて時間方向にもずらすんですよね？
その通り！ リープフロッグ（蛙跳び）方式と呼ばれていて、時刻 $n$ で $\mathbf{E}$ を更新したら、次に時刻 $n+\tfrac{1}{2}$ で $\mathbf{H}$ を更新する。また時刻 $n+1$ で $\mathbf{E}$ を更新…という繰り返しだ。例えば $E_x$ の更新式は1次元で書くと：
ここで $C_a, C_b$ は材料の誘電率 $\varepsilon$ と導電率 $\sigma$ から決まる更新係数だよ。各セルごとに掛け算と足し算だけだから、GPU並列化との相性が極めていい。
セルが直方体ってことは、曲面をきれいに表現できないですよね？
鋭い指摘だね。曲面は「階段近似（staircasing）」になってしまうのがFDTD法の弱点の一つだ。円形パッチアンテナの縁を階段で切ると共振周波数がずれることもある。対策としては、局所的にメッシュを細かくするサブグリッディングや、Dey-Mittra法のような曲面補正テクニックがあるよ。でも球体レドームみたいに全体が曲面の構造だと、正直FEMのほうが楽な場合が多い。
Courant安定条件（CFL条件）
陽解法ということは、時間刻みに制限があるんですよね？ 構造解析の陽解法でもCFL条件ってありましたけど…
まさに同じ考え方だよ。FDTD法ではCourant安定条件（CFL条件）を満たす必要がある。3次元の場合：
ここで $c$ は光速（真空中なら $3 \times 10^8$ m/s）。等方メッシュ $\Delta x = \Delta y = \Delta z = \Delta$ なら：
$$\Delta t \leq \frac{\Delta}{c\sqrt{3}}$$この条件を1 bitでも超えると、解が指数関数的に発散して使い物にならなくなるよ。
実務だと、メッシュサイズはどうやって決めるんですか？ 波長に対してどのくらい細かくすればいいんでしょう？
古典的な目安は最短波長の1/10以下。つまり最高周波数 $f_{\max}$ に対して：
例えば28 GHzのミリ波5Gアンテナなら、$\lambda_{\min} \approx 10.7$ mm だから $\Delta \leq 1.07$ mm くらい。精度を上げたいなら $\lambda/20$ にする場合もある。ただしメッシュを半分にすると3次元だと要素数が8倍、Courant条件で時間ステップ数も2倍になるから、計算コストは約16倍に跳ね上がる。ここのバランス感覚が電磁解析エンジニアの腕の見せ所だね。
Courant数を限界ギリギリに設定しても大丈夫ですか？
理論的にはギリギリでも安定だけど、実務ではCourant数を0.9〜0.95くらいに設定するのが一般的だ。浮動小数点の丸め誤差で微妙に超えてしまうリスクを避けるためだよ。あと、PML境界や分散性媒質を入れると実効的な安定限界が少し下がることもあるから、少しマージンを取るのが安全だね。
PML吸収境界条件
FDTD法って計算領域が有限ですよね。領域の端で電磁波が反射しちゃわないんですか？
まさにそこが大問題で、1994年にBerenger が発明したPML（Perfectly Matched Layer）が革命的な解決策になった。計算領域の外周に8〜16セルくらいの「仮想吸収層」を設けて、そこに入った電磁波を減衰させて消してしまうんだ。
普通に損失媒質（導体とか）を置くのとは違うんですか？
全然違う！ 普通の損失媒質を置くと、波がその境界でインピーダンス不整合を起こして反射してしまう。PMLのすごいところは、入射角にも周波数にもよらず理論上の反射係数がゼロになること。物理的にはあり得ない異方性をもった仮想媒質なんだけど、Maxwell方程式上は完全に整合するように設計されている。
今の主流はUPML（Uniaxial PML）やCPML（Convolutional PML）という改良版で、特にCPMLはエバネッセント波の吸収性能も高くて実装もしやすいから、多くのFDTDコードで標準採用されているよ。
PMLの層数って増やせば増やすほどいいんですか？
理論上はそうだけど、実用上は8〜16層あればPML表面での反射を $-60$ dB〜$-80$ dB以下に抑えられる。それ以上増やしてもメモリと計算時間の無駄になるだけだね。むしろ重要なのはPMLの導電率プロファイル（多項式グレーディングの次数）のチューニングで、一般的には3次〜4次の多項式で導電率を層の奥に向かって徐々に増やすのがベストプラクティスだ。急に大きくすると離散化誤差で逆に反射が増えてしまうよ。
広帯域パルス解析
FDTD法の強みは「1回の計算で広帯域」と聞きましたけど、具体的にどうやるんですか？
励振源としてガウスパルスや変調ガウスパルスを使うんだ。例えばこんな波形：
$f_0$ が中心周波数、$\tau$ がパルス幅で、$\tau$ を短くするほど周波数帯域が広がる。これを入力ポートから打ち込んで時間発展を計算し、入射波と反射波・透過波の時間波形を記録する。計算が終わったらFFTで周波数領域に変換すれば、Sパラメータ（$S_{11}$, $S_{21}$ など）が一気に全周波数で求まるよ。
FEMの周波数領域解法だと、1周波数ずつ計算するんですよね？ それと比べるとめちゃくちゃ効率的じゃないですか。
そう、例えばUWB（Ultra-Wide Band）アンテナで3〜12 GHzを100ポイント評価したいとする。周波数領域FEMだと100回の連立方程式求解が必要だけど、FDTD法なら1回のシミュレーションで全部カバーできる。5Gミリ波や車載レーダー（77 GHz帯）のように広い帯域が必要な用途では、FDTD法の計算効率が圧倒的に有利だね。
逆に、高Q共振器やフィルタのような非常に狭い帯域の構造だと、時間領域の波形がなかなか減衰しなくて計算ステップ数が膨大になるから、そういう場合はFEMのほうが効率的だよ。
アンテナの遠方界パターンもFDTD法で求められるんですか？ 遠くまでメッシュを引っ張るのは無理ですよね…
いい質問だね。近傍-遠方変換（NTFF: Near-to-Far-Field Transformation）を使うんだ。アンテナを囲む仮想閉曲面上の近傍界データをFDTDで計算して、それを等価電磁流に変換してから解析的に遠方界を求める。だから遠方までメッシュを延ばす必要はないよ。放射パターン、指向性利得、偏波特性なんかを効率的に評価できるんだ。
商用ソルバー比較：HFSS vs CST
実務ではCST Studio SuiteやAnsys HFSSがメジャーですよね。FDTD法との関係を教えてください。
整理するとこうなるよ：
| 項目 | CST Time Domain Solver | Ansys HFSS |
|---|---|---|
| 手法 | FIT（≒FDTD） | FEM（周波数領域） |
| メッシュ | 六面体（構造格子） | 四面体（非構造格子） |
| 得意 | 広帯域・過渡・大規模 | 高Q共振・曲面形状 |
| 帯域効率 | 1回で全帯域 | 周波数ごとに求解 |
| GPU対応 | 非常に強い | 対応あり（限定的） |
| 用途例 | 5Gアンテナ、EMC、レーダー | 導波管フィルタ、コネクタ |
CSTの時間領域ソルバーはFIT（Finite Integration Technique）ベースだけど、六面体メッシュ上ではFDTD法と数学的にほぼ等価だ。だからCSTでFDTD的な広帯域解析の恩恵がフルに受けられるよ。
じゃあ、HFSSが不要かというとそうでもないんですね？
全然そんなことない！ 導波管フィルタやキャビティ共振器みたいなQ値が数千〜数万の構造だと、FDTD法では波が減衰するまで何十万ステップも回す必要があって非効率。HFSSなら共振周波数付近の数点を解くだけで済む。あと、HFSSは適応メッシュ細分化（Adaptive Meshing）が非常に優秀で、初心者でもメッシュ設定をあまり気にせず精度の高い結果が得られるという実務上の利点がある。
最近はHFSS側にもTransient Solverがあるし、CSTにもFEM Frequency Domain Solverがあるから、両方のツールが互いの領域に攻め込んでいる状況だね。使い分けのポイントは「帯域幅が広いか」「形状が曲面主体か」「過渡現象を見たいか」の3軸で判断するといいよ。
実務Tips・よくある落とし穴
FDTD法を使うときに、実務で「やらかしがち」な失敗ってありますか？
よくある落とし穴をいくつか挙げるよ：
- PMLと構造の距離不足 — 解析対象とPML境界の間に最低 $\lambda/4$ 以上の空間を確保しないと、エバネッセント波がPMLに入り込んで反射が増える。特にパッチアンテナの接地面端部は要注意。
- メッシュの非等方比が大きすぎる — $\Delta x : \Delta z = 1 : 10$ みたいな扁平セルを使うと、Courant条件で $\Delta t$ が極端に小さくなり計算時間が爆増する。アスペクト比は3以下に抑えるのが目安。
- シミュレーション打ち切りが早すぎる — 時間波形が十分減衰する前に計算を止めると、FFT結果にリップルが乗る。信号のピーク値から $-30$ dB以下まで減衰したのを確認してから打ち切ること。
- 分散媒質の実装ミス — 周波数依存の誘電率（Debyeモデル、Drude モデル等）を使うとき、更新式の導出を間違えると因果律が破れて非物理的な結果になる。補助微分方程式（ADE）法で導入するのが安全だよ。
GPU計算が向いているという話でしたけど、実際どのくらい速くなるんですか？
FDTD法はセルごとの更新が完全に独立だから、GPUの超並列アーキテクチャと相性抜群だ。CSTのGPUアクセラレーションだと、CPUだけの場合と比較して10倍〜50倍の高速化が報告されている。例えば車体全体のEMC解析（10億セル規模）でも、GPU 4枚で数時間で終わるような時代になっているよ。学術系のオープンソースFDTDコード（gprMax、MEEP、openEMS等）もCUDA/OpenCL対応が進んでいるから、商用ツールがなくても試せるよ。
まとめ
ここまでの話をまとめると、FDTD法のキーポイントは何になりますか？
5つにまとめるとこうだ：
- Yeeセル：電場と磁場を空間的に半格子ずらし、Maxwell方程式の回転を中心差分で忠実に離散化
- リープフロッグ時間進行：陽解法で行列求解不要。GPU並列化に最適
- Courant安定条件：$\Delta t \leq \Delta/(c\sqrt{3})$。破ると即発散。実務では0.9〜0.95倍で設定
- PML吸収境界：開放空間を有限領域で模擬。CPMLが現在の標準
- 広帯域パルス解析：ガウスパルス + FFT で1回の計算から全周波数特性を取得
FDTD法はアルゴリズムが直感的でデバッグしやすく、時間発展のアニメーションで物理現象を「目で見て理解」できるのも大きな魅力だよ。電磁波がアンテナからどう放射されるか、散乱体にぶつかってどう回折するか――それを動画として観察できる手法は他にあまりないからね。
CAE用語の正確な理解は、チーム内のコミュニケーションの基盤です。 — Project NovaSolverは実務者の学習支援も視野に入れています。
FDTD法の実務で感じる課題を教えてください
Project NovaSolverは、CAEエンジニアが日々直面する課題——メッシュサイズの決定、PMLパラメータの調整、計算コストの削減——の解決を目指しています。あなたの実務経験が、より良いツール開発の原動力になります。
お問い合わせ（準備中）関連トピック
なった
詳しく
報告