摘要
本文提出一數值方法於自由液面流 (free surface flow)。使得接觸角 (contact angle) 對自由液面形狀的影響得以被考慮。文中利用Lee and Sheu [13] 計算不可壓縮黏性自由液面流之架構進行噴泉效應流場 (fountain effect flow) 計算。使用extended weighting function scheme與NAPPLE法則配合SIS求解流場。利用調和函數 (harmonic function) 外差並內差出自由液面的遷移速度，然後追蹤液面上的標記以決定自由液面的形狀。為了使自由液面滿足動態接觸角，本文引用Lee and Lee [16] 所提出的解析幾何方法 (analytic geometry method)，依據自由液面上的力平衡方程式與自由液面液體側壓力，由自由液面接和點修正自由液面，使液面於固體邊界滿足動態接觸角。過程中，採用一natural cubic spline smoothing method [24] 以準確計算力平衡方程式中的表面張力項。最後將此一數值方法用於噴泉效應流場的計算。

參考文獻
[1] Rose, W., 1961, “Fluid-fluid interface in steady motion,” Nature, Vol.191, pp.242-243.
[2] Rose, W. and Heins, R. W., 1962, “Moving interfaces and contact angle rate-dependency,” Journal of Colloid Science, Vol.17, pp.39-48.
[3] Brown, C. E., Jones, T. J. and Neustadter, E. L., 1980, “Interfacial flow during immiscible displacement,” Journal of Colloid and Interface Science, Vol.76, pp.582-586.
[4] Dussan, E. B., 1977, “Immiscible liquid displacement in a capillary tube: the moving contact line,” AIChE Journal , Vol.23, pp.131-133.
[5] Coyle, D. J., Blake, J. W., and Macosko, C. W., 1987, “The kinematics of fountain flow in mold filling,” AIChE Journal, Vol.33, pp.1168-1177.
[6] Mavridis, H., Hrymak, A. N., and Vlachopoulos, J., 1986, “Finite element simulation of fountain flow in injection molding,” Polymer Engineering and Science, Vol.26, pp.449-454.
[7] Lee, H. S., 1997, “Finite element analysis for the flow characteristics along the thickness direction in injection molding,” Polymer Engineering and Science, Vol.37, pp.559-567.
[8] Gogos, C. G., Huang, C. F., and Schmidt, L. R., 1986, “The process of cavity filling including the fountain flow in injection molding,” Polymer Engineering and Science, Vol.26, pp.1457-1466.
[9] Kamal, M. R., Chu, E., and Lafleur, P. G., 1986, “Computer simulation of injection mold filling for viscoelastic melts with fountain flow,” Polymer Engineering and Science, Vol.26, pp.190-196.
[10] Kamal, M. R., Goyal, K. S., and Chu, E., 1988, “Simulation of injection mold filling of viscoelastic polymer with fountain flow,” AIChE Journal, Vol.34, pp.94-106.
[11] Behrens, R. A., Crochet, M. J., Denson, C. D., and Metzner, A. B., 1987, “Transient free-surface flows: motion of a fluid advancing in a tube,” AIChE Journal, Vol.33, pp.1178.
[12] Kim, J. M., Ahn, K. H., Lee, S. J., and Lee, S. J., 2001, “Numerical simulation of moving free surface problems in polymer processing using volume-of-fluid method,” Polymer Engineering and Science, Vol.41, pp.858-866.
[13] Lee, S. L., and Sheu, S. R., 2001, “A new numerical formulation for incompressible viscous free surface flow without smearing the free surface,” International Journal of Heat and Mass Transfer, Vol.44, pp.1831-1848.
[14] Ngan, C. G. and Dussan V., E. B., 1984, “The moving contact line with a advancing contact angle,” Phys. Fluids, Vol.27, pp.2785-2787.
[15] Mavridis, H., Hrymak, A. N., and Vlachopoulos, J., 1988, “Transient free-surface flows in injection mold filling,” AIChE Journal, Vol.34, pp.403-410.
[16] Lee, S. L. and Lee, H. D., 2003, “Evolution of the liquid meniscus in a capillary-force-dominant flow,” The International Symposium on Transport Phenomena.
[17] Sarpkaya, T., 1996, “Vorticity, free surface and surfactants,” Annual Review of Fluid Mechanics, Vol.28, pp.83-128.
[18] Tsai, W. T., and Yue, D. K. P., 1996, “Computation of nonlinear free-surface flows,” Annual Review of Fluid Mechanics, Vol.28, pp.249-278.
[19] Lee, S. L., and Tzong, R. Y., 1992, “An enthalpy formulation for phase change problems with a large thermal diffusivity jump across the interface,” International Journal of Heat and Mass Transfer, Vol.35, pp.2705-2716.
[20] Lee, S. L., 1989, “Weighting function scheme and its application on multidimensional conservation equation,” International Journal of Heat and Mass Transfer, Vol.32, pp.2065-2073.
[21] 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.
[22] Lee, S. L., 1989, “A strongly-implicit solver for two-dimensional elliptic differential equations,” Numerical Heat Transfer, Vol.16, pp.161-178.
[23] Kistler, S. F., 1993, “Hydrodynamics of wetting,” in Wettability, Berg, J. C., ed., Chapter 6, pp.311-429, Marcel Dekker, New York.
[24] Lancaster, P. and Salkauskas, K., 1986, K., Curve and Surface Fitting: an introduction, Chapter 4, pp.101-111, Academic Press, London.