CFDメッシュ品質計算機 戻る
流体解析

CFDメッシュ品質指標計算機

アスペクト比・スキューネス・直交性をリアルタイム計算。OpenFOAM・ANSYS Fluent・Star-CCM+の許容範囲と比較してメッシュ品質を即座に評価。

ラプラシアン平滑化によるメッシュ品質の改善(リアルタイム)

一時停止中はアニメーションが静止します。再生ボタンで自動アニメーションを再開できます。

セル数
平均スキューネス
最大アスペクト比
最小直交品質
品質しきい値未満の割合
平滑化反復
理論・主要公式

このアニメーションは、歪んだ2Dメッシュにラプラシアン平滑化を反復適用し、各セルを品質(緑=良好・赤=不良)で色分けして表示します。反復が進むと平均スキューネスが低下する様子を確認できます。

等角スキューネス:

$$S = \frac{\theta_{max}- \theta_{equi}}{180 - \theta_{equi}}$$

$\theta_{equi}$: 理想等角(四辺形=90°)。$S\in[0,1]$(0=理想)

アスペクト比: $AR = L_{max}/ L_{min}$

直交品質: $Q_{orth} = 1 - S$

ラプラシアン平滑化: $\mathbf{x}_i \leftarrow \mathbf{x}_i + \alpha\,(\bar{\mathbf{x}}_{nb} - \mathbf{x}_i)$(隣接ノードの平均へ移動)

CFDメッシュ品質指標とは

🙋
メッシュ品質って、何でそんなに重要なんですか?ただの計算の「下準備」じゃないんですか?
🎓
大まかに言うと、メッシュの「形の悪さ」がそのまま計算結果の「信用度」を決めるんだ。例えば、自動車の空力解析で、ボンネット上のメッシュが極端に歪んでいると、そこでの圧力分布がガタガタになって、実際の流れを全然再現できなくなるよ。このシミュレーターのレーダーチャートで、5つの指標が一気に赤くなったら、それは「このメッシュで計算しても意味ないかも」という警告だね。
🙋
え、そうなんですか!じゃあ、上の「スキューネス」のスライダーを動かすと、具体的にどんなメッシュになるんですか?
🎓
スキューネスが大きい(例えば0.9に近い)メッシュは、四角形が菱形や、三角形が細長くつぶれた形になる。実務で多いのは、複雑な形状を無理やりメッシュで埋めた時にできる「変形要素」だ。こういう要素があると、Navier-Stokes方程式を離散化する時の「勾配計算」が大きく狂って、圧力や速度が振動する原因になる。ツールで値を変えると、レーダーチャートのスパイクが伸びて、視覚的に「まずい」とわかるようになってるよ。
🙋
「アスペクト比」と「膨張比」って似てますけど、何が違うんですか?両方とも「大きさの違い」ですよね?
🎓
良い質問だ!アスペクト比は「1つのセル内での」最大辺と最小辺の比だ。例えば、境界層を表現する細長い六面体メッシュは、意図的にAR=1000とかになる。一方、膨張比は「隣り合うセル同士の体積の急激な変化」を表す。これが大きいと、流れ場の物理量(温度や圧力)がセル境界で不連続になり、計算が発散しやすくなる。シミュレーターで両方を別々に動かして、レーダーチャートの反応を見比べてみると、その違いがよく分かるはずだ。

よくある質問

アスペクト比は、メッシュ要素の最長辺の長さを最短辺の長さで割った値です。1に近いほど理想的な正方形や立方体に近く、値が大きいほど細長い要素となります。OpenFOAMでは一般的に10以下、ANSYS Fluentでは5以下が推奨されます。
スキューネスが高いと、要素内の角度が理想形状から大きくずれるため、数値計算の離散化誤差が増大します。特に勾配計算や拡散項の精度が低下し、解の発散や非物理的な振動を引き起こす原因となります。目安として0.8以上は修正が必要です。
本ツールはOpenFOAM、ANSYS Fluent、Star-CCM+の3つの主要ソルバーに対応しています。各ソルバーが推奨するアスペクト比、スキューネス、直交性の許容範囲と自動比較し、メッシュ品質を「良好」「注意」「要修正」の3段階で評価します。
直交性が悪い場合は、メッシュのスナップ処理やレイヤー設定を見直してください。具体的には、境界層メッシュの成長率を小さくする、スムージング反復回数を増やす、または六面体メッシュの場合はブロック分割を調整してセル面が流れ方向に直交するよう改善します。

実世界での応用

自動車・航空機の外部空力解析:車体や翼表面の境界層メッシュでは、高アスペクト比の細長い要素が使われます。しかし、ボディ角部などではスキューネスが大きくなりがちで、ここをチェックしないと抗力係数(Cd値)の予測精度が大幅に低下します。

電子機器の熱流体解析:ヒートシンクのフィン間やICチップ周りは、幾何形状が複雑でメッシュ品質が悪化しやすい領域です。膨張比が急激に変化するメッシュでは、熱伝達率の計算が不安定になり、過大または過小評価されるリスクがあります。

化学反応を伴う燃焼解析:反応が起こる火炎面近傍では、濃度や温度の勾配が非常に急峻です。ここでの直交性の悪いメッシュ(スキューネス大)は、種々の輸送方程式の離散化を歪め、反応速度や火炎伝播速度の計算を誤らせます。

CAEソフトウェアとの連携:OpenFOAMのcheckMesh、ANSYS FluentのMesh Quality Report、Star-CCM+のMesh Diagnosticsは、本ツールと同様の指標を計算・警告します。高品質メッシュは、ソルバーの収束性を向上させ、より信頼性の高い結果を短時間で得るための必須条件です。

よくある誤解と注意点

まず、「レーダーチャートが全部緑なら完璧」という誤解があります。実は、解析の種類や領域によっては、ある指標を意図的に悪くすることがあります。例えば、航空機の翼まわりの境界層解析では、壁面に垂直方向に極端に細長いメッシュ(アスペクト比1000以上)を使います。これは物理現象を正しく捉えるための「戦略的な悪化」です。ツールの許容範囲はあくまで一般的な目安。最終的には「何を計算したいか」で判断しましょう。

次に、「ツールでチェックしたから大丈夫」という過信。このツールはメッシュの「幾何学的な品質」を見ていますが、それだけでは不十分な場合があります。例えば、衝撃波が発生する超音速流れでは、流れ方向に対して十分に細かいメッシュ解像度(これはツールの指標では測れない)がなければ、衝撃波の位置がずれてしまいます。幾何学的品質と物理的解像度、両方のバランスが肝心です。

最後に、局所的な悪い要素の見落とし。全体の平均値は良くても、流れの重要な部分(例えば、自動車のドアミラー後方の剥離点)にひとつだけ「スキューネス0.9」の悪い要素が混じっていると、そこから誤差が増幅され、全体の流れ場をダメにすることがあります。ツールで指標を計算した後は、数値が悪い要素が「どこに」存在するかを必ず可視化して確認する習慣をつけましょう。

使い方ガイド

  1. メッシュ統計情報からアスペクト比(AR値)を入力。スライダー範囲は1~100で、境界層ではより大きな値を別途確認します
  2. スキューネス(0~1)を設定。0に近いほど良好で、OpenFOAMは0.85未満、Fluentは0.90未満を推奨
  3. 非直交角θ(0~89°)を指定します。cos θは出力値で、角度が小さいほど直交性は良好です
  4. スムースネス比と膨張比ER(成長率)を入力し、総合スコア100点満点で判定

具体的な計算例

良好な翼型周辺メッシュ例:アスペクト比AR=20、スキューネス=0.35、非直交角θ=25°、スムースネス比=1.2、膨張比ER=1.2。ページの重み付けでは各サブスコアがAR 70、スキューネス88、直交性85、スムースネス100、ER 100となり、総合スコア85点で「優良」と判定されます。

実務での注意点