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

JOURNAL OF COMPUTATIONAL PHYSICS

卷:406

Article

Yu, Fei¹  Guo, Zhenlin¹  Lowengrub, John^1,2 

[1] Univ Calif Irvine, Dept Math, Irvine, CA 92697 USA

[2] Univ Calif Irvine, Dept Biomed Engn, Irvine, CA 92697 USA

关键词: Partial differential equations;    Complex geometry;    Smoothed boundary method;    Matched asymptotic analysis;    Multigrid methods;    Adaptive mesh refinement;   

DOI  :  10.1016/j.jcp.2019.109174

来源: Elsevier

【 摘 要 】

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 c, the parameter that characterizes the width of the region over which the characteristic function is smoothed. Typically, epsilon proportional to 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 <= p <= 2 in the L-infinity 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. (C) 2019 Published by Elsevier Inc.

【 授权许可】

Free