Loading…
A mixed hybrid finite element method for the Helmholtz equation
In this paper, a high-order hybrid mixed finite element method for solving the two-dimensional Helmholtz equation with high wave numbers is presented. The novelty of this method is in using a discrete eigenfunction basis for solving the hybrid problem. Such a basis allows an inexpensive elimination...
Saved in:
Published in: | Journal of modern optics 2011-03, Vol.58 (5-6), p.424-437 |
---|---|
Main Authors: | , , |
Format: | Article |
Language: | English |
Subjects: | |
Citations: | Items that this one cites Items that cite this one |
Online Access: | Get full text |
Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
Summary: | In this paper, a high-order hybrid mixed finite element method for solving the two-dimensional Helmholtz equation with high wave numbers is presented. The novelty of this method is in using a discrete eigenfunction basis for solving the hybrid problem. Such a basis allows an inexpensive elimination of the inner degrees of freedom, which considerably reduces the size of the resulting linear system. On a rectangular grid with hanging nodes, the eigenfunction basis is constructed by solving a one-dimensional eigenvalue problem for each pair of edge length and polynomial order in the mesh. The eigenvalue problem can be solved efficiently for polynomial orders up to one thousand. Together with the reduced size of the linear system, this makes it possible to work with very high order basis functions, and consequently high frequency waves can be resolved on a coarse mesh. The hp-refinement is used to obtain an accurate solution for a minimal number of degrees of freedom. The effectiveness of our approach is demonstrated with numerical examples using polynomials of the degree up to one thousand. |
---|---|
ISSN: | 0950-0340 1362-3044 |
DOI: | 10.1080/09500340.2010.527067 |