A conservative level set method for two phase flow
2005 (English)In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 210, no 1, 225-246 p.Article in journal (Refereed) Published
A conservative method of level set type for moving interfaces in divergence free velocity fields is presented. The interface is represented implicitly by the 0.5 level set of a function Phi being a smeared out Heaviside function, i.e., a function being zero on one side of the interface and one on the other. In a transition layer of finite, constant thickness Phi goes smoothly from zero to one. The interface is moved implicitly by the advection of Phi, which is split into two steps. First Phi is advected using a standard numerical method. Then an intermediate step is performed to make sure that the smooth profile of Phi and the thickness of the transition layer is preserved. Both these steps are performed using conservative second order approximations and thus conserving integral Phi. In this way good conservation of the area bounded by the 0.5 contour of Phi is obtained.
Numerical tests shows up to second order accuracy and very good conservation of the area bounded by the interface.
The method was also coupled to a Navier-Stokes solver for incompressible two phase flow with surface tension. Results with and without topological changes are presented.
Place, publisher, year, edition, pages
2005. Vol. 210, no 1, 225-246 p.
level set method, free boundary, two phase flow, incompressible flow, conservative method
IdentifiersURN: urn:nbn:se:kth:diva-5362DOI: 10.1016/j.jcp.2005.04.007ISI: 000231265500009ScopusID: 2-s2.0-23044452940OAI: oai:DiVA.org:kth-5362DiVA: diva2:9707
QC 201011222006-02-142006-02-142010-11-22Bibliographically approved