Higher-order accurate diffuse-domain methods for partial differential equations with Dirichlet boundary conditions in complex, evolving geometries

Abstract

The diffuse-domain, or smoothed boundary, method is an attractive approach for solving partial differential equations in complex geometries because of its simplicity and flexibility. In this method the complex geometry is embedded into a larger, regular domain. The original PDE is reformulated using a smoothed characteristic function of the complex domain and source terms are introduced to approximate the boundary conditions. The reformulated equation, which is independent of the dimension and domain geometry, can be solved by standard numerical methods and the same solver can be used for any domain geometry. A challenge is making the method higher-order accurate. For Dirichlet boundary conditions, which we focus on here, current implementations demonstrate a wide range in their accuracy but generally the methods yield at best first order accuracy in $\epsilon$, the parameter that characterizes the width of the region over which the characteristic function is smoothed. Typically, $\epsilon \propto h$, the grid size. Here, we analyze the diffuse-domain PDEs using matched asymptotic expansions and explain the observed behaviors. Our analysis also identifies simple modifications to the diffuse-domain PDEs that yield higher-order accuracy in $\epsilon$, e.g., $O({\epsilon }^2)$ in the $L^2$ norm and $O({\epsilon }^p)$ with $1.5\le p\le 2$ in the $L^{\infty }$ norm. Our analytic results are confirmed numerically in stationary and moving domains where the level set method is used to capture the dynamics of the domain boundary and to construct the smoothed characteristic function.

References in corpus (9)

• Crawling and turning in a minimal reaction-diffusion cell motility model: coupling cell shape and biochemistry
• Periodic migration in a physical model of cells on micropatterns
• Diffuse interface models of locally inextensible vesicles in a viscous fluid
• Immersed Boundary Smooth Extension (IBSE): A high-order method for solving incompressible flows in arbitrary smooth domains
• An efficient mass-preserving interface-correction level set/ghost fluid method for droplet suspensions under depletion forces
• An Immersed Interface Method for Discrete Surfaces
• A dimensionally split Cartesian cut cell method for hyperbolic conservation laws
• An adaptive fast multipole accelerated Poisson solver for complex geometries
• Adaptive Boundary Element Methods

Cited by in corpus (2)

• Improving convergence of volume penalised fluid-solid interactions
• A consistent and conservative model and its scheme for $N$-phase-$M$-component incompressible flows