×

Diffuse interface model for incompressible two-phase flows with large density ratios. (English) Zbl 1388.76403

Summary: We investigate the applicability of an incompressible diffuse interface model for two-phase incompressible fluid flows with large viscosity and density contrasts. Diffuse-interface models have been used previously primarily for density-matched fluids, and it remains unclear to what extent such models can be used for fluids of different density, thereby potentially limiting the application of these models. In this paper, the convective Cahn-Hilliard equation and the condition that the velocity field is divergence-free are derived from the conservation law of mass of binary mixtures in a straightforward way, for fluids with large density and viscosity ratios. Differences in the equations of motion with a previously derived quasi-incompressible model are shown to result from the respective assumptions made regarding the relationship between the diffuse fluxes of two species. The convergence properties of the model are investigated for cases with large density ratio. Quantitative comparisons are made with results from previous studies to validate the model and its numerical implementation. Tests show that the variation in volume during the computation is of the order of machine accuracy, which is consistent with our use of a conservative discretization scheme (finite volume methods) for the Cahn-Hilliard equation. Results of the method are compared with previous work for the change in topology of rising bubbles and Rayleigh-Taylor instability. Additional results are presented for head-on droplet collision and the onset of droplet entrainment in stratified flows.

MSC:

76T10 Liquid-gas two-phase flows, bubbly flows

Software:

SOLA-VOF
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Hirt, C. W.; Nichols, B. D., Volume of fluid (VOF) method for the dynamics of free boundaries, J. Comput. Phys., 39, 201-225 (1981) · Zbl 0462.76020
[2] B.D. Nichols, C.W. Hirt, R.S. Hotchkiss, SOLA-VOF: A solution algorithm for transient fluid flow with multiple free boundaries, Los Alamos National Lab Report LA-8355, 1980.; B.D. Nichols, C.W. Hirt, R.S. Hotchkiss, SOLA-VOF: A solution algorithm for transient fluid flow with multiple free boundaries, Los Alamos National Lab Report LA-8355, 1980.
[3] Sussman, M.; Smereka, P.; Osher, S., A level set approach for computing solutions to incompressible two-phase flow, J. Comput. Phys., 114, 146-159 (1994) · Zbl 0808.76077
[4] Adalsteinsson, D.; Sethian, J., A fast level set method for propagating interfaces, J. Comput. Phys., 118, 269-277 (1995) · Zbl 0823.65137
[5] Tryggvason, G.; Bunner, B.; Esmaeeli, A., A front-tracking method for the computations of multiphase flow, J. Comput. Phys., 169, 708-759 (2001) · Zbl 1047.76574
[6] Jacqmin, D., Calculation of two-phase Navier-Stokes flows using phase-field modelling, J. Comput. Phys., 155, 96-127 (1999) · Zbl 0966.76060
[7] Anderson, D. M.; McFadden, G. B.; A Wheeler, A., Diffuse-interface methods in fluid mechanics, Ann. Rev. Fluid Mech., 30, 139-165 (1998) · Zbl 1398.76051
[8] Badalassi, V. E.; Ceniceros, H. D.; Banerjee, S., Computation of multiphase systems with phase field models, J. Comput. Phys., 190, 371-397 (2003) · Zbl 1076.76517
[9] J.D. van der Waals, Verhandel. Konink. Akad. Weten. AmsterdamSect vol. 1 No. 8 (Dutch) 1893, Translation of J.D. van der Waals (The thermodynamic theory of capillarity under the hypothesis of a continuous density variation) J. Stat. Phys. (1979).; J.D. van der Waals, Verhandel. Konink. Akad. Weten. AmsterdamSect vol. 1 No. 8 (Dutch) 1893, Translation of J.D. van der Waals (The thermodynamic theory of capillarity under the hypothesis of a continuous density variation) J. Stat. Phys. (1979).
[10] Biben, T.; Kassner, K.; Misbah, C., Phase-field approach to three-dimensional vesicle dynamics, Phys. Rev. E, 72, 041921 (2005)
[11] Folch, R.; Casademunt, J.; Hernández-Machado, A.; Ramirez-Piscina, L., Phase-field model for Hele-Shaw flows with arbitrary viscosity contrast. I. Theoretical approach, Phys. Rev. E, 60, 1724 (1999)
[12] Yue, P.; Feng, J. J.; Liu, C.; Shen, J., A diffuse-interface method for simulating two-phase flows of complex fluids, J. Fluid Mech., 515, 293-317 (2004) · Zbl 1130.76437
[13] Seppecher, P., Moving contact lines in the Cahn-Hilliard theory, Int. J. Engng. Sci., 34, 977-992 (1996) · Zbl 0899.76042
[14] Jacqmin, D., Contact-line dynamics of a diffuse fluid interface, J. Fluid Mech., 402, 57-88 (2000) · Zbl 0984.76084
[15] Gurtin, M. E.; Polignone, D.; Vinals, J., Two-phase binary fluids and immiscible fluids described by an order parameter, Math. Models Meth. Appl. Sci., 6, 815-831 (1996) · Zbl 0857.76008
[16] Ding, H.; Spelt, P. D.M., Inertial effects in droplet spreading: a comparison between diffuse interface and level-set simulation, J. Fluid Mech., 576, 287-296 (2007) · Zbl 1125.76024
[17] Kim, J., A continuous surface tension force formulation for diffuse-interface models, J. Comput. Phys., 204, 784-804 (2005) · Zbl 1329.76103
[18] Jacqmin, D., Onset of wetting failure in liquid-liquid systems, J. Fluid Mech., 517, 209-228 (2004) · Zbl 1131.76316
[19] Antanovskii, L. K., A phase field model of capillary, Phys. Fluids, 7, 4, 747 (1995) · Zbl 1039.76502
[20] Lowengrub, J.; Truskinovsky, L., Quasi-incompressible Cahn-Hilliard fluids, Proc. R. Soc. London Ser. A, 454, 2617-2654 (1998) · Zbl 0927.76007
[21] Boyer, F., A theoretical and numerical model for the study of incompressible mixture flows, Comp. Fluids, 31, 41-68 (2002) · Zbl 1057.76060
[22] Cahn, J. W.; Hilliard, J. E., Free energy of a nonuniform system I, J. Chem. Phys., 28, 258 (1958) · Zbl 1431.35066
[23] Cahn, J. W.; Hilliard, J. E., Free energy of a nonuniform system III, J. Chem. Phys., 31, 688 (1959)
[24] Jacqmin, D., An energy approach to the continuum surface tension method, AIAA96-0858, (Proceedings of the 34th Aerospace Sciences Meeting and Exhibit (1996), American Institute of Aeronautics and Astronautics: American Institute of Aeronautics and Astronautics Reno)
[25] Liu, X. D.; Osher, S.; Chan, T., Weighted essentially non-oscillatory schemes, J. Comput. Phys., 115, 200-212 (1994) · Zbl 0811.65076
[26] Tryggvason, G., Numerical simulations of the Rayleigh-Taylor instability, J. Comput. Phys., 75, 253-282 (1988) · Zbl 0638.76056
[27] Guermond, J. L.; Quartapelle, L., A projection FEM for variable density incompressible flows, J. Comput. Phys., 165, 167-188 (2000) · Zbl 0994.76051
[28] Daly, B. J., Numerical study of two fluid Rayleigh-Taylor instabilities, Phys. Fluids, 10, 297 (1967) · Zbl 0147.45806
[29] He, X. Y.; Chen, S. Y.; Zhang, R. Y., A lattice Boltzmann scheme for incompressible multiphase flow and its application in simulation of Rayleigh-Taylor instability, J. Comput. Phys., 152, 642-663 (1999) · Zbl 0954.76076
[30] Bonometti, T.; Magnaudet, J., Transition from spherical cap to toroidal bubbles, Phys. Fluids, 18, 052102 (2006)
[31] Sussman, M.; Smereka, P., Axisymmetric free boundary problems, J. Fluid Mech., 341, 269-294 (1997) · Zbl 0892.76090
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.