PML — CAE用語解説

良い質問だ。PMLは「Perfectly Matched Layer」の略で、計算領域の境界で発生する「非物理的な反射」と「計算領域内の物理的な波動」とを整合させる、つまり反射をゼロに近づけるための人工的な層だ。例えば、アンテナから放射された電磁波が解析領域の端でいきなり壁にぶつかると、現実には存在しない反射波が発生して結果が汚染される。PMLはその反射を吸収するスポンジのような役割をする。

その通り、根本的に異なる。通常の減衰材ではインピーダンスが不整合で反射が生じる。PMLの核心は、複素数座標系への「座標変換」にある。波動方程式において、実数座標を複素数に拡張することで、波動がPML内を伝搬するにつれて指数関数的に減衰する解を得る。この変換が「整合」を生み出す。例えば、1次元の波動方程式で、座標xを

実装レベルでは、複素数座標変換の効果を等価な「複素誘電率・透磁率」または「複素屈折率」を持つ媒質としてモデル化する。つまり、PML層には「損失」と「分散」の両方の特性を持たせた仮想的な材料を定義している。Ansys HFSSでは「PML境界条件」として、ユーザーが層の厚さ（例えば、中心波長の10%から20%）と減衰プロファイル（多項式や幾何級数）を設定する。背後でこの等価媒質のパラメータが自動計算されている。

必須だ。PMLは計算領域の外側に「層」として追加される独立した領域だから、そこにも当然メッシュが必要になる。FEMでは、PML領域の要素に先ほどの複素材料特性を割り当てる。重要なのは、PML内のメッシュ品質を高く保つこと。特に減衰が急峻に変化するプロファイルを使う場合、粗いメッシュでは減衰が正しく表現されず、反射低減性能が劣化する。経験則として、PML層の厚さ方向に少なくとも3〜5要素は欲しい。

良い出発点はある。厚さは扱う波長λに対して決める。電磁波解析ではλ/4からλが一般的だ。減衰係数σの最大値は、PML出口での波動の振幅が十分小さくなるように設定する。目安として、PML内を1波長伝搬した時の減衰量が60dB（振幅で1/1000）以上になるようにσ_maxを選ぶ。COMSOL Multiphysicsの「吸収境界条件」では、多項式次数やスケーリング係数をユーザーが調整できるが、多くの場合デフォルト設定（例えば、多項式次数3、理論反射率-60dB）で十分だ。

鋭い指摘だ。その通り、PML領域の要素は複素数の材料定数を持つため、全体の係数行列も複素数になる。実数行列と比べて、記憶容量は約2倍、計算時間も増加する。しかし、代替手段である「吸収境界条件（ABC）」は反射低減性能が劣る（特に斜入射や低周波数）。トレードオフだ。大規模問題では、PMLを必要な方向だけに適用する（例えば放射問題では放射方向のみ）ことでコストを抑制する。

まず、PMLが本当に必要かだ。閉じたキャビティや導波管のように、境界が完全導体で物理的に反射がある場合はPML不要だ。PMLが必要なのは「開放境界」を模擬する時、つまりアンテナ放射、散乱断面積、音響放射などの問題だ。次に、PMLを置く位置。少なくとも光源や散乱体からλ/2以上離す。近すぎると近接場のエバネッセント波をPMLが正しく扱えず、結果が不正確になる。

二つの現実的な方法がある。第一は「参照解との比較」だ。PMLの代わりに非常に大きな計算領域（反射が無視できるほど）を使って解を求め、それを正解としてPMLを使った小領域の解と比較する。第二は「残留反射のモニタリング」だ。Ansys HFSSなら「フィールドオーバーライド」でPML手前の面を定義し、その面上での後方への電力流（ポインティングベクトル）を計算し、入射電力に対する比を見る。実用的には-50dB以下を目指す。

主に三つだ。1) 共振周波数がずれる：閉じた系ではないのに、特定周波数で鋭い共振ピークが現れる。これは境界での反射による定在波の証拠。2) 放射パターンが歪む：アンテナの指向性が、理論や測定と比べて左右非対称になる。3) 収束が悪い：周波数掃引でSパラメータがギザギザしたり、メッシュを細かくしても結果が安定しない。これらは全て、PMLの反射低減性能が不十分で、計算領域内で波動が「閉じ込められて」いる可能性を示唆する。

コンセプトは同じだが、インターフェースと適用範囲に違いがある。Ansys HFSSでは、3DオブジェクトとしてPMLを定義する（「Create」→「Boundary」→「PML」）。一方、COMSOLでは「吸収境界条件」の一種として「PML」を選択し、それを適用する「ドメイン」を指定する。CST Studio Suiteでは「境界条件」の設定で「PML」を選択し、層の数やプロファイルを細かく制御できる。また、CSTはFIT（有限積分法）ベースなので、PMLの実装がFEMベースのソフトと内部的に異なる点には注意だ。

使われる。音響（圧力音場）や弾性波（固体中の振動）の放射・散乱問題だ。COMSOLの「音響モジュール」や「構造力学モジュール」には専用の「PML」機能がある。ただし、支配方程式が異なるので、内部で定義される等価な複素材料特性は電磁気の場合と違う。例えば音響では、複素数密度と複素数体積弾性率が定義される。Ansys Mechanicalでも「無限要素」が同様の目的で使われるが、数学的基礎はPMLとは異なる場合がある。

商用ソフトに比べて難易度は高い。ユーザー自身がPML領域を定義し、その領域に適切な複素材料特性や座標変換に基づく弱形式を実装する必要がある。例えば、FEniCSを使った電磁波シミュレーションでは、弱形式にPMLの複素座標変換のヤコビアンを明示的に組み込むコードを書く。OpenFOAMの音響ソルバー「acousticFoam」にはPML機能が実装されているものの、その設定ファイル（`fvOptions`）を正しく記述する知識が必要だ。教育目的やカスタマイズ性が求められる研究では強力だが、実務での導入障壁は高い。

まず疑うべきは「数値的不安定性」だ。PMLの減衰係数σ_maxを極端に大きく設定しすぎると、PML層内で波動が急激に減衰する。この急峻な変化をメッシュが追従できず、離散化誤差が大きくなり、ソルバーが発散することがある。対策は、σ_maxを小さくする（例えばデフォルト値から始める）、PMLの厚さを増やす（同じ減衰量をより緩やかな勾配で達成する）、そしてPML層内のメッシュを細かくする、の3点だ。Ansysのログに「matrix is ill-conditioned」といった警告が出ていないか確認しよう。

その通り、PMLは波長に対して厚さを設計するので、低周波数では相対的にPMLが薄くなり、吸収が不十分になる。具体的な対策は二つ。第一に「PMLの厚さを波長比ではなく絶対長で確保する」。例えば、最低周波数1MHz（空気中で波長300m）でも現実的な厚さ0.5mを確保する。第二に「周波数依存の減衰係数プロファイルを使う」。COMSOLなどでは、減衰係数σを周波数に比例させる（

注意が必要だ。古典的なPMLの理論は等方性の線形媒質を前提としている。異方性材料に面したPMLは、理論的に完全な整合が達成できず、多少の反射が生じる。商用ソフト（Ansys, COMSOL）は内部で対応を試みているが、性能は等方性の場合より劣る可能性がある。非線形媒質が主役の問題（例えば、高強度光による非線形光学効果）では、PMLの適用は非常に難しく、通常は使わない。その場合は計算領域を物理的に大きく取るか、他の開放境界手法を検討することになる。実装前にソフトウェアのマニュアルで対応状況を必ず確認すべきだ。