Efficient MFS Algorithms for Problems in Thermoelasticity

摘要

We propose efficient fast Fourier transform (FFT)-based algorithms using the method of fundamental solutions (MFS) for the numerical solution of certain problems in planar thermoelasticity. In particular, we consider problems in domains possessing radial symmetry, namely disks and annuli and it is shown that the MFS matrices arising in such problems possess circulant or block-circulant structures. The solution of the resulting systems is facilitated by appropriately diagonalizing these matrices using FFTs. Numerical experiments demonstrating the applicability of these algorithms are also presented.