導波管モード解析

矩形導波管ではTE₁₀が基本モードだ。カットオフ周波数以上でのみ伝搬できる。レーダーや衛星通信で使うKu帯（12〜18 GHz）の標準導波管WR-62ではTE₁₀のカットオフが9.49 GHzになる。

TEは Transverse Electric（横電界）で、伝搬方向の電界成分 $E_z = 0$ のモード。TMは Transverse Magnetic（横磁界）で $H_z = 0$ のモードだ。同軸ケーブルのようにTEMモード（$E_z = H_z = 0$）が存在する構造もあるが、単一導体の導波管ではTEMモードは存在しない。これは導波管内部の電界がかならず横方向に「定在波」を形成する必要があるからだ。

ざっくり言うと、導波管の壁で電磁波が反射を繰り返しながらジグザグに進んでいくイメージだ。壁の間隔で「ちょうど整数個の半波長が収まる」条件を満たすモードだけが存在できる。TE_mnの m と n は、それぞれ幅方向と高さ方向に何個の半波長が入るかを表す。TE₁₀なら「幅方向に半波長1個、高さ方向に0個」だ。

そうだ。カットオフ周波数 $f_c$ 未満では電磁波はエバネッセント波になって、指数関数的に減衰する。導波管はハイパスフィルタとして機能するんだ。

例えばXバンド（8〜12 GHz）で使うWR-90導波管は $a = 22.86$ mm、$b = 10.16$ mmだ。TE₁₀のカットオフは $f_{c,10} = 3 \times 10^8 / (2 \times 0.02286) = 6.56$ GHz。次のモードTE₂₀のカットオフは $13.12$ GHz。だから6.56〜13.12 GHzの帯域ではTE₁₀だけが伝搬する「シングルモード帯域」になる。設計では通常、カットオフの1.25倍〜1.9倍の範囲で使うんだ。

違う。導波管内の位相速度は自由空間の光速より速い。ただし情報やエネルギーが伝わる群速度は光速より遅い。相対論に矛盾しないよ。

そういうことだ。だから実務では、カットオフの1.25倍以上の周波数で使う。それ未満だと分散が大きすぎて、パルス信号が歪んでしまう。レーダーのパルス圧縮などでは群速度の周波数依存性を正確に把握しておかないと、距離精度がガタ落ちになる。

まさにそこがマイクロ波回路設計のキモだ。導波管-同軸変換（例えばプローブカプラ）では、プローブの挿入長や位置を調整して広帯域でVSWR（電圧定在波比）を1.2以下に抑えるのが実務の目標。CAEではSパラメータ解析で反射特性S₁₁を周波数掃引して、-20 dB以下であればまず合格だ。

TE₀₁モードは壁面の電流が周方向にしか流れないんだ。周波数が上がるほど壁面の電流密度が下がって損失が減る。1970年代にAT&Tが長距離通信用にTE₀₁円形導波管を研究していたが、光ファイバに取って代わられた。でも今、テラヘルツ帯の伝送路として再び注目されている。

マクスウェル方程式から出発して、導波管断面上の2次元ヘルムホルツ方程式に帰着させるんだ。$z$方向に $e^{-j\beta z}$ で伝搬する仮定を置くと、断面内の成分だけの固有値問題になる。

電磁場解析で一番重要なポイントだ。普通の節点ベースの要素（スカラー要素）で電界ベクトルを解くと、物理的に存在しない偽のモード——スプリアスモード——が大量に出てくる。これはマクスウェル方程式の $\nabla \cdot \mathbf{E} = 0$ の条件を要素間で保証できないからだ。

辺要素（Nedelec要素）は「辺に沿った接線成分」を自由度として持つ。これにより電界の接線成分の連続性は保証しつつ、法線成分の不連続（異種媒質界面で物理的に正しい挙動）も自然に表現できる。結果として $\nabla \cdot \mathbf{E} = 0$ が近似的に満たされ、スプリアスモードが排除される。現代の商用ソフト（HFSS、CST FEM、COMSOL RF Module）はすべて辺要素ベースだ。

その通り。固有値 $k_c^2$ からカットオフ周波数 $f_c = c k_c / (2\pi)$ が得られ、固有ベクトルが対応するモードの電界分布を表す。HFSSのEigenmode Solverは、このタイプの問題を直接解く。

FDTD（時間領域差分法）も使える。導波管内でパルスを励振し、時間応答をFFTして共振周波数を特定する方法だ。CSTのTransient Solverがこのアプローチを取る。FEMと比べて広帯域特性を一度の計算で得られるのが利点だが、高Q構造では長い時間ステップが必要になるのが欠点だ。

基本的なワークフローは以下の4ステップだ。矩形導波管なら最初の解析は30分もかからない。

大きく2つ。まず最高周波数の波長 $\lambda$ に対して要素サイズを $\lambda/5$ 以下にすること。二次要素なら $\lambda/3$ でもOKだ。次に、壁面近傍で電界が急変する（表皮効果）ので、損失計算する場合はスキンデプス $\delta$ の3分の1以下のメッシュ層が壁面に必要になる。

基本的には優秀だが、盲信は禁物だ。初期メッシュがあまりにも粗いと、適応精密化が正しい方向に進まないことがある。最初の手動メッシュで $\lambda/3$ 程度のシードを入れておくのが現場の知恵だ。あと、適応精密化の収束基準は $\Delta S$ だけでなく、求めたいカットオフ周波数やQ値の変化率も別途確認すべきだ。

TE₁₀モードは幅方向（$x$方向）に $\sin(\pi x / a)$ の分布を持つ。$x = a/2$ の中央面で電界がピークになるので、この面にPMC対称条件を置けば半分のモデルで計算できる。ただし対称条件を使うと、その対称性を持たないモード（TE₂₀など）は結果に出てこないので要注意だ。全モードの網羅的探索が目的なら対称条件は使わないほうがいい。

高周波電磁界解析のデファクトスタンダードはAnsys HFSSだ。元々はAnsoft社が開発した製品で、2008年にAnsysが買収。適応メッシュ精密化（Adaptive Mesh Refinement）が最大の特徴で、ユーザーがメッシュサイズを手動で追い込む必要がほとんどない。

CST（現Dassault Systemes SIMULIA）の最大の特徴はFDTDベースのTime Domain Solverだ。パルス励振→FFTで広帯域のSパラメータを1回の計算で得られる。もちろんFEM Frequency Domain SolverやEigenmode Solverも搭載している。HFSSがFEM一本勝負なのに対し、CSTは問題に応じて手法を選べる「マルチソルバー」思想だ。

COMSOLの強みは、電磁場解析に熱解析や構造解析を簡単に連成できるところだ。例えば、大電力導波管で壁面の発熱→熱膨張→断面寸法の変化→カットオフ周波数のシフト、という連成問題を1つのモデルで解ける。衛星用の高出力TWTAの導波管設計では、このような熱-電磁連成が必要になる場面がある。

OpenEMSを推奨する。MATLABのインターフェースで矩形導波管のモード解析をやれば、カットオフ周波数や電界分布の理解が深まる。論文や卒論ではOpenEMS + Pythonの組み合わせで十分な品質の結果が出せる。就職後に商用ソフトに移行する際も、電磁界解析の原理が分かっていれば習得は早い。

典型的なスプリアスモードだ。節点ベースの要素を使っているか、辺要素でも定式化に問題がある場合に発生する。$f = 0$ 付近のモードは静電界モードに対応する非物理的な解だ。

メッシュ収束性の確認は必須だ。3水準以上のメッシュ密度で $f_c$ を計算し、変化率が0.1%以下になったら収束とみなす。リチャードソン外挿を使えば、収束値の推定精度をさらに上げられる。

ポートの定義ミスが最も多い原因だ。Wave Portの積分線が正しい方向を向いていないか、ポートのサイズが導波管断面と一致していないケースだ。HFSSではポートのフィールド分布を「Port Field Display」で可視化して、TE₁₀モードのパターン（半波長の正弦波）が出ていることを必ず確認しよう。

以下のチェックリストは最低限やるべきことだ。1つでもNGなら結果を信用してはいけない。