本研究以相位場法配合熱流模式模擬固化過程中氣泡生成之機制。模擬軟體為COMSOL之兩相流模組，其中以溫度和相位場函數區分固、液及氣三相。方程式包含動量、質量、能量以及濃度守恆方程式。模擬結果顯示氣孔在固化過程中會受入口流速及溫度冷卻效應之影響，造成氣孔形狀之差異。

This study applies the phase- field method to simulate the pore shapes in solid, and explore mechanisms responsible for pore formation. The simulation method is based on two-dimensional, two-phase flow module in the COMSOL software.
Conservation equations of mass, momentum, energy and concentration are solved in the entire domain by incorporating with temperature to distinguish solid and liquid phases, and phase field function to distinguish liquid and gas phases. The computed results show that pore shape in solid is strongly influenced by the heat convection and boundary velocity.

論文審定書 i
中文摘要 ii
Abstract iii
目錄 iv
圖目錄 vi
符號說明 viii
下標符號說明 xi
第一章 緒論 1
1-1研究背景 1
1-2相位場法(PFM)及二相流(Two phase flow) 1
1-3 研究內容簡介與架構 2
第二章 模型設定與理論之分析 3
2-1研究方法、進行步驟及執行進度規劃: 3
2-2研究模擬之流程圖 4
2-3模型架構 5
2-3-1 模型架構設定 5
2-3-2 網格分布與設定 5
2-3-3 初始值與邊界設定 7
2-3-4 流體性質 8
2-4模組之統御方程式 9
2-4-1 相位場法方程式 9
2-4-2質量及動量守恆方程式 11
2-4-3 能量方程式 15
2-4-4 濃度方程式 16
第三章 結果與討論 17
3-1 基本性質 17
3-2 氣泡拉長 19
3-2-1 密度圖 19
3-2-2 濃度圖 21
3-2-3 溫度圖 23
3-3 入口流速加大 25
3-3-1 密度圖 25
3-3-2 濃度圖 27
3-4 固化速度差異 29
3-4-1密度圖 29
3-5 網格驗證 31
第四章 結論與未來展望 33
參考文獻 34

[1] S. Kou, Welding Metallurgy. Wiley, New York, 1987.
[2] M. C. Flemings, Solidification Processing, McGraw-Hill, New York, 1974
[3] Y. Sun and C. Beckermann, 2010,“Phase-field modeling of bubble growth and flowin a Hele-shaw cell”, J. Heat and Mass Transfer, 2969-2978.
[4] David C. Venerus and Nadia Yala,1997,“Transport Analysis of Diffusion-Induced Bubble Growth and Collapse in Viscous
[5] Tanai L. Marin,“Solidification of a Liquid Metal Droplet Impinging on a Cold Surface”,Excerpt from the Proceedings of the COMSOL Users Conference 2006 Boston.
[6] V.R.Voller and C.Prakash,“ A fixed grid numerical Modelling Methodology for convection-diffusion mushy region phase-change problem”,Jourmal of Heat and Mass Transfer,30(8),1709-1719(1987).
[7] Luiz C. Wrobel and M. H. Aliabadi, 2003, The boundary element methods. Wiley, UK.
[8] Vittorio Cristini, Jerzy Bławzdziewicz, and Michael Loewenberg, 1998, “Drop breakup in three-dimensional viscous flows”, Phsics fluids, Vol. 10, pp.1781-1783
[9] Howard H. Hu, N. A. Patankar and M. Y. Zhu, 2000, “Direct Numerical Simulations of Fluid–Solid Systems Using the Arbitrary Lagrangian–Eulerian Technique”, Journal of Computational Physics 169, pp.427–462
[10] S. Ramaswamy and L.G. Leal, 1998, “ The deformation of a viscoelastic drop subjected to steady uniaxial extensional flow of a Newtonian fluid”, J. Non-Newtonian Fluid Mech., 85, pp.127-163
[11] Stanley Osher and Nikos Paragios, 2003, Geometric level set methods in imaging, vision, and graphics. Springer-Verlag. New York.
[12] Ruben Scardovelli and Stephane Zaleski,1999,“Direct numerical simulation of free-surface and interfacial flow”,Annu. Rev. Fluid Mech.,Vol.31,567-603.
[13] Y. Sun and C. Beckermann,2007,“Sharp interface tracking using the phase-field equation”, Journal of Computational Physics 220,pp.626-653.
[14] Cahn, J. W. and Hilliard, J. E, 1958, “Free energy of a nonuniform system. I. Interfacial free energy”, J. Chem. Phys. 28, pp.258-267.
[15] F. Kong, H. Zhang and G. Wang,2008,“Numerical Simulation of Transient Multiphase Field during Hybrid Plasma-Laser Deposition Manufacturing”J. Heat Transfer, Vol.130, NO.112101,pp.1-7.a
[16] Shyamprasad Karagadde, Suresh Sundarraj, Pradip Dutta,2012, “A model for growth and engulfment of gas microporosity during aluminum alloy solidification process”,Computational Materials Science 65 ,pp.383-394.
[17] comsol1998–2008, “Chemical Engineering MODULE”,BOILING WATER ,pp.155-156.
[18] PengtaoYue, JamesJ. Feng, Chun Liu and Jie Shen,2004, “A diffuse-interface method for simulating two-phase flows of complex fluids”, J. Fluid Mech. , Vol. 515,pp.293-317.
[19] S. Karagadde, S. Sundarrij,P.Dutta,2012: “A model for growth and engulfment of gas microporosity during aluminum alloy solidification process”Computational Meterials Science 65,pp.383-394
 
[20] Christopher J. Forster, Marc K. Smith ,“The Transient Modeling of Single-Bubble Nucleate Boiling in a Sub-Cooled Liquid Using an ALE Moving Mesh”, 2011 COMSOL Conference in Boston