The aim of this study is to calculate two fluids flow and develop a robust, accurate and simple scheme by Volume-of-Fluid method. It is based on the solution of an indicator equation which is chosen to be the volume fraction. After a simple test and compared with the well-known HRIC, CICSAM and STACS, a composite scheme which uses the combination of SUPERBEE and MUSCL is adopted. To avoid the divergence and maintain high accuracy in the calculation of velocity field, the convective fluxes of momentum equation uses Van Leer scheme. The numerical results are demonstrated and show that this CFD methodology is well performed.
In addition, because of the investigations of two-fluids flow have many applications, some further improvements of this CFD methodology can be made in the future work. First, the two-fluids simulation can be applied in three-dimensional or axial symmetry for real cases.
Besides, the analysis of heat transfer with two immiscible fluids is popular subject in industry, such as the soldering of IC package process, electroplate engineering and metal casting. The heat transfer calculation can be acceded into the numerical simulated system.
Reference
[1] Jyi-Tyan Yeh, “A VOF-FEM and coupled inkjet simulation”, 2001 ASME Fluid Engineering Division Summer Meeting, New Orleans, Louisiana, 2001
[2] Kamisuki, S. et al., “A low power, small, electrostatically-driven commercial inkjet head”, Proc. IEEE MEMS Workshop, Heidelberg, Germany, 25, 1998
[3] Liou et al., “Three-dimensional simulations of the droplet formulation during the inkjet printing process”, Int. Comm. Heat Mass Transfer, Vol.29, No.8, pp. 1009-1118,2002.”
[4] Hirt, C. W. and Nichols, B.D., “Volume of fluid (VOF) method for the dynamics of free boundaries”, J. Comput. Phys., Vol. 39, pp. 201-225, 1981.
[5] Feziger, J. H. and Peric, M., “Computational methods for fluid dynamics”, Heidelberg:
Springer, 1996.
[ 6 ] Daly, B. J., “A technique for including surface tension effects in hydrodynamic calculations,” J. Comput. Phys., Vol. 4, pp.97-117, 1969.
[7] Nichols, B. D. and Hirt, C. W., “Calculating three-dimensional free surface flows in the vicinity of submerged and exposed structures,” J. Comput. Phys., Vol. 12, pp, 234-246, 1973.
[8] Osher, S. and Sethian, J. A., “Front propagating with curvature-dependant speed:
algorithms based on Hamilton-Jacobi formulations,” J. Comput. Phys., Vol.79, pp.12-49, 1988.
[9] Harlow, F. H. and Welch, J. E., “Numerical calculation of time-dependent viscous incompressible flow of fluid with free surface”, Phys. Fluids, Vol. 8, No.12, pp.
2182-2189, 1965.
[10] Chen, S., Johnson, D. B. and Raad, P. E., “Velocity boundary conditions for the simulation of free surface fluid flow”, J. Comput. Phys., Vol.116, No.2, pp. 262-276, 1995.
[11] Noh, W. F. and Woodward, P. R., “SLIC (Simple Line Interface Method)”, Lecture Note in Physics, Vol. 59, pp. 330-340, 1976
[12] Youngs, D. L., “Time-dependent multi-materiel flow with large fluid distortion”, in K. W.
Morton and M. J. Baines (eds), Numerical Methods for Fluid Dynamics, Academic, New York, pp. 273-285, 1982.
[13] Rider, W. J. and Kothe, D. B., “Reconstructing volume tracking,” J. Comput. Phys., Vol.
141, pp. 112-152, 1998.
[14] Ramshaw, J. D. and Trapp, J. A. ”A numerical technique for low-speed homogeneous two-phase flow with sharp interfaces”, J. Comput. Phys., Vol. 21, pp. 438-453, 1976.
[15] Shyy, W. , Thakur, S. and Wright, J., “Second-order upwind and central difference schemes for recirculating flow computation”, AIAA Journal, vol.30, pp. 923-932, 1992.
[16] Leonard, B. P., “A stable and accurate convective modeling procedure based on quadratic interpolation”, Computer Methods in Applied Mechanics and Engineering, Vol.19, pp.56-98, 1979.
[17] Zalesak, S. T. , “Fully multi-dimensional flux corrected transport algorithms for fluid flow”, J. Comput. Phys., Vol.31, pp.335-362, 1979.
[18] Rudman, M., “Volume-tracking methods for interfacial flow calculations”, Int. J. Numer.
Methods Fluids. Vol. 24, pp.671-691, 1997.
[19] Leonard, B. P., “Simple high-accuracy resolution program for convective modeling of discontinuities”, In. J. Numer. Methods Fluids, Vol. 8, pp1291-1318, 1988
[20] Harten, A. “High resolution schemes for hyperbolic conservation laws”, J. Comput.
Phys., vol.49, pp. 357-393, 1983.
[21] Van Leer, B., “Towards the ultimate conservative difference scheme. . Monotonicity Ⅱ and conservation combined in a second-order scheme”, J. Comput. Phys., Vol.14, pp.361-370 1974.
[22] Sweby, P. K., “High resolution schemes using flux limiters for hyperbolic conservation laws”, SIAM J. Numer. Analysis., Vol.21, No.5, pp. 995-1011, 1984.
[23] Gaskell, P. H. and Lau, A. K. C., ”Curvature-compensated convective transport: SMART, a new boundedness-preserving transport algorithm.” Int. J. Numer. Meth. Fluids, vol. 8, pp. 617-641, 1988.
[24] Jasak, H., Weller, H. G. and Gosman, A. D., “High Resolution NVD Differencing Scheme for Arbitrary Unstructured Meshes”, Int. J. for Numer. Meth. Fluids, vol. 31, pp.431-449, 1999.
[25] Roe, P. L., “Large scale computations in fluid mechanics, Part2”, In Lectures in Applied
Mathematics, Vol. 22, pp. 163-193, 1985.
[26] Darwish, M. S., “A new high-resolution scheme based on the normalized variable formulation”, Numer. Heat Trans, Part B, Vol. 24, pp. 353-371, 1993.
[27] Van Leer, B., “Towards the Ultimate Conservative Difference Scheme V. A Second-Order Sequel to Godunov's Method”, J. Comput. Phys., vol. 135, No.2, pp.
229-248, 1997.
[28] Muzaferija, S., Peric, M., Sames, P. and Schellin, T., “A Two-Fluid Navier-Stokes Solver to Simulate Water Entry”, Twenty-Second Symp. On Naval Hydrodynamics, Washington, D.C., pp.638-649, 1998.
[29] Darwish, M. and Moukalled, F., “Convective schemes for capturing interfaces of free-surface flows on unstructured grids.”, Numer. Heat Trans., Part B, 49, pp.19-42, 2006 [30] Ubbink, O., “Numerical prediction of two fluid systems with sharp interfaces,” thesis
submitted from Department of Mechanical Engineering Imperial College of Science, Technology & Medicine, January 1997.
[31] Issa, R. I., “Solution of the implicitly discretised fluid flow equations by operator-splitting,” J. Comput. Phys, Vol.62, pp. 40-65, 1985.
[32] Brackbill et al., “A continuum method for modeling surface tension”, J. Comput. Phys, vol.100, pp. 335-354, 1992.
[33] Batchelor, G. K., ”An introduction to fluid dynamics,” Cambridge: Cambridge University Press, 1967.
[34] Adamson, A. W., “Physical chemistry of surfaces,” New York: John Wiley & Sons, 1976.
[35] Orlanski, "A simple boundary condition for unbounded hyperbolic flows", J.
Computational Physics, vol.21,pp.251-269,1976
[36] Ubbink, O. and Issa, R. I., “A method for capturing sharp fluid interfaces on arbitrary meshes,” J. Comput. Phys., Vol. 153, pp.26-50, 1999.
[37] Lafaurie, B., et al., “Modeling merging and fragmentation in multiphase flows with SURFER,” J. Comput. Phys., Vol. 113, pp.134-147,1994
[38] Puckett, E. G., et al., “A high-order projection method for tracking fluid interfaces in
variable density incompressible flows,” J. Comput. Phys., Vol. 130, pp. 269-282, 1997.
[39] Raad, P. E., Chen, S. and Johnson, D. B., “The introduction of micro cells to treat pressure in free surface fluid flow problem,” J. Fluids Engineering, Vol. 117, pp.683-690, 1995.
[40] Martin, J. C. and Moyce, W. J., “An experimental study of the collapse of liquid columns on a rigid horizontal plane,” Philos. Trans. Roy. Soc. London, Vol. A244, pp.312-324, 1952.
[41] Koshizuka, S., Tamako, H. and Oka, Y., “A particle method for incompressible viscous flow with fluid fragmentation,” Comput. Fluid Dynamics JOURNAL, Vol.4(1), pp.29-46, 1995.
[42] Stoker, J. J., “Water waves,” New York: John Wiley & Sons, 1958.