クエット流れ（平行平板間せん断流） — ソルバー比較

5セル/高さ方向のメッシュでの速度プロファイル比較だ。

全ソルバーで機械精度レベルの一致。COMSOLは二次有限要素で線形場を厳密に表現するから丸め誤差も極小になる。

壁面近傍の速度補間方法がソルバーによって異なるため、壁隣接セルでのグラディエント計算に微小な差が出ることがある。特にOpenFOAMのcorrected snGrad と uncorrected snGrad で壁面フラックスの精度が変わる。非直交メッシュでは差が拡大するから、skewness > 0 のメッシュでの比較テストも有用だ。

クエット流れの場合、div(phi,U) のスキームは upwind でも linearUpwind でも limited linear でも結果は同じ（線形場に対して全て厳密）。差が出るのは非直交メッシュでの laplacian(nu,U) のスキーム選択だ。

• Gauss linear corrected: 非直交補正あり。推奨
• Gauss linear uncorrected: 非直交補正なし。直交メッシュでは問題ないが非直交で精度低下
• Gauss linear limited 0.5: 非直交補正を制限。安定性重視

その通り。クエット流れ + 直交メッシュは全てのスキームが同じ結果になる。これこそがCode Verificationの起点として優れている理由だ。ここから非直交メッシュや非線形問題に進むことで、各設定の影響を段階的に切り分けられる。

領域分割の境界でのプロセッサ間通信が正しく実装されていれば、逐次実行と完全に一致するはずだ。ただし浮動小数点の加算順序が変わることで丸め誤差レベル（$10^{-15}$オーダー）の差は出る。

クエット流れを2分割・4分割・8分割で実行し、結果が$10^{-10}$以内で一致することを確認する。これは並列化の回帰テストとしても有用で、MPIライブラリのバージョンアップ時に走らせるべきだ。