Computing Reviews
Today's Issue Hot Topics Search Browse Recommended My Account Log In
Review Help
A fast direct solver for the biharmonic problem in a rectangular grid
Ben-Artzi M., Croisille J., Fishelov D. SIAM Journal on Scientific Computing31 (1):303-333,2008.Type:Article
Date Reviewed: Mar 26 2009

The biharmonic problem comes up in important application areas such as fluid dynamics. Methods to solve it have become more complex as computers have improved in speed and storage capacity. The methods are based on a finite-difference approximation to the given partial differential equation (PDE), leading to a system of linear algebraic equations.

This paper presents a method from Stephenson [1] that combines a second-order approximation of the solution with a fourth-order approximation of the gradient, using a nine-point stencil. Ben-Artzi et al. modify the lower order part so that the entire approximation is fourth order. Having defined the basic discrete operators, they construct a matrix representation from simple difference operators such as the symmetric tridiagonal matrix usually associated with the standard second difference operator, namely 2 on the diagonal and -1 on the off-diagonal. The modified Stephenson operator is then written as the sum of a matrix constructed from this standard matrix and a lower-rank matrix. Kronecker products are used to construct two-dimensional operators from the one-dimensional difference matrices. The resulting system is solved with a combination of the fast Fourier transform and the Sherman-Morrison formula. The cost of the solution is O(N²log N), where N is the number of grid points in one direction. Numerical results from a Fortran90 implementation of the method are given, together with error estimates.

Reviewer:  Charles R. Crawford Review #: CR136627 (0911-1060)
1) Stephenson, J. W. Single cell discretizations of order two and four for biharmonic problems. J. Comput. Phys. 55, (1984), 65–80.
Bookmark and Share
Would you recommend this review?
Other reviews under "Linear Systems (Direct And Iterative Methods)": Date
A node-addition model for symbolic factorization
Law K., Fenives S. ACM Transactions on Mathematical Software 12(1): 37-50, 1986. Type: Article
Apr 1 1987
A feature selection Newton method for support vector machine classification
Fung G., Mangasarian O. Computational Optimization and Applications 28(2): 185-202, 2004. Type: Article
Jun 29 2005
Computing the MDMT decomposition
Kaufman L. ACM Transactions on Mathematical Software 21(4): 476-489, 1995. Type: Article
Jul 1 1996

E-Mail This Printer-Friendly
Send Your Comments
Contact Us
Reproduction in whole or in part without permission is prohibited.   Copyright 2004™
Terms of Use
| Privacy Policy