利用一基於流體體積(VOF) 之介面追蹤法求解薄膜沸騰

本研究使用單一流體模型搭配CISIT 介面追蹤法求解雙相流之熱傳與質傳現象。其中所用的離散方法為有限體積法 (finite volume method)，而CISIT 介面追蹤法則是利用求解流體體積分率的傳輸方程式來追蹤介面之位置。利用介面處的質量與能量跳越條件求出相變化所造成的質量變化，並將其加入連續方程式中的源項作為修正。將動量方程式與連續方程式進行耦合，所得速度場將滿足介面處的質量跳躍函數。為了驗證相變化模型的準確性先求解具有理論解之相變化問題，發現溫度場之準確性與介面溫度梯度之計算都會影響模擬相變化的準確性。最後模擬一薄膜沸騰的氣泡生成，可觀察到氣泡的生成與成長過程中對於熱傳及質傳變化的影響。

關鍵字

雙相流；相變化；介面追蹤法；薄膜沸騰

並列摘要

A single-fluid model with a VOF based interface tracking method (CISIT) is used for simulation of two-fluid flows with heat and mass transfer in this study. The finite volume method is used as the discretiztion method. The motion of the interface is tracking by the solution of a transport equation for the volume fraction. By using the interface position obtained from the interface tracking method the heat flux which drives the phase change on interface can be calculated. The mass change across the phase boundary is computed by taking the mass and energy jump conditions at the interface. It is added as a source term in the continuity equation. This continuity equation is used in the pressure-velocity coupling algorithm. After the velocity field obtained from solving the momentum equation the jump condition for the continuity is satisfied. To assess the accuracy of this mode two test problems for which analytical solutions are available, are considered. It is shown that the precision of the temperature profile and the length used for calculating temperature gradients at the interface would affect the accuracy of simulation. Finally, a simulation of film boiling on a horizontal plate is presented. It shows that how the bubble growth affects the heat and mass transfer of the flow.

並列關鍵字

two-phase flow ； phase change ； interface tracking method ； film boiling

參考文獻

[1] Heyes, R. E., “Numerical simulation of mold filling in reaction injection molding”, Polym. Eng. Sci., v31, pp.842-848, 1991

[2] Muzaferija,S., Peric, M., “Computation of free-surface flows using the finite-volume method and moving girds”, Numerical Heat Transfer., Part B, v32, p.p. 369-384, 1997

[3] Trggvasson, G., “A front tracking method for the computations of multiphase flow”, J. Comput. Phys., v169, p.p. 708-759, 2001

[4] Osher, S., Sethian, J. A., “Fronts propagating with curvature-dependent speed: algorithms based on Hamilton Jacobi formulations”, J. Comput. Phys., v79, p.p. 12-49, 1988

[5] Li, X.L., “Study of three-dimensional Rayleigh-Taylor instability compressible fluids through level set method and parallel computation”, Phys. Fluids A, v5, p.p. 1904-1913, 1993

國際替代計量

利用一基於流體體積(VOF) 之介面追蹤法求解薄膜沸騰

全文下載

主題瀏覽