モーメント法（MoM）— 積分方程式ベースの電磁界解析手法
理論と物理
概要 — MoMとは何か
モーメント法ってFEMやFDTDと何が違うんですか？ どれも電磁界を解く手法ですよね？
いい質問だ。MoMは積分方程式ベースで、導体表面の電流分布だけを未知量にする。体積メッシュが不要で、開領域問題に圧倒的に強い。ワイヤーアンテナやRCS（レーダー断面積）計算ではFEMより桁違いに効率的だよ。
体積メッシュ不要って、めちゃくちゃ楽じゃないですか！
そう、それがMoMの最大の武器だ。FEMでは構造物の周囲の空気まで含めて3D体積メッシュを切らなきゃいけない。しかも開放空間を扱うには吸収境界条件（PML等）が必要。MoMでは導体表面だけを離散化すれば、グリーン関数が無限遠の放射条件を自動的に満たしてくれる。NEC2はフリーのMoMコードで、アマチュア無線家にも広く使われているよ。
NEC2ってあのアマチュア無線の人がよく使っているやつですよね？ プロの現場でも使えるんですか？
NEC2はワイヤーアンテナに特化した薄線近似のMoMコードで、ダイポールアンテナや八木アンテナの設計では今でも第一選択肢だ。ただし面要素が扱えないから、パッチアンテナや筐体を含む複雑な形状にはAltair FekoやAnsys HFSSのMoMソルバーを使う。本記事ではMoMの理論的基礎から実装まで、体系的に見ていこう。
FEM・FDTDとの本質的な違い
電磁界の数値解析には大きく3つのアプローチがある。それぞれの定式化と得意領域を理解することが、適切な手法選択の鍵となる。
| 手法 | 定式化 | 離散化対象 | 行列の性質 | 得意分野 |
|---|---|---|---|---|
| MoM | 積分方程式 | 導体表面のみ | 密行列（小規模） | アンテナ放射、RCS、開領域 |
| FEM | 微分方程式（弱形式） | 体積全体 | 疎行列（大規模） | 複雑形状、不均質媒質、共振器 |
| FDTD | 微分方程式（差分） | 体積全体（直交格子） | 行列不要（陽解法） | 広帯域過渡応答、非線形 |
支配方程式 — EFIE（電界積分方程式）
MoMの数式を見てみたいです。どんな方程式がベースなんですか？
基本となるのはEFIE（Electric Field Integral Equation、電界積分方程式）だ。完全導体（PEC）表面上では接線方向の全電界がゼロになるという境界条件から導かれる。
完全導体表面 $S$ 上で、入射電界 $\mathbf{E}^{inc}$ と表面電流 $\mathbf{J}(\mathbf{r}')$ の関係は次のEFIEで記述される：
ここでグリーン関数 $G(\mathbf{r}, \mathbf{r}')$ は自由空間中の点源応答を表す：
うわ、結構複雑ですね…。ざっくり言うとどういう意味ですか？
ざっくり言うと「導体上のあらゆる場所の電流が、グリーン関数という"伝播の重み"を介して、互いに電磁的に影響し合う」という式だ。左辺が入射波（外から来る電磁波）、右辺がその電流から生まれる散乱界。両者が導体上で釣り合う電流分布 $\mathbf{J}$ を求めるのがMoMの仕事だよ。
各項の物理的意味
- $\hat{n} \times \mathbf{E}^{inc}$：入射電界の接線成分。外部から導体に入射する電磁波が、表面電流を誘起する「ソース」に相当する。例えば、レーダー波がステルス機に当たるときの入射場。
- $j\omega \mu \int G \cdot \mathbf{J} \, dS'$：ベクトルポテンシャル $\mathbf{A}$ に由来する項。電流の時間変動が生み出す磁界成分で、アンテナの放射に直結する。
- $\frac{1}{j\omega\varepsilon} \nabla \int G \, \nabla' \cdot \mathbf{J} \, dS'$：スカラーポテンシャル $\phi$ に由来する項。電荷蓄積（$\nabla' \cdot \mathbf{J}$）が生み出す電界成分。低周波ではこの項が支配的になる。
- $G(\mathbf{r}, \mathbf{r}')$（グリーン関数）：「$\mathbf{r}'$ にある点電流源が $\mathbf{r}$ に作る電界」を表す核関数。$e^{-jkR}/4\pi R$ という形は球面波の伝播そのもの。距離 $R$ に反比例して減衰し、位相 $kR$ だけ遅延する。
仮定条件と適用限界
- 完全導体（PEC）仮定：表面の有限導電率を無視。マイクロ波帯の金属導体ではほぼ正確だが、光周波数や薄膜ではインピーダンス境界条件への拡張が必要
- 自由空間グリーン関数：均質な無限空間を前提。地面反射がある場合はSommerfeld積分や像理論で修正
- 表面等価定理：構造物の体積電流を表面電流に置き換えられるのはPECまたは均質誘電体に限られる
- 周波数領域：基本のMoMは単一周波数での解析。広帯域特性は周波数掃引またはAWE（漸近波形評価）で取得
次元解析と単位系
| 変数 | SI単位 | 備考 |
|---|---|---|
| 表面電流密度 $\mathbf{J}$ | A/m | 体積電流密度 [A/m²] と区別すること |
| グリーン関数 $G$ | 1/m | $e^{-jkR}/(4\pi R)$ の次元 |
| 波数 $k$ | rad/m | $k = 2\pi/\lambda = \omega\sqrt{\mu\varepsilon}$ |
| インピーダンス行列要素 $Z_{mn}$ | Ω（オーム） | 基底関数の積分結果として得られる |
| 電界 $\mathbf{E}$ | V/m | 入射波の電界強度 |
Galerkin法による離散化
この積分方程式をコンピュータで解くにはどうするんですか？
未知の表面電流 $\mathbf{J}$ を基底関数（basis function）$\mathbf{f}_n$ の線形結合で近似する。これがMoMの核心だ。
未知電流を $N$ 個の基底関数で展開する：
$I_n$ は各基底関数の係数（未知数）。この展開をEFIEに代入し、テスト関数 $\mathbf{f}_m$ との内積（テスティング）を取る。Galerkin法ではテスト関数と基底関数に同じ関数を用いる：
ここで $\mathcal{L}$ はEFIEの積分演算子、$\langle \cdot, \cdot \rangle$ は内積を表す。
なるほど、FEMのGalerkin法と同じ考え方なんですね！
そう、重み付き残差法という意味ではFEMと同じ枠組みだ。違いは演算子が微分ではなく積分であること。だからMoMでは"Method of Weighted Residuals applied to Integral Equations"とも呼ばれる。Harringtonが1968年の著書で体系化したのが始まりだよ。
インピーダンス行列 [Z][I]=[V]
Galerkin法のテスティングにより、連続的な積分方程式は以下の行列方程式に変換される：
各成分は次のように定義される：
$[Z]$ はインピーダンス行列、$[I]$ は電流係数ベクトル、$[V]$ は励振ベクトルってことですね。FEMの $[K]\{u\} = \{F\}$ と構造が似てる！
その通り。ただし決定的な違いが1つある。FEMの剛性マトリクス $[K]$ は疎行列（隣接要素間しか結合がない）だが、MoMのインピーダンス行列 $[Z]$ は密行列（フル行列）だ。グリーン関数を通じて全ての基底関数が互いに結合するからね。
ということは、未知数が $N$ 個あったらメモリが $O(N^2)$ 必要ってことですか？
鋭いね。直接法で解くと $O(N^2)$ のメモリと $O(N^3)$ の計算時間が必要。だからMoMは「小規模だけど正確」な手法として長年使われてきた。これを打破したのがMLFMM（多重レベル高速多重極法）で、$O(N \log N)$ まで計算量を削減できる。詳しくは数値解法のセクションで説明するよ。
基底関数の選択
基底関数 $\mathbf{f}_n$ って具体的にどんな形をしているんですか？
代表的なものを表にまとめよう。問題の形状に応じて使い分けるんだ。
| 基底関数 | 形状 | 適用対象 | 特徴 |
|---|---|---|---|
| パルス関数 | 区間内一定値 | ワイヤー（1D） | 最も単純。NEC2で使用 |
| 三角関数 | 区間内線形補間 | ワイヤー（1D） | 電流の連続性を保証 |
| RWG関数 | 三角形ペア上のベクトル | 面要素（2D） | Rao-Wilton-Glisson。面MoMの標準 |
| 高次基底関数 | 多項式 | 面要素（2D） | WIPL-Dで採用。少ない要素で高精度 |
RWG関数が面のMoMでは定番なんですね。FEMでいう辺要素みたいなものですか？
いい比較だ。辺要素が接線成分の連続性を保証するのと同様に、RWG関数は隣接三角形ペアの共有辺を通じて法線方向の電流連続性を保証する。1982年にRao・Wilton・Glissonが提案して以来、面MoMの事実上の標準だ。例えばAltair FekoもRWG基底関数をベースにしているよ。
Harringtonが切り拓いた積分方程式の世界
モーメント法の体系はRoger Harrington（Syracuse大）が1968年の著書『Field Computation by Moment Methods』で確立した。表面電流を基底関数で展開し、テスト関数との内積でインピーダンス行列を構築するアプローチは、それまで解析的に扱えなかった任意形状アンテナの放射・散乱問題を数値的に解く道を開いた。この著書はたった110ページほどの薄い本だが、電磁界計算の歴史を変えた一冊として今も読み継がれている。Harrington自身は「学生が理解できない理論は、理論が悪い」と常に語っていたそうだ。
数値解法と実装
密行列問題とその対処
さっき密行列で $O(N^2)$ メモリが必要と聞きましたが、具体的にどれくらいの規模が限界なんですか？
目安として、$N = 10{,}000$ の未知数なら $[Z]$ は $10^4 \times 10^4 = 10^8$ の複素数要素で、メモリは約1.6 GB。$N = 100{,}000$ なら160 GB。LU分解による直接解法だと $O(N^3)$ の演算量だから、$N$ が増えるとあっという間に計算不能になる。
| 未知数 $N$ | 行列サイズ | メモリ（複素倍精度） | LU分解時間（概算） |
|---|---|---|---|
| 1,000 | $10^6$ | 16 MB | 数秒 |
| 10,000 | $10^8$ | 1.6 GB | 数分 |
| 100,000 | $10^{10}$ | 160 GB | 数日 |
| 1,000,000 | $10^{12}$ | 16 TB | 直接法では不可能 |
16 TB って...一般的なワークステーションでは無理ですよね。大きい問題はどうするんですか？
そこで登場するのがMLFMMだ。行列を陽に作らず、行列ベクトル積を $O(N \log N)$ で近似計算する。反復解法（GMRES等）と組み合わせれば、百万未知数クラスの問題もワークステーションで解ける。
MLFMM（多重レベル高速多重極法）
MLFMM（Multi-Level Fast Multipole Method）は、遠方の基底関数グループ間の相互作用を多重極展開で近似する手法だ。
多重極展開って…何か直感的なイメージはありますか？
例えるなら「遠くの都市の人口」を調べるとき、個人ごとに数えるのではなく「東京都は約1400万人」とまとめて扱うイメージだ。近くの人（隣接セル）は正確に計算し、遠くのグループはまとめて近似する。この「近くは正確、遠くはまとめる」階層構造がMLFMMの本質だよ。
| 解法 | メモリ | 行列ベクトル積 | 全体計算量 |
|---|---|---|---|
| 直接法（LU分解） | $O(N^2)$ | — | $O(N^3)$ |
| 反復法（行列を陽に保持） | $O(N^2)$ | $O(N^2)$ | $O(N^2 \cdot n_{iter})$ |
| MLFMM + 反復法 | $O(N \log N)$ | $O(N \log N)$ | $O(N \log N \cdot n_{iter})$ |
EFIE・MFIE・CFIEの使い分け
EFIEのほかにMFIEとかCFIEもあるって聞いたんですが、何が違うんですか？
3つの積分方程式にはそれぞれ特徴がある：
- EFIE（Electric Field IE）：開放面・ワイヤー構造に適用可能。ただし閉じた導体では内部共振周波数で解が一意に定まらない
- MFIE（Magnetic Field IE）：閉じた導体のみに適用可能。内部共振問題はEFIEと同様
- CFIE（Combined Field IE）：EFIEとMFIEの線形結合 $\alpha \cdot \text{EFIE} + (1-\alpha) \cdot \text{MFIE}$。内部共振を完全に排除でき、MLFMMとの組み合わせで反復収束も良好。$\alpha = 0.5$ が一般的
閉じた導体のRCS計算なら迷わずCFIEってことですね。
その通り。Altair FekoでもCFIEがデフォルト設定だ。ただしワイヤーアンテナのように開いた構造ではMFIEが使えないからEFIE一択になる。FekoやNEC2はこの使い分けを自動的に判断してくれるよ。
特異積分の処理
グリーン関数って $|\mathbf{r} - \mathbf{r}'| \to 0$ で発散しますよね。自分自身の要素を積分するときどうするんですか？
いいところに気づいた！ これが「特異積分（singular integral）」問題で、MoM実装の最も難しい部分の1つだ。主な対処法は：
- 特異点抽出法：$1/R$ の特異部分を解析的に積分し、残りの滑らかな部分を数値積分
- Duffy変換：三角形要素上の自己項を座標変換で非特異化
- 適応的ガウス積分：観測点と積分領域の距離に応じて積分点数を増やす
実務では、近傍要素（$R < 0.5\lambda$）は高精度積分、遠方要素は低次ガウス積分またはMLFMMで処理する。
周波数領域と時間領域の使い分け
基本のMoMは周波数領域の手法で、単一周波数での定常応答を計算する。「ラジオの特定のチャンネルに合わせる」感覚だ。広帯域特性が必要な場合は周波数掃引を行うが、周波数点ごとにインピーダンス行列を再構築する必要がある（計算コスト大）。これを回避するためにAWE（漸近波形評価）やモデル次数低減が使われる。一方、時間領域MoM（TD-MoM）は「全チャンネルを同時録画」するようなもので、超広帯域パルス応答の解析に向いているが、安定性の問題があり実装が難しい。
実践ガイド
MoM解析の実務フロー
実際にMoMで解析するとき、どんなステップを踏むんですか？
典型的なMoM解析フローを説明しよう：
- 形状モデリング：ワイヤー構造はセグメント分割、面構造はCADからの三角形メッシュ作成
- 励振の設定：給電点（voltage gap）、入射平面波、近傍源など
- メッシュの妥当性確認：セグメント長・三角形辺長が $\lambda/10$ 以下であることを確認
- ソルバー設定：EFIE/CFIE選択、MLFMM使用の有無、周波数掃引範囲
- 解析実行：インピーダンス行列構築 → 連立方程式求解 → 電流分布取得
- 後処理：遠方界パターン、ゲイン、入力インピーダンス、RCSの計算
表面メッシュのベストプラクティス
$\lambda/10$ のルールは分かりましたが、他に注意すべきポイントはありますか？
MoMのメッシュはFEMとは異なるポイントがある：
- 辺長の上限：最高解析周波数での $\lambda/10$ が目安。$\lambda/7$ まで緩められる場合もあるが、精度に注意
- エッジ部：導体のエッジ（端部）には電流が集中するため、エッジ近傍はメッシュを細かくする
- アスペクト比：三角形のアスペクト比は5以下を推奨。細長い三角形はRWG関数の精度を劣化させる
- ワイヤーと面の接合：ワイヤーセグメントの端点を面メッシュの節点に一致させる
- 給電ギャップ：NEC2のVoltage Sourceでは給電セグメントの長さが入力インピーダンスに影響する。$\lambda/1000$ 程度が推奨
ワイヤーモデリングの注意点
NEC2でダイポールアンテナを解析するとき、ありがちな失敗ってありますか？
NEC2のワイヤーモデリングは「薄線近似」が前提で、ワイヤー半径 $a$ がセグメント長 $\Delta$ より十分小さくなければならない。具体的には：
- $a / \Delta < 1$：これが破れると薄線近似が崩壊。$a / \Delta < 0.01$ が理想
- 隣接ワイヤーの間隔：2本のワイヤーが近接する場合、間隔がワイヤー半径の5倍以上必要
- セグメント数の等分割：1本のワイヤーを等間隔に分割するのが基本。電流変化が大きい部分だけ密にする不等分割も可能
- 接合角度：ワイヤー同士の接合角度が90度未満だと電流の不連続が生じやすい
なるほど、FEMのメッシュ品質とは全然違う指標なんですね。ワイヤーの半径とセグメント長の比率って意外と見落としそうです。
よくある失敗と対策
| 症状 | 原因 | 対策 |
|---|---|---|
| 入力インピーダンスが異常値 | 給電セグメントが長すぎる / メッシュが粗い | 給電近傍のメッシュを細密化 |
| 内部共振で解が発散 | 閉じた導体にEFIEのみ使用 | CFIEに切替（$\alpha = 0.5$） |
| MLFMM反復が収束しない | 多重極展開の打切り次数が不足 | 展開次数を増加、前処理を変更 |
| パターンの対称性が崩れる | メッシュの非対称性 | 対称メッシュの再生成 |
| 低周波でコンディションが悪化 | EFIEの低周波分解（low-frequency breakdown） | ループ-スター基底関数を使用 |
初心者が陥りやすい落とし穴 — 「周波数を上げたら未知数が爆発した」
MoMでは表面メッシュの辺長を波長の1/10以下に保つ必要がある。つまり周波数を2倍にすると波長が半分になり、2D表面メッシュの要素数は約4倍に増える。10 GHzでは十分だったモデルが、30 GHzにしたら未知数が約9倍、メモリが約81倍（密行列の場合）になるのだ。高周波問題でMoMを使う場合はMLFMMの利用が必須。それでも限界を超える場合は物理光学近似（PO）やGTDとのハイブリッド手法を検討しよう。
境界条件の考え方 — MoMには「遠方境界」の設定がない
FEMやFDTDでは空間を有限領域で打ち切り、吸収境界条件（PML等）を設定する必要がある。PMLの距離や層数が不適切だと反射が生じて結果がおかしくなる。MoMではグリーン関数が無限空間の放射条件を自動的に満たすため、遠方境界を設定する必要がない。これがアンテナやRCSのような開放問題でMoMが好まれる最大の理由だ。
アンテナと筐体の影響 — MoM-POハイブリッド法
車載アンテナや航空機アンテナを設計するとき、アンテナ単体ではなく「大型導体筐体（車体・機体）と組み合わせた状態」でのパターンが重要だ。筐体の回折・反射がアンテナゲインを大幅に変化させる。MoM（特にMoM-PO混合法）は、アンテナ周辺の細部をMoMで精密に扱い、遠方の大型筐体を物理光学近似（PO）で近似することで、巨大モデルを現実的な計算時間で解析できる。Altair FekoのHYBRIDソルバがこのアプローチの代表実装として広く使われている。
ソフトウェア比較
NEC2 — フリーのMoMコード
まずはNEC2の詳しい話を聞きたいです。どんな経緯で生まれたんですか？
NEC（Numerical Electromagnetics Code）はLawrence Livermore国立研究所がNASAの資金で1970年代に開発した。NEC-2は1981年にリリースされてパブリックドメインになった。もう40年以上前のコードだが、ワイヤーアンテナ解析の精度は今でも折り紙付きだ。
- 定式化：EFIE（ワイヤー）+ MFIE（面）の混合
- 基底関数：ワイヤー上の三角関数（正弦・余弦の組み合わせ）
- 得意分野：ダイポール、八木、ログペリオディック、ヘリカル等のワイヤーアンテナ
- 制限：面パッチは平面のみ（曲面不可）、誘電体非対応
- ライセンス：パブリックドメイン（NEC-2）。NEC-4はLLNLのライセンスが必要
アマチュア無線の人たちはGUIフロントエンドを使ってるんですよね？
そう。代表的なのは4nec2（Windows、無料）とEZNEC（有料、より直感的なUI）だ。NECのカードデッキ入力（GW/EX/FR等のコマンド形式）を知らなくてもGUIでアンテナを設計できる。アマチュア無線家がアンテナの指向性パターンやSWRを最適化するのに日常的に使っているよ。
Altair Feko
プロの現場ではどんなツールが使われていますか？
MoM専用ソルバーとして業界で最も広く使われているのがAltair Fekoだ。南アフリカのステレンボッシュ大で開発が始まり、EM Software & Systems社を経て、2014年にAltairが買収した。
- MoMソルバー：RWG基底関数ベース。EFIE/MFIE/CFIE対応
- MLFMM：大規模問題向け高速ソルバー内蔵
- ハイブリッド手法：MoM/MLFMM/PO/RL-GO/UTDを自動混合。業界最多の混合解法
- 得意分野：車載アンテナ配置、航空機RCS、EMC/EMI、反射鏡アンテナ
- 特徴：CMA（特性モード解析）に対応。アンテナの基本モードを物理的に理解できる
その他のMoM対応ツール
| ツール名 | 開発元 | MoM実装 | 特徴 |
|---|---|---|---|
| WIPL-D | WIPL-D社（セルビア） | 高次基底関数MoM | 要素数を大幅に削減でき効率的 |
| Ansys HFSS（IE-Solver） | Ansys Inc. | MoM（Integral Equation Solver） | FEMソルバーとの切替・連成が容易 |
| CST Studio Suite | Dassault Systemes SIMULIA | MoM（Integral Equation Solver） | FDTD/FEM/MoMの切替が可能 |
| COMSOL（RF Module） | COMSOL AB | BEM（MoM相当） | マルチフィジクス連成が強み |
| 概念（CONCEPT-II） | ハンブルク工科大 | MoM | EMC解析向け、学術用途 |
機能比較マトリクス
結局、MoMで解析するならどれを選べばいいですか？
用途別に整理しよう。
| 用途 | 推奨ツール | 理由 |
|---|---|---|
| ワイヤーアンテナの設計 | NEC2（4nec2） | 無料・高精度・豊富な実績 |
| 車載/航空機アンテナ配置 | Altair Feko | 大型筐体のハイブリッド解析 |
| RCS計算 | Altair Feko / WIPL-D | CFIE + MLFMM の成熟度 |
| FEMとの連成が必要 | Ansys HFSS（IE-Solver） | 同一環境でFEM↔MoM切替 |
| EMC解析（全波+ケーブル） | CST Studio Suite | ケーブルハーネスのMoM解析 |
| 学術研究・カスタム実装 | 自作コード / NEC2改造 | アルゴリズム検証の自由度 |
なるほど、MoMと一口に言っても用途で全然ツールが違うんですね。
その通り。「MoMに強い」と謳っていても、FekoのようにMoM専業として磨かれたツールと、HFSSやCSTのように「FEMやFDTDの補助としてMoMも使える」ツールでは、大規模問題での性能や安定性がかなり違う。ベンチマーク問題で比較してから選定することを勧めるよ。
NEC2とアマチュア無線コミュニティ
NEC2はパブリックドメインとなったことで、世界中のアマチュア無線家（ham radio operator）に爆発的に広まった。CQ誌やARRL（米国無線中継連盟）の記事では「NEC2で最適化した八木アンテナ」の設計例が無数に紹介されている。特に4nec2というWindowsフロントエンドが登場してから、テキストベースのカードデッキ入力から解放され、一気にユーザー層が広がった。プロのアンテナエンジニアの中にも「最初のアンテナ設計はNEC2だった」という人は多い。40年以上前のコードが今も現役というのは、Harringtonの理論とMoMの普遍性の証だろう。
先端技術
ハイブリッド手法（MoM-PO, MoM-FEM）
MoMの先端研究ってどんな方向に進んでいるんですか？
最も実用的で成果が出ているのがハイブリッド手法だ。異なる解法の長所を組み合わせることで、単一手法では扱えない問題を解く。
- MoM-PO：アンテナ周辺の精密部分をMoM、大型筐体を物理光学（PO）で近似。航空機やRCS解析で標準的
- MoM-FEM：導体部分をMoM、誘電体・吸収体をFEMで解析。レドーム付きアンテナの定番
- MoM-UTD/GTD：遠方のエッジ回折を幾何光学（GTD/UTD）で処理。計算コストをさらに削減
FekoのHYBRIDソルバーがまさにそれですね！
そう。Fekoは1つのモデル内でMoM/MLFMM/PO/RL-GO/UTD/FEMを混在させられる。ユーザーが各領域の解法を指定すると、境界面で自動的に結合条件を処理してくれる。例えば航空機の翼端に載せたアンテナの解析では、アンテナ本体（MoM）+ 翼表面（MLFMM）+ 機体遠方部（PO）という3層の混合が典型的だ。
広帯域MoMとモデル次数低減
周波数掃引が遅いという問題はどう解決されていますか？
いくつかのアプローチがある：
- AWE（Asymptotic Waveform Evaluation）：ある周波数点でのインピーダンス行列とその微分から、近傍周波数の応答をPade近似で外挿。周波数点数を10〜100倍削減
- CBFM（Characteristic Basis Function Method）：大型構造を部分領域に分割し、各部分の特性基底関数を事前計算。周波数掃引時の計算量を大幅に削減
- CMA（Characteristic Mode Analysis）：構造固有のモードを周波数非依存で計算。アンテナの物理的挙動の理解と最適設計に活用
機械学習との融合
最近のAI/MLブームはMoMにも影響していますか？
注目されている研究分野がいくつかある：
- サロゲートモデル：MoMの入出力関係をニューラルネットワークで学習。パラメトリック最適化の反復計算を高速化
- 前処理の学習：MLFMMの反復解法に最適な前処理行列をデータ駆動で構築
- メッシュ最適化：電流分布の予測に基づく適応的メッシュ細密化
ただし、物理法則に基づくMoMの精度はMLで代替できるものではない。あくまで「MoMを加速する道具」としてMLを活用する方向が主流だ。
トラブルシューティング
内部共振問題
閉じた導体のMoM解析で結果がおかしくなることがあると聞いたのですが…
症状：特定の周波数でインピーダンス行列の条件数が急激に悪化し、解が発散する
原因：閉じた導体のEFIE（またはMFIE）では、導体内部の共振周波数で非一意性が生じる。これは物理的な共振ではなく、数学的な定式化に起因するアーティファクト。
対策：
- CFIEに切り替える（$\alpha \cdot \text{EFIE} + (1-\alpha) \cdot \text{MFIE}$、$\alpha = 0.5$推奨）
- CFIEは内部共振を完全に除去できる
- 開放構造（ワイヤー）では内部共振は発生しない
MLFMM収束不良
MLFMMで反復が収束しないとき、何をチェックすべきですか？
症状：GMRES/BiCGStab等の反復解法が指定回数内に収束しない
考えられる原因と対策：
- 多重極展開の次数不足 → Fekoでは "Number of multipoles per box" を増加
- EFIEの使用 → CFIEに切替。CFIEは行列のスペクトル特性が良好で収束が速い
- 前処理の不適切 → ブロック対角前処理やSAI（Sparse Approximate Inverse）を試す
- メッシュ品質 → 極端に小さい/大きい三角形がないか確認
- 反復回数の上限 → デフォルト値が小さい場合は増加（500〜1000）
メモリ不足（密行列）
先生、MoMで「Out of Memory」が出たんですが…
症状：インピーダンス行列の確保でメモリ不足
対策（優先順）：
- MLFMMを有効にする：行列を陽に保持しなくなるため、メモリが $O(N^2) → O(N \log N)$ に削減
- メッシュを粗くする：$\lambda/10$ → $\lambda/7$ に緩和（精度低下に注意）
- ハイブリッド手法：遠方部分をPO近似に切替え、MoM未知数を削減
- Out-of-core解法：行列をディスクに退避（速度低下あり）
- 並列計算：MPI分散メモリで複数ノードに行列を分散
MoMの全体像がつかめました！ 積分方程式ベースでFEMとは全然違うアプローチなのに、行列方程式に帰着させるところは同じなんですね。まずはNEC2でワイヤーアンテナから試してみます。
いいね！ NEC2で半波長ダイポールを解析して、理論値（入力インピーダンス 73+j42.5 Ω）と比較するのが最初の一歩としてお勧めだ。セグメント数を変えて収束性も確認してみてくれ。分からないことがあったらいつでも聞いてくれ。
「解析が合わない」と思ったら
- まず深呼吸 — 焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る — 例えば半波長ダイポールのような既知の問題で、同じ設定が正しく動くか確認
- 1つだけ変えて再実行 — メッシュ密度、積分方程式の種類、周波数のいずれか1つだけ変更して影響を切り分ける
- 物理に立ち返る — 電流分布に非物理的なスパイクがないか可視化して確認。対称構造の電流が対称でなければメッシュの問題
なった
詳しく
報告