誘電体共振器の電磁界シミュレーション

ざっくり言うと、高誘電率セラミック（$\varepsilon_r = 20 \sim 90$）の内部に電磁波が閉じ込められて共振するデバイスだ。金属壁で囲む空洞共振器とは違って、誘電体と空気の誘電率差で電磁波を「内部反射」させるんだよ。

うん。光ファイバーと同じ全反射の原理だよ。$\varepsilon_r$ が大きいほど波長が短くなる——$\lambda = \lambda_0 / \sqrt{\varepsilon_r}$ だから、同じ周波数でも共振器を小さくできる。例えば $\varepsilon_r = 36$ のセラミックなら、金属空洞の $1/6$ のサイズで済む。5G基地局のフィルタや発振器の周波数安定化に使われているんだ。

無負荷Q値（$Q_u$）で1万〜10万と極めて高い。金属空洞は導体損失で $Q$ が数千程度に留まるけど、誘電体は $\tan\delta$ が $10^{-4}$ 以下の低損失材料を使えるからね。だから基地局フィルタで帯域外の信号をスパッと切れるわけだ。

金属空洞では純粋なTEモード・TMモードが存在するけど、DRは開放構造だから軸方向に電場も磁場も残る「ハイブリッドモード」になるんだ。HEモードは磁場成分が支配的、EHモードは電場成分が支配的なハイブリッドモードだよ。

最も広く使われるのは TE_01δ モードだ。軸方向（z方向）に電場成分 $E_z$ を持たない最低次のモードで、電界が円周方向に分布する。添字の $\delta$ は、軸方向の閉じ込めが不完全——つまり電磁場が誘電体の外にも若干しみ出すことを示しているんだ。Q値が最も高く、他モードとの周波数分離も良いからフィルタ設計の定番だよ。

デュアルモードフィルタでは HE_11δ モードの二重縮退（2つの直交偏波）を利用して、1つのDRで2つの共振を実現することもある。これで素子数を半分にできるから、衛星通信の入出力フィルタで重宝されているよ。

円柱形DRの場合、厳密解は存在しないけど、Courtney のモデルや Kajfez の近似式がよく使われる。TE_01δ モードの代表的な近似式はこれだ：

$a$ は円柱の半径、$h$ は高さだ。例えば、$\varepsilon_r = 36$、$a = 5\,\text{mm}$、$h = 4\,\text{mm}$ のDRだと $f_0 \approx 4.8\,\text{GHz}$ になる。実務ではこの式で初期値を推定して、あとはFEMで精密に追い込むという流れだね。

Q値は「蓄えたエネルギー」と「1周期あたりの損失」の比だ。式で書くと：

DRの場合、損失は3つに分けられる。誘電体損失 $Q_d$、導体損失 $Q_c$（金属筐体がある場合）、放射損失 $Q_r$ だ：

誘電体損失が支配的なことが多くて、$Q_d = 1/\tan\delta$ で近似できる。$\tan\delta = 10^{-4}$ の材料なら $Q_d \approx 10{,}000$ だ。金属シールド筐体に入れると $Q_c$ は数万程度、$Q_r$ は筐体内なら事実上無限大と見なせる。だから $Q_u \approx Q_d$ になることが多いんだよ。

温度係数 $\tau_f$（単位: ppm/K）で温度変化による周波数シフトを評価する：

ここで $\tau_\varepsilon$ は誘電率の温度係数、$\alpha_L$ は線膨張係数だ。実用的なDR材料では、組成を調整して $\tau_f \approx 0\,\text{ppm/K}$ に追い込む。例えばBa-Zn-Ta系（BZT）は $\tau_f = 0 \pm 2\,\text{ppm/K}$ を実現していて、$-40°$C〜$+60°$Cの環境変化でも周波数ずれが 0.01% 以下に収まるんだ。

そう。COMSOLなどのマルチフィジクスツールで熱-電磁場連成解析をやると、筐体内の温度分布に起因する局所的な $\varepsilon_r$ 変動が共振周波数にどう影響するかまで評価できる。特に高出力フィルタでは自己発熱の効果が無視できないから、連成解析が必須になってきているよ。

近似式は孤立した単純円柱DRにしか使えない。実製品では金属筐体の中に入っていたり、結合孔でDR同士が繋がっていたり、支持台や調整ネジが付いていたりする。こういう複雑な形状の共振周波数やQ値を正確に出すには、マクスウェル方程式をフルウェーブで解く必要があるんだ。

時間調和場を仮定して $e^{j\omega t}$ を分離すると、電場に関するヘルムホルツ型の固有値方程式が得られる：

ここで $k_0 = \omega/c$ が固有値になる。辺要素で離散化すると、一般化固有値問題になる：

$[S]$ はカール-カール行列、$[T]$ は質量行列に相当する。$k_0^2$ を固有値として解くと、各固有値から $f_n = k_0 c / (2\pi)$ で共振周波数が求まる。Q値は固有値を複素数に拡張して、$Q = \text{Re}(k_0) / (2 \cdot \text{Im}(k_0))$ で計算するんだ。

これは電磁場FEMの最重要ポイントだ。節点要素（スカラー要素）でベクトル場を近似すると、$\nabla \cdot \mathbf{E} = 0$ の発散条件を自動的に満たせない。その結果、物理的に存在しない「スプリアスモード」が大量に固有値に混入して、本物のモードを見分けるのが困難になるんだ。

辺要素（Nedelec要素）は自由度を辺上に配置して、電場の接線成分の連続性だけを保証する。法線成分は不連続でもよい——これが異なる誘電率の界面で電場の法線成分が不連続になるという物理を正しく再現するんだ。結果として発散条件が自動的に満たされ、スプリアスモードが排除される。

ポイントは3つだ。まず、誘電体と空気の界面を十分に細かく切ること。ここで電場が急激に変化するからね。目安は誘電体内の波長 $\lambda_d = \lambda_0 / \sqrt{\varepsilon_r}$ の $1/10$ 以下。次に、金属壁近傍——表皮深さ $\delta_s$ の $1/3$ 以下の要素サイズが必要。最後に、空気領域の外側境界は誘電体から少なくとも $\lambda_0 / 4$ 以上離すこと。

大きく影響する。空気領域を狭くしすぎると、しみ出した電磁場が外側境界で反射して、共振周波数がずれてしまう。特に $\varepsilon_r$ が低いDR（$\varepsilon_r < 30$）では電磁場のしみ出しが大きいから、空気領域を広めに取る必要があるんだ。

主に4つだ。金属壁は PEC（完全電気導体）、外側の開放境界には PML（完全整合層） または ABC（吸収境界条件）、DRの軸対称性を利用するなら 対称面の活用（PMCまたはPEC面）。対称性を使えば計算規模を $1/4$ や $1/8$ に削減できるから、必ず検討すべきだよ。

ABCは1次近似だから入射角によっては反射が残る。PMLは人工的な吸収層で、理論上はどんな角度でも無反射。DRの固有値解析では電磁場のしみ出し方向が事前に分からないから、PMLの方が安全だよ。HFSSはデフォルトで放射境界条件を使うけど、正確なQ値を出すにはPMLに切り替えた方がいい場合もある。

典型的な解析フローはこうだ：

カタログ値は「室温・特定周波数」での代表値だから、注意が必要だ。実際には $\varepsilon_r$ と $\tan\delta$ は周波数依存性と温度依存性を持っている。例えば Ba(Mg_1/3Ta_2/3)O₃（BMT）系材料は1GHzと10GHzで $\tan\delta$ が倍近く異なることがある。精密な解析をするなら、実測データをテーブルで入力するのがベストだ。

いい質問だ。フィルタ設計では、まず各DRの $f_0$ と $Q$ を個別に求めて、次にDR間の結合係数 $k_{ij}$ を計算する。結合係数は、2つのDRを配置したときの2つの固有周波数 $f_1, f_2$ から次の式で求められる：

この $k_{ij}$ と各DRの $Q$ を結合行列に入れて、フィルタ理論（Chebyshev特性やelliptic特性）に基づいてSパラメータを合成する。最終的にフルウェーブ3D解析で全体のSパラメータを検証するんだ。HFSSのDriven Mode解析やCSTのSパラメータ解析がここで活躍するよ。

業界標準はAnsys HFSSだろうね。FEMベースの固有値ソルバー（Eigenmode Solver）があって、自動適応メッシュ（Adaptive Meshing）で共振周波数を自動的に収束させてくれる。DR解析では Eigenmode 解析で $f_0$ と $Q$ を求めて、その後 Driven Modal 解析でフィルタ全体のSパラメータを検証する、という2段階ワークフローが定番だ。

3つある。(1) 空気領域（Airbox）はDR半径の3〜5倍に設定、(2) 金属壁のQ値も評価したい場合は Finite Conductivity 境界条件を使う（Perfect E だと導体損失がゼロになる）、(3) Adaptive Meshing の収束基準を $\Delta f < 0.01\%$ に設定すること。デフォルトの $\Delta S < 0.02$ はDR解析では粗すぎることがあるよ。

CSTはFIT（有限積分法）やFDTDベースで、時間領域から広帯域のSパラメータを一発で取れるのが強みだ。Eigenmode Solverも搭載していて、Jacobi-Davidson法でDRの共振モードを解ける。特にDRフィルタの広帯域特性を見たいときや、スプリアスモードの全体像を把握したいときに時間領域解析は便利だよ。FITは自動的に $\nabla \cdot \mathbf{B} = 0$ を満たすから、スプリアスモードの心配も少ない。

COMSOLの強みはマルチフィジクス連成だ。DRの自己発熱による温度上昇が $\varepsilon_r$ に与える影響——つまり電磁場-熱連成解析をシームレスにできる。高出力フィルタや基地局のパワーアンプ段のDR設計では、この連成解析がないと $\tau_f$ の影響を正確に評価できないんだ。RF Moduleの Eigenfrequency Study で $f_0$ と $Q$ を求められるよ。

使えるよ。ただし課題もある。ミリ波帯ではDRのサイズが $1\sim2\,\text{mm}$ と極小になるから、加工精度が共振周波数に直結するんだ。寸法公差 $\pm 10\,\mu\text{m}$ で $f_0$ が数百MHz動く。だからCAEでの公差解析（Monte Carloシミュレーション）が設計段階で不可欠になってきている。

テープキャスティングやLTCC（低温同時焼成セラミック）技術で対応している。最近ではセラミック積層造形（3Dプリンティング）による複雑形状DRの研究も進んでいるよ。

かなり活発だ。フルウェーブ3D解析は1ケースあたり数時間かかるから、パラメトリック最適化に何百回も回すのは現実的じゃない。そこでサロゲートモデル（ニューラルネットワーク、ガウス過程回帰など）をFEMデータで学習させて、最適化ループを高速化するアプローチが主流になりつつある。HFSSのDesign of Experiments + Neural Network Modelが代表例だね。

まさにその通り。従来のプレス成形+焼結では作れなかった「空洞入りDR」や「勾配誘電率構造」が積層造形で実現できる可能性がある。トポロジー最適化とCAEを組み合わせて、「最も低損失で最も小型な形状」を自動探索する研究が進行中だ。ただしセラミックの3Dプリントは焼成後の収縮（15〜20%）を正確に予測する必要があるから、まだ研究段階のテーマが多い。

まずは要素タイプを確認しよう。節点要素を使っていたらそれが原因だ。辺要素に切り替えればほぼ解消する。

Q値の不一致はDR解析で最もよくある問題だ。原因を系統的にチェックしよう：

周波数ずれの原因を大きい順に並べるとこうなる：

だからこそ実務では「$\varepsilon_r$ を実測して、その値でシミュレーションを再実行する」というフローが鉄則なんだ。Hakki-Coleman法やCylindrical Cavity法でDRの $\varepsilon_r$ と $\tan\delta$ を実測してからCAEに入力するのがベストプラクティスだよ。