Article ID | Journal | Published Year | Pages | File Type |
---|---|---|---|---|
9507149 | Applied Mathematics and Computation | 2005 | 17 Pages |
Abstract
We develop a simple and efficient FFT-based fast direct solver for the biharmonic equation on a disk. The biharmonic equation is split into a coupled system of harmonic problems. We first use the truncated Fourier series expansion to derive a set of coupled singular ODEs, then we solve those singular equations by second-order finite difference discretizations. Using a radial grid with shifting a half mesh away from the origin, we can handle the coordinate singularity easily without pole conditions. The Sherman-Morrison formula is then applied to solve the resultant linear system in a cost-efficient way. The computational complexity of the method consists of O(MNÂ log2Â N) arithmetic operations for MÂ ÃÂ N grid points. The numerical accuracy check and some applications to the incompressible Navier-Stokes flows inside a disk are conducted.
Related Topics
Physical Sciences and Engineering
Mathematics
Applied Mathematics
Authors
Ming-Chih Lai, Hsi-Chi Liu,