Numerical solution of partial differential equations with Powell-Sabin splines

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.

[1]  Paul Dierckx,et al.  Surface fitting using convex Powell-Sabin splines , 1994 .

[2]  M. Rivara.  Selective refinement/derefinement algorithms for sequences of nested triangulations , 1989 .

[3]  P. Oswald,et al.  Hierarchical conforming finite element methods for the biharmonic equation , 1992 .

[4]  Joris Windmolders.  Powell-Sabin splines for computer aided geometric design , 2003 .

[5]  Philippe G. Ciarlet,et al.  The finite element method for elliptic problems , 2002, Classics in applied mathematics.

[6]  Malcolm A. Sabin,et al.  Piecewise Quadratic Approximations on Triangles , 1977, TOMS.

[7]  Charles K. Chui,et al.  Multivariate vertex splines and finite elements , 1990 .

[8]  Adhemar Bultheel,et al.  Powell-Sabin spline Wavelets , 2004, Int. J. Wavelets Multiresolution Inf. Process..

[9]  Paul Dierckx,et al.  Curve and surface fitting with splines , 1994, Monographs on numerical analysis.

[10]  C. Chui,et al.  Bivariate C 1 Quadratic Finite Elements and Vertex Splines , 1990 .

[11]  Paul Dierckx,et al.  On calculating normalized Powell-Sabin B-splines , 1997, Comput. Aided Geom. Des..

[12]  C. D. Boor,et al.  On Calculating B-splines , 1972 .

[13]  Gerald Farin,et al.  Triangular Bernstein-Bézier patches , 1986, Comput. Aided Geom. Des..

[14]  M. Rivara.  Algorithms for refining triangular grids suitable for adaptive and multigrid techniques , 1984 .

[15]  Adhemar Bultheel,et al.  On the choice of the PS-triangles , 2003 .

[16]  Larry L. Schumaker,et al.  Smooth Macro-Elements Based on Powell–Sabin Triangle Splits , 2002, Adv. Comput. Math..

[17]  Adhemar Bultheel,et al.  Uniform Powell--Sabin spline wavelets , 2003 .

[18]  Paul Dierckx,et al.  Algorithms for surface fitting using Powell-Sabin splines , 1992 .