不规则区域上泊松方程的快速求解器OA

Based on the Fourier spectral method,a fast solver was established for the numerical solution of the Poisson equation on the irregular domain.Considering the computation on the unit disc,by the polar coordinate transform and Fourier series expansion,a series of two point boundary value problems satisified by Fourier coefficients were obtained.We constructed the Chebyshev polynomial approximation of the right hand side of these boundary value problems and derived the generalized expressions of Fourier coefficients,solved the undetermined coefficients in the expressions using the boundary conditions.For the computation on the irregular domain,we embedded the irregular domain in the unit disc.With the help of the solving process on the unit disc,combined the generalized expressions of Fourier coefficients and discrete Fourier inverse transform,we constructed the system of linear equations which was satisfied by the boundary conditions and the undetermined coefficients.The advantage of the new scheme was by using the FFT and IFFT algorithm,it could improve the computational efficiency,reduced the scale of computation and attained the exponential convergence accuracy for the numerial solution on the irregular domain.These implied that the proposed scheme provided an effective way for the high precision numerical computation of the Poisson equation.