高周波FEM
理論と物理
高周波電磁界のFEM
先生、高周波の電磁界をFEMで解くときの支配方程式は？
ヘルムホルツ方程式（周波数領域）：
$k_0 = \omega/c = 2\pi f/c$: 自由空間の波数。エッジ要素（Nédélec要素）で離散化し、複素対称行列を解く。
HFSSがこの手法ですね。
そう。Ansys HFSSは周波数領域FEMの代表格。四面体メッシュ+2次エッジ要素で複雑形状を高精度に解析できる。適応メッシュ機能でSパラメータの収束まで自動反復する。
まとめ
高周波FEMの「数値分散」問題——波長が短いほど誤差が積み重なる
高周波FEMには「数値分散」という本質的な問題があります。FEM近似の中を伝搬する電磁波の位相速度が、真の光速からずれてしまう現象です。この誤差はメッシュが粗いほど、また周波数が高いほど大きくなります（位相誤差は伝搬距離に比例して蓄積）。対策は①メッシュを細かくする（要素辺長を波長の1/6〜1/10以下）、②高次要素（2次・3次多項式）を使う、③DG-FEM（不連続Galerkin有限要素法）など数値分散の少ない手法を採用する——の3つです。ミリ波帯（30〜300 GHz）になるとこの問題が設計精度に直結します。
各項の物理的意味
- 電場項 $\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial t$：ファラデーの電磁誘導法則。時間変動する磁束密度が起電力を生じさせる。【日常の例】自転車のダイナモ（発電機）は、磁石を回転させることで近くのコイルに電圧が発生する——磁場が時間的に変化すると電場が誘起されるというこの法則の直接的応用。IHクッキングヒーターも同じ原理で、高周波磁場の変化が鍋底に渦電流を誘起し、ジュール熱で加熱する。
- 磁場項 $\nabla \times \mathbf{H} = \mathbf{J} + \partial \mathbf{D}/\partial t$：アンペア-マクスウェルの法則。電流と変位電流が磁場を生成する。【日常の例】電線に電流を流すと周囲に磁場が生じる——これがアンペアの法則。電磁石はこの原理で動作し、コイルに電流を流して強力な磁場を作る。スマートフォンのスピーカーも、電流→磁場→振動板の力というこの法則の応用。高周波（GHz帯のアンテナ等）では変位電流 $\partial D/\partial t$ が無視できなくなり、電磁波の放射を記述する。
- ガウスの法則 $\nabla \cdot \mathbf{D} = \rho_v$：電荷が電束の発散源であることを示す。【日常の例】下敷きで髪の毛をこすると静電気で髪が逆立つ——帯電した下敷き（電荷）から電気力線が放射状に広がり、軽い髪の毛に力を及ぼす。コンデンサ（キャパシタ）の設計では、電極間の電場分布をこの法則で計算する。ESD（静電気放電）対策もガウスの法則に基づく電場解析が基盤。
- 磁束保存 $\nabla \cdot \mathbf{B} = 0$：磁気単極子が存在しないことを表す。【日常の例】棒磁石を半分に割っても、N極だけ・S極だけの磁石は作れない——必ずN極とS極がペアで存在する。これは磁力線が「始点も終点もない閉じたループ」を描くことを意味する。数値解析では、この条件を満たすためにベクトルポテンシャル $\mathbf{B} = \nabla \times \mathbf{A}$ という定式化を用い、磁束保存を自動的に保証する。
仮定条件と適用限界
- 線形材料仮定：透磁率・誘電率が磁場・電場強度に依存しない（飽和領域では非線形B-Hカーブが必要）
- 準静的近似（低周波）：変位電流項を無視可能（$\omega \varepsilon \ll \sigma$）。渦電流解析で一般的
- 2D仮定（断面解析）：電流方向が一様で、端部効果を無視できる場合に有効
- 等方性仮定：異方性材料（珪素鋼板の圧延方向等）では方向別の特性定義が必要
- 適用外ケース：プラズマ（電離気体）、超伝導体、非線形光学材料では追加の構成則が必要
数値解法と実装
FEMの実装
高周波FEMの実装上の注意点を教えてください。
1. メッシュ: 波長$\lambda$あたり5〜10要素（2次要素の場合）
2. ポート設定: 導波路ポートにモードパターンを設定しSパラメータを抽出
3. 吸収境界: PML（Perfectly Matched Layer）で開放空間をモデル化
4. 行列解法: 大規模問題では反復法（GMRES+前処理）。HFSSのDomain Decomposition Method（DDM）で領域分割並列
周波数スイープはどうしますか？
HFSSのFast Frequency Sweepは、数点の直接解から有理関数近似で帯域全体のSパラを補間する。各周波数で解き直す必要がないため非常に効率的。
まとめ
辺要素（Nedelec要素）——なぜ高周波FEMで必須なのか
高周波電磁界FEMに通常のノード要素（ラグランジュ要素）を使うと「偽解（スプリアスモード）」という物理的に存在しないニセの固有モードが現れます。これを防ぐために開発されたのが辺要素（Nedelec要素、Whitney要素とも呼ぶ）で、自由度を節点ではなく要素の辺に割り当て、∇·E=0（電場の発散ゼロ）条件を自動的に満足させます。HFSSをはじめ主要な高周波FEMツールはほぼ全て辺要素ベース。「なぜHFSSはNedelec要素を使うのか？」と聞かれたら、「偽解を消すため」が一番シンプルな答えです。
辺要素（Nedelec要素）
電磁場解析に特化した要素。接線成分の連続性を自動的に保証し、スプリアスモードを排除。3D高周波解析の標準。
節点要素
スカラーポテンシャル定式化に使用。静磁場のスカラーポテンシャル法や静電場解析で有効。
FEM vs BEM（境界要素法）
FEM: 非線形材料・非均質媒質に対応。BEM: 無限領域（開領域問題）を自然に扱える。ハイブリッドFEM-BEMも有効。
非線形収束（磁気飽和）
B-Hカーブの非線形性をニュートン・ラフソン法で処理。残差基準: $||R||/||R_0|| < 10^{-4}$が一般的。
周波数領域解析
時間高調波仮定により定常問題に帰着。複素数演算が必要だが、広帯域特性は時間領域解析で取得。
時間領域の時間刻み
最高周波数成分の1/20以下の時間刻みが必要。暗黙的時間積分ではより大きな刻みも可能だが精度に注意。
周波数領域と時間領域の使い分け
周波数領域解析は「ラジオの特定の周波数に合わせる」ようなもの——1つの周波数での応答を効率的に計算できる。時間領域解析は「全チャンネルを同時に録画する」ようなもの——あらゆる周波数成分を含む過渡現象を再現できるが計算コストが高い。
実践ガイド
実務での適用
導波路フィルタ、RF-MEMS、パッケージ内の電磁界解析、コネクタのSパラメータ抽出が代表的。
実務チェックリスト
HFSSのアダプティブメッシュ——「デルタS」の正体
ANSYS HFSSを使うと解析設定に「最大デルタS」という収束判定パラメータがあります。これは連続する2回のメッシュ細分化サイクル間でSパラメータ（反射係数・透過係数）の変化量がこの値以下になったら「収束した」と判定する基準値です。デフォルトの0.02（2%）は多くのケースで十分ですが、フィルタの帯域外減衰特性やアンテナの広帯域特性を精密評価したい場合は0.005以下に設定することが推奨されます。ただし厳しくすればするほどパス数が増えて計算時間が伸びるので、要求精度に合った値の選択が実務の肝です。
解析フローのたとえ
モータの電磁界解析は「ギターの調律」に近い感覚です。弦の太さ（コイル巻数）とブリッジの位置（磁石配置）を調整して、最も美しい音色（効率の良いトルク特性）を引き出す。1つのパラメータを変えると全体のバランスが変わる——だからパラメトリックスタディが重要なんです。
初心者が陥りやすい落とし穴
「空気領域？ なんで空気をメッシュで切るの？」——初めて電磁界解析に触れた人がほぼ全員抱く疑問です。答えは「磁力線は鉄心の外にも広がるから」。解析領域を鉄心ぎりぎりにすると、行き場を失った磁束が壁に「ぶつかって」反射し、実際にはありえない磁束集中が起きます。部屋が狭すぎてボールが壁に跳ね返りまくる状態を想像してみてください。
境界条件の考え方
遠方の境界条件って地味ですが超重要です。「ここから先は無限に広がる空間」ということを数値的に表現する必要がある。設定を間違えると、まるで「見えない壁」があるかのように磁束が跳ね返されてしまいます。
ソフトウェア比較
ツール
| ツール | 特徴 |
|---|---|
| Ansys HFSS | 高周波FEMの業界標準。適応メッシュ+DDM |
| COMSOL RF | マルチフィジックス連成。熱-RF連成が容易 |
| CST Studio Suite | FEM Solver（周波数領域）も搭載 |
| Altair FEKO | FEM+MoMハイブリッド。大規模問題 |
高周波FEMツール——ANSYS HFSS vs COMSOL RF Module
高周波FEMの2大ツールはANSYS HFSSとCOMSOL RF Moduleだ。HFSSは設計自動化（Optimetrics）・ネットワークアナライザ連携・アダプティブメッシュが業界標準として定着しており、HPC並列計算で大規模問題にも対応する。COMSOL RF Moduleは多物理連成（熱・構造）が得意で、誘電体加熱や熱放射との連成解析が直感的にできる。2023年にはCSTがHPCクラウド対応を強化し、大規模アンテナアレイ解析での採用が増えた。選定は連成物理場の必要性とCAD連携の優先度で決まることが多い。
選定で最も重要な3つの問い
- 「何を解くか」：高周波FEMに必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端技術
hp-FEMと機械学習——次世代高周波FEMの最前線
高周波FEMの最先端では、要素サイズ（h）と多項式次数（p）の両方を自動適応させるhp-FEM（hp-適応有限要素法）が注目されています。従来のh-適応（メッシュ細分化のみ）に比べ、指数関数的な収束速度を実現できる可能性があり、特定の高周波問題では必要な自由度を1/10以下に削減できた事例もあります。さらに最近は機械学習でメッシュ生成を最適化する研究も活発化しており、「どこを細かくすべきか」をニューラルネットが予測するアプローチが論文で多数報告されています。
トラブルシューティング
トラブル
「メッシュを細かくしたら答えが変わった」——FEM収束の落とし穴
高周波FEM解析で「メッシュを細かくするほど正確になる」とは限らない。電気的に小さな要素に分割しすぎると行列条件数が悪化し、反復ソルバが収束しなくなる現象が起きる。経験則として1波長あたり少なくとも10要素（λ/10）を確保しつつ、高曲率部だけ細分するAdaptive Meshingが有効だ。HFSSは「Δ S収束基準」で自動的に再メッシュ化するが、この収束誤差設定が緩すぎると見た目は収束でも誤差が残る。解析周波数と要素サイズのバランスを可視化したコンバージェンスプロットの確認が必須だ。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——高周波FEMの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
なった
詳しく
報告