Abstract

A new algorithm for the volume tracking of interfaces in two dimensions is presented. The algorithm is based upon a well-defined, second-order geometric solution of a volume evolution equation. The method utilizes local discrete material volume and velocity data to track interfaces of arbitrarily complex topology. A linearity-preserving, piecewise linear interface geometry approximation ensures that solutions generated retain second-order spatial accuracy. Second-order temporal accuracy is achieved by virtue of a multidimensional unsplit time integration scheme. We detail our geometrically based solution method, in which material volume fluxes are computed systematically with a set of simple geometric tasks. We then interrogate the method by testing its ability to track interfaces through large, controlled topology changes, whereby an initially simple interface configuration is subjected to vortical flows. Numerical results for these strenuous test problems provide evidence for the algorithm's improved solution quality and accuracy.

Links and resources

Tags