Source author record

J. Thomas Beale

J. Thomas Beale appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA math.AP Numerical Analysis

Catalog footprint

What is connected

4works

3topics

2close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2020arXiv

Solving Partial Differential Equations on Closed Surfaces with Planar Cartesian Grids

We present a general purpose method for solving partial differential equations on a closed surface, based on a technique for discretizing the surface introduced by Wenjun Ying and Wei-Cheng Wang [J. Comput. Phys. 252 (2013), pp. 606-624] which uses projections on coordinate planes. Assuming it is given as a level set, the surface is represented by a set of points at which it intersects the intervals between grid points in a three-dimensional grid. They are designated as primary or secondary. Discrete functions on the surface have independent values at primary points, with values at secondary points determined by an equilibration process. Each primary point and its neighbors have projections to regular grid points in a coordinate plane where the equilibration is done and finite differences are computed. The solution of a p.d.e. can be reduced to standard methods on Cartesian grids in the coordinate planes, with the equilibration allowing seamless transition from one system to another. We observe second order accuracy in examples with a variety of equations, including surface diffusion determined by the Laplace-Beltrami operator and the shallow water equations on a sphere.

preprint2018arXiv

Solution of the Dirichlet problem by a finite difference analog of the boundary integral equation

Several important problems in partial differential equations can be formulated as integral equations. Often the integral operator defines the solution of an elliptic problem with specified jump conditions at an interface. In principle the integral equation can be solved by replacing the integral operator with a finite difference calculation on a regular grid. A practical method of this type has been developed by the second author. In this paper we prove the validity of a simplified version of this method for the Dirichlet problem in a general domain in $R^2$ or $R^3$. Given a boundary value, we solve for a discrete version of the density of the double layer potential using a low order interface method. It produces the Shortley-Weller solution for the unknown harmonic function with accuracy $O(h^2)$. We prove the unique solvability for the density, with bounds in norms based on the energy or Dirichlet norm, using techniques which mimic those of exact potentials. The analysis reveals that this crude method maintains much of the mathematical structure of the classical integral equation. Examples are included.

preprint2016arXiv

A Simple Method for Computing Singular or Nearly Singular Integrals on Closed Surfaces

We present a simple, accurate method for computing singular or nearly singular integrals on a smooth, closed surface, such as layer potentials for harmonic functions evaluated at points on or near the surface. The integral is computed with a regularized kernel and corrections are added for regularization and discretization, which are found from analysis near the singular point. The surface integrals are computed from a new quadrature rule using surface points which project onto grid points in coordinate planes. The method does not require coordinate charts on the surface or special treatment of the singularity other than the corrections. The accuracy is about $O(h^3)$, where $h$ is the spacing in the background grid, uniformly with respect to the point of evaluation, on or near the surface. Improved accuracy is obtained for points on the surface. The treecode of Duan and Krasny for Ewald summation is used to perform sums. Numerical examples are presented with a variety of surfaces.

preprint2015arXiv

Uniform error estimates for Navier-Stokes flow with an exact moving boundary using the immersed interface method

We prove that uniform accuracy of almost second order can be achieved with a finite difference method applied to Navier-Stokes flow at low Reynolds number with a moving boundary, or interface, creating jumps in the velocity gradient and pressure. Difference operators are corrected to $O(h)$ near the interface using the immersed interface method, adding terms related to the jumps, on a regular grid with spacing $h$ and periodic boundary conditions. The force at the interface is assumed known within an error tolerance; errors in the interface location are not taken into account. The error in velocity is shown to be uniformly $O(h^2|\log{h}|^2)$, even at grid points near the interface, and, up to a constant, the pressure has error $O(h^2|\log{h}|^3)$. The proof uses estimates for finite difference versions of Poisson and diffusion equations which exhibit a gain in regularity in maximum norm.