Spectrally accurate solutions to inhomogeneous elliptic PDE in smooth geometries using function intension
arXiv:2203.01798 · doi:10.1016/j.jcp.2022.111594
Abstract
We present a spectrally accurate embedded boundary method for solving linear, inhomogeneous, elliptic partial differential equations (PDE) in general smooth geometries, focusing in this manuscript on the Poisson, modified Helmholtz, and Stokes equations. Unlike several recently proposed methods which rely on function extension, we propose a method which instead utilizes function `intension', or the smooth truncation of known function values. Similar to those methods based on extension, once the inhomogeneity is truncated we may solve the PDE using any of the many simple, fast, and robust solvers that have been developed for regular grids on simple domains. Function intension is inherently stable, as are all steps in the proposed solution method, and can be used on domains which do not readily admit extensions. We pay a price in exchange for improved stability and flexibility: in addition to solving the PDE on the regular domain, we must additionally (1) solve the PDE on a small auxiliary domain that is fitted to the boundary, and (2) ensure consistency of the solution across the interface between this auxiliary domain and the rest of the physical domain. We show how these tasks may be accomplished efficiently (in both the asymptotic and practical sense), and compare convergence to several recent high-order embedded boundary schemes.
References in corpus (12)
- SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
- Quadrature by Expansion: A New Method for the Evaluation of Layer Potentials
- Immersed Boundary Smooth Extension: A high-order method for solving PDE on arbitrary smooth domains using Fourier spectral methods
- A fast platform for simulating flexible fiber suspensions applied to cell mechanics
- Immersed Boundary Smooth Extension (IBSE): A high-order method for solving incompressible flows in arbitrary smooth domains
- A scalable computational platform for particulate Stokes suspensions
- An adaptive fast multipole accelerated Poisson solver for complex geometries
- A stable and accurate scheme for solving the Stefan problem coupled with natural convection using the Immersed Boundary Smooth Extension method
- Fourier continuation method for incompressible fluids with boundaries
- Partition of Unity Extension of Functions on Complex Domains
- A Fast Integral Equation Method for the Two-Dimensional Navier-Stokes Equations
- Hydrodynamics of a Semipermeable Vesicle Under Flow and Confinement