A robust solver for elliptic PDEs in 3D complex geometries

02/11/2020
by   Matthew J. Morse, et al.
0

We develop a boundary integral equation solver for elliptic partial differential equations on complex 3D geometries. Our method is high-order accurate with optimal O(N) complexity and robustly handles complex geometries. A key component is our singular and near-singular layer potential evaluation scheme, hedgehog : a simple extrapolation of the solution along a line to the boundary. We present a series of geometry-processing algorithms required for hedgehog to run efficiently with accuracy guarantees on arbitrary geometries and an adaptive upsampling scheme based on a iteration-free heuristic for quadrature error that incorporates surface curvature. We validate the accuracy and performance with a series of numerical tests and compare our approach to a competing local evaluation method.

READ FULL TEXT

page 22

page 27

page 42

research
08/30/2019

Solution of Stokes flow in complex nonsmooth 2D geometries via a linear-scaling high-order adaptive integral equation scheme

We present a fast, high-order accurate and adaptive boundary integral sc...
research
03/31/2021

Analytical computation of boundary integrals for the Helmholtz equation in three dimensions

A key issue in the solution of partial differential equations via integr...
research
08/20/2019

A Fast Integral Equation Method for the Two-Dimensional Navier-Stokes Equations

The integral equation approach to partial differential equations (PDEs) ...
research
09/18/2021

Quadrature by fundamental solutions: kernel-independent layer potential evaluation for large collections of simple objects

Well-conditioned boundary integral methods for the solution of elliptic ...
research
11/03/2021

A reduced order Schwarz method for nonlinear multiscale elliptic equations based on two-layer neural networks

Neural networks are powerful tools for approximating high dimensional da...

Please sign up or login with your details

Forgot password? Click here to reset