検証 — CAE用語解説

いい着眼点です。主な誤差は、(1)離散化誤差、(2)数値（丸め）誤差、(3)反復（打ち切り）誤差です。メッシュを細かくすると(1)は減りますが、要素数が増えると(2)が累積され、(3)の反復回数も増える可能性があります。例えば、Ansys Mechanicalで非線形解析を行う際、力の収束基準をデフォルトの0.5%から0.1%に厳しくすると、反復誤差は減りますが計算時間は増加します。このトレードオフを理解することが重要です。

残差は支配方程式（例えばナビエ-ストークス方程式）の「満たされなさ」を示すので、反復誤差の指標です。ただし、残差が10^-6まで下がっても、離散化誤差が大きければ全体の精度は低いままです。実務では、残差の収束と並行して、モニタポイント（例えば翼の抗力係数Cd）の値が定常状態に達したかも確認します。FluentやSTAR-CCM+では、このモニタ値を履歴としてプロットする機能があります。

離散化スキームは離散化誤差の「次数」を決定します。1次風上では誤差がメッシュサイズhに比例（O(h)）ですが、2次風上ではhの2乗に比例（O(h^2)）します。つまり、メッシュを2倍細かくすると、1次スキームでは誤差が約1/2に、2次スキームでは約1/4になります。検証のメッシュ収束性調査では、高次スキームを使うほど少ないメッシュ数で高精度が得られる傾向があります。ただし、解の勾配が急峻な領域では発振する可能性があり、その場合は数値粘性を追加するなどの対策が必要です。

まず、関心ある出力量（Quantity of Interest: QoI）を明確にします。「強度が知りたい」ではなく「ボルト締結部の第一主応力の最大値」などと具体化します。次に、そのQoIに対して最も影響が大きいと思われるパラメータ（通常はメッシュサイズ）を選び、パラメータスタディを計画します。例えば、基本メッシュサイズを5mmとし、それを1.5倍、0.75倍、0.5倍にした4ケースを実行し、結果の変化をグラフ化します。

その通りで、局所細分化が鍵です。最初の解析で、応力集中部や流速・温度の勾配が大きい領域を特定します。Abaqus/CAEの「メッシュ制御」機能やAnsys Meshingの「Face Sizing」を使い、関心領域だけメッシュサイズを1/2、1/4と段階的に細かくします。コーナーや穴の縁など、理論的に特異点となる場所はメッシュを細かくしても応力が発散するので、その場合は応力集中係数などの別の評価指標に切り替える必要があります。

業界によって様々ですが、自動車業界では多くの場合、内部のチェックリスト様式に従います。少なくとも以下の項目を記録すべきです：使用ソフトウェアとバージョン（例：Abaqus 2023）、要素タイプと積分点数（例：C3D8R：8節点六面体低減積分）、メッシュ統計（要素数、節点数）、収束判定基準と結果（例：最大ミーゼス応力の相対変化が2%以下）、そして最終的に採用したメッシュケースの根拠。これらはトレーサビリティのために不可欠です。

はい、近年は高度化しています。Ansys Workbenchには「Convergence」ツールがあり、指定したQoIに対してh-refinementまたはp-refinement（要素次数向上）を自動実行し、結果の収束をグラフで示してくれます。COMSOL Multiphysicsでは、「メッシュシーケンス」研究ステップを設定でき、一連の異なるメッシュでパラメータスイープのように解析を連続実行できます。Abaqus/CAEには標準ではそのような自動ツールは少ないですが、Pythonスクリプトでメッシュサイズを変えてバッチ実行するマクロを組むことが一般的です。

非常に重要です。悪い品質の要素は、離散化誤差を著しく増大させ、場合によってはソルバーが発散する原因になります。各ソフトは強力なチェック機能を持ちます。例えば、Siemens NX Nastranの「Element Quality Check」では、Jacobian比、ワーピング量、アスペクト比などを一括評価し、しきい値を超える要素を赤色でハイライトします。一般的な目安として、アスペクト比は10以下、スキュー角は60度以下を目標とします。ただし、どうしても悪い要素が生じる領域では、その影響を評価するための感度解析（良いメッシュの領域と結果を比較）が必要です。

主に2点です。1点目は「コード検証」の情報が限定的であることです。商用ソフト（Ansys, Dassaultなど）は巨額の開発費をかけ、ベンチマークテストでコード検証を重ねていますが、オープンソースではその履歴が不明確な場合があります。2点目は、デフォルト設定の信頼性です。OpenFOAMは非常に柔軟ですが、その分、離散化スキームやソルバー設定をユーザーが正しく選択する責任があります。したがって、既知の厳密解やベンチマーク問題（例：NASAの後流円柱問題）で自分自身でソルバー設定を検証する作業が必須になります。

これは「メッシュロック」または「数値的不安定性」の兆候かもしれません。メッシュロックは、例えば完全積分の1次要素（C3D8）で体積不可圧材料（ポアソン比0.5に近い）を扱う時に発生し、要素が過剰に剛性を示します。対策は低減積分要素（C3D8R）やハイブリッド要素の使用です。数値的不安定性は、CFDで高次風上差分スキームを使い、メッシュが細かすぎる場合に発生することがあります。この場合、ソルバーの緩和係数を調整するか、スキームを若干1次寄りに変更する必要があります。

問題です。これは、粗いメッシュでは物理現象を捉えるのに必要な空間分解能が足りていないことを示しています。特に非線形問題（接触や材料非線形）で顕著です。例えば、Abaqusでの接触解析で、粗いメッシュでは接触面が「跳び越えて」しまい、過大な貫入が発生して発散することがあります。検証の観点では、収束した細かいメッシュの結果を「より正確」と見なせますが、その計算コストが現実的でない場合、中間的なメッシュで安定して解が得られるかを確認し、そのメッシュを採用する判断をします。

これはよくある、そして重要な事象です。変位は全体的な挙動で、比較的粗いメッシュでも収束しやすい。一方、応力は微分値（変位の勾配）なので、特に応力集中部では非常に細かいメッシュが必要になります。判断基準は、設計判断に用いる「最もクリティカルなQoI」です。強度設計が目的なら、最大応力の収束を優先します。その際、全体のメッシュを一律細かくするのではなく、応力集中部のみを局所細分化するアプローチが計算効率上、必須となります。最終レポートには、「変位は5mmメッシュで収束、応力は2mmメッシュで収束を確認したため、応力評価には2mmメッシュを使用」と明記します。