メッシュ独立性検証シミュレーター

Q: メッシュ独立性検証（Mesh Independence Study）とは何ですか？

CFDやFEMでは、メッシュを段階的に細かくしていくと解（応力、温度、抗力など）が真値に向かって収束していきます。隣り合う段階の解の相対変化が許容値（典型的には1%）以下になった時点で「これ以上細かくしても結果が変わらない」と判断し、その段階をメッシュ独立達成段階と呼びます。この検証なしに発表する解析結果は、結論がメッシュサイズに依存している可能性があり、V&V（検証と妥当性確認）の観点から不十分とみなされます。

Q: 観測収束次数 p はどう計算しますか？

メッシュサイズを r 倍（通常 r=2、つまり代表寸法を半分）で順次細分化した場合、観測収束次数は p = ln(ε₃₂ / ε₄₃) / ln(r) で求めます。ここで ε_{k+1,k} は隣接解の相対変化で、最も細かい3つの解から計算した比を使います。理想的には p ≈ 設計次数（1次精度スキームなら 1、2次精度スキームなら 2）になります。p が設計値より大きく低い場合、メッシュがまだ漸近収束領域に入っていないか、解に不連続性・境界層の未解像があります。

Q: Richardson補外（外挿）で何が得られますか？

Richardson補外は、有限のメッシュ解列から「無限に細かいメッシュでの推定真値」を逆算する手法です。式は u_ext = u_n + (u_n − u_{n−1}) / (r^p − 1) で、有限のメッシュ計算からメッシュ依存誤差を除去した参照値を作れます。実機解析でメッシュ独立な解が得られていなくても、この補外値と最も細かいメッシュ解の差を「離散化誤差」として定量化できるため、論文や報告書でV&Vのエビデンスとして提示します。GCI（Grid Convergence Index）の基礎にもなっています。

Q: なぜ毎回最低3段階のメッシュで計算する必要があるのですか？

観測収束次数 p を求めるには、隣接2段階の相対変化（ε₃₂ と ε₄₃）の比が必要で、その計算に最低3段階（つまり u_粗・u_中・u_細）の解が要ります。さらに「補外値の信頼性」を見るには5段階あると安心です。手間とコストを惜しんで2段階だけで「ほぼ同じだから独立」と結論づける現場が今も多いですが、それは偶然2点が近かっただけの可能性を排除できず、メッシュ独立の証明としては不十分です。

メッシュ独立性検証シミュレーター

CFDや有限要素解析で必須の「メッシュ独立性」を、5段階のグリッドで一気に評価するツールです。粗→最微細のメッシュで得た解を入力するだけで、隣接段階の相対変化・独立達成段階・観測収束次数 p・Richardson補外値が瞬時に算出され、V&V（検証と妥当性確認）のエビデンスとして使えます。

パラメータ設定

粗メッシュ解 u₁

第1段階（最も粗いメッシュ）で得た代表解

中メッシュ解 u₂

代表寸法を半分にした第2段階の解

細メッシュ解 u₃

さらに半分（基準寸法 ÷4）の第3段階

微細メッシュ解 u₄

基準寸法 ÷8 の第4段階

最微細メッシュ解 u₅

基準寸法 ÷16 の第5段階（最も精細）

収束許容値 ε_tol

相対変化がこの値未満になった段階で独立と判定

計算結果

—

ε₂₁ (%)

—

ε₃₂ (%)

—

ε₄₃ (%)

—

ε₅₄ (%)

—

メッシュ独立達成段階

—

観測収束次数 p

—

5段階メッシュ細分化の可視化

左から右へメッシュが段階的に細分化される様子をアニメーション表示。各段階の上に解 u_k と隣接相対変化 ε(%)を表示。独立達成段階には緑のマーカー、未達なら赤のヒントを表示します。

解 u_k のメッシュ段階収束

相対変化 ε(%) と許容値

理論・主要公式

$$\varepsilon_{k+1,k}=\frac{|u_{k+1}-u_k|}{|u_{k+1}|}\times 100\%,\quad p=\frac{\ln(\varepsilon_{32}/\varepsilon_{43})}{\ln r},\quad u_{ext}=u_n+\frac{u_n-u_{n-1}}{r^{p}-1}$$

通常 r=2 で順次倍細分化する。観測 p は理想的には設計次数（1次・2次など）に近づき、Richardson補外 u_ext はメッシュ無限細分時の推定真値を与える。

$$\text{独立判定: } \varepsilon_{k+1,k} < \varepsilon_{tol}\ \text{(典型 1\%)}$$

隣接段階の相対変化が許容値を下回った段階で「メッシュ独立達成」と判断。許容値は対象量と要求精度により 0.1〜2% の範囲で選ぶ。

$$\text{推定離散化誤差: }e_n \approx \frac{u_n-u_{ext}}{u_{ext}}\times 100\%$$

最微細メッシュ解 u_n と補外値 u_ext の差から、現メッシュの離散化誤差を見積もる。報告書では補外値・観測 p・GCIをセットで提示する。

メッシュ独立性検証 (Mesh Independence Study) とは

🙋

CFDの結果を上司に出したら「これ、メッシュ独立性は取った？」って必ず聞かれます。要するに何を見ればいいんですか？

🎓

よく聞かれる質問だね。CFDも有限要素解析も、要素を細かくするほど真の解に近づいていく性質がある。逆に言うと「今のメッシュサイズで出た数字が、もう少し細かくしても変わらないか？」を確かめないと、その結果はメッシュ依存の偶然値かもしれない。だから粗→中→細→…と何段階かメッシュを変えて同じ問題を解き、注目する量（最大応力、抗力係数、出口温度など）がメッシュサイズに対してどう変化するかを見るんだ。それが「メッシュ独立性検証」。

🙋

じゃあ何段階くらいやればいいんですか？2段階で「差が小さいからOK」って報告書もよく見ます。

🎓

2段階だと「たまたま近かった」のか「本当に独立した」のか区別がつかない。最低3段階、できれば5段階欲しい。理由は、Richardson補外で観測収束次数 p を逆算するのに最低3段階の解（最微細3点）が必要だからだ。p が分かれば「無限に細かいメッシュでの推定値」も逆算できる。このツールでは5段階入れてもらって、ε₂₁〜ε₅₄の相対変化と、最微細3点から p を出している。

🙋

デフォルト値で計算したら ε₅₄ がまだ3%もあって、verdict が赤の「未達」になりました。なんで？

🎓

いいケース観察だね。u=[50,75,90,95,98] だと隣接の差が33%→17%→5.3%→3.1% と、まだ下がり方が緩い。1%の許容値に対してまだ届いていない。これは「漸近収束領域にまだ入っていない」サインで、もう1〜2段階細分化するか、メッシュの細分化方法（一様か境界層集中か）を見直す必要がある。観測 p も0.78と1未満になっていて、理論的な2次精度（p≈2）には遠い。境界層やキャプチャすべき急勾配領域がまだ解像できていない典型例だよ。

🙋

じゃあ u₅ を99とか99.5に変えてみたら、verdictが緑になりますか？

🎓

試してみて。u₅=99.5 にすると ε₅₄=|99.5-95|/99.5×100=4.52%だから、まだ届かない。u₅=95.5 だと0.52%で1%未満、独立達成（段階5）になる。重要なのは「数字を捏造して達成にする」ことじゃなく、「実機の細分化結果がそのカーブを描いている」ことを示すことだ。実務では、観測 p と Richardson補外値も併せて報告して、「現メッシュの離散化誤差は概ね2%以内」のように定量的に書くのが正しいV&Vの作法だよ。

🙋

最後にひとつ。境界層のある問題（翼まわりとか）でメッシュ独立を取るコツはありますか？

🎓

境界層問題は要注意だ。全域を一律に細分化しても、壁面近くの y+ が変わらなければ意味がない。コツは「一様倍率での全体細分化」と「境界層メッシュの y+ を独立にコントロール」の両軸でスタディすること。例えば全体は r=2 で順次倍細分化、それと別に第1層 y+ を1→0.5→0.2と振る2軸スタディだ。本ツールは1パラメータ系を扱うが、現場では2軸スタディの「全体細分化」軸の収束カーブをここに投入するイメージで使うといい。

よくある質問

CFDやFEMでは、メッシュを段階的に細かくしていくと解（応力、温度、抗力など）が真値に向かって収束していきます。隣り合う段階の解の相対変化が許容値（典型的には1%）以下になった時点で「これ以上細かくしても結果が変わらない」と判断し、その段階をメッシュ独立達成段階と呼びます。この検証なしに発表する解析結果は、結論がメッシュサイズに依存している可能性があり、V&V（検証と妥当性確認）の観点から不十分とみなされます。

メッシュサイズを r 倍（通常 r=2、つまり代表寸法を半分）で順次細分化した場合、観測収束次数は p = ln(ε₃₂ / ε₄₃) / ln(r) で求めます。ここで ε_{k+1,k} は隣接解の相対変化で、最も細かい3つの解から計算した比を使います。理想的には p ≈ 設計次数（1次精度スキームなら 1、2次精度スキームなら 2）になります。p が設計値より大きく低い場合、メッシュがまだ漸近収束領域に入っていないか、解に不連続性・境界層の未解像があります。

Richardson補外は、有限のメッシュ解列から「無限に細かいメッシュでの推定真値」を逆算する手法です。式は u_ext = u_n + (u_n − u_{n−1}) / (r^p − 1) で、有限のメッシュ計算からメッシュ依存誤差を除去した参照値を作れます。実機解析でメッシュ独立な解が得られていなくても、この補外値と最も細かいメッシュ解の差を「離散化誤差」として定量化できるため、論文や報告書でV&Vのエビデンスとして提示します。GCI（Grid Convergence Index）の基礎にもなっています。

観測収束次数 p を求めるには、隣接2段階の相対変化（ε₃₂ と ε₄₃）の比が必要で、その計算に最低3段階（つまり u_粗・u_中・u_細）の解が要ります。さらに「補外値の信頼性」を見るには5段階あると安心です。手間とコストを惜しんで2段階だけで「ほぼ同じだから独立」と結論づける現場が今も多いですが、それは偶然2点が近かっただけの可能性を排除できず、メッシュ独立の証明としては不十分です。

実世界での応用

航空機・自動車の空力解析：翼まわりや車両周りのCFD（RANS / LES）では、抗力係数 Cd や揚力係数 Cl がメッシュサイズに敏感に変わります。「Cd の変化が1%以下に収まる段階」を独立達成と定め、そのメッシュで設計検討する、というのがエアロチームの標準フローです。観測 p や Richardson補外値はSAEなどの論文投稿で必ず添付資料として要求されます。

原子炉・圧力容器のFEM応力解析：溶接ノズルの応力集中など局所最大応力が支配する問題では、メッシュ独立性の証明が規制当局への提出条件です。ASME BPVCコードでも「最大応力点で隣接細分化に対する相対変化が許容値以下」が要件として明記されており、本ツールのような相対変化追跡が解析報告書の必須セクションとなります。

射出成形・鋳造の熱流体解析：金型内充填や凝固解析では、ウェルドラインの位置やゲート近傍の流速ピークがメッシュに依存します。粗いメッシュで作った成形条件は実機でハズすことが多く、独立性検証なしの結果を量産設計に直接使うのは禁物です。3段階以上の細分化スタディを最低限のV&V手順として組み込みます。

大学教育・自習：CFD/FEM初学者がよくつまずくのが「結果の信頼性」です。本ツールは実際に解析を回さなくても、教科書例題の収束カーブを入れて p や補外値を直感的に体験できるため、講義や演習の補助に使えます。教員からも「学生がメッシュ独立の意味をすぐに掴める」教材としての利用例があります。

よくある誤解と注意点

まず最大の落とし穴が、「全体メッシュ数だけ報告して、独立性検証を省略する」こと。「総要素500万」と書いても、それが粗いか細かいかは問題依存で、独立性の証明にはなりません。学会発表や論文で必ず聞かれるのが「収束履歴」と「観測 p」です。例えば翼の Cd 解析で「メッシュ200万→500万→1000万→2000万→4000万で Cd が 0.0312→0.0298→0.0291→0.0289→0.0288」のように細分化ステップとセットで提示するのが最低ラインです。本ツールの5段階フォーマットは、まさにこの提出形式に対応しています。

次に、「観測 p が設計次数と一致するはず」と思い込むこと。2次精度スキームを使っていても、メッシュがまだ漸近収束領域に入っていなければ p は 0.5〜1.5 とばらつきます。逆に p が極端に大きい（>3）場合は、解の振動や乱数的なノイズで偶然 ε が近づいたケースが多く、信頼できません。理想的な p は「収束領域に十分入った3点で、設計次数の ±50% の範囲」が目安です。範囲外ならメッシュをさらに細分化するか、収束履歴そのものを疑ってください。

最後に、「補外値 u_ext を真値だと信じてしまう」こと。Richardson補外はあくまで「現メッシュ列を直線的に外挿した推定値」であり、本当の真値ではありません。特に乱流モデルの近似誤差や境界条件の不確定性はこの補外で取り除けません。報告書には「メッシュ離散化誤差は u_ext を基準として推定〜%」と書き、「真値」とは書かない慎重さが求められます。GCI（Grid Convergence Index）では補外誤差にさらに安全係数 1.25 を掛けるのが標準実務です。

メッシュ独立性検証 (Mesh Independence Study) とは

よくある質問

実世界での応用

よくある誤解と注意点

使い方ガイド

具体的な計算例

実務での注意点