Reinitialization of the Level-Set Function in 3d Simulation of Moving Contact Lines

Thumbnail Image



Journal Title

Journal ISSN

Volume Title

Repository Usage Stats


Citation Stats


The level set method is one of the most successful methods for the simulation of multi-phase flows. To keep the level set function close the signed distance function, the level set function is constantly reinitialized by solving a Hamilton-Jacobi type of equation during the simulation. When the fluid interface intersects with a solid wall, a moving contact line forms and the reinitialization of the level set function requires a boundary condition in certain regions on the wall. In this work, we propose to use the dynamic contact angle, which is extended from the contact line, as the boundary condition for the reinitialization of the level set function. The reinitialization equation and the equation for the normal extension of the dynamic contact angle form a coupled system and are solved simultaneously. The extension equation is solved on the wall and it provides the boundary condition for the reinitialization equation; the level set function provides the directions along which the contact angle is extended from the contact line. The coupled system is solved using the 3rd order TVD Runge-Kutta method and the Godunov scheme. The Godunov scheme automatically identifies the regions where the angle condition needs to be imposed. The numerical method is illustrated by examples in three dimensions.





Published Version (Please cite this version)


Publication Info

Xu, S, and W Ren (2016). Reinitialization of the Level-Set Function in 3d Simulation of Moving Contact Lines. Communications in Computational Physics, 20(5). pp. 1163–1182. 10.4208/cicp.210815.180316a Retrieved from

This is constructed from limited available data and may be imprecise. To cite this article, please review & use the official citation provided by the journal.



Shixin Xu

Assistant Professor of Mathematics at Duke Kunshan University

Shixin Xu is an Assistant Professor of Mathematics.   His research interests are machine learning and data-driven models for diseases,  multiscale modeling of complex fluids, Neurovascular coupling, homogenization theory, and numerical analysis.  The current projects he is working on are

  • image data-based for the prediction of hemorrhagic transformation in acute ischemic stroke,
  • electrodynamics modeling of saltatory conduction along a myelinated axon
  • electrochemical modeling
  • fluid-structure interaction with mass transportation and reaction

Unless otherwise indicated, scholarly articles published by Duke faculty members are made available here with a CC-BY-NC (Creative Commons Attribution Non-Commercial) license, as enabled by the Duke Open Access Policy. If you wish to use the materials in ways not already permitted under CC-BY-NC, please consult the copyright owner. Other materials are made available here through the author’s grant of a non-exclusive license to make their work openly accessible.