An Improved Coupled Level Set and Continuous Moment-of-Fluid Method for Simulating Multiphase Flows with Phase Change

Abstract

An improved algorithm for computing multiphase flows is presented in which the multimaterial Moment-of-Fluid (MOF) algorithm for multiphase flows, initially described by Li et al. (2015), is enhanced addressing existing MOF difficulties in computing solutions to problems in which surface tension forces are crucial for understanding salient flow mechanisms. The Continuous MOF (CMOF) method is motivated in this article. The CMOF reconstruction method inherently removes the “checkerboard instability” that persists when using the MOF method on surface tension driven multiphase (multimaterial) flows. The CMOF reconstruction algorithm is accelerated by coupling the CMOF method to the level set method and coupling the CMOF method to a decision tree machine learning (ML) algorithm. Multiphase flow examples are shown in the two-dimensional (2D), three-dimensional (3D) axisymmetric “RZ”, and 3D coordinate systems. Examples include two material and three material multiphase flows: bubble formation, the impingement of a liquid jet on a gas bubble in a cryogenic fuel tank, freezing, and liquid lens dynamics.

Keywords

Moment-of-Fluid (MOF) / Surface tension / Two phase flow / Phase change / Deforming boundaries with change(s) in topology / Two-dimensional (2D) / Three-dimensional (3D) axisymmetric / 3D

References

[1.] Ahn HT, Shashkov M. Multi-material interface reconstruction on generalized polyhedral meshes. J. Comput. Phys., 2007, 226(2): 2096-2132

[2.] Ahn HT, Shashkov M. Adaptive moment-of-fluid method. J. Comput. Phys., 2009, 228(8): 2792-2821

[3.] Ancellin M, Després B, Jaouen S. Extension of generic two-component vof interface advection schemes to an arbitrary number of components. J Comput Phys., 2022, 473

[4.] Arienti M, Sussman M. An embedded level set method for sharp-interface multiphase simulations of diesel injectors. Int. J. Multiph. Flow, 2014, 59: 1-14

[5.] Asuri Mukundan A, Ménard T, Brändle de Motta JC, Berlemont A. A hybrid moment of fluid-level set framework for simulating primary atomization. J. Comput. Phys., 2022, 451: 110864

[6.] Bentz, M., Knoll, R., Hasan, M., Lin, C. Low-g fluid mixing: further results from the tank pressure control experiment. In: 29th Joint Propulsion Conference and Exhibit, AIAA-93-2423. AIAA (1993)

[7.] Bentz, M., Meserole, J., Knoll, R.: Jet mixing in low gravity-results of the tank pressure control experiment. In: 28th Joint Propulsion Conference and Exhibit, p. 3060 (1992)

[8.] Bonhomme R, Magnaudet J, Duval F, Piar B. Inertial dynamics of air bubbles crossing a horizontal fluid-fluid interface. J. Fluid Mech., 2012, 707: 405-443

[9.] Brackbill JU, Kothe DB, Zemach C. A continuum method for modeling surface tension. J. Comput. Phys., 1992, 100(2): 335-354

[10.] Breiman, L., Friedman, J., Olshen, R., Stone, C.: Classification and Regression Trees. Wadsworth & Brooks/Cole Advanced Books & Software, Monterey, CA (1984)

[11.] Caboussat A, Francois MM, Glowinski R, Kothe DB, Sicilian JM. A numerical method for interface reconstruction of triple points within a volume tracking algorithm. Math. Comput. Model., 2008, 48(11): 1957-1971

[12.] Colella P, Woodward PR. The piecewise parabolic method (ppm) for gas-dynamical simulations. J. Comput. Phys., 1984, 54(1): 174-201

[13.] Cummins SJ, Francois MM, Kothe DB. Estimating curvature from volume fractions. Comput Struct, 2005, 83(6–7): 425-434

[14.] De Gennes P-G, Brochard-Wyart F, Quéré D. . Capillarity and wetting phenomena: drops, bubbles, pearls, waves, 2013 New York Springer Science & Business Media

[15.] Dyadechko, V., Shashkov, M.: Moment-of-fluid interface reconstruction. Los Alamos report LA-UR-05-7571 (2005)

[16.] Dyadechko V, Shashkov M. Reconstruction of multi-material interfaces from moment data. J. Comput. Phys., 2008, 227(11): 5361-5384

[17.] Enright D, Fedkiw R, Ferziger J, Mitchell I. A hybrid particle level set method for improved interface capturing. J. Comput. Phys., 2002, 183(1): 83-116

[18.] Francois MM, Cummins SJ, Dendy ED, Kothe DB, Sicilian JM, Williams MW. A balanced-force algorithm for continuous and sharp interfacial surface tension models within a volume tracking framework. J. Comput. Phys., 2006, 213(1): 141-173

[19.] Gibou F, Fedkiw RP, Cheng L-T, Kang M. A second-order-accurate symmetric discretization of the Poisson equation on irregular domains. J. Comput. Phys., 2002, 176(1): 205-227

[20.] Glimm J, Isaacson E, Marchesin D, McBryan O. Front tracking for hyperbolic systems. Adv. Appl. Math., 1981, 2(1): 91-119

[21.] Godunov, S.: Different methods for shock waves. PhD Dissertation. Moscow State University (1954)

[22.] Godunov SK, Bohachevsky I. Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics. Matematičeskij sbornik, 1959, 47(3): 271-306

[23.] Harten A. High resolution schemes for hyperbolic conservation laws. J. Comput. Phys., 1997, 135(2): 260-278

[24.] Helsby F, Tuson K. Behaviour of air bubbles in aqueous solutions. Research, 1955, 8: 270-275

[25.] Hirt CW, Nichols BD. Volume of fluid (vof) method for the dynamics of free boundaries. J. Comput. Phys., 1981, 39(1): 201-225

[26.] Hu H, Jin Z. An icing physics study by using lifetime-based molecular tagging thermometry technique. Int. J. Multiph. Flow, 2010, 36(8): 672-681

[27.] Huang Z, Lin G, Ardekani AM. Consistent, essentially conservative and balanced-force phase-field method to model incompressible two-phase flows. J. Comput. Phys., 2020, 406: 109192

[28.] Huang Z, Lin G, Ardekani AM. A consistent and conservative phase-field model for thermo-gas-liquid-solid flows including liquid-solid phase change. J. Comput. Phys., 2022, 449: 110795

[29.] Jemison M, Loch E, Sussman M, Shashkov M, Arienti M, Ohta M, Wang Y. A coupled level set-moment of fluid method for incompressible two-phase flows. J. Sci. Comput., 2013, 54(2–3): 454-491

[30.] Jemison M, Sussman M, Arienti M. Compressible, multiphase semi-implicit method with moment of fluid interface representation. J. Comput. Phys., 2014, 279: 182-217

[31.] Kim J. Phase field computations for ternary fluid flows. Comput. Methods Appl. Mech. Eng., 2007, 196(45): 4779-4788

[32.] Kucharik M, Garimella RV, Schofield SP, Shashkov MJ. A comparative study of interface reconstruction methods for multi-material ale simulations. J. Comput. Phys., 2010, 229(7): 2432-2452

[33.] Li G, Lian Y, Guo Y, Jemison M, Sussman M, Helms T, Arienti M. Incompressible multiphase flow and encapsulation simulations using the moment-of-fluid method. Int. J. Numer. Meth. Fluids, 2015, 79(9): 456-490

[34.] Liu Y, Sussman M, Lian Y, Hussaini MY, Vahab M, Shoele K. A novel supermesh method for computing solutions to the multi-material stefan problem with complex deforming interfaces and microstructure. J. Sci. Comput., 2022, 91(1): 1-40

[35.] Lyu S, Wang K, Zhang Z, Pedrono A, Sun C, Legendre D. A hybrid vof-ibm method for the simulation of freezing liquid films and freezing drops. J. Comput. Phys., 2021, 432: 110160

[36.] Markstein G. Interaction of flow pulsations and flame propagation. J Aeronaut Sci, 1951, 18(6): 428-429

[37.] Miao F, Wu B, Sun Z, Peng C. Calibration method of the laser beam based on liquid lens for 3d precise measurement. Measurement, 2021, 178: 109358

[38.] Milcent T, Lemoine A. Moment-of-fluid analytic reconstruction on 3d rectangular hexahedrons. J. Comput. Phys., 2020, 409: 109346

[39.] Ohta M, Kikuchi D, Yoshida Y, Sussman M. Robust numerical analysis of the dynamic bubble formation process in a viscous liquid. Int. J. Multiph. Flow, 2011, 37(9): 1059-1071

[40.] Olsson E, Kreiss G. A conservative level set method for two phase flow. J. Comput. Phys., 2005, 210(1): 225-246

[41.] Osher S, Sethian JA. Fronts propagating with curvature-dependent speed: algorithms based on hamilton-jacobi formulations. J. Comput. Phys., 1988, 79(1): 12-49

[42.] Pathak A, Raessi M. A three-dimensional volume-of-fluid method for reconstructing and advecting three-material interfaces forming contact lines. J. Comput. Phys., 2016, 307: 550-573

[43.] Pei C, Vahab M, Sussman M, Hussaini MY. A hierarchical space-time spectral element and moment-of-fluid method for improved capturing of vortical structures in incompressible multi-phase/multi-material flows. J. Sci. Comput., 2019, 81(3): 1527-1566

[44.] Qiu R, Huang R, Xiao Y, Wang J, Zhang Z, Yue J, Zeng Z, Wang Y. Physics-informed neural networks for phase-field method in two-phase flow. Phys. Fluids, 2022, 34(5): 052109

[45.] Remmerswaal RA, Veldman AE. Parabolic interface reconstruction for 2d volume of fluid methods. J. Comput. Phys., 2022, 469: 111473

[46.] Salas MD. Shock fitting method for complicated two-dimensional supersonic flows. AIAA J., 1976, 14(5): 583-588

[47.] Saurel R, Abgrall R. A multiphase godunov method for compressible multifluid and multiphase flows. J. Comput. Phys., 1999, 150(2): 425-467

[48.] Schofield SP, Garimella RV, Francois MM, Loubère R. Material order-independent interface reconstruction using power diagrams. Int. J. Numer. Meth. Fluids, 2008, 56(6): 643

[49.] Schofield SP, Garimella RV, Francois MM, Loubère R. A second-order accurate material-order-independent interface reconstruction technique for multi-material flow simulations. J. Comput. Phys., 2009, 228(3): 731-745

[50.] Shetabivash H, Dolatabadi A, Paraschivoiu M. A multiple level-set approach for modelling containerless freezing process. J. Comput. Phys., 2020, 415: 109527

[51.] Shin S, Juric D. A hybrid interface method for three-dimensional multiphase flows based on front tracking and level set techniques. Int. J. Numer. Meth. Fluids, 2009, 60(7): 753-778

[52.] Shu C-W, Osher S. Efficient implementation of essentially non-oscillatory shock-capturing schemes. J. Comput. Phys., 1988, 77(2): 439-471

[53.] Sijoy C, Chaturvedi S. Volume-of-fluid algorithm with different modified dynamic material ordering methods and their comparisons. J. Comput. Phys., 2010, 229(10): 3848-3863

[54.] Smith KA, Solis FJ, Chopp D. A projection method for motion of triple junctions by level sets. Interfac Free Bound, 2002, 4(3): 263-276

[55.] Starinshak DP, Karni S, Roe PL. A new level set model for multimaterial flows. J. Comput. Phys., 2014, 262: 1-16

[56.] Sussman M. A second order coupled level set and volume-of-fluid method for computing growth and collapse of vapor bubbles. J. Comput. Phys., 2003, 187(1): 110-136

[57.] Sussman M, Ohta M. A stable and efficient method for treating surface tension in incompressible two-phase flow. SIAM J. Sci. Comput., 2009, 31(4): 2447-2471

[58.] Sussman M, Puckett EG. A coupled level set and volume-of-fluid method for computing 3d and axisymmetric incompressible two-phase flows. J. Comput. Phys., 2000, 162(2): 301-337

[59.] Sussman M, Smereka P, Osher S. A level set approach for computing solutions to incompressible two-phase flow. J. Comput. Phys., 1994, 114(1): 146-159

[60.] Tatebe, O.: The multigrid preconditioned conjugate gradient method. In: the Sixth Copper Mountain Conference on Multigrid Methods, Part 2. NASA, Copper Mountain (1993)

[61.] Unverdi SO, Tryggvason G. A front-tracking method for viscous, incompressible, multi-fluid flows. J. Comput. Phys., 1992, 100(1): 25-37

[62.] Vahab, M., Pei, C., Hussaini, M.Y., Sussman, M., Lian, Y.: An adaptive coupled level set and moment-of-fluid method for simulating droplet impact and solidification on solid surfaces with application to aircraft icing. In: 54th AIAA Aerospace Sciences Meeting, p. 1340 (2016)

[63.] Vahab M, Sussman M, Shoele K. Fluid-structure interaction of thin flexible bodies in multi-material multi-phase systems. J. Comput. Phys., 2021, 429: 110008

[64.] Van Leer B. Towards the ultimate conservative difference scheme. v. a second-order sequel to godunov’s method. J. Comput. Phys., 1979, 32: 101-136

[65.] Vu TV, Tryggvason G, Homma S, Wells JC. Numerical investigations of drop solidification on a cold plate in the presence of volume change. Int. J. Multiph. Flow, 2015, 76: 73-85

[66.] Welch SW, Wilson J. A volume of fluid based method for fluid flows with phase change. J. Comput. Phys., 2000, 160(2): 662-682

[67.] Zalesak ST. Fully multidimensional flux-corrected transport algorithms for fluids. J. Comput. Phys., 1979, 31(3): 335-362

Funding

National Aeronautics and Space Administration(80NSSC20K0352)