本研究使用相位場法並結合熱流模式以模擬固化過程中氣孔生成之機制，及探討影響氣孔形狀的原因。模擬軟體為COMSOL之二維兩相流模組，其中以溫度和相位場之變數區分固、液、氣三相。制衡方程式包含動量、質量、能量以及濃度守恆方程式。模擬結果顯示為氣泡在固化過程中會因上冷卻效應之影響，造成氣孔形狀之差異。

This study use thermal and heat flow mode cooperating phase- field method to simulate 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 the temperature to distinguish solid and liquid phases, phase field function to distinguish liquid and gas phase. The computed results show that pore shape in solid is strongly influenced by cooling effect.

論文審定書 i
中文摘要 ii
Abstract iii
目錄 iv
圖目錄 vi
符號說明 vii
下標符號說明 ix
第一章緒論 1
1-1 研究背景及動機 1
1-2. 相位場法(PFM)及二相流(Two phase flow) 2
1-3. 本論文研究內容簡介 3
第二章 系統模型設定與理論之分析 4
2-1.研究方法、進行步驟及執行進度規劃: 4
2-2. 研究模擬之流程圖 5
2-3.模型架構 6
2-3-1.模型幾合 6
2-3-2.模型網格設置 7
2-3-3.邊界條件設定 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.模擬結果 18
3-2-1.氣泡拉長 18
3-2-2.固化速度差異 20
3-2-3.熱毛細力 22
3-3. 模擬結果驗證 24
3-3-1.網格加密及氣孔寬度 24
3-3-2.邊界拉長 27
第四章 討論 29
參考文獻 30

[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] 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.
[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] 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 ,2011, “The Transient Modeling of Single-Bubble Nucleate Boiling in a Sub-Cooled Liquid Using an ALE Moving Mesh” ,COMSOL Conference,Boston