Solving pdes in complex geometries

Abstract

Abstract. 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 confirm 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.Key words. Partial differential equations, phase field approximation, complex geometry, diffuse interface, adaptive finite element methods, adaptive finite difference methods, multigrid methods.

Extended smoothed boundary method for solving partial differential equations with general boundary conditions on complex boundaries

“…However, the method was only applicable to no-flux boundary conditions, and no further extensions to other types of equations or boundary conditions have been reported. Recently, Lowengrub and coworkers [27,28,29,30,31,32,33] developed an alternative formulation for solving partial differential equations with various boundary conditions, based on asymptotic analyses commonly conducted in phase field modeling, which is different from the general derivation of the smoothed boundary method presented in this paper. Although such an implementation for imposing boundary conditions differs from the 'formal' practice suggested by Cahn [17], it dramatically simplifies the formulation, provides a justification of the method, and increases the applicability of the approach.…”

Continuum modelling of semiconductor heteroepitaxy: an applied perspective

“…To avoid solving the mechanical equilibrium equation with the free-surface boundary condition (see Section 2.2) in a time varying domain, a diffuse-interface approach can be conveniently used [141,157,158]. Thus, the problem is extended to the whole domain by introducing material properties dependent on ϕ:…”

A depth semi-averaged model for coastal dynamics

“…To make the scheme well-balanced, we implemented the surface gradient of Zhou et al 19 , which prescribes the reconstruction of the free-surface signal η rather than the total water Neumann conditions have been enforced by using a diffuse boundary approach, as described in Li X. et al 22 . In brief, the Poisson equation in (17) is solved over a domain D δ that is obtained by diffusing the original domain D in the normal outer direction over a reference smoothing length δ.…”