A new volume conservation enforcement method for PLIC reconstruction in general convex grids. (English) Zbl 1349.76627
Summary: A comprehensive study is made of methods for resolving the volume conservation enforcement problem in the PLIC reconstruction of an interface in general 3D convex grids. Different procedures to bracket the solution when solving the problem using previous standard methods are analyzed in detail. A new interpolation bracketing procedure and an improved analytical method to find the interface plane constant are proposed. These techniques are combined in a new method to enforce volume conservation, which does not require the sequential polyhedra truncation operations typically used in standard methods. The new methods have been implemented into existing geometrical routines described in López and Hernández , which are further improved by using more efficient formulae to compute areas and volumes of general convex 2 and 3D polytopes. Different tests using regular and irregular cell geometries are carried out to demonstrate the robustness and substantial improvement in computational efficiency of the proposed techniques, which increase the computation speed of the mentioned routines by up to 3 times for the 3D problems considered in this work.
MSC:
| 76M25 | Other numerical methods (fluid mechanics) (MSC2010) |
| 76T30 | Three or more component flows |
Keywords:
volume of fluid method; PLIC reconstruction; volume conservation enforcement; general convex grids; unstructured grids; solution bracketing; volume truncationReferences:
| [1] | Ahn, H. T.; Shashkov, M., Multi-material interface reconstruction on generalized polyhedral meshes (2007), Los Alamos National Laboratory, Tech. Rep. LA-UR-07-0656 · Zbl 1388.76232 |
| [2] | Anbarlooei, H. R.; Mazaheri, K., ‘Moment of fluid’ interface reconstruction method in axisymmetric coordinates, Int. J. Numer. Methods Biomed. Eng., 27, 1640-1651 (2011) · Zbl 1469.76082 |
| [3] | Barth, T. J.; Frederickson, P. O., Higher-order solution of the Euler equations on unstructured grids using quadratic reconstruction, (28th AIAA Aerosp. Sci. Meeting (1990)) |
| [4] | Brent, R. P., Algorithms for Minimization Without Derivatives (1973), Prentice-Hall: Prentice-Hall Englewood Cliffs, New Jersey, 195 pp · Zbl 0245.65032 |
| [5] | Cummins, S. J.; François, M. M.; Kothe, D. B., Estimating curvature from volume fractions, Compos. Struct., 83, 425-434 (2005) |
| [6] | Diot, S.; François, M. M.; Dendy, E. D., An interface reconstruction method based on analytical formulae for 2D planar and axisymmetric arbitrary convex cells, J. Comput. Phys., 275, 53-64 (2014) · Zbl 1349.76600 |
| [7] | Diot, S.; François, M. M., An interface reconstruction method based on an analytical formula for 3D arbitrary convex cells, J. Comput. Phys., 305, 63-74 (2016) · Zbl 1349.76601 |
| [8] | Doob, J. L., Stochastic Processes (1967), Wiley: Wiley New York |
| [9] | Dyadechko, V.; Shashkov, M., Moment-of-fluid interface reconstruction, 1-33 (2005), LANL Report, LA-UR-05-7571 |
| [10] | Gelder, A. V., Efficient computation of polygon area and polyhedron volume, (Paeth, A. W., Graphics Gems V (1995), Academic Press: Academic Press New York), 35-41 |
| [11] | Goldman, R. N., Area of planar polygons and volume of polyhedra, (Arvo, J., Graphics Gems II (1991), Academic Press: Academic Press Boston), 170-171 |
| [12] | Gueyffier, D.; Li, J.; Nadim, A.; Scardovelli, R.; Zaleski, S., Volume-of-fluid interface tracking with smoothed surface stress methods for three-dimensional flows, J. Comput. Phys., 152, 423-456 (1999) · Zbl 0954.76063 |
| [13] | Harvie, D. J.E.; Fletcher, D. F., A new volume of fluid advection algorithm: the stream scheme, J. Comput. Phys., 162, 1-32 (2000) · Zbl 0964.76068 |
| [14] | Li, J., Calcul d’interface affine par morceaux, C. R. Acad. Sci. Paris, Sér. IIb (Paris), 320, 391-396 (1995) · Zbl 0826.76065 |
| [15] | López, J.; Hernández, J., Analytical and geometrical tools for 3D volume of fluid methods in general grids, J. Comput. Phys., 227, 5939-5948 (2008) · Zbl 1142.76042 |
| [16] | López, J.; Zanzi, C.; Gómez, P.; Zamora, R.; Faura, F.; Hernández, J., An improved height function technique for computing interface curvature from volume fractions, Comput. Methods Appl. Mech. Eng., 198, 2555-2564 (2009) · Zbl 1228.76099 |
| [17] | Peterson, W. W., Addressing for random storage, IBM J. Res. Dev., 1, 131-132 (1957) |
| [18] | Rider, W. J.; Kothe, D. B., Reconstructing volume tracking, J. Comput. Phys., 141, 112-152 (1998) · Zbl 0933.76069 |
| [19] | Rycroft, C. H., Voro++: a three-dimensional Voronoi cell library in C++ (2009) |
| [20] | Rycroft, C. H., Multiscale modeling in granular flow (2007), Massachusetts Institute of Technology, Ph.D. thesis |
| [21] | Scardovelli, R.; Zaleski, S., Direct numerical simulation of free-surface and interfacial flow, Annu. Rev. Fluid Mech., 31, 567-603 (1999) |
| [22] | Scardovelli, R.; Zaleski, S., Analytical relations connecting linear interfaces and volume fractions in rectangular grids, J. Comput. Phys., 164, 228-237 (2000) · Zbl 0993.76067 |
| [23] | Schneider, P. J.; Eberly, D. H., Geometric Tools for Computer Graphics (2003), Morgan Kaufmann Publishers: Morgan Kaufmann Publishers London |
| [24] | Sunday, D., Fast polygon area and Newell normal computation, J. Graph. Tools, 7, 9-13 (2002) |
| [25] | Yang, X.; James, A. J., Analytical relations for reconstructing piecewise linear interfaces in triangular and tetrahedral grids, J. Comput. Phys., 214, 41-54 (2006) · Zbl 1137.76456 |
| [26] | The code and pseudo-code are available at the web site |
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.
