Calculation of the interface curvature and normal vector with the level-set method. (English) Zbl 1290.76101
Summary: This article addresses the use of the level-set method for capturing the interface between two fluids. One of the advantages of the level-set method is that the curvature and the normal vector of the interface can be readily calculated from the level-set function. However, in cases where the level-set method is used to capture topological changes, the standard discretization techniques for the curvature and the normal vector do not work properly. This is because they are affected by the discontinuities of the signed-distance function half-way between two interfaces. This article addresses the calculation of normal vectors and curvatures with the level-set method for such cases. It presents a discretization scheme based on the geometry-aware curvature discretization by P. Macklin and J. Lowengrub [J. Comput. Phys. 215, No. 2, 392–401 (2006; Zbl 1089.92024)]. As the present scheme is independent of the ghost-fluid method, it becomes more generally applicable, and it can be implemented into an existing level-set code more easily than Macklin and Lowengrub’s scheme [Zbl 1089.92024]. The present scheme is compared with the second-order central-difference scheme and with Macklin and Lowengrub’s scheme [Zbl 1089.92024], first for a case with no flow, then for a case where two drops collide in a 2D shear flow, and finally for a case where two drops collide in an axisymmetric flow. In the latter two cases, the Navier-Stokes equations for incompressible two-phase flow are solved. The article also gives a comparison of the calculation of normal vectors with the direction difference scheme presented by P. Macklin and J. Lowengrub [J. Comput. Phys. 203, No. 1, 191–220 (2005; Zbl 1067.65111)] and with the present discretization scheme. The results show that the present discretization scheme yields more robust calculations of the curvature than the second-order central difference scheme in areas where topological changes are imminent. The present scheme compares well to Macklin and Lowengrub’s method [Zbl 1089.92024]. The results also demonstrate that the direction difference scheme [Zbl 1067.65111] is not always sufficient to accurately calculate the normal vectors.
MSC:
| 76M20 | Finite difference methods applied to problems in fluid mechanics |
| 76T99 | Multiphase and multicomponent flows |
Keywords:
level-set method; curvature discretization; normal-vector discretization; curve-fitting discretization scheme; finite differences; ghost-fluid methodSoftware:
pchipReferences:
| [1] | Macklin, P.; Lowengrub, J., An improved geometry-aware curvature discretization for level set methods: application to tumor growth, J Comput Phys, 215, 392-401 (2006) · Zbl 1089.92024 |
| [2] | Macklin, P.; Lowengrub, J., Evolving interfaces via gradients of geometry-dependent interior Poisson problems: application to tumor growth, J Comput Phys, 203, 191-220 (2005) · Zbl 1067.65111 |
| [3] | Osher, S.; Sethian, J. A., Fronts propagating with curvature dependent speed: algorithms based on Hamilton-Jacobi formulations, J Comput Phys, 79, 12-49 (1988) · Zbl 0659.65132 |
| [4] | Sethian, J. A.; Smereka, P., Level set methods for fluid interfaces, Ann Rev Fluid Mech, 35, 341-372 (2003) · Zbl 1041.76057 |
| [5] | Macklin, P.; Lowengrub, J. S., A new ghost cell/level set method for moving boundary problems: application to tumor growth, J Sci Comput, 35, 266-299 (2008) · Zbl 1203.65144 |
| [6] | Salac, D.; Lu, W., A local semi-implicit level-set method for interface motion, J Sci Comput, 35, 330-349 (2008) · Zbl 1203.65147 |
| [7] | Wang, Z.; Tong, A. Y., A sharp surface tension modeling method for two-phase incompressible interfacial flows, Int J Numer Methods Fluids, 64, 709-732 (2010) · Zbl 1426.76438 |
| [8] | Herrmann, M., A balanced force refined level set grid method for two-phase flows on unstructured flow solver grids, J Comput Phys, 227, 2674-2706 (2008) · Zbl 1388.76252 |
| [9] | Marchandise, E.; Geuzaine, P.; Chevaugeon, N.; Remacle, J.-F., A stabilized finite element method using a discontinuous level set approach for the computation of bubble dynamics, J Comput Phys, 225, 949-974 (2007) · Zbl 1118.76040 |
| [10] | Desjardins, O.; Moureau, V.; Pitsch, H., An accurate conservative level set/ghost fluid method for simulating turbulent atomization, J Comput Phys, 227, 8395-8416 (2008) · Zbl 1256.76051 |
| [12] | Xu, J.-J.; Li, Z.; Lowengrub, J.; Zhao, H.-K., A level set method for interfacial flows with surfactants, J Comput Phys, 212, 2, 590-616 (2006) · Zbl 1161.76548 |
| [14] | Teigen, K. E.; Munkejord, S. T., Influence of surfactant on drop deformation in an electric field, Phys Fluids, 22, 11, 112104 (2010) |
| [15] | Brackbill, J. U.; Kothe, D. B.; Zemach, C., A continuum method for modeling surface tension, J Comput Phys, 100, 335-354 (1992) · Zbl 0775.76110 |
| [16] | Kang, M.; Fedkiw, R. P.; Liu, X.-D., A boundary condition capturing method for multiphase incompressible flow, J Sci Comput, 15, 3, 323-360 (2000) · Zbl 1049.76046 |
| [18] | Zhao, H.-K.; Chan, T.; Merriman, B.; Osher, S., A variational level set approach to multiphase motion, J Comput Phys, 127, 179-195 (1996) · Zbl 0860.65050 |
| [19] | 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 |
| [20] | Gottlieb, S.; Shu, C. W.; Tadmor, E., Strong stability-preserving high-order time discretization methods, SIAM Rev, 43, 89-112 (2001) · Zbl 0967.65098 |
| [21] | Adalsteinsson, D.; Sethian, J. A., A fast level set method for propagating interfaces, J Comput Phys, 118, 269-277 (1995) · Zbl 0823.65137 |
| [22] | Liu, X.-D.; Fedkiw, R. P.; Kang, M., A boundary condition capturing method for Poisson’s equation on irregular domains, J Comput Phys, 160, 151-178 (2000) · Zbl 0958.65105 |
| [23] | Osher, S.; Fedkiw, R. P., The level-set method and dynamic implicit surfaces (2003), Springer · Zbl 1026.76001 |
| [24] | Lorensen, W. E.; Cline, H. E., Marching cubes: a high resolution 3d surface construction algorithm, Comput Graphics, 21, 4, 163-169 (1987) |
| [25] | Prautzsch, H.; Boehm, W.; Paluszny, M., B(́e)zier and B-spline techniques (2002), Springer · Zbl 1033.65008 |
| [26] | Fritsch, F. N.; Carlson, R. E., Monotone piecewise cubic interpolation, SIAM J Numer Anal, 17, 2, 238-246 (1980) · Zbl 0423.65011 |
| [27] | Waerden, B.; Artin, E.; Noether, E., Algebra, vol. 1 (2003), Springer-Verlag, <http://books.google.com/books?id=XDN8yR8R1OUC> · Zbl 1032.00002 |
| [28] | Lang, S., Algebra. Algebra, Graduate texts in mathematics (2002), Springer, <http://books.google.com/books?id=Fge-BwqhqIYC> · Zbl 0984.00001 |
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.
