A surface marker algorithm coupled to an area-preserving marker redistribution method for three-dimensional interface tracking Eugenio Aulisa, Sandro Manservisi 1 , Ruben Scardovelli * INFM-BO and DIENCA, Lab. di Montecuccolino, Via dei Colli 16, 40136 Bologna, Italy Received 21 May 2003; received in revised form 10 December 2003; accepted 11 December 2003 Available online 23 January 2004 Abstract In this work we present a new interface tracking algorithm based on surface markers to reconstruct and advect interfaces in three-dimensional space. The algorithm is coupled to a marker redistribution method, which is area- preserving along the interface. The interface is described by a set of closed lines which define a coarse Lagrangian quadrangular mesh in the fixed computational grid. Fixed markers are located where the lines cross and are maintained during the whole simulation. At each time step the interface is first reconstructed and then all markers are advected by following streamlines. In the reconstruction step, new markers are determined by computing the intersections of the interface lines with the grid cell faces and by adding area conservation markers between fixed and grid intersection markers. Intersection and conservation markers defined at the previous time step are discarded. The method maintains a smooth geometrical description of the interface for both two-dimensional and three-dimensional tests with accurate volume conservation even in very challenging situations where the fluid bodies progressively deform and stretch de- veloping localized regions with very high curvature and thin fluid filaments. The method compares favorably with front capturing methods, such as volume-of-fluid and level set, and other hybrid techniques, such as the particle level set method. Ó 2004 Elsevier Inc. All rights reserved. PACS: 65C20; 68U20 Keywords: Multiphase flows; Interface reconstruction and advection; Front tracking and front capturing; Hybrid methods www.elsevier.com/locate/jcp Journal of Computational Physics 197 (2004) 555–584 * Corresponding author. Tel.: +39-51-644-1720; fax: +39-51-644-1747. E-mail address: raus@mail.ing.unibo.it (R. Scardovelli). 1 Present address: Department of Mathematics and Statistics, Texas Tech, Lubbock, TX, USA. 0021-9991/$ - see front matter Ó 2004 Elsevier Inc. All rights reserved. doi:10.1016/j.jcp.2003.12.009