Communications in Mathematical Sciences

Volume 13 (2015)

Analysis of the diffuse-domain method for solving PDEs in complex geometries

Pages: 1473 – 1500

DOI: https://dx.doi.org/10.4310/CMS.2015.v13.n6.a6

Authors

Karl Yngve Lervåg (Department of Energy and Process Engineering, Norwegian University of Science and Technology, Trondheim, Norway; and SINTEF Energy Research, Trondheim, Norway)

John Lowengrub (Department of Mathematics, University of California at Irvine)

Abstract

In recent work, Li et al. [Commun. Math. Sci., 7, 81-107, 2009] developed a diffusedomain method (DDM) for solving partial differential equations in complex, dynamic geometries with Dirichlet, Neumann, and Robin boundary conditions. The diffuse-domain method uses an implicit representation of the geometry where the sharp boundary is replaced by a diffuse layer with thickness $\epsilon$ that is typically proportional to the minimum grid size. The original equations are reformulated on a larger regular domain and the boundary conditions are incorporated via singular source terms. The resulting equations can be solved with standard finite difference and finite element software packages. Here, we present a matched asymptotic analysis of general diffuse-domain methods for Neumann and Robin boundary conditions. Our analysis shows that for certain choices of the boundary condition approximations, the DDM is second-order accurate in $\epsilon$. However, for other choices the DDM is only first-order accurate. This helps to explain why the choice of boundary-condition approximation is important for rapid global convergence and high accuracy. Our analysis also suggests correction terms that may be added to yield more accurate diffuse-domain methods. Simple modifications of first-order boundary condition approximations are proposed to achieve asymptotically second-order accurate schemes. Our analytic results are confirmed numerically in the $L^2$ and $L^{\infty}$ norms for selected test problems.

Keywords

numerical solution of partial differential equations, phase-field approximation, implicit geometry representation, matched asymptotic analysis

2010 Mathematics Subject Classification

35B40, 35K51, 35K52, 35K57, 65Mxx

Published 13 May 2015