メッシュ収束・Richardson外挿・GCIとは
🙋
メッシュ収束って何ですか？CAEでメッシュを細かくすればするほど結果が変わるって聞いたんですけど、どこまで細かくすればいいんですか？
🎓
大まかに言うと、メッシュを細かくしても答えがほとんど変わらなくなった状態を「収束した」って言うんだ。例えば自動車のボンネットの変形を計算する時、粗いメッシュだと変形量が10mm、細かくすると12mm、もっと細かくしても12.1mm…となれば、12mm付近で収束していると考えられるよ。このシミュレーターの「次元」パラメータを1D, 2D, 3Dに変えてみると、必要なメッシュ数の増え方が大きく異なるのが体感できる。
🙋
え、そうなんですか！でも、毎回非常に細かいメッシュで計算するのは時間がかかりすぎますよね。細かくする前にある程度答えを予測できないんですか？
🎓
そこで使うのがRichardson外挿法だ。3段階の異なる粗さのメッシュで計算した結果から、理論上の「無限に細かいメッシュの答え」を推定する方法なんだ。このツールで、粗・中・細の3つのメッシュで得た「注目する値」（例えば最大応力）を入力してみて。すると、推定された「真の解」が計算されるよ。実務では、この推定値を使って細かいメッシュ計算を省略できるケースもある。
🙋
なるほど！でも、その推定がどれくらい信用できるか、判断する指標はあるんですか？「収束してます」って言うための数字的な基準が欲しいです。
🎓
そのためにGCI（格子収束指数）がある。これは「推定値と現在の解の間の相対誤差が、メッシュをさらに細かくしたら何%改善される見込みか」を示す指標だ。このツールでは、計算されたGCIの値が2%未満なら緑、5%未満なら橙、それ以上なら赤で表示される。現場で多いのは、GCI_fineが3%を切ることを目標にメッシュを決めることだね。パラメータを動かして、GCIの色が変わる様子を確認してみよう。
物理モデルと主要な数式
メッシュサイズ$h$と数値解$\phi$の関係は、理論的な精度次数$p$を用いて以下のようにモデル化されます。これがRichardson外挿の基礎です。
$$\phi = \phi_{exact}+ g h^p + \cdots$$
ここで、$\phi_{exact}$は真の解（連続体の解）、$g$は定数、$h$は代表メッシュサイズです。$p$は離散化スキームの理論的精度次数（例えば中心差分なら2次）ですが、実際の問題では複雑なため、3つの解から計算される「見かけの精度次数」が使われます。
3つの異なるメッシュ（細：1, 中：2, 粗：3）から、見かけの精度次数$p$と外挿解$\phi_{ext}^{21}$、GCIを計算します。
$$ p = \frac{1}{\ln(r_{21})}\left| \ln \left| \frac{\epsilon_{32}}{\epsilon_{21}}\right| + q(p) \right| $$
$$ \phi_{ext}^{21}= \phi_1 + \frac{\phi_1 - \phi_2}{r_{21}^p - 1}$$
$$ GCI_{fine}^{21}= F_s \frac{|\epsilon_{21}|}{r_{21}^p - 1}$$
$r_{21}=h_2/h_1$は細化比、$\epsilon_{21}=\phi_2-\phi_1$は解の差、$F_s$は安全係数（通常1.25）です。$q(p)$は$p$が未知の場合の補正項。ツールはこれらの式に基づき、入力された3つの解から自動的に$p$, 外挿解, GCIを計算し、収束性を評価します。
よくある質問
粗・中・細の3段階以上で、メッシュサイズの比（r）を一定（例：r=2）に設定するのが基本です。例えば、粗メッシュの要素数をN、中を2N、細を4Nとします。比が一定でないと、見かけの精度次数pの計算精度が低下します。
いいえ。外挿解は、解が漸近領域（メッシュが十分細かい領域）にある場合にのみ信頼できます。GCIが小さく、見かけの精度次数pが理論値に近いことを確認してください。漸近領域外では外挿解が真の解から大きく乖離する可能性があります。
一般的な目安として、GCIが5%未満であれば解はメッシュに依存していないと見なせます。ただし、解析の目的や要求精度によります。また、粗・中・細の各メッシュ間でGCIが一貫して減少していることも重要です。
主な原因として、(1)メッシュが粗すぎて漸近領域に達していない、(2)メッシュ品質が悪い（歪みやアスペクト比不良）、(3)境界条件や物理モデルに不連続点がある、(4)数値スキームが理論通りに動作していない、などが考えられます。まずはメッシュをさらに細かくして確認してください。
実世界での応用
構造強度解析（自動車・航空機）：部品の最大応力や変形量を評価する際、メッシュ依存性を排除した信頼できる結果を得るために必須の手順です。特に安全率が小さい軽量化設計では、GCIを用いた定量的な収束確認が設計判断の根拠となります。
流体解析（CFD）：翼型の揚力・抗力係数や、熱交換器の伝熱量・圧力損失の予測精度を高めるために広く利用されます。複雑な乱流モデルでは理論次数$p$が不明なため、計算される見かけの次数から収束性を判断します。
電磁界解析：アンテナの放射パターンや共振周波数、電子デバイスの寄生容量など、メッシュの切り方に敏感な量を評価する際に、Richardson外挿によって計算コストを抑えつつ高精度な解を推定します。
学術研究・論文執筆：数値計算結果の信頼性を示すための重要な検証項目として、ほぼ必須の手法です。再現性のある研究を行うため、使用したメッシュサイズとともにGCIの値を報告することが推奨されています。
よくある誤解と注意点
まず、「メッシュを細かくすれば必ず正解に近づく」とは限らないという点を押さえましょう。例えば、モデルの幾何学的な特異点（鋭い角や接触点）付近では、メッシュを細かくしても応力が発散し続け、「収束」しないことがあります。この場合は、理論的に応力が無限大になるので、CAEの目的を「最大応力値の特定」から「その周辺の応力分布の把握」に切り替えるなど、解釈を変える必要があります。
次に、3つのメッシュは「系統的に」細かくすることが鉄則です。例えば、全体の要素サイズを単純に1mm、0.7mm、0.5mmとするのではなく、細化比 $r$ を一定（推奨は1.3以上）に保ち、1mm、0.77mm（1/1.3）、0.59mm（1/1.3^2）のように設定します。でないと、見かけの精度次数 $p$ の計算が不安定になり、外挿結果が信用できなくなります。
また、GCIの値だけを盲信しないでください。GCI_fineが2%を切っていても、注目している物理量が「たまたま」収束しているだけかもしれません。例えば変形量は収束していても、最大応力はまだ収束していない、というケースはよくあります。必ず複数の重要な出力変数について個別に収束性を確認しましょう。