1 メッシュ分割の概観
1.1 用語と目的
メッシュ分割とは、連続的な空間領域を多数の小領域（要素）に分割し、境界や内部の幾何形状を離散的に近似する作業である。これにより、偏微分方程式などの支配方程式を、要素ごとの未知量として扱える形に変換するための基盤が得られる。
目的は、(1) 計算結果の精度を確保すること、(2) 計算資源（計算時間・メモリ）を許容範囲に収めること、(3) 解法の安定性を保つこと、の三点を同時に満たすように分解能と配置を設計する点にある。実務では、目的変数（例：応力、速度、温度、濃度）に対して必要な情報量が確保されているかを中心に判断する。
1.2 数値解析との関係
1.2.1 有限要素法での役割
有限要素法では、メッシュは形状表現と近似関数（基底）を定義する土台になる。要素ごとに形状を定義し、節点や要素内の積分により弱形式が構成されるため、要素の幾何精度と分割のきめ細かさがそのまま近似の品質に結びつく。特に境界付近や勾配の大きい領域では、要素サイズと曲率追従性が誤差に直結しやすい。
1.2.2 有限体積法での役割
有限体積法では、要素は制御体として解釈され、保存則が要素境界を介した流束の釣り合いとして離散化される。したがって、メッシュは幾何情報だけでなく、面積、法線方向、面間距離などの幾何量を通じて流束計算の精度を左右する。格子がねじれたり極端に歪んだりすると、流束評価が不適切になりやすく、結果として物理量の保存や整合性が崩れる場合がある。
1.2.3 有限差分法での位置づけ
有限差分法では、典型的には規則格子（構造格子）上で微分演算子を格子点間差分に置き換えるが、複雑形状では座標変換や格子生成により準構造的な配置が導入されることもある。メッシュは差分係数や離散化誤差の源になるため、格子幅の偏りや境界近傍の取り扱いが誤差・安定性に影響する。特に段階的な格子幅変更が大きい場合、誤差が局所に蓄積しやすい。
1.3 メッシュの種類の全体像
メッシュは、分割形状と格子の接続構造に基づき大別できる。代表的には、四面体や六面体などの多面体要素で構成される体メッシュ、境界形状に対応した表面メッシュ、そしてこれらを連携させた一連の手順がある。要素配置の規則性により、構造格子（格子線が連続して規則的に並ぶ）と非構造格子（要素が多様な形状で接続される）が区別される。
構造格子はデータ参照が容易で計算効率に優れることが多い。一方、非構造格子は複雑形状の追従性に強みがある。近年は、境界適合性と自動生成の両立を目指し、ハイブリッドな構成（例：内部は非構造、境界近傍は層状に整える）もよく用いられる。
2 メッシュの生成手法
2.1 自動生成
2.1.1 形状ベース分割
形状ベース分割では、CAD形状や三角形などの表現から計算領域を抽出し、その幾何情報を基に要素を生成する。最初に境界（表面）を分割し、次に内部（体）へ拡張するのが典型的な流れである。自動化では、形状の欠損（穴、自己交差、表面の反転）や冗長なエッジがそのまま生成失敗の原因になり得るため、前処理として幾何を整える作業が重要になる。
2.1.1.1 表面メッシュと体メッシュの連携
表面メッシュと体メッシュの連携は、境界の近似精度を確保するための中核工程である。表面側で規定した節点配置や分割密度が、体側の要素サイズや方向付けに反映されるため、両者の不整合は境界直近の誤差増大や流束の不自然な挙動につながる。連携が適切に行われると、境界法線方向や曲率変化に沿って要素が自然に配置されやすい。
2.1.2 格子生成（構造格子）
構造格子の自動生成では、格子の生成ルール（格子線の方向、分割数、伸長率など）を指定し、領域全体を規則配置で埋める。入口や壁面など特定方向に対して格子を整える用途では、方向制御が可能な点が利点になる。ただし複雑な曲面や入り組んだ空間では、規則格子だけで高品質を得るのが難しく、非構造部との併用が検討される。
2.1.3 非構造格子生成
非構造格子生成では、対象形状をメッシュ用の要素（例：三角形、四面体）で充填する。代表的手法として、まず表面を三角形に分割し、その後に体へ拡張することで内部要素を作る方法がある。生成では、節点密度、要素サイズ分布、曲率追従の目標、ならびに許容される品質基準が使われる。品質は生成時に評価され、基準を満たさない要素は局所的に再生成・改善される。
2.2 手動・半自動生成
2.2.1 境界形状の指定
手動・半自動生成では、境界形状への適合度を高めるために、重要境界の選定と分割条件の指定を行う。曲面の近似をより滑らかにするために、分割密度や許容誤差を制御し、要素が境界の形状変化に追従するよう調整する。境界に沿ったメッシュは、境界条件の適用精度を左右するため、誤差の大きい領域を優先して整える運用が多い。
2.2.2 局所細分化の設計
局所細分化では、解析結果が急激に変化する可能性のある領域に対して、要素サイズを小さくする方針を立てる。典型例として、応力集中点、接触周辺、境界層、衝撃や混合が起きやすい領域などが挙げられる。設計では、細分化の範囲と遷移（粗い側から細い側への変化）を適切に設定し、急な格子幅変化による数値誤差の増幅を抑える。
2.3 リメッシング（再分割）
2.3.1 メッシュ変形時の対応
リメッシングの一形態として、変形問題ではメッシュを動かしながら計算する場合がある。このとき、要素が歪みすぎると品質基準を逸脱し、離散化の信頼性が低下する。そこで、所定の品質指標が閾値を超える段階で局所再生成や節点移動を行う。完全再分割まで行わない場合でも、品質の劣化を検知する仕組みが実務では重要になる。
2.3.2 アダプティブ分割
アダプティブ分割では、計算結果の誤差推定や指標に基づいて、必要な場所だけを再細分化する。誤差が大きい領域に計算資源を集中させるため、同等の精度をより少ない要素数で達成しやすい。実装上は、誤差指標の定義、再分割のルール（過度な細分化の抑制、粗化の条件）、再計算に伴う安定性の確保が設計ポイントとなる。
3 メッシュ品質と評価指標
3.1 要素品質
3.1.1 歪み（スキュー）とアスペクト比
要素品質は、形状の理想からのずれとして評価されることが多い。歪み（スキュー）は要素の角度や形状が理想配置からどれだけ外れているかを示し、アスペクト比は寸法間の比がどれほど極端かを表す。一般に、これらの値が悪化すると、勾配の評価や積分の精度が低下し、さらに数値解法によっては収束性にも影響する。
3.1.2 要素サイズの整合性
要素サイズの整合性は、隣接要素間の大きさの変化が過大でないことを意味する。粗い要素から急に細い要素へ切り替わると、離散化誤差の分布が不均一になり、局所的な残差が過大になることがある。適切な遷移を設けることで、誤差の急増を抑えつつ全体の要素数を増やしすぎないバランスが得られる。
3.1.3 角度・曲率の表現
曲面が多い問題では、角度の表現と曲率追従が品質を左右する。表面メッシュが十分に細かくない場合、境界は見かけ上滑らかでも離散的には段差や局所的な誤差として現れる。曲率の変化が速い箇所では、要素を増やすだけでなく、要素配置が曲率方向に沿うような工夫が必要になることがある。
3.2 トポロジーと整合性
3.2.1 接続性（隣接関係）
接続性は、要素がどのように隣接し、節点や辺、面を介してどう繋がるかを指す。有限体積法では特に、面を共有する隣接関係が流束評価の前提になるため、誤った接続は物理量の保存を壊す原因になる。有限要素法でも、整合した接続は形状関数の組み立てや積分領域の定義に影響する。
3.2.2 縫い目・ギャップの扱い
複雑形状では、分割時に境界が完全に一致せず、縫い目（分割境界）や微小ギャップが残ることがある。ギャップや重なりがあると、境界条件の適用先が曖昧になったり、物理的には連続であるはずの領域が分断されたりする。したがって、許容寸法を含めて「連結すべきか」「切り分けるべきか」を事前に方針化し、必要に応じてジオメトリ修復やメッシュ結合を行う。
3.3 メッシュ品質が与える影響
3.3.1 精度への影響
メッシュ品質の低下は近似誤差の増大に繋がる。要素の形が不適切であれば、勾配や曲率に関する情報が離散的に歪み、結果として応答変数が局所的に過大あるいは過小に評価される。特に誤差が境界近傍や特異性の近くに集中する場合、品質の差が顕著に表れやすい。
3.3.2 計算コストへの影響
要素数の増加は計算量とメモリ使用量を増やし、反復解法では線形代数処理が支配的になることが多い。品質改善のために全体を細かくする代わりに、必要な場所のみを局所的に解像する設計が望ましい。加えて、品質劣化が原因で反復回数が増えると、要素数が同程度でも計算コストが跳ねる場合がある。
3.3.3 数値安定性への影響
安定性は、離散化の性質と解法の相性で左右される。歪みが大きいメッシュでは、離散演算子が想定より条件が悪化し、解の発散や振動の原因になり得る。さらに、強い勾配や急峻な境界条件の取り扱いと重なると、品質の改善が安定性に直結することがある。したがって、メッシュ品質の閾値は解析設定と同時に整合させる必要がある。
4 メッシュ設計の実務指針
4.1 目的に応じた解像度設計
4.1.1 応力集中・境界層への配慮
応力集中が予想される箇所では、幾何学的特徴（段差、曲率の急変、孔周り）と力学応答の変化が重なるため、要素サイズの細分化が有効になる。境界層のように速度や温度が短い距離で変化する領域では、壁近傍での解像度確保が重要である。単に密にするだけでなく、壁面からの距離に応じた増加率や方向性が設計の鍵になる。
1.1 物性変化や勾配の大きい領域
物性が空間的に変化する問題、あるいは温度・濃度・速度勾配が大きくなる領域では、解像度が不足すると局所的な物理過程が平均化され、誤差が系統的に残りやすい。細分化の対象は必ずしも全領域ではなく、支配方程式の特性や期待される支配領域に応じて選ぶことが望ましい。設計段階では予備計算や既知の傾向を参照し、細分化の範囲を合理化する。
4.2 境界条件とメッシュの整合
4.2.1 形状と境界の適合度
境界条件の適用は、メッシュ上で定義される境界の幾何に依存する。たとえば壁面での速度条件、断熱境界での熱条件、接触条件の扱いなどは、境界形状の近似誤差に影響される。適合度が不足すると、物理的に正しい境界であるはずの位置が実質的にずれて見え、誤差の原因となる。
4.2.2 物理量の取り扱い（離散化）
離散化における物理量の定義は、メッシュの配置と密接に連動する。有限要素法では節点自由度と積分手順が、有限体積法では制御体境界の流束表現が鍵になる。境界や内部に対して、物理量の取り扱いが整合しているかを確認し、例えば保存則が成立する形で流束が計算されるか、または境界近傍での勾配評価が破綻していないかを点検する。
4.3 メッシュ独立性の検証
4.3.1 収束傾向の確認
メッシュ独立性は、要素数を増減したときに主要な結果が大きく変わらないことを意味する。実務では複数段階の分割で同一条件を再計算し、誤差が減少する傾向、あるいは指標が安定する兆候を確認する。理想的には単調性や加速された収束が期待されるが、境界近傍の解像や離散化方式の影響で必ずしも単純にならないため、複数指標で判断する。
4.3.2 代表量（応力・流量等）の選定
独立性の判定に用いる代表量は、目的とする現象を反映している必要がある。応力解析では最大応力、平均応力、特定位置での応答などが選ばれることが多い。流体解析では流量、圧力差、壁面摩擦に相当する量、場合によっては渦度のような指標が用いられる。代表量の選び方が不適切だと、メッシュを増やしても判定が曖昧になり、実際の誤差を見落とす可能性がある。
4.4 よくある問題と対策
4.4.1 解析が発散する原因
発散は、メッシュ品質の不十分さ、境界条件の整合性不足、離散化設定の不適切さなど複数要因で起こり得る。特に歪みや極端なサイズ比があると、離散演算子の条件が悪化し、反復計算で解が発散することがある。まずメッシュ品質指標を点検し、その後に解法パラメータ（緩和、時間刻み、収束基準）や境界近傍の設定を確認すると、原因の切り分けが進む。
4.4.2 不自然な結果（ロッキング等）
ロッキングは接続や幾何関係の不整合、あるいは要素の変形により数値的に望ましくない自由度挙動が生じるときに観測されることがある。有限要素法では要素の形状や制約の与え方が影響し、有限体積法では不適切な面法線やサイズ比が不自然な分布につながる場合がある。対策としては、該当領域の要素改善（再分割や品質向上）、境界条件の再定義、あるいはモデル化の見直しが行われる。
4.4.3 メッシュ改善の手順
改善の手順は、観測された異常を起点に局所的に介入する方針が一般的である。まず不自然な残差分布や勾配の過大箇所を特定し、次にその周辺の要素サイズと形状品質を確認する。必要に応じて、細分化範囲を拡大する、遷移率を緩和する、曲面追従のために表面メッシュを修正する、といった具体策を段階的に適用する。最後に、再計算して代表量の変化と収束挙動を確認し、改善が有効であることを検証する。