Paper Info
Title
Parallel solution of high frequency Helmholtz equations using high order finite difference schemes.
Abstract
High frequency Helmholtz equations pose a challenging computational task, especially at high frequencies. This work examines the parallel solution of such problems, using 2nd, 4th and 6th order finite difference schemes. The examples include problems with known analytic solutions, enabling error evaluation of the different schemes on various grids, including problems on L-shaped domains. The resulting linear systems are solved with the block-parallel CARP-CG algorithm which succeeded in lowering the relative residual, indicating that it is a robust and reliable parallel solver of the resulting linear systems. However, lowering the error of the solution to reasonable levels with practical mesh sizes is obtained only with the higher order schemes. These results corroborate the known limitations of the 2nd order scheme at modeling the Helmholtz equation at high frequencies, and they indicate that CARP-CG can also be used effectively with high order finite difference schemes. Furthermore, the parallel scalability of CARP-CG improves when the wave number is increased (on a fixed grid), or when, for a fixed wave number, the grid size is decreased. (C) 2012 Elsevier Inc. All rights reserved.
Year
DOI
Venue
2012
10.1016/j.amc.2012.04.052
Applied Mathematics and Computation
Keywords
Field
DocType
CARP-CG,Helmholtz equation,Heterogeneous domain,High frequency,High order scheme,L-shaped domain,Marmousi,Parallel processing
Residual,Mathematical optimization,Linear system,Finite difference,Mathematical analysis,Helmholtz equation,Solver,Flux limiter,Mathematics,Grid,Scalability
Journal
Volume
Issue
ISSN
218
21
0096-3003
Citations 
PageRank 
References 
5
0.58
9
Authors
2
Name
Order
Citations
PageRank
Dan Gordon121021.44
Rachel Gordon218317.97