• norsk
    • English
  • English 
    • norsk
    • English
  • Login
View Item 
  •   Home
  • Øvrige samlinger
  • Publikasjoner fra CRIStin - NTNU
  • View Item
  •   Home
  • Øvrige samlinger
  • Publikasjoner fra CRIStin - NTNU
  • View Item
JavaScript is disabled for your browser. Some features of this site may not work without it.

Calculation of the interface curvature and normal vector with the level-set method

Lervåg, Karl Yngve; Müller, Bernhard; Munkejord, Svend Tollak
Journal article
Submitted version
Thumbnail
View/Open
curvatures.pdf (3.095Mb)
URI
http://hdl.handle.net/11250/2468881
Date
2013
Metadata
Show full item record
Collections
  • Institutt for energi og prosessteknikk [2690]
  • Publikasjoner fra CRIStin - NTNU [20842]
Original version
Computers & Fluids. 2013, 84 (Sept), 218-230.   10.1016/j.compfluid.2013.06.004
Abstract
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 Macklin and Lowengrub [1]. 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 [1]. The present scheme is compared with the second-order central-difference scheme and with Macklin and Lowengrub’s scheme [1], 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 Macklin and Lowengrub in [2] 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 [1]. The results also demonstrate that the direction difference scheme [2] is not always sufficient to accurately calculate the normal vectors.
Publisher
Elsevier
Journal
Computers & Fluids

Contact Us | Send Feedback

Privacy policy
DSpace software copyright © 2002-2019  DuraSpace

Service from  Unit
 

 

Browse

ArchiveCommunities & CollectionsBy Issue DateAuthorsTitlesSubjectsDocument TypesJournalsThis CollectionBy Issue DateAuthorsTitlesSubjectsDocument TypesJournals

My Account

Login

Statistics

View Usage Statistics

Contact Us | Send Feedback

Privacy policy
DSpace software copyright © 2002-2019  DuraSpace

Service from  Unit