マクスウェル方程式 — CAE用語解説
マクスウェル方程式
先生、マクスウェル方程式って電磁気学の全てですよね?CAEとの関係を教えてください。
理論と物理
基本概念と支配方程式
マクスウェル方程式って、結局何を表しているんですか?電磁気学の全てをまとめたもの、というのは聞くのですが、具体的にどういう関係を記述しているのかイメージが湧きません。
良い質問だ。4つの方程式は、電場
「アンペール・マクスウェルの法則」には「変位電流」という項が追加されていますよね。なぜ電流ではない「変位電流」が必要なんですか?
それがマクスウェルの最大の洞察だ。従来のアンペールの法則だけでは、コンデンサが充電される際の回路の磁場を説明できない。導線中の伝導電流は途切れているが、コンデンサ極板間の電場の時間変化
支配方程式は積分形と微分形がありますが、CAEソフトウェアではどちらを使って計算しているんですか?
数値計算の基礎となるのはほぼ微分形だ。例えば、有限要素法(FEM)で解くのは、微分形のマクスウェル方程式から導出される弱形式だ。具体的には、波動方程式
数値解法と実装
FEMによる離散化とソルバー
電磁界解析で「エッジ要素」を使うと聞きました。通常のFEMの節点要素と何が違うのですか?
本質的な違いだ。節点要素は節点でのスカラー値(例えば電位)を定義するが、電場や磁場のようなベクトル場には不向きな場合がある。エッジ要素は、要素のエッジに沿ったベクトル場の接線成分を未知数とする。これにより、材料界面での法線成分の不連続性を自然に表現でき、また発散のない場(磁束密度
周波数領域と時間領域の解法(FEMとFDTD)は、どうやって使い分けるんですか?
用途と求解対象による。FEM(周波数領域)は、共振器の固有モードや特定周波数(例:5Gの28GHz帯)での定常応答を求めるのに適している。一方、FDTD(時間領域)は、広帯域の過渡応答、例えば雷サージがケーブルを伝わる様子や、パルスレーダーの波形の伝搬を解くのに向く。メーカーで言えば、Ansys HFSSはFEM、CST Studio SuiteはFDTDが強みだ。ただし、HFSSにも時間領域ソルバー「Transient」は存在する。
「PML(完全整合層)」という境界条件をよく見ます。これは何をしているんですか?
PMLは、計算領域の端で電磁波を反射せずに吸収する(減衰させる)仮想的な層だ。無限遠方まで広がる場を有限の計算領域で扱うための必須技術だ。原理は、座標を複素数に拡張し、波動が層に入ると指数関数的に減衰するようにする。実務では、使用周波数に応じて層の厚さを調整する。例えば10GHzの解析では、波長の約4分の1、つまり7.5mm程度のPML厚さが目安になる。設定を誤ると反射率が悪化し、Sパラメータの結果が信頼できなくなる。
実践ガイド
ワークフローとチェックリスト
アンテナのSパラメータをシミュレーションするとき、最初に何を確認すべきですか?
まず「励振」の設定だ。給電点(ポート)のモードの定義が正しいか。例えば、同軸コネクタを使う場合は、TEMモードのインピーダンスを50Ωに設定する。次に「メッシュ」。波長に対して十分細かいか。経験則として、最高周波数の波長の1/10以下のメッシュサイズが必要だ。2.4GHzなら空気中での波長は約125mmなので、少なくとも12.5mm以下の要素サイズを目標にする。Ansys HFSSでは「Lambda Refinement」設定でこれを自動化できる。
材料特性で「比誘電率」と「誘電正接」を入力しますが、周波数によって値が変わると聞きました。どう扱えばいいですか?
その通りで、これが実務上の大きな落とし穴だ。FR-4基板の比誘電率は、1MHzでは約4.7だが、10GHzでは4.3程度に低下する。シミュレーションの中心周波数に対応したデータシートの値を使うこと。広帯域解析の場合は、デバイバー(周波数分散)モデルが必要だ。多くのソフトウェアは「Debye」や「Djordjevic-Sarkar」モデルを用意している。メーカーの測定データがなければ、IPC-2141Aなどの規格で定義されている代表値を使うこともあるが、高周波設計ではそれでは不十分だ。
結果の収束をどう判断すればいいですか?「Adaptive Meshing」が終わったらそれで終了?
適応メッシュ分割は第一歩に過ぎない。収束の判断は、求めたい物理量の変化を見る。例えばSパラメータなら、メッシュステップ(パス)ごとのS11の値の変化が、設定した収束条件(デフォルトで0.02以下)を満たしているか。さらに、エネルギー誤差(Ansysでは「Delta S」)が十分小さいか(例えば0.01以下)を確認する。最終的には、メッシュ数を2倍にした「収束解析」を別途行い、結果が変わらないことを確認するのが確実な方法だ。
ソフトウェア比較
Ansys / Abaqus / COMSOL等
Ansys HFSSとCOMSOLの「RFモジュール」は、どちらも電磁波を解くと思いますが、根本的な違いは何ですか?
アーキテクチャと強みの分野が異なる。HFSSは「エッジ要素FEM」と「適応型メッシュ分割」を核とした、高周波・アンテナ・アイソレーション解析の業界標準だ。特にSパラメータの精度と信頼性で評価が高い。一方、COMSOL Multiphysicsは「弱形式」をユーザーが直接編集できる柔軟性が最大の特徴で、電磁界と熱、構造、流体などの連成解析が容易だ。例えば、誘導加熱(電磁界+熱伝導)やMEMS(静電+構造)を一つの環境でシームレスに解ける。単体の電磁波解析速度はHFSSに分があるが、マルチフィジックスではCOMSOLが有利だ。
Abaqusにも電磁気解析機能はあるのですか?
Abaqus/Standardには「Low-frequency Electromagnetic」分析機能がある。これはマクスウェル方程式のうち、変位電流項
オープンソースの電磁界解析ソルバーは使えるものがありますか?
研究や学習用途では有力な選択肢がある。「Elmer FEM」は熱流体構造と電磁気を解けるマルチフィジックスソルバーだ。「OpenFOAM」にも電磁気拡張モジュール(foam-extendのelectromagnetics)がある。ただし、商用ソフトのような高度なメッシャーやGUI、技術サポートはない。また、産業界で事実上の標準である「VHDL-AMS」や「IBIS」モデルとの互換性も限られる。実務で使うには、自前で前処理・後処理を整備するか、有償のフロントエンド(例えばSalome Platform)と組み合わせる必要がある。
トラブルシューティング
よくあるエラーと対策
シミュレーション結果のS11が現実の測定値と比べて全く異なります。どこから調査すべきですか?
まず「ポートの校正面」を疑え。ソフトウェア上でのポート位置(Reference Plane)と、ベクトルネットワークアナライザ(VNA)で校正した位置が一致しているか。次に「損失」。現実の導体には表面粗さによる損失がある。例えば、HFSSで銅の導電率を5.8e7 S/m(理想値)のまま使っていないか?実測値を入力するか、「Surface Roughness」モデルを適用する。最後に「周囲環境」。シミュレーションは理想的なフリースペースだが、実測ではアンテナ測定用チャンバーを使っていない限り、床や壁からの反射が影響する。
メッシュを細かくすると「メモリ不足」エラーで計算が止まります。どう対処すればいいですか?
まず、本当にその細かさが必要か見直す。局所的なメッシュ制御を使いすぎていないか。例えば、エッジや面に対して一律に細かいメッシュを適用するのではなく、電場や磁場が集中する領域(導体エッジ、誘電体界面、給電点近傍)に限定する。次にソルバー設定。Ansys HFSSの「Direct Solver」はメモリ消費が大きい。「Iterative Solver」に切り替えるとメモリ使用量が大幅に削減できる場合があるが、収束性は問題によって異なる。根本的には、64GBや128GB以上の大容量RAMを搭載したワークステーションの導入を検討せざるを得ない。
「共振周波数」のシミュレーション値がデータシートと数%ずれます。この誤差は許容範囲ですか?
用途による。一般のフィルタやアンテナ設計では1%以内を目標とする。数%の誤差は無視できない。原因は主に三つだ。1. 材料定数:比誘電率の実測値とシミュレーション入力値の差。2. 構造寸法:データシートの図面は概略値であることが多く、特に曲げ部分の半径や導体厚みが不正確。3. メッシュとソルバー誤差。対策としては、既知の共振周波数を持つ簡単な構造(例えば、半波長ダイポールアンテナ)でソフトウェアの設定を「校正」し、誤差の傾向を把握しておくことが有効だ。その後、実デバイスの解析では、その誤差分を考慮して設計マージンを見積もる。
放射境界条件(Radiation Boundary)を設定したのに、領域端で場がゼロになっていません。これは正常ですか?
「放射境界」は完全な吸収体ではないので、ある程度の反射は残る。特に、境界が放射源に近すぎる場合や、境界面に対して電磁波が浅い角度(Grazing Angle)で入射する場合に反射が大きくなる。正常かどうかは、境界での反射率をモニターするか、境界までの距離を1.5倍に広げて結果が変わらないか確認する。HFSSの「Radiation Boundary」の理論的反射率は、垂直入射で約1%程度だ。より高性能なPMLを使うか、境界までの距離を最低でもλ/4(中心周波数の波長の1/4)以上確保するのが原則だ。
関連トピック
なった
詳しく
報告