[論文レビュー] Linearly decoupled energy-stable numerical methods for multi-component two-phase compressible flow

Jisheng Kou, Shuyu Sun+1 more|arXiv (Cornell University)|Dec 6, 2017

Computational Fluid Dynamics and Aerodynamics参考文献 2被引用数 3

ひとこと要約

本稿では、ペング・ロビンソン状態方程式（例：Peng-Robinson）を用いた多成分二相圧縮流体のための、線形で分離可能かつエネルギー安定な数値スキームを2つ提案する。中間速度と成分別スカラー補助変数（SAV）アプローチを導入することで、速度とモル密度の結合を解消し、各時間ステップで線形系のみを解くことにより、無条件的なエネルギー散逸を保証し、複雑な流体系の効率的で安定なシミュレーションを可能にする。

ABSTRACT

In this paper, for the first time we propose two linear, decoupled, energy-stable numerical schemes for multi-component two-phase compressible flow with a realistic equation of state (e.g. Peng-Robinson equation of state). The methods are constructed based on the scalar auxiliary variable (SAV) approaches for Helmholtz free energy and the intermediate velocities that are designed to decouple the tight relationship between velocity and molar densities. The intermediate velocities are also involved in the discrete momentum equation to ensure the consistency with the mass balance equations. Moreover, we propose a component-wise SAV approach for a multi-component fluid, which requires solving a sequence of linear, separate mass balance equations. We prove that the methods preserve the unconditional energy-dissipation feature. Numerical results are presented to verify the effectiveness of the proposed methods.

研究の動機と目的

• ペング・ロビンソンを含む現実的な状態方程式を用いた多成分二相圧縮流体のための効率的で、熱力学的に整合性があり、圧力方程式の構築を回避する数値法の開発。
• 質量保存および運動量保存方程式における速度とモル密度の強い結合性が数値スキームの複雑化を引き起こす問題の克服。
• 無条件的に離散エネルギー散逸法則を保持する線形で分離可能なスキームの構築により、数値的安定性を保証。
• 多成分流体における計算を簡素化しつつエネルギー安定性を維持できる成分別 SAV アプローチの導入。
• 離散的設定において、運動量方程式と質量保存方程式の整合性を保つために、中間速度を離散的定式化に組み込むこと。

提案手法

• 非線形ヘルムホルツ自由エネルギー密度を扱うためにスカラー補助変数（SAV）アプローチを用い、エネルギー安定スキームを線形系に変換する。
• 質量保存方程式の対流項における速度とモル密度の密接な結合性を解消するために、2つの中間速度を導入する。
• 中間速度を組み込んだ離散的運動量方程式を設計し、質量保存方程式と整合性を保ちつつ、運動エネルギーの変化を保存する。
• 各成分ごとに独立した線形質量保存方程式を解けるよう、成分別 SAV 定式化を提案し、計算効率を著しく向上させる。
• 圧力方程式を必要とせず、熱力学的に整合性のあるモデル化を可能にする、NVTに基づくモデルフレームワーク（モル数、体積、温度を主変数とする）を採用。
• 離散エネルギー関数を構築し、それが時間ステップごとに減少することを示すことにより、無条件的なエネルギー散逸を証明する。

実験結果

リサーチクエスチョン

• RQ1 現実的な状態方程式（例：ペング・ロビンソン）を用いた多成分二相圧縮流体のための、エネルギー安定的で線形的かつ分離可能な数値スキームをどのように構築できるか？
• RQ2 質量保存方程式の対流項における速度とモル密度の結合性を、エネルギー安定性を保ちつつ効果的に解消する方法は何か？
• RQ3 離散的設定において、質量保存方程式と運動量方程式の両者と整合性を持つように中間速度を設計できるか？
• RQ4 計算効率を向上させるために、SAV アプローチを多成分系に成分別に拡張する有効な方法は何か？
• RQ5 成分別 SAV アプローチを用いることで、圧縮性多成分流体に対する数値スキームの安定性および精度にどのような影響を与えるか？

主な発見

• 提示されたスキームは、各時間ステップで減少する離散エネルギー関数を用いて厳密に証明された通り、無条件的にエネルギー散逸的である。
• 各時間ステップで非線形結合ソルバーに比べて著しく計算効率が向上する、単一の線形方程式系の連鎖を解くことのみを要する。
• 成分別 SAV アプローチにより、ヘルムホルツ自由エネルギーを個々の成分寄与に分解することで、多成分流体の効率的取り扱いが可能になる。
• 数値結果により、長時間シミュレーションにおいても、流れのクイバー（渦）や速度大きさの等高線などの複雑な流れのパターンを効果的にシミュレートできることを示している。
• 中間速度の使用により、運動量方程式と質量保存方程式の整合性が保たれ、離散系の物理的整合性が維持されている。
• 実際のペング・ロビンソンパラメータを用いた3成分系（メタン、ペンタン、ドデカン）において、さまざまな時間ステップでスキームのロバスト性と安定性が検証された。