On a High Order Numerical Method for Solving Partial Differential Equations in Complex Geometries

  • Authors:
  • Knut S. Eckhoff

  • Affiliations:
  • Department of Mathematics, University of Bergen, Johs. Brunsgt. 12, N-5008 Bergen, Norway

  • Venue:
  • Journal of Scientific Computing
  • Year:
  • 1997

Quantified Score

Hi-index 0.01

Visualization

Abstract

Initial-boundary value problems for hyperbolic and parabolic partial differential equations with Diriclilet boundary conditions are considered by the method of lines approach in an arbitrarily given domain D in 2-D or 3-D. With D embedded in a rectangular domain, a new high order method for the space discretization problem is constructed in D by employing a Fourier collocation method in a uniform Cartesian system of gridpoints. Singularities are systematically removed by utilizing properties of the Bernoulli polynomials. Theoretical estimates for the accuracy of the method are established. The estimates are confirmed by numerical experiments for simple approximation problems.