A high-precision calculation method for interface normal and curvature on an unstructured grid

• We develop a new height function method which works on an unstructured grid.
• Geometric relationships are used to calculate interface normal and curvature.
• The calculation accuracy of interface normal shows roughly second-order accuracy.
• The numerical error in interface curvature shows the monotonic decreasing trends with cell size.

In the volume-of-fluid algorithm, the calculations of the interface normal and curvature are crucially important for accurately simulating interfacial flows. However, few methods have been proposed for the high-precision interface calculation on an unstructured grid. In this paper, the authors develop a height function method that works appropriately on an unstructured grid. In the process, the definition of the height function is discussed, and the high-precision calculation method of the interface normal is developed to meet the necessary condition for a second-order method. This new method has highly reduced computational cost compared with a conventional high-precision method because the interface normal calculation is completed by solving relatively simple algebraic equations. The curvature calculation method is also discussed and the approximated quadric curve of an interface is employed to calculate the curvature. Following a basic verification, the developed height function method is shown to successfully provide superior calculation accuracy and highly reduced computational cost compared with conventional calculation methods in terms of the interface normal and curvature. In addition, the height function method succeeds in calculating accurately the slotted-disk revolution problem and the oscillating drop on unstructured grids. Therefore, the developed height function method is confirmed to be an efficient technique for the high-precision numerical simulation of interfacial flows on an unstructured grid.