×

Mixed element FEM level set method for numerical simulation of immiscible fluids. (English) Zbl 1426.76270

Summary: A new realization of a finite element level set method for simulation of immiscible fluid flows is introduced and validated on numerical benchmarks. The new method involves a mixed discretization of the dependent variables, discretizing the flow variables with non-conforming Rannacher – Turek finite elements while using a simple first order conforming discretization of the level set field. A three step segregated solution approach is employed, first a discrete projection method is used to decouple and compute the velocity and pressure separately, after which the level set field can be computed independently.
The developed method is tested and validated on a static bubble test case and on a numerical rising bubble test case for which a very accurate benchmark solution has been established. The new approach is also compared against two commercial simulation codes, Ansys Fluent and Comsol Multiphysics, which shows that the developed method is a magnitude or more accurate and at the same time significantly faster than state of the art commercial codes.

MSC:

76M10 Finite element methods applied to problems in fluid mechanics
76M25 Other numerical methods (fluid mechanics) (MSC2010)
76T10 Liquid-gas two-phase flows, bubbly flows
76D05 Navier-Stokes equations for incompressible viscous fluids

Software:

FEATFLOW; COMSOL
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] 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
[2] Bristeau, M. O.; Glowinski, R.; Periaux, J., Numerical methods for the Navier-Stokes equations. Applications to the simulation of compressible and incompressible viscous flows, Comput. Phys. Rep., 6, 73-187 (1987)
[3] Brooks, A. N.; Hughes, T. J.R., Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations, Comput. Methods Appl. Mech. Eng., 32, 199-259 (1982) · Zbl 0497.76041
[4] Hackbusch, W., Multi-Grid Methods and Applications. Multi-Grid Methods and Applications, Springer Series in Computational Mathematics (1985), Springer-Verlag, ISBN:0-387-12761-5
[5] Harlow, F. H.; Welch, J. E., Numerical calculation of time-dependent viscous incompressible flow of fluid with free surface, Phys. Fluids, 8, 2182-2189 (1965) · Zbl 1180.76043
[6] 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
[7] Hinton, E.; Rock, T.; Zienkiewicz, O. C., A note on mass lumping and related processes in the finite element method, Earthq. Eng. Struct. D., 4, 245-249 (1976)
[8] Hughes, T. J.R.; Brooks, A., A multidimensional upwind scheme with no crosswind diffusion in Finite Element Methods for Convection Dominated Flows, (Hughes, T. J.R., AMD, Vol. 34 (1979), ASME: ASME New York), 19-35
[9] Hysing, S.; Turek, S., The Eikonal equation: Numerical efficiency vs. algorithmic complexity on quadrilateral grids, (Proceedings of Algoritmy (2005), Slovak University of Bratislava: Slovak University of Bratislava Bratislava), 22-31, ISBN: 80-227-2192-1 · Zbl 1180.65108
[10] Hysing, S., A new implicit surface tension implementation for interfacial flows, Int. J. Num. Meth. Fluids, 51, 659-672 (2006) · Zbl 1158.76350
[11] S. Hysing, Numerical Simulation of Immiscible Fluids with FEM Level Set Techniques, PhD Thesis, University of Dortmund, Institute of Applied Mathematics, Dortmund, 2007.; S. Hysing, Numerical Simulation of Immiscible Fluids with FEM Level Set Techniques, PhD Thesis, University of Dortmund, Institute of Applied Mathematics, Dortmund, 2007. · Zbl 1196.76003
[12] Hysing, S.; Turek, S.; Kuzmin, D.; Parolini, N.; Burman, E.; Ganesan, S.; Tobiska, L., Quantitative benchmark computations of two-dimensional bubble dynamics, Int. J. Num. Meth. Fluids, 60, 1259-1288 (2009) · Zbl 1273.76276
[13] Rising bubble benchmark reference data: http://www.featflow.de/en/benchmarks/cfdbenchmarking/bubble.html; Rising bubble benchmark reference data: http://www.featflow.de/en/benchmarks/cfdbenchmarking/bubble.html
[14] Kuzmin, D., A high-resolution finite element scheme for convection-dominated transport, Communications in Numerical Methods in Engineering, 16, 215-223 (2000) · Zbl 0948.65099
[15] Kuzmin, D., On the design of general-purpose flux limiters for finite element schemes. I. Scalar convection, J. Comput. Phys., 219, 513-531 (2006) · Zbl 1189.76342
[16] Kuzmin, D.; Turek, S., High-resolution FEM-TVD schemes based on a fully multidimensional flux limiter, J. Comput. Phys., 198, 131-158 (2004) · Zbl 1107.76352
[17] Kuzmin, D.; Löhner, R.; Turek, S., Flux-Corrected Transport: Principles, Algorithms, and Applications. Flux-Corrected Transport: Principles, Algorithms, and Applications, Scientific Computation (2005), Springer, ISBN: 978-3-540-23730-3
[18] Osher, S. J.; Fedkiw, R. P., Level Set Methods and Dynamic Implicit Surfaces, (Applied Mathematical Sciences, vol. 153 (2002), Springer), ISBN: 978-0-387-95482-0 · Zbl 1026.76001
[19] Osher, S. J.; Sethian, S. A., Fronts propagating with curvature-dependent speed: Algorithms based on Hamilton-Jacobi formulations, J. Comput. Phys., 79, 12-49 (1988) · Zbl 0659.65132
[20] Peskin, C. S., Numerical analysis of blood flow in the heart, J. Comput. Phys., 25, 220-252 (1977) · Zbl 0403.76100
[21] W.J. Rider, D.B. Kothe, Stretching and tearing interface tracking methods, in the proceeding of the 12th AIAA Computational Fluid Dynamics Conference, AIAA 95-1717, 1995.; W.J. Rider, D.B. Kothe, Stretching and tearing interface tracking methods, in the proceeding of the 12th AIAA Computational Fluid Dynamics Conference, AIAA 95-1717, 1995.
[22] Sethian, J. A., A fast marching level set method for monotonically advancing fronts, Proc. Natl. Acad. Sci., 93, 1591-1595 (1996) · Zbl 0852.65055
[23] Turek, S.; Rannacher, R., A simple nonconforming quadrilateral Stokes element, Num. Meth. Part. D.E., 8, 97-111 (1992) · Zbl 0742.76051
[24] Turek, S., On discrete projection methods for the incompressible Navier-Stokes equations: an algorithmical approach, Comput. Methods Appl. Mech. Eng., 143, 271-288 (1997) · Zbl 0898.76069
[25] Turek, S., Efficient Solvers for Incompressible Flow Problems, An Algorithmic and Computational Approach, (Lecture Notes in Computational Science and Engineering, Vol. 6 (1999), Springer-Verlag), ISBN: 3-540-65433-X · Zbl 0930.76002
[26] Unverdi, S. O.; Tryggvason, G., A front-tracking method for viscous, incompressible, multi-fluid flows, J. Comput. Phys., 100, 25-37 (1992) · Zbl 0758.76047
[27] Zhang, Z., Polynomial preserving gradient recovery and a posteriori estimate for bilinear element on irregular quadrilaterals, Int. J. Num. Anal. Mod., 1, 1, 1-24 (2004) · Zbl 1071.65145
[28] Zienkiewicz, O. C.; Zhu, J. Z., The superconvergent patch recovery and a posteriori error estimators. part 1: The recovery technique, Int. J. Num. Meth. Eng., 33, 1331-1364 (1992) · Zbl 0769.73084
[29] FeatFlow CFD benchmark repository: http://www.featflow.de/en/benchmarks/cfdbenchmarking.html; FeatFlow CFD benchmark repository: http://www.featflow.de/en/benchmarks/cfdbenchmarking.html
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.