Abstract
We develop a High Order Fourier solver for nonseparable, selfadjoint elliptic equations with variable (diffusion) coefficients. The solution of an auxiliary constant coefficient equation, serves in a transformation of the dependent variable. There results a ”modified Helmholtz” elliptic equation with almost constant coefficients. The small deviations from constancy are treated as correction terms. We developed a highly accurate, fast, Fourier-spectral algorithm to solve such constant coefficient equations. A small number of correction steps is required in order to achieve very high accuracy. This is achieved by optimization of the coefficients in the auxiliary equation. For given coefficients the approximation error becomes smaller as the domain decreases. A highly parallelizable hierarchical procedure allows a decomposition into smaller sub-domains where the solution is efficiently computed. This step is followed by hierarchical matching to reconstruct the global solution. Numerical experiments illustrate the high accuracy of the approach even at coarse resolutions.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
Preview
Unable to display preview. Download preview PDF.
References
Averbuch, A., Vozovoi, L., Israeli, M.: On a Fast Direct Elliptic Solver by a Modified Fourier Method. Numerical Algorithms 15, 287–313 (1997)
Skölermo, G.: A Fourier method for numerical solution of Poisson’s equation. Mathematics of Computation 29(131), 697–711 (1975)
Concus, P., Golub, G.H.: Use of fast direct methods for the efficient numerical solution of nonseparable elliptic equations. SIAM J. Numer. Anal. 10(6), 1103–1120 (1973)
Israeli, M., Braverman, E., Averbuch, A.: A hierarchical domain decomposition method with low communication overhead. In: Domain decomposition methods in science and engineering, Lyon, 2000, pp. 395-403. Theory Eng. Appl. Comput. Methods, Internat. Center Numer. Methods Eng (CIMNE), Barcelona (2002)
Author information
Authors and Affiliations
Editor information
Rights and permissions
Copyright information
© 2007 Springer Berlin Heidelberg
About this paper
Cite this paper
Israeli, M., Sherman, A. (2007). High Order Fourier-Spectral Solutions to Self Adjoint Elliptic Equations. In: Daydé, M., Palma, J.M.L.M., Coutinho, Á.L.G.A., Pacitti, E., Lopes, J.C. (eds) High Performance Computing for Computational Science - VECPAR 2006. VECPAR 2006. Lecture Notes in Computer Science, vol 4395. Springer, Berlin, Heidelberg. https://doi.org/10.1007/978-3-540-71351-7_29
Download citation
DOI: https://doi.org/10.1007/978-3-540-71351-7_29
Publisher Name: Springer, Berlin, Heidelberg
Print ISBN: 978-3-540-71350-0
Online ISBN: 978-3-540-71351-7
eBook Packages: Computer ScienceComputer Science (R0)