質量流量 — CAE用語解説

その理解で正しいです。定義式は

その通りです。気体では圧力・温度で密度が大きく変動するため、境界条件は質量流量で与える方が物理的に安定します。Ansys FluentやSiemens STAR-CCM+では、インレット境界条件として「Mass Flow Inlet」を直接指定できます。例えば、自動車エンジンの吸気ポート解析では、バルブリフトに応じた瞬時質量流量プロファイルを入力します。逆に、送風機の出口など、圧力がほぼ一定の場合は「Pressure Outlet」を設定し、質量流量は結果として算出させるのが一般的です。

微分形の質量保存則（連続の方程式）は

根本的な違いは、未知数として扱う変数です。速度境界条件では、速度ベクトルを固定値とし、圧力と密度をソルバーが調整して質量流量を決めます。一方、質量流量境界条件では、

いいえ、単純な算術平均ではありません。対流項の離散化スキームに依存します。例えば、1次風上差分なら、上流側セルの値を使って

大きく変わります。Pressure-Basedソルバー（低速～亜音速流れ向け）は、圧力補正法（SIMPLE, PISO）を用い、連続の式を圧力補正方程式の形で満たします。質量流量はこの補正プロセスの中心変数です。一方、Density-Basedソルバー（超音速流れ向け）は、密度、運動量、エネルギーの保存式をカップリングして解きます。こちらも質量保存は厳密に扱われますが、解法のアプローチが異なります。例えば、Fluentでマッハ数0.3以上の圧縮性流れを解く時、Density-Basedソルバーを選択すると、質量流量境界条件がより直接的に扱われ、衝撃波の解像度が向上することがあります。

実務では、プロジェクトや業界の慣習に合わせます。自動車業界ではg/sやkg/hがよく使われます。例えば、エンジン吸気量は数十から数百kg/hのオーダーです。航空宇宙ではlb/s（ポンド毎秒）も見かけます。重要なのは、ソフトウェア内で単位系を一貫させること。Ansys Workbenchでは「単位系」を設定すると、入力フィールドにその単位が表示されるので便利です。また、値が極端に小さい/大きいとソルバーのスケーリング問題を引き起こすことがあるので、例えば1e-5 kg/sではなく10 mg/sと入力するなど、オーダーを1前後に調整する工夫も有効です。

主要なCAEソフトには専用のモニター機能やレポート機能があります。Ansys Fluentなら、「Reports」→「Fluxes」で、選択した境界面（インレット、アウトレット、壁など）を通る質量流量の合計を算出できます。すべての流入面の合計と流出面の合計が、収束後の定常解析ではほぼ一致し（例えば差が1e-6 kg/s以下）、非定常解析ではその差が計算領域内の質量変化率と一致していることを確認します。また、ソルバーの収束履歴に「Net Mass Flow Rate」や「Mass Imbalance」のモニターを追加して、反復計算中にゼロに近づいていくことをリアルタイムで確認するのが標準的な手順です。

設定順序自体に大きな影響はありませんが、物理的に「過制約」や「未制約」にならないように注意が必要です。典型的な過制約は、すべての開口部に質量流量を設定すること。これでは圧力が決まらず、ソルバーが発散します。少なくとも1か所は圧力基準（Pressure Inlet/Outlet）が必要です。逆に、すべて圧力条件だと質量流量が決まらない未制約にはなりませんが、現実の駆動力が再現できない場合があります。例えば、換気扇のある部屋の解析では、ファン面に質量流量（または体積流量）を、外部開口部に圧力条件を設定するのが現実的です。

基本的な概念は同じですが、インターフェースと高度な機能に違いがあります。Fluentでは、境界条件ダイアログで「Mass Flow Rate」を直接入力し、「Direction Specification Method」で流入方向を定義します。STAR-CCM+では、「Mass Flow」を選択し、値と共に「Flow Direction」を「Normal to Boundary」や「Components」で指定します。STAR-CCM+の特徴は、「Field Function」を使って、位置や時間の関数として質量流量を柔軟に定義できる点です。例えば、ピストン運動に連動した正弦波状の質量流量プロファイルを、GUI上で数式として直接記述できます。FluentでもUDF（ユーザ定義関数）で同様のことは可能ですが、プログラミング知識が必要です。

COMSOLでは、物理場インターフェースごとに流入条件の名称が異なりますが、本質は同じです。「CFDモジュール」や「微視的流れ」では「流入」（Inlet）境界に「質量流量」を指定できます。「非等温流」や「反応流」を扱う場合も同様です。COMSOLの強みは、これに「固体力学」や「電磁界」を直接連成させられる点です。例えば、MEMSのマイクロポンプ解析で、ピエゾ駆動による膜変形（固体力学）で生じる体積変化を、流路内の質量流量源項として自動的に連成させることができます。この場合、ユーザが質量流量を直接指定するのではなく、連成物理から間接的に決まります。

OpenFOAMでは、境界条件は「0/」ディレクトリ以下の各場（U, p, Tなど）のファイルで、パッチごとに指定します。質量流量を直接指定する境界タイプは「flowRateInletVelocity」です。この設定では、`massFlowRate` または `volumetricFlowRate` のキーワードに値を書き、ソルバーが境界での速度を内部で調整します。例えば、`inlet { type flowRateInletVelocity; massFlowRate 0.01; // kg/s value uniform (0 0 0); }` と記述します。商用ソフトに比べて設定ファイルを直接編集する必要があり、初期学習コストは高いですが、カスタマイズ性は極めて高く、特殊な質量流量プロファイルもプログラミングで自由に実装できます。

「逆流」警告自体は、計算の初期段階や流れが複雑な領域ではよく起こり、必ずしもエラーではありません。しかし、定常計算で収束後も継続する場合は問題です。原因は、設定した質量流量に対して出口の圧力が高すぎる、または流路内に大きな障害物があって圧力損失が想定より大きいなどです。対策は、(1) 初期条件を現実に近づける（例えば、インレット近傍のセルに適切な速度を設定）、(2) ソルバー設定で「Pressure Outlet」の「Backflow Direction Specification Method」を「Normal to Boundary」から「From Neighboring Cell」などに変更する（Fluentの場合）、(3) 非定常計算として解く、などがあります。根本的には境界条件の設定が物理的に妥当か再検討が必要です。

はい、質量流量の値自体は正しくても、関連する他の設定が原因の場合があります。主な原因と対策は以下の通りです。
1. **メッシュ品質**: インレット面のメッシュが極端に粗い、または歪んでいる。界面での流量計算が不安定になります。面メッシュの品質を改善しましょう。
2. **密度の設定**: 材料の密度が間違っている、または理想気体則を使う場合に操作条件（圧力、温度）が非現実的。例えば、温度を300Kと設定すべきところを3Kと入力すると、密度が100倍になり、実質的な質量流量が想定より大幅に増加します。
3. **ソルバースキーム**: 高次精度スキームと大きな質量流量の組み合わせで発振することがあります。まずは1次風上差分で安定化させ、収束後にスキームを上げる「スキームの緩和」が有効です。
4. **タイムステップ（非定常）**: タイムステップが大きすぎて、1ステップで流入する質量がセルの容量を超える（Courant数が大きすぎる）。タイムステップを小さくします。

メッシュ依存性も一因ですが、より直接的な原因は「数値的な流れの分配」にあります。合流部の複雑な渦や逆流によって、設定通りの流量が各枝に分配されないことが起こり得ます。対策としては、
1. **合流部のメッシュを細かくする**: 特に、各流路が衝突するコーナー部分の解像度を上げます。
2. **初期条件を工夫する**: すべてのインレットに一様な速度（質量流量に比例した値）を初期条件として与え、計算の出だしを安定させます。
3. **ソルバー設定を調整する**: 圧力-速度連成のソルバー（SIMPLE, PISO）の緩和係数を小さく（保守的に）して、各反復での流量調整を穏やかにします。
4. **現実的な境界条件に変更**: 可能であれば、マニホールドの上流側の共通ダクトまでモデルに含め、そこに一つの質量流量を設定し、下流の分配は流れの抵抗（圧力損失）によって自然に決まるようにモデル化する方が、現実に近く、数値的にも安定することが多いです。