Smoothed Boundary Method for Solving Partial Differential Equations with General Boundary Conditions on Complex Boundaries

Abstract

In this article, we describe an approach for solving partial differential equations with general boundary conditions imposed on arbitrarily shaped boundaries. A function that has a prescribed value on the domain in which a differential equation is valid and smoothly but rapidly varying values on the boundary where boundary conditions are imposed is used to modify the original differential equations. The mathematical derivations are straight forward, and generically applicable to a wide variety of partial differential equations. To demonstrate the general applicability of the approach, we provide four examples: (1) the diffusion equation with both Neumann and Dirichlet boundary conditions, (2) the diffusion equation with surface diffusion, (3) the mechanical equilibrium equation, and (4) the equation for phase transformation with additional boundaries. The solutions for a few of these cases are validated against corresponding analytical and semi-analytical solutions. The potential of the approach is demonstrated with five applications: surface-reaction diffusion kinetics with a complex geometry, Kirkendall-effect-induced deformation, thermal stress in a complex geometry, phase transformations affected by substrate surfaces, and a self-propelling droplet.Comment: A better smooth algorithm has been developed and tested, will soon replace Eq. 58 in page 16. We have also developed a level-set moving boundary SBM method, and it will replace the Navier-Stokes-Cahn-Hilliard type domain parameter tracking method in Section 5.