  A hybrid Eulerian-Lagrangian method is proposed to simulate passive scalar
transport on arbitrary shape interface. In this method, interface deformation
is tracked by an Eulerian method while the transport of the passive scalar on
the material interface is solved by a single-layer Lagrangian particle method.
To avoid particle clustering, a novel remeshing approach is proposed. This
remeshing method can resample particles, adjust the position of particles by a
relaxation process, and transfer mass from pre-existing particles to resampled
particles via a redistribution process, which preserves mass both globally and
locally. Computational costs are controlled by an adaptive remeshing strategy.
Accuracy is assessed by a series of test cases.
