A matrix-free Legendre spectral method for initial-boundary value problems

Electronic Transactions on Numerical Analysis

Bernd Brumm¹ 

关键词: spectral methods;    matrix-free methods;    method of lines;    stability;    computational wave propagation;    boundary conditions;   

DOI  :  

学科分类：数学（综合）

来源: Kent State University * Institute of Computational Mathematics

【 摘 要 】

We present a Legendre spectral method for initial-boundary value problems with variable coefficients and of arbitrary dimensionality, where the computational work in each time step scales linearly with the number of unknowns. Boundary conditions are enforced weakly, allowing for stable solutions of many classes of problems. Working in coefficient space, derivatives can be evaluated recursively in linear time. We show how also the action of variable coefficients can be implemented without transforming back to coordinate space using a recursive, linearly scaling matrix-free algorithm, under the assumption that the coefficients vary on a much longer scale than the solution. We also prove that spectral accuracy is preserved for smooth solutions. Numerical results for the wave equation in two and three dimensions corroborate the theoretical predictions.

【 授权许可】

Unknown