[論文レビュー] A model for unsteady mixed flows in non uniform closed water pipes and a well-balanced finite volume scheme

Christian Bourdarias, Mehmet Ersoy+1 more|arXiv (Cornell University)|Nov 29, 2008

Computational Fluid Dynamics and Aerodynamics参考文献 25被引用数 5

ひとこと要約

本稿では、非一様な閉管内の非定常混合流に適した新しいPFSモデルを提案する。自由水面流と圧力流の力学を、源項を含む統一的定式化によって統合する。曲線座標系と漸近解析を用いて、VFRoeスキームから派生させた、静水状態を正確に保持するWell-balancedな有限体積スキームを導出する。このスキームは、自由水面流と圧力流の間の遷移、特に水ハンマーおよび圧力低下現象を正確に捉えることができる。

ABSTRACT

We present the derivation of a new unidirectional model for We present the derivation of a new unidirectional model for unsteady mixed flows in non uniform closed water pipes. We introduce a local reference frame to take into account the local perturbation caused by the changes of section and slope. Then an asymptotic analysis is performed to obtain a model for free surface flows and another one for pressurized flows. By coupling these models through the transition points by the use of a common set of variables and a suitable pressure law, we obtain a simple formulation called PFS-model close to the shallow water equations with source terms. It takes into account the changes of section and the slope variation in a continuous way through transition points. Transition point between the two types of flows is treated as a free boundary associated to a discontinuity of the gradient of pressure. The numerical simulation is performed by making use of a Roe-like finite volume scheme that we adapted to take into account geometrical source terms in the convection matrix. Finally some numerical tests are presented.

研究の動機と目的

• 非一様な閉管内の非定常混合流を記述する統一的な数学的モデルの構築を目的とする。自由水面流および圧力流の両流動状態を含む。
• プレイスマンスロット法や非保存的Allievi方程式といった既存モデルの限界を克服し、共通の変数セットを用いて自由水面流と圧力流を自然に結合できるようにすることを目的とする。
• 曲線座標変換と漸近解析を用いて、局所的な断面および勾配の変化に起因する幾何的効果を連続的にモデルに組み込むこと。
• 複雑な幾何形状と源項を有する非一様管においても、静水状態の定常状態を正確に保持するWell-balancedな有限体積スキームの設計を目的とする。
• 非一様管（収縮、一様、拡大）における一時的境界条件下での数値的妥当性を確認する。

提案手法

• パイプの曲率と勾配を考慮するため、局所的な曲線座標系（Serret-Frenet枠）を用いて、3次元非圧縮性Euler方程式から自由水面（FS）モデルを導出する。
• 漸近解析を適用し、断面および勾配の変化に起因する局所的効果を捉える一次近似を得る。
• 同様の導出手法を3次元圧縮性Euler方程式に適用し、水の体積弾性率と音速を考慮した圧力流（P）モデルを導出する。
• 共通の保存変数セットと連続的遷移を可能にする圧力則を用いて、PおよびFSモデルを統合し、単一のPFSモデルを構築する。
• VFRoe有限体積スキームを改良し、幾何的源項を対流行列に組み込むことで、静水状態の正確なバランスを達成する。
• 数値フラックスにおける源項を処理するために界面上向き補間（interfacial upwinding）を用い、一時的シミュレーションにおける安定性と精度を確保する。

実験結果

リサーチクエスチョン

• RQ1 断面および勾配が変化する非一様な閉管内における自由水面流と圧力流を、連続的かつ一貫したモデルで記述する方法は何か？
• RQ2 局所的な幾何的変化（断面および勾配の変化）は、自由水面流と圧力流の間の遷移ダイナミクスにどのような影響を及えるか？
• RQ3 複雑な源項を有する非一様管においても、静水状態の定常状態を正確に保持する有限体積スキームを構築可能か？
• RQ4 提案されたPFSモデルは、異なる径の管における水ハンマーや圧力低下といった一時的現象を、どの程度正確に捉えることができるか？
• RQ5 古典的なRoe型スキームは、この文脈でどのような限界を示し、対流行列の修正によってどのように克服できるか？

主な発見

• PFSモデルは、一貫した圧力則を用いた統一的定式化により、自由水面流と圧力流を効果的に結合し、流動状態の連続的遷移を可能にする。
• 数値シミュレーションでは、上流の圧力頭の急上昇が圧力波を引き起こすことが確認され、収縮管では約19秒、一様管では約24秒、拡大管では約28秒に圧力低下が発生する。
• 変数$ A / extbf{S} $の変化が、予測通りの時期に圧力低下状態の発生を示し、モデルの臨界流動遷移の検出能力を裏付けている。
• 古典的なRoe行列近似では、非一様管内では静水状態が正確に保持されないが、対応する対流行列の修正により、再びWell-balancingが回復される。
• 一次元空間離散化による滑らか化効果が、特に不連続付近で顕著に現れるが、全体的な物理的一致性は維持されている。
• モデルは水ハンマーおよび圧力波の伝播を正確に捉えており、異なる管の幾何形状において一貫性のある結果を示しているが、非一様系の実験的妥当性は限定的である。