本研究是利用數值方法模擬在常溫常壓下的純水溶液底部板面施加二氧化碳通量，當物理性質在不變的情況下水溶液中的二氧化碳濃度達到過飽和時，則氣泡表面因為二氧化碳濃度梯及氣泡底部灌入二氧化碳通量的影響，氣泡在板面上某特定區域內成長之後脫離，並且在氣泡成長過程中氣泡是以固定接觸角和氣泡形狀模式的方式在板面上滑移。
在固定氣泡高度和氣泡接觸角下會有無限多組氣泡形狀模式，在數值計算中我們皆用同1種氣泡形狀模式，以避免在模擬氣泡成長過程中，氣泡在任何固定高度下其外型皆以不同之氣泡形狀模式成長。
氣泡底部的接觸面積和氣泡的體積是影響氣泡脫離板面的關鍵，
氣泡成長過成中可以得知氣泡底部的接觸面積和氣泡的體積會逐漸增加後在減小，氣泡底部的接觸面積達到最大時，氣泡的底部面積會開始快速內縮，在底部面積內縮期間氣泡體積卻維持一定，在相同氣泡體積下的氣泡形狀模式並不唯一，氣泡將固定體積下改變氣泡形狀模式，並且將氣泡的體積除以底面積後我們發現隨著氣泡的高度增加其斜率是以倍數增加，底部接觸面積最大值時的斜率是小於體積的最大值，故我們以底部接觸面積最大值時作為氣泡脫離板面的時機。
在氣泡成長數顆之後水溶液內的二氧化碳濃是不斷的升高，第1顆氣泡成長是需要比較長的時間，這是因為水溶液內的二氣化碳濃度並未完全受到水溶液底部二氧化碳通量的影響，但在成長數顆後單1氣泡成長時間是成線性增加。
在計算過程中改變水溶液半徑大小和底部二氧化碳通量大小來討論成長20顆後每1顆氣泡所花費的成長時間，若將水溶液半徑減小和增加水溶液底部二氧化碳通量皆有助於縮短單1氣泡成長時間。

[1] Bisperink, and C, C. J., Prins, A., 1994,” Bubble Growth in Carbonated Liquids,” Colloids Surf. A: Physicochem. Eng. Aspects, 85(2–3), 237,
[2] Frank, F. C., 1950, “ Radially Symmetric Phase Growth Controlled by Diffusion,” proceedings of the royal society of London series a-mathematical and physical sciences 201 : 568,
[3] Van Voorst Vader, F., and Erkens, F., and Van den Tempel, M., 1964,” Measurement of Dilatational Surface Properties,” Transactions of the faraday society 60 : 1170,
[4] Mori, B. K., Douglas, Baines W., 2001,” Bubble Departure from Cavities,” International Journal of Heat and Mass Transfer , Vol.44 , pp.771-783 ,
[5] Divinis, N., Karapantsios, T. D., Kostoglou, M., Panoutsos, C. S., Bontozoglou, V., Michels, A.C., Sneep, M.C., De, Bruijn R., Lotz, H., 2004,” Bubbles Growing in Supersaturated Solutions at Reduced Gravity,” AIChE J., 50:2369-2382,
[6] 江俊緯, 2004,” 靜態液滴與氣泡形狀之數值計算,” 國立清華大學碩士論文,
[7] Lee, S. L., Lee, H. D., 2004,” Evolution of The Liquid Meniscus in Capillary-Force-Dominant Flow,” Proceedings of International Symposium on Advances in Computational Heat Transfer, April 19-24, Norway, Paper Number: CHT-04-151,
[8] 楊宇宏, 2005,” 解壓對二氧化碳氣泡於過飽和水溶液中之消長現象研究,” 國立清華大學碩士論文,
[9] Cha, Y. S., Henry, R. E., 1981,” Bubble Growth During Decompression of A Liquild,”Transactions of the ASME, Vol. 103, pp.56-60,
[10] 劉建鴻,2006,”二氧化碳氣泡於多孔性壁面成長之研究,”國立清華大學碩士論文,
[11] Peter, G. T. Fogg, William, Gerrard., 1991,” Solubility of Gases in Liquids : A Critical Evaluation of Gas/Liquid Systems in Theory and Practice,”Wiley,New York,
[12] Lee, S. L., 1989,” Weighting Function Scheme and its Application on Multidimensional Conservation Equations,”Int. J. Heat Mass Transfer, 32, pp. 2065-2073,
[13] Lee, S. L., and Tzong, R. Y., 1992,” Artificial Pressure for Pressure-Linked Equation,”International Journal of Heat and Mass Transfer, Vol.35, pp.2705-2716,
[14] Lee, S. L., 1989,” A Strongly-Implicit Solver for Two-Dimensional Elliptic Differential Equations,”Numerical Heat Transfer, Vol.16, pp.161-178,
[15] Gerlach, D., Biswas, G.., Durst, F., Kolobaric, V., 2005,”Quasi-Static Bubble Formation on Submerged Orifices,” International Journal of Heat and Mass Transfer,pp.425-438,
[16] Lee, S. L., and Sheu, S. R., 2001, “A New Numerical Formulation for Incompressible Viscous Free Surface Flow without Smearing the Free Surface,” Int. J. Heat Mass Transfer, 44, pp. 1837-1848.
[17] Lee, S. L., and Sheu, S. R., 2003, “Filling Process in An Open Tank,” ASME J. of Fluids Engineering, 125, pp. 1016-1021.
[18] Lee, S. L., and Tzong, R. Y., 1991, “An Enthalpy Formulation for Phase Change Problems,” Int. J. Heat Mass Transfer, 34, pp. 1491-1502.
[19] 彭琪鈺, 2004,”氣泡於圓柱容器中上升之數值模擬研究,” 國立清華大學碩士論文,
[20] Hey, M. J., Hilton, A. M., 1994,” The Formation and Growth of Carbon Dioxide Gas Bubbles from Supersaturated Aqueous Solutions” Food Chemisitry, 349-357,