Indian Journal of Science and Technology
Year: 2019, Volume: 12, Issue: 20, Pages: 1-8
A. G. Shaikh1*, A. H. Sheikh1, Asif Ali2 and Salman Zeb3
1Department of Mathematics, Quaid-e-Awam University of Engineering, Science and Technology, Nawabshah Sindh, Pakistan; [email protected], [email protected]
2Mehran University of Engineering and Technology, Jamshoro, Pakistan; [email protected]
3Department of Mathematics, Malakand University, Lower Dir Khyber Pakhtunkhwa, Pakistan; [email protected]
*Author for correspondence
A. G. Shaikh
Department of Mathematics, Quaid-e-Awam University of Engineering, Science and Technology, Nawabshah Sindh, Pakistan.
Email: [email protected]
Objectives: The study reviews the iterative solvers for solving Helmholtz equation. It also surveys the preconditioners for Krylov subspace methods for discrete Helmholtz equation. Numerical experiments based conclusions are drawn, after a thorough comparison of various preconditioners. Methods/Statistical Analysis: Five point finite difference scheme for two-dimensional Helmholtz equation is employed to obtain discrete analogue of Helmholtz equation. Krylov subspace method for is considered for iterations. Matrix based preconditioner ILU with variety of fill-in and operator based preconditioner complex shifted Laplace preconditioner (CSLP) are tested. Different shifts are proposed and experimented in CSLP. Optimization of shifts is recommended on basis of results. Spectral analysis is performed to validate results. Findings: For small wave number problem, ILU type preconditioned performs better and is fairly comparable with CSLP. However, increasing wave number makes unaffordable in terms of computational cost. The pure imaginary shift (0,1) and complex shift (1, / 4 π ) costsless time and number of iterations compared to other choices of shift. These choices are recommended when preconditioner is inverted directly. Approximate inversion may affect iterations, but will save solve time for large problems. Extensive spectral analysis, validating results, highlights this work. Application/Improvements: Dependency upon wave number is decreased with this recommended choice of shifts in CSLP. The number of iterations increases at very low rate when wave number is increased at rate of 2.
Keywords: Helmholtz Equation, Iterative Methods, Krylov Subspace Methods, Preconditioners, Spectral Analysis
Subscribe now for latest articles and news.