Gómez, Pablo; Zanzi, Claudio; López, Joaquín; Hernández, Julio Simulation of high density ratio interfacial flows on cell vertex/edge-based staggered octree grids with second-order discretization at irregular nodes. (English) Zbl 1416.76037 J. Comput. Phys. 376, 478-507 (2019). Summary: A numerical code has been developed for the simulation of unsteady incompressible interfacial flows with large density ratios, based on discretizing the conservation equations on a rectangular adaptive grid with a graded octree data structure, in which the pressure and velocity components are stored at the cell vertices and edges, respectively. With this arrangement, which is novel for octree grids, node alignment simplifies the Poisson equation discretization at nodes common to cells with different refinement levels (irregular nodes), while the staggered storage of variables avoids the pressure-velocity coupling difficulties associated with collocated grids. Three different discretization approaches at irregular nodes are proposed: second- and first-order schemes, and an efficient scheme based on a linear interpolation from the surrounding nodes. A grid refinement test in two dimensions, and 3D deformation and static bubble tests were carried out to assess the accuracy and efficiency of the proposed discretization methods at irregular nodes, the performance of the different schemes used to solve the level set transport equation and the capability of the numerical code to reduce spurious currents. The results of the tests are discussed and compared with results available in the literature. Finally, the ability of the code to accurately simulate the complex phenomena involved in the impact of a water drop on a free surface is demonstrated by thoroughly comparing numerical and experimental results. Cited in 4 Documents MSC: 76D05 Navier-Stokes equations for incompressible viscous fluids 76M20 Finite difference methods applied to problems in fluid mechanics 76D45 Capillarity (surface tension) for incompressible viscous fluids 76T99 Multiphase and multicomponent flows Keywords:high density ratio interfacial flows; graded octree grids; cell vertex/edge staggered variable arrangement; projection method; level set method; drop impact on a pool Software:PROST; VOFTools; SPARSKIT; Gerris PDFBibTeX XMLCite \textit{P. Gómez} et al., J. Comput. Phys. 376, 478--507 (2019; Zbl 1416.76037) Full Text: DOI References: [1] Berger, M. J.; Oliger, J., Adaptive mesh refinement for hyperbolic partial differential equations, J. Comput. Phys., 53, 484-512 (1984) · Zbl 0536.65071 [2] Berger, M. J.; Colella, P., Local adaptive mesh refinement for shock hydrodynamics, J. Comput. Phys., 82, 64-84 (1989) · Zbl 0665.76070 [3] Almgren, A.; Bell, J.; Colella, P.; Howell, L.; Welcome, M., A conservative adaptive projection method for the variable density incompressible Navier-Stokes equations, J. Comput. Phys., 142, 1-46 (1998) · Zbl 0933.76055 [4] Sussman, M.; Almgren, A.; Bell, J.; Colella, P.; Howell, L.; Welcome, M., An adaptive level set approach for incompressible two-phase flows, J. Comput. Phys., 148, 81-124 (1999) · Zbl 0930.76068 [5] English, R. E.; Qiu, L.; Yu, Y.; Fedkiw, R., An adaptive discretization of incompressible flow using a multitude of moving Cartesian grids, J. Comput. Phys., 254, 107-154 (2013) · Zbl 1349.76456 [6] Samet, H., The Design and Analysis of Spatial Data Structures (1989), Addison-Wesley: Addison-Wesley Reading, MA [7] Strain, J., Fast tree-based redistancing for level set computations, J. Comput. Phys., 152, 664-686 (1999) · Zbl 0944.65020 [8] Popinet, S., Gerris: a tree-based adaptive solver for the incompressible Euler equations in complex geometries, J. Comput. Phys., 190, 572-600 (2003) · Zbl 1076.76002 [9] Min, C.; Gibou, F., A second order accurate projection method for the incompressible Navier-Stokes equation on non-graded adaptive grids, J. Comput. Phys., 219, 912-929 (2006) · Zbl 1330.76096 [10] Losasso, F.; Fedkiw, R.; Osher, S., Spatially adaptive techniques for level set methods and incompressible flow, Comput. Fluids, 35, 995-1010 (2006) · Zbl 1177.76295 [11] Gómez, P.; Hernández, J.; López, J., On the reinitialization procedure in a narrow-band locally-refined level set method for interfacial flows, Int. J. Numer. Methods Eng., 63, 10, 1478-1512 (2005) · Zbl 1086.76559 [12] Min, C.; Gibou, F., A second order accurate level set method on non-graded adaptive Cartesian grids, J. Comput. Phys., 225, 300-321 (2007) · Zbl 1122.65077 [13] Min, C.; Gibou, F.; Ceniceros, H., A supra-convergent finite difference scheme for the variable coefficient Poisson equation on non-graded grids, J. Comput. Phys., 218, 123-140 (2006) · Zbl 1106.65093 [14] Losasso, F.; Gibou, F.; Fedkiw, R., Simulating water and smoke with an octree data structure, ACM Trans. Graph. (SIGGRAPH Proc.), 457-462 (2004) [15] Manteuffel, T.; White, A., The numerical solution of second-order boundary value problems on nonuniform meshes, Math. Comput., 47, 511-535 (1986) · Zbl 0635.65092 [16] Batty, C., A cell-centred finite volume method for the Poisson problem on non-graded quadtrees with second order accurate gradients, J. Comput. Phys., 331, 49-72 (2017) · Zbl 1380.65399 [17] Guittet, A.; Theillard, M.; Gibou, F., A stable projection method for the incompressible Navier-Stokes equations on arbitrary geometries and adaptive quad/octrees, J. Comput. Phys., 292, 215-238 (2015) · Zbl 1349.76336 [18] Ni, M.-J., Consistent projection methods for variable density incompressible Navier-Stokes equations with continuous surface forces on a rectangular collocated mesh, J. Comput. Phys., 228, 6938-6956 (2009) · Zbl 1261.76034 [19] López, J.; Hernández, J.; Gómez, P.; Faura, F., A new volume conservation enforcement method for PLIC reconstruction in general convex grids, J. Comput. Phys., 316, 338-359 (2016) · Zbl 1349.76627 [20] López, J.; Hernández, J.; Gómez, P.; Faura, F., VOFTools - a software package of calculation tools for volume of fluid methods using general convex grids, Comput. Phys. Commun., 223C, 45-54 (2018) · Zbl 07693323 [21] Chang, Y.; Hou, T.; Merriman, B.; Osher, S., A level set formulation of Eulerian interface capturing methods for incompressible fluid flows, J. Comput. Phys., 124, 449-464 (1996) · Zbl 0847.76048 [22] 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 [23] López, J.; Hernández, J.; Gómez, P.; Faura, F., An improved PLIC-VOF method for tracking thin fluid structures in incompressible two-phase flows, J. Comput. Phys., 208, 1, 51-74 (2005) · Zbl 1115.76359 [24] Hernández, J.; López, J.; Gómez, P.; Zanzi, C.; Faura, F., A new volume of fluid method in three dimensions, part I: multidimensional advection method with face-matched flux polyhedra, Int. J. Numer. Methods Fluids, 58, 897-921 (2008) · Zbl 1151.76552 [25] Sussman, M.; Fatemi, E.; Smereka, P.; Osher, S., An improved level set method for incompressible two-phase flows, Comput. Fluids, 27, 663-680 (1998) · Zbl 0967.76078 [26] Gómez, P.; Zanzi, C.; Palacios, J.; Hernández, J.; López, J., An adaptive-grid projection method for high density ratio interfacial flows, (ASME 2013 International Mechanical Engineering Congress and Exposition. ASME 2013 International Mechanical Engineering Congress and Exposition, San Diego, California (2013)) [27] Ni, M.-J.; Abdou, M.; Komori, S., A variable-density projection method for interfacial flows, Numer. Heat Transf., Part B, Fundam., 44, 553-574 (2003) [28] Sussman, M.; Smith, K. M.; Hussaini, M. Y.; Ohta, M.; Zhi-Wei, R., A sharp interface method for incompressible two-phase flows, J. Comput. Phys., 221, 469-505 (2007) · Zbl 1194.76219 [29] Almgren, A.; Bell, J.; Szymczak, W., A numerical method for the incompressible Navier-Stokes equations based on an approximate projection, SIAM J. Sci. Comput., 17, 358-369 (1996) · Zbl 0845.76055 [30] Brown, D.; Cortez, R.; Minion, M., Accurate projection methods for the incompressible Navier-Stokes equations, J. Comput. Phys., 168, 464-499 (2001) · Zbl 1153.76339 [31] Francois, M. M.; Cummins, S. J.; Dendy, E. D.; Kothe, D. B.; Sicilian, J. M.; Williams, M. W., A balanced-force algorithm for continuous and sharp interfacial surface tension models within a volume tracking framework, J. Comput. Phys., 213, 141-173 (2006) · Zbl 1137.76465 [32] Renardy, Y.; Renardy, M., PROST: a parabolic reconstruction of surface tension for the volume-of-fluid method, J. Comput. Phys., 183, 400-421 (2002) · Zbl 1057.76569 [33] Lee, B.; Min, C., An energy-stable method for solving the incompressible Navier-Stokes equations with non-slip boundary condition, J. Comput. Phys., 360, 104-119 (2018) · Zbl 1391.76483 [34] Puckett, E.; Almgren, A.; Bell, J.; Marcus, D.; Rider, W., A high order projection method for tracking fluid interfaces in variable density incompressible flows, J. Comput. Phys., 130, 269-282 (1997) · Zbl 0872.76065 [35] Colella, P., Multidimensional upwind methods for hyperbolic conservation laws, J. Comput. Phys., 87, 171-200 (1990) · Zbl 0694.65041 [36] Colella, P., A direct Eulerian MUSCL scheme for gas dynamics, SIAM J. Sci. Stat. Comput., 6, 104-117 (1985) · Zbl 0562.76072 [37] Popinet, S., Numerical models of surface tension, Annu. Rev. Fluid Mech., 50, 1, 49-75 (2018) · Zbl 1384.76016 [38] Hysing, S., A new implicit surface tension implementation for interfacial flows, Int. J. Numer. Methods Fluids, 51, 6, 659-672 (2006) · Zbl 1158.76350 [39] Sussman, M.; Ohta, M., A stable and efficient method for treating surface tension in incompressible two-phase flow, SIAM J. Sci. Comput., 31, 4, 2447-2471 (2009) · Zbl 1387.76107 [40] Denner, F.; van Wachem, B. G., Numerical time-step restrictions as a result of capillary waves, J. Comput. Phys., 285, 24-40 (2015) · Zbl 1351.76028 [41] Liu, X.-D.; Osher, S.; Chan, T., Weighted essentially non-oscillatory schemes, J. Comput. Phys., 115, 1, 200-212 (1994) · Zbl 0811.65076 [42] Jiang, G. S.; Peng, D., Weighted ENO schemes for Hamilton-Jacobi equations, SIAM J. Sci. Comput., 21, 6, 2126-2143 (2000) · Zbl 0957.35014 [43] Sussman, M.; Fatemi, E., An efficient, interface-preserving level set redistancing algorithm and its application to interfacial incompressible fluid flow, SIAM J. Sci. Comput., 20, 4, 1165-1191 (1999) · Zbl 0958.76070 [44] Shu, C.-W., Total-variation-diminishing time discretizations, SIAM J. Sci. Stat. Comput., 9, 6, 1073-1084 (1988) · Zbl 0662.65081 [45] Saad, Y., SPARSKIT: A Basic Tool Kit for Sparse Matrix Computations - Version 2 (1994) [46] Minion, M. L., A projection method for locally refined grids, J. Comput. Phys., 127, 1, 158-178 (1996) · Zbl 0859.76047 [47] Enright, D.; Fedkiw, R.; Ferziger, J.; Mitchell, I., A hybrid particle level set method for improved interface capturing, J. Comput. Phys., 183, 83-116 (2002) · Zbl 1021.76044 [48] Solomenko, Z.; Spelt, P. D.M.; Náraigh, L. O.; Alix, P., Mass conservation and reduction of parasitic interfacial waves in level-set methods for the numerical simulation of two-phase flows: a comparative study, Int. J. Multiph. Flow, 95, 235-256 (2017) [49] Popinet, S., Accurate adaptive solver for surface tension driven flow, J. Comput. Phys., 228, 5838-5866 (2009) · Zbl 1280.76020 [50] Abadie, T.; Austin, J.; Legendre, D., On the combined effects of surface tension force calculation and interface advection on spurious currents within volume of fluid and level set frameworks, J. Comput. Phys., 297, 611-636 (2015) · Zbl 1349.76422 [51] Meland, R.; Gran, I.; Olsen, R.; Munkejord, S., Reduction of parasitic currents in level-set calculations with a consistent discretization of the surface-tension force for the CSF model, (16th Australasian Fluid Mechanics Conference. 16th Australasian Fluid Mechanics Conference, Gold Coast, Queensland, Australia (2007)), 862-865 [52] Yokoi, K., A density-scaled continuum surface force model within a balanced force formulation, J. Comput. Phys., 278, 221-228 (2014) · Zbl 1349.76410 [53] Palacios, J.; Hernández, J.; Gómez, P.; Zanzi, C.; López, J., Experimental study of splashing patterns and the splashing/deposition threshold in drop impacts onto dry smooth solid surfaces, Exp. Therm. Fluid Sci., 44, 571-582 (2012) [54] Yu, J.-D.; Sakai, S.; Sethian, J., A coupled level set projection method applied to ink jet simulation, Interfaces Free Bound., 5, 459-482 (2003) · Zbl 1042.35059 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.