uu.seUppsala University Publications
Change search
ReferencesLink to record
Permanent link

Direct link
Domain Decomposition and Preconditioned Iterative Methods for the Helmholtz Equation
Uppsala University, Disciplinary Domain of Science and Technology, Mathematics and Computer Science, Department of Information Technology, Division of Scientific Computing. Uppsala University, Disciplinary Domain of Science and Technology, Mathematics and Computer Science, Department of Information Technology, Numerical Analysis. (ANLA)
2000 (English)Doctoral thesis, comprehensive summary (Other academic)
Abstract [en]

New preconditioned iterative solution methods for the Helmholtz equation are constructed. To evaluate the performance of the methods, two-dimensional problems with a waveguide geometry are used as model problems. Second- and fourth-order accurate finite difference discretizations are used. By introducing a domain decomposition framework, problems where the physical domain consists of a number of layers of different materials can be solved. The new algorithms are shown to perform well compared with standard methods.

First, a discretization of the Helmholtz equation in a curvilinear waveguide with smoothly varying material properties is studied. Nonlocal radiation boundary conditions are constructed for the artificial in- and outflow boundaries. A preconditioner that can be applied using fast transform methods is constructed. Experiments show that the total arithmetic complexity is much less with the fast transform preconditioner than when using a standard symmetric successive over-relaxation preconditioner. Compared with band Gaussian elimination, the gain is large both in arithmetic complexity and memory requirements.

Next, an application with layers of different materials is considered. It is shown that nonlocal radiation boundary conditions can be constructed also in this case. A domain decomposition formulation is employed, where the fast transform preconditioner is used as a subdomain preconditioner. The performance is excellent compared with band Gaussian elimination.

The multilayer solver is used for an investigation of the properties of the solutions of a number of underwater acoustics problems. The quality of the solutions obtained by a widely used approximation, the parabolic wave equation, is evaluated. For some problems the agreement is excellent, whereas in other cases the approximation cannot replace the full Helmholtz solution.

A parallel version of the solver is implemented, where the algorithms are modified for increased parallel performance. With the parallel code larger problems can be solved in less computational time.

Place, publisher, year, edition, pages
Uppsala: Acta Universitatis Upsaliensis , 2000. , 22 p.
Comprehensive Summaries of Uppsala Dissertations from the Faculty of Science and Technology, ISSN 1104-232X ; 523
Keyword [en]
Helmholtz equation, domain decomposition, preconditioned iterative method, fast transforms, parallelization
National Category
Computational Mathematics
Research subject
Numerical Analysis
URN: urn:nbn:se:uu:diva-441ISBN: 91-554-4689-2OAI: oai:DiVA.org:uu-441DiVA: diva2:164867
Public defence
2000-05-19, Room 2347, Polacksbacken, Uppsala University, Uppsala, 13:15 (English)
Available from: 2000-04-28 Created: 2000-04-28 Last updated: 2011-10-26Bibliographically approved

Open Access in DiVA

No full text
Buy this publication >>

Search in DiVA

By author/editor
Larsson, Elisabeth
By organisation
Division of Scientific ComputingNumerical Analysis
Computational Mathematics

Search outside of DiVA

GoogleGoogle Scholar
The number of downloads is the sum of all downloads of full texts. It may include eg previous versions that are now no longer available

Total: 840 hits
ReferencesLink to record
Permanent link

Direct link