体積流量 — CAE用語解説

良い質問だ。CAEでは、扱う流体の種類と「どの状態での体積か」が重要になる。非圧縮性流体の場合は単純で、密度ρが一定だから、体積流量Qは断面積Aと流速vの積で定義される：

その通り。例えば、Ansys Fluentで圧縮性流れを解く時、境界条件として「Mass Flow Rate」を設定する方が一般的だ。なぜなら質量は保存されるが、体積は状態で変化するからだ。もし「Volume Flow Rate」を設定するなら、その値が参照する温度・圧力を明確にしなければならない。航空機エンジンのシミュレーションでは、国際標準大気（ISA）の海面状態（15℃, 101325 Pa）に換算した標準体積流量（Normalized Volume Flow Rate）を使うことが多い。

連続の式（質量保存則）の積分形はこれだ：

「流速（Velocity Inlet）」を設定するのは最も直接的だが、全体の流量が計算結果に依存してしまう。一方、「体積流量（Volume Flow Rate）」や「質量流量（Mass Flow Rate）」を設定する場合、ソルバーは内部でそれを満たすような流速分布を逆算する。例えば、FVMでは、流量指定境界は通常「ニューマン型境界条件」として扱われる。ソルバーは境界セル面の流速を調整する係数を反復計算の中で更新し、指定された流量値に収束させる。OpenFOAMの`flowRateInletVelocity`境界条件がその典型例だ。

する。特に非定常計算や、流入方向が複雑な形状の場合は注意が必要だ。対策として、まず「流速指定」で定常解を求め、その結果の流速分布を初期場として「流量指定」の計算に引き継ぐワークフローが使われる。また、Ansys CFXでは「Mass Flow Rate」を設定する際、流入方向の指定方法（Normal to Boundary, Direction Vector, From Subdomain）を適切に選ばないと、非物理的な流れが発生することがある。

全ての流入・流出境界の正味の流量、つまり「ネット流量」がゼロに収束することを確認する。非定常なら時間平均値だ。具体的には、ソフトのモニタ機能で「Net Mass Flow Rate」や「Net Volume Flow Rate」をプロットし、その値が十分小さな値（例えば最大流量の0.1%以下）に振動しているかを見る。Abaqus/CFDでは、`History Output`で`CFNM`（質量流量）をリクエストできる。これがゼロから大きく外れていたら、質量保存が満たされておらず、計算結果は信用できない。

まず、カタログ値がどの状態を指すかを確認する。ポンプの性能曲線は通常、常温（20℃）の水のような特定の流体に対する「体積流量」で記載されている。これをそのままCFDの入力値にする。第一歩は、その流量値と吸込口の概算断面積から、平均流速を計算し、レイノルズ数を見積もることだ。例えば、流量Q=0.01 m³/s、管径D=0.1mなら、平均流速v≈1.27 m/s、水ならRe≈127,000となり、乱流であることがわかる。この情報が後々の乱流モデル選択やメッシュ設計の基礎になる。

非常に良い指標になる。特に、圧力損失（ΔP）が重要な場合、異なるメッシュ密度で計算した「流量-圧力損失特性」を比較する。例えば、3段階のメッシュ（粗い、中程度、細かい）で計算し、出口の質量流量が入口設定値に対してどれだけバランスしているか（ネット流量）、および目標とする圧力損失の値がメッシュを細かくしても±2%以内に収束しているかを確認する。自動車の排気管シミュレーションでは、バックプレッシャー（圧力損失）がエンジン性能に直結するため、この検証は必須だ。

実験データがなくても、簡易的な理論式や経験式と比較する「サニティチェック」はできる。例えば、円管内の層流（ハーゲン・ポアズイユ流れ）なら、理論的な流量と圧力損失の関係は

Fluentでは、境界条件タイプを「Velocity Inlet」や「Mass-Flow Inlet」などに設定した後、詳細設定で「Volume Flow Rate」を指定できる。ただし、これは非圧縮流れが前提だ。選択基準は明確だ：非圧縮流れ（液体や低速気体）→ 体積流量、圧縮性流れ（高速気体、温度変化大）→ 質量流量。例えば、家庭用給湯管の水流解析なら体積流量（例：15 L/min）でいい。一方、ターボチャージャー内部の流れでは、入口・出口で密度が大きく変わるので、質量流量（kg/s）を設定する。Fluentの「Mass-Flow Inlet」は、指定した質量流量を満たすようにソルバーが入口の総圧を調整する仕組みだ。

COMSOLの「流入」境界は物理場に依存するが、流体（CFDモジュール）では主に3つ：「流速」「圧力」「質量流量/体積流量」。体積流量を設定するには、「流入」境界を選択し、「境界条件」を「体積流量」に設定する。ここで注意すべきは、COMSOLは「流入方向」を別途指定する必要があることだ。デフォルトは「境界に垂直」だが、斜め流入の場合は「速度場」を指定する。また、COMSOLはマルチフィジックス結合が得意なので、例えば「微粒子の体積流量」を別の物理インターフェースで定義し、流体と連成させるような複雑な設定も可能だ。

根本的な数値手法（FVM）は同じだが、設定の「手動度」が高い。OpenFOAMでは、`0/`ディレクトリの境界条件ファイル（例：`U`, `p`）を直接編集する。体積流量を設定するには、`flowRateInletVelocity` という境界条件タイプを`U`ファイルで指定し、`volumetricFlowRate` の値を`constant/`ディレクトリなどで定義する。Salome-MECA（Code_Aster GUI）では、コマンド（`AFFE_CHAR_MECA`）を使って境界条件を定義する。いずれも、商用ソフトのGUI操作に比べて柔軟性は高いが、設定ミスのリスクも高く、理論的な理解が必須となる。産業界では再現性と効率から商用ソフトが使われることが多いが、研究開発ではOpenFOAMの柔軟性が活きる。

大いに関係ある。「Negative volume」はメッシュの極端な歪みが原因だが、その一因として、設定した流量に対し流路が細すぎて流速が異常に高くなり、セル内の圧力勾配が大きすぎることが挙げられる。まず、設定した流量から推定される平均流速を計算せよ。例えば、直径5mmの管に100 m³/hという非現実的な流量を設定していないか？ 一般的な目安として、管内流速は水で2-3 m/s、空気で10-30 m/s程度が妥当だ。これを大幅に超える設定は、数値的不安定性を招く。

以下のチェックリストを順に実行せよ。

考えられる原因は二つ。第一に、**流れのレジーム変化**だ。例えば、流量を上げていって臨界レイノルズ数を超えると、層流から乱流に遷移する。使用した乱流モデル（例えば、層流しか解けない設定や、標準k-εモデル）がその遷移を適切に捉えられていない可能性がある。第二に、**圧縮性効果の無視**だ。気体の場合、流速がマッハ数0.3を超えるあたりから圧縮性を考慮する必要がある。非圧縮ソルバーで高速流れを解こうとすると、エネルギー保存が破綻し、非物理的な結果が出る。対策は、流量を少しずつ段階的に増加させ、各段階の解を初期値として用いる「継続計算」を行い、どこで不安定化するかを特定することだ。

流体の密度次第で、影響は軽微から致命的まで幅広い。例えば、水（ρ≈1000 kg/m³）の解析で、体積流量1 m³/sを設定すべきところを、うっかり質量流量1 kg/sと設定してしまった場合、実際の流量は1/1000の0.001 m³/sになってしまう。流速や圧力損失は流量の2乗に比例する場合が多いので、結果は完全に無意味なものになる。逆に、低速空気（ρ≈1.2 kg/m³）なら、誤差は2割程度で気づかないかもしれない。発見したら、密度値と設定値を照合し、