核心概念
本文提出了一種基於彈性速度的動力學模型,用於求解多組分歐拉方程,並著重於確保數值解的正性和準確捕捉穩態接觸不連續性。
文獻資訊
Shashi Shekhar Roya, S. V. Raghurama Raob. (2024). A Kinetic Scheme Based On Positivity Preservation For Multi-component Euler Equations. arXiv:2411.00285v1
研究目標
本研究旨在開發一種新的數值方法,用於求解多組分可壓縮流動的控制方程——多組分歐拉方程。其目標是解決現有數值方法在處理多組分流動時遇到的兩個主要挑戰:確保每個組分的密度保持正值,以及消除跨接觸不連續面的非物理壓力振盪。
方法
提出了一種基於彈性速度的動力學模型,並在向量動力學框架下實現。
該模型在 1D 中使用兩個速度,在 2D 中使用三個速度,這些速度與單元界面對齊,以確保宏觀法向通量局部一維化。
通過設定速度的大小來滿足一階精度下的正性守恆條件,並採用類似 CFL 的時間步長限制。
為了準確捕捉穩態接觸不連續性,對速度的定義進行了修改。
使用 Chakravarthy-Osher 類型的通量限制方法和強穩定性守恆 Runge-Kutta (SSPRK) 方法將基本數值方案擴展到三階精度。
主要發現
所提出的數值方案不需要計算 Roe 平均值,也不強烈依賴於底層的特征結構。
通過求解一系列 1D 和 2D 基准測試案例,包括衝擊波-氣泡交互作用測試案例,證明了該方案的穩健性、熵守恆性和準確性。
主要結論
本文提出的基於彈性速度的動力學模型為求解多組分歐拉方程提供了一種穩健且準確的數值方法。
該方案能夠有效地保持每個組分的密度和總壓力的正性,並能準確捕捉穩態接觸不連續性。
通過使用高階通量限制方法和時間積分方案,該方案可以實現高階精度。
意義
本研究為多組分可壓縮流動的數值模擬提供了一種有價值的工具,並為進一步開發更先進的數值方法奠定了基礎。
局限性和未來研究方向
本文僅考慮了理想氣體混合物,未來可以將該方法擴展到非理想氣體混合物。
未來的研究可以集中於對高階精度方案進行嚴格的正性分析。
可以進一步探索該方法在更複雜的流動問題中的應用,例如湍流和化學反應流動。