Solving pdes in complex geometries

We extend previous work and present a general approach for solving partial differential equations in complex, stationary, or moving geometries with Dirichlet, Neumann, and Robin boundary conditions. Using an implicit representation of the geometry through an auxilliary phase field function, which replaces the sharp boundary of the domain with a diffuse layer (e.g. diffuse domain), the equation is reformulated on a larger regular domain. The resulting partial differential equation is of the same order as the original equation, with additional lower order terms to approximate the boundary conditions. The reformulated equation can be solved by standard numerical techniques. We use the method of matched asymptotic expansions to show that solutions of the reformulated equations converge to those of the original equations. We provide numerical simulations which confim this analysis. We also present applications of the method to growing domains and complex three-dimensional structures and we discuss applications to cell biology and heteroepitaxy.

Keywords

Partial differential equations, phase field approximation, complex geometry, diffuse interface, adaptive finite element methods, adaptive finite difference methods, multigrid methods

2010 Mathematics Subject Classification

35B40, 35K57, 65Mxx, 82C24

Full Text (PDF format)

Published 1 January 2009