We develop a numerical method for realizing mean curvature motion of interfaces separating multiple phases, whose volumes are preserved throughout time. The foundation of the method is a thresholding algorithm of the Bence-Merriman-Osher type. The original algorithm is reformulated in a vector setting, which allows for a natural inclusion of constraints, even in the multiphase case. Moreover, a new method for overcoming the inaccuracy of thresholding methods on non-adaptive grids is designed, since this inaccuracy becomes especially prominent in volume-preserving motions. Formal analysis of the method and numerical tests are presented. © 2013 Elsevier B.V. All rights reserved.
Svadlenka, K., Ginder, E., & Omata, S. (2014). A variational method for multiphase volume-preserving interface motions. Journal of Computational and Applied Mathematics, 257, 157–179. https://doi.org/10.1016/j.cam.2013.08.027