Fourth-Order Interface Tracking in Two Dimensions Via an Improved Polygonal Area Mapping Method

Qinghai Zhang,Aaron Fogelson

DOI: https://doi.org/10.1137/140951886

IF: 2.968

2014-01-01

SIAM Journal on Scientific Computing

Abstract:We present an improved PAM (iPAM) methodas the first fourth-order interface tracking methodwhose convergence rates areindependent of $C^1$ (derivative) discontinuities of the interface.As an improved versionof the polygonal area mapping (PAM) method[Q. Zhang and P. L.-F. Liu, J. Comput. Phys. ,227 (2008), pp. 4063--4088],the accuracy of the iPAM method is achievedvia (i) augmenting the abstract data structure of PAMto faithfully represent multiple components of material regionswithin a single cell,(ii) removing restrictive assumptions of PAM,(iii) adjusting the volume of represented cell material regionsvia polygon ear removal,and (iv)maintaining a relation $(h_L=r_hh^{\alpha})$between the Eulerian grid size $h$and the Lagrangian length scale that measuresthe distance between adjacent interface markers.Unlike volume-of-fluid methods and level-set methods,merging and separation of the tracked materialis not automatically handled by the iPAM method.Instead, flexible subgrid resolutions are providedso that algorithmic behaviorsof interface merging and separation can be determinedfrom the specific physics of the problem under study.The iPAM method is a product of interdisciplinary researchof different fieldssuch as ordinary differential equations,computational geometry,and general topology.Its superior accuracy, efficiency, and versatilityover previous interface tracking methodsare demonstrated by a set of benchmark tests.In particular,for the vortex shear test on the $128\times 128$ grid,the fourth-order iPAM methodcould be hundreds of timesmore efficient than state-of-the-art volume-of-fluid methods.