Chapter 13
収束 — その結果を信じてよい条件
この章のゴール.
「収束した」という言葉の 3 つの異なる意味を区別できること。 メッシュ収束の確認手順を、自分の手で実行できるようになること。 そしてエネルギー保存とパッシブ性という、結果を疑うための検算を持つこと。
この章で使う既出の用語（定義は各リンク先）. FDTD（01 章 5 節）、実測（01 章 3 節）、領域（01 章 5 節）、PEC（03 章 4 節）、階段近似（05 章 6 節）、非一様格子（06 章 2 節）、メッシュ収束（09 章 5 節）、急変（10 章 1 節）、PML（11 章 1 節）、ディエンベッド（12 章 5 節）、ポート（12 章 1 節）
1. 「収束」の 3 つの意味
この言葉は混同されやすい。ソルバのログに「Converged」と出ても、それがどの意味かで価値が全く違う。
| # | 種類 | 何が収束したか | 意味 |
|---|---|---|---|
| ① | 反復解法の収束 | 連立方程式の残差が閾値以下になった | 「その離散化された方程式を正しく解けた」だけ |
| ② | 時間発展の収束（FDTD） | エネルギーが十分減衰した | 「計算を打ち切ってよい」だけ |
| ③ | メッシュ収束 | メッシュを細かくしても答えが変わらなくなった | 「離散化誤差が小さい」——これが本命 |
①と③を混同するのが最大の落とし穴である.
①は「1000 セルの粗いメッシュで作った方程式を、機械が厳密に解いた」という報告にすぎない。 粗いメッシュの答えを、極めて正確に計算しただけである。 物理的に正しいかどうかは何も保証していない。
ソルバが「収束しました」と出しても、③を自分で確認していないなら、その結果は信じてはいけない。
2. メッシュ収束の確認手順
唯一の正しい方法は、メッシュを変えて答えが変わらないことを確かめることである。
| 段 | 処理 | ||
|---|---|---|---|
| 1 | 基準のメッシュで解く。評価する量を 1 つ決める（共振周波数、$` | S_{21} | `$ @特定周波数、利得など） |
| 2 | セルサイズを 1.3〜1.5 倍細かくして解き直す | ||
| 3 | 評価量の変化を見る | ||
| 4 | 変化が許容値以下（例: 共振周波数で 0.5 % 未満、$` | S | `$ で 0.02 未満）になるまで 2〜3 を繰り返す |
| 5 | 変化が小さくなった 1 つ前のメッシュを採用する（それ以上細かくする価値がない） |
「最後に細かくした結果」ではなく「変化しなくなった点」を採用する.
収束曲線が平坦になったところが真値の近傍である。 平坦になっていないなら、まだ真値が分からない—— どれだけ「答えらしく」見えても、それは途中経過にすぎない。
外挿という上級テクニック
収束が 2 次（誤差 \(\propto h^2\)、04 章）と分かっていれば、 2 つのメッシュ結果からリチャードソン外挿で真値を推定できる:
粗い 2 点から、細かいメッシュで解いた場合の答えを予測できる。 外挿値と実際の細メッシュ結果が一致すれば、収束次数が想定通りであることの確認にもなる （一致しなければ、階段近似などで次数が落ちている疑いがある——06 章）。
メッシュ収束を確認する
メッシュを段階的に細かくしたときの評価量の変化を見て、 「どこで打ち切るか」「まだ収束していないのはどれか」を判断する練習ができる。
3. FDTD の時間発展の打ち切り
時間領域ではいつ計算をやめるかが問題になる。
判定は残留エネルギーで行う。ソルバは領域内の全エネルギー
を監視し、ピーク値から十分下がったら終了する。
| 打ち切り基準 | 用途 |
|---|---|
| −30 dB（0.1 %） | 粗い評価、傾向を見るだけ |
| −40 dB | 標準的な設定 |
| −50 〜 −60 dB | 高精度、高 Q 構造、低いサイドローブ（アンテナが主方向以外へ漏らす放射、15 章）や深いノッチを見たいとき |
早く切りすぎると何が起きるか——「打ち切り誤差」.
時間波形をエネルギーが残ったまま打ち切ると、FFT（16 章）では波形を突然ゼロにしたことになる。 これは矩形窓を掛けたのと同じで、スペクトルにリップル（ギブス現象）が乗る。
症状は特徴的で、S パラメータが本来滑らかなはずの帯域で細かく波打つ。 「S パラメータがギザギザする」ときは、まずこれを疑う。 対策は「もっと長く計算する」が正道で、 窓関数を掛けて誤魔化す方法もあるが分解能が落ちる（DSP シリーズ 04 章と同じトレードオフ）。
4. 物理法則による検算
数値解が物理的に破綻していないかは、保存則で確かめられる。これは非常に強力な検算である。
① エネルギー保存
S パラメータで書けば、無損失・無放射の 2 ポート回路なら
これを超えたら明確な異常である（後述のパッシブ性）。 損失や放射があれば 1 より小さくなるが、その差が「どこへ行ったか」を説明できなければならない。
② パッシブ性（受動性）
電源を含まない構造では、出ていく電力が入る電力を超えることはあり得ない:
1 を超えたら計算が間違っている。（ただしポートの参照インピーダンスが実数でない場合など、 定義次第で見かけ上超えることはある。まず定義を確認する。）
③ 相反性
非相反素子（磁化フェライトなど、03 章）を含まなければ
\(S_{21} \ne S_{12}\) になったら、メッシュが非対称かポート設定がおかしい。 対称構造なのに \(S_{11} \ne S_{22}\) になる場合も同様である。
④ 因果性
時間領域応答が、入力が入る前に反応していないこと。 周波数領域の結果から時間領域に変換したとき（16 章）、 \(t<0\) に応答が現れるなら、周波数点が足りないか帯域が不足している。
これらの検算は 1 分で終わり、致命的な誤りを高い確率で捕まえる.
S パラメータを出したら、まずパッシブ性と相反性を見る習慣をつけたい。 「\(|S_{11}|^2+|S_{21}|^2 = 1.4\)」のような結果は、 どれだけ綺麗なグラフに見えても、ポート定義かメッシュに深刻な問題があることを意味する。
5. 収束しない・おかしいときの切り分け
| 症状 | 疑うもの |
|---|---|
| FDTD が発散する（NaN＝Not a Number、数値として意味をなさない値。無限大も同様） | クーラン数の設定（06 章）／不安定な材料モデル／PML の実装／非一様格子の急変 |
| エネルギーが減らない | 高 Q 共振（正常。長時間必要）／PML が効いていない／PEC で囲まれた閉領域になっている |
| 反復解法が収束しない | 要素品質が悪い（10 章）／前処理の選択／周波数が共振点に極めて近い |
| S パラメータがギザギザ | 時間打ち切りが早い（本章 3 節）／掃引の補間が外れている（07 章） |
| メッシュを細かくすると答えが動き続ける | まだ収束していない／特異点（無限に細かくしても収束しないエッジ）／モデルの不整合 |
| パッシブ性を破る | ポート定義／参照インピーダンス／ディエンベッドの誤り（12 章） |
| 共振周波数が実測とずれる | 材料定数（03 章）／メッシュ不足／給電のインダクタンス（12 章）／製造公差 |
6. 信頼できる結果の条件 — チェックリスト
- ☐ メッシュ収束を確認した（評価量を決めて、細分化しても変わらないこと）
- ☐ 時間領域ならエネルギーが −40 dB 以上減衰するまで計算した
- ☐ パッシブ性（\(\sum|S|^2 \le 1\)）を満たしている
- ☐ 相反性（\(S_{ij}=S_{ji}\)）を満たしている（非相反素子がない場合）
- ☐ 対称な構造なら、結果も対称である
- ☐ 境界までの距離を変えても結果が変わらない（11 章）
- ☐ ポートのサイズ・位置を変えても結果が変わらない（12 章）
- ☐ 既知の解（解析解、実測、文献値）と桁が合っている
- ☐ できれば別手法でも解いた（09 章）
上 4 つは必ず、残りは重要な結果なら——これが実務の基準になる。
7. この章のまとめ
| ポイント | 内容 |
|---|---|
| 3 つの収束 | 反復解法・時間発展・メッシュ。本命は 3 番目 |
| メッシュ収束 | 評価量を決めて細分化し、変化しなくなった点を採る |
| 外挿 | 2 次収束を仮定して真値を推定。収束次数の確認にもなる |
| FDTD の打ち切り | 残留エネルギー −40 dB が標準。早いとS パラメータがギザギザになる |
| 検算 | パッシブ性・相反性・エネルギー保存・因果性。1 分で致命傷を見つける |