Code_Saturne — CAE用語解説

本質は「圧縮性・非圧縮性の定常・非定常ナビエ-ストークス方程式」を解く汎用CFDソルバーです。最大の特徴は、フランス電力公社(EDF)が自社の原子力プラントやタービン、建屋の換気解析のために開発した「産業向けオープンソース」ソフトウェアである点です。商用ソフトのようにGUIが充実しているわけではなく、CやFortranでユーザー定義関数を書いて物理モデルを拡張する文化が根付いています。

デフォルトの求解アルゴリズムは、非圧縮性または弱圧縮性流れを想定した「圧力補正法」です。連続の式と運動量式は以下の通りです。

いえ、非常に豊富です。産業用途を強く意識しており、RANSモデルでは標準k-ε、RNG k-ε、k-ω SST、低レイノルズ数対応のv2-f (Elliptic Blending)モデルが使えます。特にv2-fモデルは、壁面近くのアニソトロピーを比較的安価に捉えられるため、原子炉内の複雑な流路解析で重宝されます。また、LES（Large Eddy Simulation）もネイティブでサポートしており、SmagorinskyモデルやDynamic Smagorinskyモデル、WALEモデルが利用可能です。

その通り、FVMがコアです。任意の形状の非構造格子（テトラ、プリズム、六面体、ピラミッド、多面体）を扱えます。離散化スキームは、対流項に対して1次風上差分から2次中心差分、または各種の高次精度スキーム（例えば、Blended Gammaスキーム）を選択可能です。実務では、数値拡散と安定性のバランスから、2次精度の中心差分とわずかな1次風上をブレンドすることが多いです。

デフォルトはSIMPLEの改良系であるSIMPLEアルゴリズムです。ただし、非定常計算ではPISO（Pressure Implicit with Splitting of Operators）アルゴリズムも選択できます。PISOは1タイムステップ内で複数回の圧力修正を行うため、大きな時間ステップを使う非定常計算で有利です。Code_Saturneでは、この圧力-速度連成ソルバーに加え、共役勾配法（CG）や代数的マルチグリッド法（AMG）を前処理とした線形ソルバーが組み込まれており、大規模問題（数億セル）でも効率的に求解できます。

はい、MPI（Message Passing Interface）による領域分割法が基本です。計算格子を複数のサブドメインに自動分割し、各サブドメインを1つのMPIプロセスが担当します。通信ライブラリにはMPIの他、大規模クラスタ向けにOpenMPIやMPICHとの組み合わせが一般的です。EDFの内部では、数千コアを使った原子炉スケールの解析が日常的に行われており、そのノウハウがソルバーのスケーラビリティに反映されています。

典型的なワークフローは「前処理→求解→後処理」に分かれます。前処理（メッシュ生成）には、Salome Platform（同じくEDF系のオープンソースCAEプラットフォーム）や、商用のANSYS Fluent Meshing、Siemens Star-CCM+が使われます。メッシュ形式はCGNSまたはMED（Salome形式）が標準的です。後処理には、ParaViewやEnSightが用いられます。Code_Saturne自体には簡易GUI「GUI」が付属しますが、本格的な設定はXML形式の設定ファイルを直接編集することが多いです。

特に重要なのは「非直交性」と「アスペクト比」です。Code_Saturneのソルバーはロバストですが、非直交性が60度を超えるセルが大量にあると、圧力修正の収束が極端に悪化したり、発散したりします。アスペクト比（セルの伸び）は1,000を超えないようにするのが無難です。また、壁面の最初のセル厚さは、使用する壁面関数（例えば、標準壁関数ならy+≈30〜300）に合うように設定する必要があります。これらは`check_mesh`ユーティリティで事前に診断できます。

残差の低下（例えば、連続の式の残差が10^-3以下）は必要条件ですが、十分条件ではありません。必ず物理量のモニタリングを併用してください。例えば、流出口の平均圧力、代表点の速度、抗力係数などが定常値に収束しているかを確認します。Code_Saturneでは`probes`や`monitoring`機能で任意の点や面の値を時系列で出力できます。産業用途では、これらモニター値の変動が初期値の0.1%以下で安定することを収束の目安とすることが多いです。

3つのオープンソース/商用ソフトを比べるとこうなります。

顕著なのは「多相流モデルの充実度」と「組み込みの最適化ツール」です。FluentにはVOF、Eulerian、Mixtureモデルに加え、破砕・合体モデルなど非常に細かいオプションがあります。Code_Saturneにも欧州原子力共同体のコードであるNEPTUNE_CFDと連携する機能はありますが、標準装備としては限定的です。また、Fluentに組み込まれたAdjoint Solverによる形状最適化や、ANSYS DesignXplorerとの連携によるDOE（実験計画法）は、Code_Saturne単体ではできません。これらはユーザーが外部スクリプトで実装する必要があります。

最大の強みは「特定の産業分野における検証済みのモデルとノウハウ」です。例えば、原子炉のコア内の熱水力解析、タービンブレードの冷却、建屋内の火災・煙流動（火災シミュレーション専用フロントエンドSalamandreと連携）、放射性エアロゾルの拡散など、EDFが実際に規制当局に提出する安全審査で使用してきた実績があります。これらのテンプレート的な設定と、大規模並列計算に対する高いスケーラビリティは、同分野では大きなアドバンテージです。

まずは「初期条件」と「境界条件」の不一致を疑います。例えば、流入速度を設定したのに流出境界が全て壁面だったり、圧力境界が重複していたりすると、ソルバーが発散します。次に「メッシュ品質」、特に体積がゼロまたは負の「潰れたセル」がないか`check_mesh`で確認します。最後に「物性値」です。密度や粘度に極端に小さい値（ほぼゼロ）や大きい値を設定していないかチェックしてください。デバッグには、`--log`オプションを詳細モード（`debug`）にして実行し、どのルーチンで落ちたかを特定します。

これはほぼ間違いなく「数値的不安定性」です。以下の順で対策します。

これは典型的な「並列効率の低下」と「通信オーバーヘッド」の問題です。原因は主に二つ。