Poisson's Equation on the US Map

us map

The region is specified by a description of the boundary, called a skeleton . The skeleton also describes any internal interfaces, in this case state lines.

mesh0a mesh0b

This initial mesh was created from the skeleton. It is colored by element diameter. All internal interfaces are incorporated into the mesh as element edges. On the right is the same mesh, but without the triangle boundaries drawn in black.

mesh1 error

We solved Poisson's equation (-Laplace u = 1) in this domain, with homogeneous Dirichlet boundary conditions. On the left is the final mesh with 100,000 vertices, created using adaptive local mesh refinement. On the right is the final a posteriori error estimate for this problem.

soln1 sonl2

Here are some perspective plots of the solution.