利用等位函数法模拟2D黏性驻波.PDF

利用等位函數法模擬 2D 黏性駐波 梁興杰 1 林家賢 2 林慶緯 3 摘要 本研究中，我們發展一數值模式，利用不可壓縮奈威爾-史托克方程式與介面運動方程式， 以求解自由液面流流場；將自由液面流轉換為固定計算網格之二相流系統，氣相與液相以等位 函數之「正」、「負」符號區分，零等位函數表示自由液面，氣相與液相之密度與黏滯係數以平 滑之等位函數描述之。氣相與液相之流場統一以壓力校正法 SIMPLE 求解；介面運動以等位函 數法捕捉。我們將本模式應用於模擬 2D 黏性駐波，以雷諾數 Re = 20 與 200 為算例，並與線性 解(Wu, et al., 2001 ）做比較。此外，我們探討Re = 200 時，不同氣相與液相密度比與黏滯係數 比對自由液面流流場，與質量守恆誤差之影響。 關鍵字：壓力校正法、等位函數法、黏性駐波 Numerical Study of Two-Dimensional Viscous Standing Wave using Level Set Method 1 2 3 Shin-Jye Liang Jia-Sian Lin Cing-Wei Lin ABSTRACT A two-dimensional numerical algorithm that couples the incompressible Navier-Stokes equations with the Hamilton-Jacobi equation for free-surface flows was developed. The free surface flow is converted into a two-phase flow system on a fixed grid in which the free surface is implicitly expressed by the zero level set function. Pressure correction method based SIMPLE was used to solve the two-phase flow filed. Level set method was employed to capture the evolution of interface. The developed model was applied to simulate 2D viscosity standing waves. The effect of density ratio and viscosity of air and liquid on the free surface flows were investigated. Computed results agreed well with the linear solution which considers the interaction of free surface and viscosity (Wu, et al., 2001). Key Words ：Pressure correction method, Level Set Method, standing wave ㄧ、緒論 法(Arbitrary Lagrangian – Eulerian Method, ALE) 自由液面流與二相流之數值研究在過去二、三 (Hirt,