Numerical solution of partial differential equations with Powell-Sabin splines

  • Authors:
  • Hendrik Speleers;Paul Dierckx;Stefan Vandewalle

  • Affiliations:
  • Department of Computer Science, Katholieke Universiteit Leuven, Celestijnenlaan 200A, B-3001 Leuven, Belgium;Department of Computer Science, Katholieke Universiteit Leuven, Celestijnenlaan 200A, B-3001 Leuven, Belgium;Department of Computer Science, Katholieke Universiteit Leuven, Celestijnenlaan 200A, B-3001 Leuven, Belgium

  • Venue:
  • Journal of Computational and Applied Mathematics
  • Year:
  • 2006

Quantified Score

Hi-index 7.29

Visualization

Abstract

Powell-Sabin splines are piecewise quadratic polynomials with global C^1-continuity. They are defined on conforming triangulations of two-dimensional domains, and admit a compact representation in a normalized B-spline basis. Recently, these splines have been used successfully in the area of computer-aided geometric design for the modelling and fitting of surfaces. In this paper, we discuss the applicability of Powell-Sabin splines for the numerical solution of partial differential equations defined on domains with polygonal boundary. A Galerkin-type PDE discretization is derived for the variable coefficient diffusion equation. Special emphasis goes to the treatment of Dirichlet and Neumann boundary conditions. Finally, an error estimator is developed and an adaptive mesh refinement strategy is proposed. We illustrate the effectiveness of the approach by means of some numerical experiments.