🙋
CFDの解析結果を見て「メッシュ依存性が…」って先輩がよく言うんですけど、結局メッシュを何分割すれば信頼できるんですか？感覚で決めてるんですかね？
🎓
いい質問だね。実は感覚で決めるべきじゃなくて、ちゃんと標準化された手法があるんだ。それが「GCI（Grid Convergence Index）」と呼ばれるもので、ASME V&V 20 という国際標準で決まっている。同じ問題を細・中・粗の3つの異なるメッシュで解いて、その差を見れば「今の細メッシュ解が真値からどれくらいズレている可能性があるか」をパーセントで出せるんだ。CFDの論文でも、ジャーナルに投稿するなら GCI の表示はほぼ必須になっているよ。
🙋
なるほど、3メッシュ必要なんですね。でもデフォルト値（f₁=100, f₂=105, f₃=115）でやると、細メッシュ GCI が 6.25% って出てます。これって「真値は100±6.25くらい」って意味ですか？
🎓
そう、まさにその解釈で合ってる。GCI 6.25% は「細メッシュ解 f₁=100 の不確かさが ±6.25% 程度」という意味で、保守的に「真値は 93.75〜106.25 の間にある可能性が高い」と読む。同時に Richardson 補外値 f_ext=95 も出ているけど、これは「メッシュを無限に細かくしたときの推定値」だ。100 と 95 の差が、まさに離散化誤差そのものだね。GCI はその誤差を 1.25 倍の安全係数で膨らませた、より保守的な上限と思えばいい。
🙋
「漸近収束チェック」っていう値も出てますね。これは何を見てるんですか？
🎓
これは「Richardson 補外を信用していいか」のセルフチェックなんだ。理論上、漸近収束領域にいれば GCI32/(r^p·GCI21) はちょうど 1 になるべき。デフォルト設定では 0.952 が出ていて、これは「ほぼ1なので漸近域内、信頼できる」という意味。逆にこれが 0.5 や 2.0 みたいに大きく外れていると、「まだメッシュが粗すぎて理論通りに収束していない」というサイン。その場合はもう一段細いメッシュで計算をやり直す必要があるんだ。
🙋
観測収束次数 p ってのが 1.000 になってます。これは2次精度スキームを使ったら 2 になるはずですよね？
🎓
鋭い！デフォルトの例だと、解の差が（f₁→f₂で5, f₂→f₃で10）と完全に2倍ずつ変化しているから、p = ln(2)/ln(2) = 1 になる。これは「実効的に1次精度で収束している」ということ。理論的に2次精度スキームを使っていても、観測 p が1になる原因はいくつかあって、(1) 境界層メッシュが粗すぎて1次精度に落ちている、(2) 限定子（limiter）が頻繁に発動している、(3) まだ漸近域に入っていない、などが典型例。実務で2次スキームを使っていて p=1 が観測されたら、まずメッシュをさらに細かくして p が 2 に近づくか確認するのが定石だよ。
🙋
リファインメント比 r って、どれくらいに設定すればいいんですか？2 で固定でもいいんでしょうか？
🎓
r = 2 は最も使われる値で、構造格子なら各方向のメッシュを半分にするだけだから実装も楽。ただし ASME V&V 20 では「r ≥ 1.3 が推奨」とあって、必ずしも 2 でなくていい。非構造格子だと体積比で 8 倍（線形寸法で 2 倍）にするのが大変なので、r = 1.5 程度で3メッシュを用意するケースも多い。逆に r を大きくしすぎる（4 以上）と、粗メッシュ解がそもそも実用的な精度に届かず、解の傾向が変わってしまうことがあるので注意。下の「GCI(%) vs r」グラフを動かしてみると、r による GCI の変化が見えるよ。
GCI (Grid Convergence Index) は、Roache (1994) が提案した離散化誤差の保守的な指標で、ASME V&V 20-2009 で標準化されています。3つの異なるメッシュ（細・中・粗）で同じ問題を解き、その解の差から「現在の細メッシュ解が真値からどれくらいズレている可能性があるか」をパーセントで表します。たとえば GCI21 = 6.25% なら、細メッシュ解の不確かさが ±6.25% 程度であることを意味します。CFDや有限要素解析の検証 (Verification) で必須のメトリクスです。
F_s = 1.25 は Roache (1998) と ASME V&V 20 で推奨されている、3つのメッシュを使う場合の標準値です。Richardson補外による誤差推定は理想的な漸近域では正確ですが、実際の CFD/FEM では完全な漸近域にいるとは限らないため、推定誤差を 1.25 倍に膨らませて保守的に評価します。2メッシュしか使えない場合は F_s = 3.0 が推奨されます。F_s を 1.0 にすると Richardson 誤差そのものになり、保守性が失われます。
漸近収束チェック GCI32 / (r^p · GCI21) は理論上 1.0 になるべき値です。0.9〜1.1 程度であれば「漸近収束領域内」とみなし、Richardson補外が信頼できると判断します。1.0 から大きく外れる場合、メッシュがまだ十分細かくない、解が単調収束していない（振動収束）、あるいは観測収束次数 p が理論次数と大きく異なる、といった理由が考えられます。この場合はさらにメッシュを細分化するか、ソルバー設定（離散化スキーム・境界条件）を見直す必要があります。
観測収束次数 p が理論次数（例えば 2次精度スキームなら p=2）と一致しないのはよくあります。原因は (1) メッシュがまだ漸近域に入っていない、(2) 境界条件や境界層・衝撃波付近で局所的に精度が落ちている、(3) 異方性メッシュやハイブリッドメッシュで実効精度が下がっている、(4) 限定子 (limiter) が働いて1次精度に落ちている、などです。p が 0.5 未満や 5 を超えると Richardson 補外は信頼できないため、メッシュ系列を見直してください。理想的には p = 1.0〜2.5 程度が健全な範囲です。
航空機・自動車空力 CFD：翼や車体まわりの抗力係数・揚力係数を求める CFD では、メッシュ依存性を排除した値を提示することが必須です。AIAA Drag Prediction Workshop や自動車メーカーの空力解析仕様書では、GCI 5% 以下を目標値とすることが多く、3〜4 メッシュレベルで計算して GCI を報告することが標準的なフローになっています。
ターボ機械・燃焼器の性能予測：蒸気タービンやガスタービン・ロケットエンジン燃焼器の効率や燃焼温度予測では、剥離点・衝撃波・反応帯の解像度が結果を大きく左右します。GCI を用いれば「効率の見積もり 92.3% に対して不確かさ ±1.5pt」のように、メッシュに由来する不確かさを明示でき、設計判断や認証申請に使えます。
原子力安全解析：NRC（米国原子力規制委員会）や ASME V&V 30 では、安全裕度の証明に GCI 評価を要求します。冷却材喪失事故 (LOCA) 時の最高被覆温度予測などで、解析の不確かさを定量化し、規制値からの余裕（マージン）を明示するために使われます。Richardson 補外で推定した「真値」と、安全係数を加味した GCI 上限値を併記するのが標準です。
構造解析（FEM）の応力集中検証：切欠き・溶接止端・ボルト穴まわりの応力集中点では、メッシュサイズで応力値が大きく変動します。最大ミーゼス応力に対して GCI を評価し、「メッシュ依存性が 3% 以下に収まった上での最大応力 320 MPa」と報告することで、過小評価・過大評価のリスクを抑えられます。特に疲労寿命予測のように応力値の精度がそのまま寿命予測に効く解析では必須の手順です。
まず最大の落とし穴が、「漸近収束領域にいないのに GCI を信用してしまう」ことです。GCI の理論は Richardson 補外を前提にしており、それが成り立つのは解が漸近的に r^p のオーダーで収束しているとき、つまり「漸近収束領域 (asymptotic range)」に入っている時だけ。漸近収束チェック GCI32/(r^p·GCI21) が 1 から大きく離れていたり、観測収束次数 p が極端な値（0.5 未満や 5 超）になっている場合、メッシュがまだ粗すぎて補外が信頼できません。この状態で「GCI=10%だから問題ない」と判断すると、実際には 30〜50% ズレていることもあり得ます。漸近域に入るまでメッシュを細かくするのが大前提です。
次に、「振動収束を単調収束と混同する」こと。GCI の式は f₁, f₂, f₃ が同じ方向（単調に増加または減少）で収束していることを暗黙に仮定しています。実際には乱流モデルの不安定性や反復計算の収束不足で、メッシュ細分化に対して解が振動することがあります。例えば f₁=100, f₂=98, f₃=102 のように真値の周りで振動している場合、|f₃-f₂|/|f₂-f₁| が異常な値になり、観測 p も意味を持ちません。プロットして振動の有無を確認し、振動収束の場合は「不確かさ ≥ max(|fi - fj|)」のように別の保守的見積もりに切り替えるべきです。本ツールでも、収束判定が「非漸近域」になる時は要警戒です。
最後に、「グローバル指標と局所指標を混同する」こと。GCI は基本的に積分量や代表点での値（抗力係数、最大温度など）に対して評価するものです。場の各点 (cell-by-cell) で GCI を出すと、ノイズや局所的な非単調性で異常値が大量に出ます。論文や報告書で GCI を提示する際は、「揚抗比 L/D の GCI が 3%」「最大ミーゼス応力 (point A) の GCI が 5%」のように、明確に評価対象の物理量を限定してください。また、複数の評価量がある場合（揚力・抗力・モーメントなど）、それぞれ独立に GCI を出して、最も保守的な値を採用するのが標準的なやり方です。