空洞共振器の電磁界シミュレーション

まさにそう。電子レンジのマグネトロンは2.45GHzの空洞共振器だ。金属の壁で完全に囲まれた空間に電磁波を閉じ込めると、特定の周波数でだけ定在波が立つ。これが「共振」で、空洞共振器の基本原理なんだ。

たくさんある。通信衛星の帯域フィルタ、粒子加速器の加速空洞、レーダーの発振器、5G基地局のフィルタ——どれも空洞共振器の応用だ。FEM解析では共振周波数とQ値を求めるのが主目的で、Q値が高いほど狭帯域で鋭い共振になる。フィルタでは高いQ値が求められるし、広帯域アンテナでは逆に低いQ値が欲しい。用途によって設計方針がまるで違うんだ。

まず一番シンプルな矩形空洞から始めよう。辺の長さが $a \times b \times d$ の矩形空洞の共振周波数は、マクスウェル方程式を境界条件（壁面で接線電場がゼロ）つきで解くと、きれいな解析解が出る。

円筒空洞（半径 $a$、高さ $d$）は少し複雑で、ベッセル関数の零点が出てくる。

覚える必要はない。表を引くか、数値計算で求める。実務で重要なのは、モード指数 $(m,n,p)$ が異なると共振周波数も電磁場パターンも変わること。低次モード（TE₁₀₁やTM₀₁₀）は場分布が単純で使いやすいけど、周波数が上がると高次モードが密集して目的のモードを見分けにくくなる。これが「モード混雑（mode crowding）」で、シミュレーションの大きな課題なんだ。

ざっくり言うと「振動の持続力」だ。ワイングラスを指で弾くと「キーン」と長く響くよね。あれはQ値が高い。段ボール箱を叩いても「ボコッ」としか鳴らない。Q値が低い。空洞共振器も同じで、壁面での電力損失が少ないほどQ値が高くなる。

壁面損失は表面抵抗 $R_s$ と壁面での接線磁場 $\mathbf{H}_t$ から計算する。

実務では「無負荷Q」と「負荷Q」を区別する。外部回路との結合を考慮しないのが無負荷Q（$Q_0$）、結合まで含めたのが負荷Q（$Q_L$）。結合係数を $\beta$ とすると次の関係がある。

その通り。反射がゼロになる最適結合状態だ。粒子加速器では臨界結合が理想だけど、フィルタ設計では意図的に過結合（$\beta > 1$）や疎結合（$\beta < 1$）を使い分ける。代表的なQ値の目安としては、通信フィルタ用の銅空洞で $Q_0 \sim 5{,}000$〜$20{,}000$、超伝導空洞（ニオブ）だと $Q_0 \sim 10^{10}$ に達する。桁がまったく違うんだ。

空洞の軸方向（たとえば $z$ 方向）に対して、電場に $z$ 成分がないのがTE（Transverse Electric）モード、磁場に $z$ 成分がないのがTM（Transverse Magnetic）モード。矩形空洞の最低次モードは通常 TE₁₀₁ で、電子レンジの基本モードもこれだ。

TE₀₁₁モードは壁面電流が接合部（角や端面と側面の境界）を横切らないんだ。だから接触抵抗の影響を受けにくい。しかも周波数が上がるほどQ値も上がるという珍しい特性を持つ。だからこのモードは周波数標準器に使われる——世界の時間を決めるセシウム原子時計にも空洞共振器が入っているんだよ。

実際の空洞共振器は、結合用のアイリス（スリット）や同軸プローブ、チューニングスクリュー、表面の微小変形がある。こういう「理想形状からのずれ」があると解析解は使えなくなる。そこでFEM（有限要素法）の出番だ。ヘルムホルツ方程式を弱形式に変換して、固有値問題として解く。

$[S]$ はcurl-curl行列（剛性行列に相当）、$[T]$ は質量行列に相当するもの。固有値 $k_0^2$ から共振周波数 $f = k_0 c / (2\pi)$ が求まり、固有ベクトル $\{e\}$ から電磁場のモードパターンが分かる。構造解析の固有振動解析と全く同じ枠組みなんだ。

いい質問だ。マクスウェル方程式には $\nabla \cdot \mathbf{D} = 0$（空洞内に電荷がない場合）という条件がある。普通の節点要素で電場をスカラー的に補間すると、この発散ゼロ条件が自動的に満たされない。その結果、物理的に存在しない「偽の共振モード」が計算結果に紛れ込む。これがスプリアスモードだ。

辺要素（Nedelec要素）は自由度を「辺」に割り当てて、電場の接線成分の連続性だけを保証する。法線成分は要素界面で不連続を許容するから、$\nabla \cdot \mathbf{D} = 0$ の条件を離散レベルで正確に取り込めるんだ。これによりスプリアスモードが排除される。現在の高周波FEM解析では辺要素が事実上の標準だ。

そうだ。2次辺要素は曲面形状のフィッティングが良く、同じメッシュ密度でも1次要素より精度が格段に高い。特に円筒空洞のような曲面を持つ形状では、1次要素だと形状近似誤差が支配的になる。HFSSのAdaptive Mesh Refinementも2次辺要素の上で動く設計になっている。

方法は大きく3つある。

設計初期は摂動法（計算が速い）、詳細設計ではSパラメータ法（結合構造の影響まで含む）が一般的だ。超伝導空洞のように極端にQ値が高い場合は、複素固有値法だと固有値の虚部が極めて小さくなって数値精度が足りなくなるから、摂動法が必須になる。

目的による。

その通り。Q値が10,000の空洞だと、エネルギーが $1/e$ に減衰するまでに約 $Q/\pi \approx 3{,}000$ 周期分の時間ステップが必要になる。10 GHzの空洞なら $3{,}000 / (10 \times 10^9) = 0.3 \,\mu\text{s}$ のシミュレーション時間。CFL条件で決まる時間刻みが数psだとすると、10万ステップ以上回さないといけない。FEMの固有値解析なら1回の行列分解で終わる。

典型的な解析フローはこうなる。

表皮深さが壁厚よりはるかに小さい（銅壁で数mm以上、表皮深さは1μm以下）から、電磁波は壁を透過しない。だからPEC近似が非常に良い精度で成り立つ。Q値計算のときだけ壁面の表面抵抗 $R_s$ を後処理で入れればいい。ただし、壁厚が表皮深さと同程度の薄膜構造（超伝導薄膜空洞など）では壁厚のモデリングが必要になるよ。

空洞共振器のメッシュで最も重要なのは「波長あたりの要素数」だ。経験則として、2次辺要素で1波長あたり6要素以上が目安。1次辺要素だと10〜15要素必要になるから、結局2次の方が効率的。

そう、HFSSのAdaptive Mesh Refinementは非常に優秀だ。初期メッシュで解いた後、誤差が大きい領域を自動的に細分化して再計算する。共振周波数の変化が収束基準（デフォルト0.02%）以下になったら完了。ただし万能ではなくて、初期メッシュがあまりに粗すぎると間違ったモードに収束する危険がある。初期メッシュの波長解像度はユーザーが責任を持って設定すべきだ。

空洞共振器解析で使う境界条件は主に以下の4種類だ。

PMCは自然界に存在しないけど、電磁場の対称性を表現するのに使う。例えば矩形空洞のTE₁₀₁モードは $xz$ 面に対して電場が対称だから、その面にPMC条件を置いて1/2モデルにできる。ただし注意点がある——対称面にPECを置くかPMCを置くかでモードが変わるんだ。PECを置くと電場が反対称なモードだけが解に出て、PMCを置くと対称なモードだけが出る。両方計算して合わせると全モードが得られる。

矩形空洞なら3つの対称面があるから、1/8モデルにできる。メッシュ要素数が1/8になると、メモリは概ね1/8、計算時間は固有値ソルバーの特性にもよるけど1/10〜1/30程度に短縮される。円筒空洞の軸対称モード（TM₀₁₀など $m=0$ のモード）は2D軸対称モデルに帰着できて、3Dの100分の1以下の計算量になる。

ただし、対称性を使うとその対称性を持たないモードは出てこない。全モードを網羅したい場合はフルモデルが必要になる。設計段階で「狙いのモードがどの対称性を持つか」を把握してからモデリングするのが基本だ。

HFSSの最大の強みはEigenmode Solverの完成度だ。2次辺要素ベースのFEMにAdaptive Mesh Refinementを組み合わせて、共振周波数を自動的に高精度で収束させてくれる。操作手順は以下の通り。

HFSSの Eigenmode結果テーブルに自動表示される。内部では壁面の表面抵抗を摂動法で評価している。壁面材料を Finite Conductivity に変更すると、表面抵抗 $R_s$ が自動計算されるんだ。注意点としては、PEC面のままだと $Q = \infty$ と出る。Q値を求めるには壁面材料を銅やアルミに変更する必要があるよ。

確かにCSTの看板はFIT（Finite Integration Technique、FDTDの変形）による時間領域解析だけど、Eigenmode Solverも搭載している。JDM（Jacobi-Davidson Method）ベースの固有値ソルバーで、FIT格子上で動作する。六面体メッシュの自動生成が速い点が利点で、形状が比較的単純な空洞なら効率的だ。

CSTが真価を発揮するのは「広帯域Sパラメータ解析」。時間領域ソルバーで短いパルスを入力して、FFTで全帯域のSパラメータを一度に取得できる。複数モードの結合特性を見たい場合や、帯域フィルタの設計には非常に便利だ。HFSSの周波数掃引より効率がいい場合が多い。

RF Moduleの Eigenfrequency Study で空洞共振器解析ができる。COMSOLの強みはマルチフィジクス連成だ。例えば粒子加速器の超伝導空洞では、電磁場の壁面損失による発熱 → 温度上昇 → 表面抵抗変化（BCS抵抗の温度依存性）→ Q値変化、というループをCOMSOL内で一貫して解ける。HFSSやCSTでは熱解析に別ソフトが必要になる。

Adaptive Mesh Refinementが弱い。HFSSのように「固有値の収束」を基準にした自動リファインメントは組み込まれていないから、メッシュ収束性の確認をユーザーが手動で行う必要がある。また、大規模モデルでの固有値ソルバーのスケーラビリティはHFSSに劣る印象だ。

できる。Elmer FEMは辺要素ベースの固有値解析に対応していて、ARPACK を使った固有値ソルバーで空洞のモードを計算できる。OpenEMSはFDTDベースで、時間領域からFFTで共振周波数を抽出する方式だ。どちらも学術研究レベルでは十分実用的だけど、商用ソフトに比べてGUIが弱く、トラブル時に自力解決が求められる。粒子加速器コミュニティではACE3P（SLAC開発）やSuperfish（LANL開発）が定番の無償ツールだ。

単体の空洞で共振周波数とQ値を求めるだけなら、HFSSが最も効率的で信頼性が高い。帯域フィルタのように結合空洞の広帯域特性が重要ならCSTの時間領域ソルバー。熱連成や構造変形の影響まで含めたい場合はCOMSOL。予算がなくて自力でコードが読めるならElmer + ARPACK。用途で使い分けるのがベストだ。

それは典型的なスプリアスモードだ。以下を確認しよう。

電場分布を可視化して $\nabla \cdot \mathbf{E}$ をプロットする。スプリアスモードでは $\nabla \cdot \mathbf{E} \neq 0$ になっている。また、メッシュを細かくしたときに周波数が大きく変動するモードはスプリアスの可能性が高い——物理的なモードはメッシュ収束性がある。

Q値の不一致は実務で最も頻繁に遭遇する問題だ。主な原因と対策をまとめる。

$\Delta$ はRMS（二乗平均平方根）表面粗さ、$\delta$ は表皮深さだ。$\Delta / \delta \ll 1$ なら影響は無視できるけど、10 GHz の銅だと $\delta \approx 0.66 \,\mu\text{m}$ で、機械加工面の表面粗さ $\Delta \approx 0.5$〜$1.5 \,\mu\text{m}$ と同オーダーだから無視できない。電解研磨すると $\Delta \approx 0.05 \,\mu\text{m}$ まで下がるから、Q値が理論値に近づく。

モード同定（mode identification）は高次モードが密集する帯域で特に厄介だ。対策は以下の通り。

固有値ソルバーの収束不良は大きく2パターンある。

パターン1：Adaptive Meshが収束しない

パターン2：固有値ソルバー自体が失敗する

鉄則は3つ。まず深呼吸——焦って設定をいじると問題が複雑化する。次に、最小再現ケースを作る——アイリスやチューニングスクリューを外した単純空洞で解析解と一致するか確認する。最後に、一度に1つだけ変えて再実行する——科学実験と同じ「対照実験」の原則だ。複数の変更を同時にすると、何が効いたか分からなくなるからね。

形状関数 $N_i$ を用いて未知量を近似:

これを数式で表すとこうなるよ。

これを数式で表すとこうなるよ。

連続体の支配方程式を離散化すると、以下の代数方程式系が得られる:

ここで $[K]$ は全体剛性マトリクス（または同等のシステムマトリクス）、$\{u\}$ は未知節点変数ベクトル、$\{F\}$ は外力ベクトルなんだ。

収束速度: 二次要素で $O(h^2)$ のオーダーで誤差が減少（滑らかな解の場合）

空洞共振器解析の実務的な解析フローと注意点を解説する。

1. 前処理 (Pre-processing)

2. 求解 (Solving)

3. 後処理 (Post-processing)

Ansoft Corporationが開発した3D高周波電磁界シミュレータ。2008年にAnsysがAnsoftを買収。

Computer Simulation Technology (ドイツ) が開発。2016年にDassault Systèmesが買収しSIMULIAに統合。

うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。

空洞共振器解析に対応する主要な商用CAEツールの機能比較と、各製品の歴史的背景を詳述する。

Ansoft Corporationが開発した3D高周波電磁界シミュレータ。2008年にAnsysがAnsoftを買収。

Computer Simulation Technology (ドイツ) が開発。2016年にDassault Systèmesが買収しSIMULIAに統合。

1986年スウェーデンで設立。MATLAB連携のFEMLABとして開始、後にCOMSOLに改名。マルチフィジックスに強み。

空洞共振器解析のツール選定においては以下を考慮:

うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。

空洞共振器解析における最新の研究動向と先進的手法を見ていこう。

症状: ソルバーが指定反復回数内に収束せず異常終了

考えられる原因:

対策:

症状: 応力/変位/温度等が物理的に非現実的な値

考えられる原因:

対策:

症状: 計算が想定時間の何倍もかかる

対策:

症状: Out of Memory エラー

対策: