Abstract
A two-step quadratic spline collocation method is formulated for the solution of the Dirichlet biharmonic problem on the unit square rewritten as a coupled system of two second-order partial differential equations. This method involves fast Fourier transforms and, in comparison to its one-step counterpart, it has the advantage of requiring the solution a symmetric positive definite Schur complement system rather than a nonsymmetric one. As a consequence, the corresponding step of the new method is performed using a preconditioned conjugate gradient method. The total cost of the method on a N × N partition of the unit square is \(O(N^{2}\log N)\). To demonstrate the optimal accuracy of the method, the results of numerical experiments are provided.
Similar content being viewed by others
References
Abushama, A.A., Bialecki, B.: Modified nodal cubic spline collocation for biharmonic equations. Numer. Algorithms 43, 331–353 (2006)
Ben-Artiz, M., Croisille, J.-P., Fishelov, J.: A fast direct solver for the biharmonic problem in a rectangular grid. SIAM J. Sci. Comput. 31, 303–333 (2008)
Bialecki, B.: A fast solver for the orthogonal spline collocation solution of the biharmonic Dirichlet problem on rectangles. J. Comput. Phys. 191, 601–621 (2003)
Bialecki, B.: A fourth order finite difference method for the Dirichlet biharmonic problem. Numer. Algorithms 61, 351–375 (2012)
Bialecki, B., Fairweather, G., Karageorghis, A.: Matrix decomposition algorithms for elliptic boundary value problems: a survey. Numer. Algorithms 56, 253–295 (2011)
Bialecki, B., Fairweather, G., Karageorghis, A., Maack, J.: A quadratic spline collocation method for the Dirichlet biharmonic problem. Numer. Algorithms 83, 165–199 (2020)
Christara, C.C.: Quadratic spline collocation methods for elliptic partial differential equations. BIT 34, 33–61 (1994)
Christara, C.C., Ng, K.S.: Fast Fourier transform solvers and preconditoner for quadratic spline collocation. BIT 42, 702–730 (2002)
Knudson, D.B.: A piecewise Hermite bicubic finite element Galerkin method for the biharmonic Dirichlet problem, Ph.D. Thesis, Colorado School of Mines, Golden, Colorado, 1997
Acknowledgements
The authors are grateful for the constructive remarks and useful suggestions of an anonymous referee which helped to improve the paper.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare no competing interests.
Additional information
Data sharing
Data sharing is not applicable to this article as no datasets were generated or analysed during the current study.
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Bialecki, B., Fairweather, G. & Karageorghis, A. An optimal two-step quadratic spline collocation method for the Dirichlet biharmonic problem. Numer Algor 91, 1115–1143 (2022). https://doi.org/10.1007/s11075-022-01294-y
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11075-022-01294-y
Keywords
- Biharmonic equation
- Quadratic spline collocation
- Fast Fourier transforms
- Preconditioned conjugate gradient method
- Optimal global convergence rates
- Superconvergence