Results 1 - 10
of
40
CONVERGENCE ANALYSIS FOR FINITE ELEMENT DISCRETIZATIONS OF THE HELMHOLTZ EQUATION WITH DIRICHLET-TO-NEUMANN BOUNDARY CONDITIONS
"... Abstract. A rigorous convergence theory for Galerkin methods for a model Helmholtz problem in R d, d ∈{1, 2, 3} is presented. General conditions on the approximation properties of the approximation space are stated that ensure quasi-optimality of the method. As an application of the general theory, ..."
Abstract
-
Cited by 21 (7 self)
- Add to MetaCart
(Show Context)
Abstract. A rigorous convergence theory for Galerkin methods for a model Helmholtz problem in R d, d ∈{1, 2, 3} is presented. General conditions on the approximation properties of the approximation space are stated that ensure quasi-optimality of the method. As an application of the general theory, a full error analysis of the classical hp-version of the finite element method (hp-FEM) is presented for the model problem where the dependence on the mesh width h, the approximation order p, andthewavenumberk is given explicitly. In particular, it is shown that quasi-optimality is obtained under the conditions that kh/p is sufficiently small and the polynomial degree p is at least O(log k). 1.
Wave-Number Explicit Convergence Analysis for Galerkin Discretizations of the Helmholtz Equation
- SIAM J. Numer. Anal
"... In this paper, we develop a new stability and convergence theory for highly indefinite elliptic partial differential equations by considering the Helmholtz equation at high wave number as our model problem. The key element in this theory is a novel k-explicit regularity theory for Helmholtz boundary ..."
Abstract
-
Cited by 20 (7 self)
- Add to MetaCart
In this paper, we develop a new stability and convergence theory for highly indefinite elliptic partial differential equations by considering the Helmholtz equation at high wave number as our model problem. The key element in this theory is a novel k-explicit regularity theory for Helmholtz boundary value problems that is based on decomposing the solution into in two parts: the first part has the H2-Sobolev regularity expected of elliptic PDEs but features k-independent regularity constants; the second part is an analytic function for which k-explicit bounds for all derivatives are given. This decomposition is worked out in detail for several types of boundary value problems including the case Robin boundary conditions in domains with analytic boundary and in convex polygons. As the most important practical application we apply our full error analysis to the classical hp-version of the finite element method (hp-FEM) where the dependence on the mesh width h, the approximation order p, and the wave number k is given explicitly. In particular, under the assumption that the solution operator for Helmholtz problems grows only polynomially in k, it is shown that quasi-optimality is obtained under the conditions that kh/p is sufficiently small and the polynomial degree p is at least O(log k).
A class of discontinuous Petrov–Galerkin methods. Part IV: the optimal test norm and time-harmonic wave propagation in 1D
, 2011
"... Abstract The phase error, or the pollution effect in the finite element solution of wave propagation problems, is a well known phenomenon that must be confronted when solving problems in the highfrequency range. This paper presents a new method with no phase errors for one-dimensional (1D) time-har ..."
Abstract
-
Cited by 17 (8 self)
- Add to MetaCart
(Show Context)
Abstract The phase error, or the pollution effect in the finite element solution of wave propagation problems, is a well known phenomenon that must be confronted when solving problems in the highfrequency range. This paper presents a new method with no phase errors for one-dimensional (1D) time-harmonic wave propagation problems using new ideas that hold promise for the multidimensional case. The method is constructed within the framework of the Discontinuous Petrov-Galerkin (DPG) method with optimal test functions. We have previously shown that such methods select solutions that are the best possible approximations in an energy norm dual to any selected test space norm. In this paper, we advance by asking what is the optimal test space norm that achieves error reduction in a given energy norm. This is answered in the specific case of the Helmholtz equation with L 2 -norm as the energy norm. We obtain uniform stability with respect to the wave number. We illustrate the method with a number of 1D numerical experiments, using discontinuous piecewise polynomial hp spaces for the trial space and its corresponding optimal test functions computed approximately and locally. A 1D theoretical stability analysis is also developed.
Absolutely stable local discontinuous Galerkin methods for the Helmholtz equation with large wave number
- Mathematics of Computation
"... Abstract. Two local discontinuous Galerkin (LDG) methods using some non-standard numerical fluxes are developed for the Helmholtz equation with the first order absorbing boundary condition in the high frequency regime. It is shown that the proposed LDG methods are absolutely stable (hence well-posed ..."
Abstract
-
Cited by 13 (0 self)
- Add to MetaCart
(Show Context)
Abstract. Two local discontinuous Galerkin (LDG) methods using some non-standard numerical fluxes are developed for the Helmholtz equation with the first order absorbing boundary condition in the high frequency regime. It is shown that the proposed LDG methods are absolutely stable (hence well-posed) with respect to both the wave number and the mesh size. Optimal order (with respect to the mesh size) error estimates are proved for all wave numbers in the preasymptotic regime. To analyze the proposed LDG methods, they are recasted and treated as (non-conforming) mixed finite element methods. The crux of the analysis is to establish a generalized inf-sup condition, which holds without any mesh constraint, for each LDG method. The generalized inf-sup conditions then easily infer the desired absolute stability of the proposed LDG methods. In return, the stability results not only guarantee the well-posedness of the LDG methods but also play a crucial role in the derivation of the error estimates. Numerical experiments, which confirm the theoretical results and compare the proposed two LDG methods, are also presented in the paper. 1.
AN EXPONENTIALLY CONVERGENT NONPOLYNOMIAL FINITE ELEMENT METHOD FOR TIME-HARMONIC SCATTERING FROM POLYGONS
"... Abstract. In recent years nonpolynomial finite element methods have received increasing attention for the efficient solution of wave problems. As with their close cousin the method of particular solutions, high efficiency comes from using solutions to the Helmholtz equation as basis functions. We pr ..."
Abstract
-
Cited by 7 (3 self)
- Add to MetaCart
(Show Context)
Abstract. In recent years nonpolynomial finite element methods have received increasing attention for the efficient solution of wave problems. As with their close cousin the method of particular solutions, high efficiency comes from using solutions to the Helmholtz equation as basis functions. We present and analyze such a method for the scattering of two-dimensional scalar waves from a polygonal domain that achieves exponential convergence purely by increasing the number of basis functions in each element. Key ingredients are the use of basis functions that capture the singularities at corners, and representing the scattered field towards infinity by a combination of fundamental solutions. The solution is obtained by minimizing a least-squares functional, which we discretize in such a way that a matrix least-squares problem is obtained. We give computable exponential bounds on the rate of convergence of the least-squares functional that are in very good agreement with the observed numerical convergence. Challenging numerical examples, including a nonconvex polygon with several corner singularities, and a cavity domain, are solved to around 10 digits of accuracy with a few seconds of CPU time. The examples are implemented concisely with MPSpack, a MATLAB toolbox for wave computations with nonpolynomial basis functions, developed by the authors. A code example is included.
HYBRIDIZING RAVIART-THOMAS ELEMENTS FOR THE HELMHOLTZ EQUATION
"... Abstract. This paper deals with the application of hybridized mixed methods for discretizing the Helmholtz problem. Starting from a mixed formulation, where the flux is considered a separate unknown, we use Raviart-Thomas finite elements to approximate the solution. We present two ways of hybridizin ..."
Abstract
-
Cited by 6 (0 self)
- Add to MetaCart
(Show Context)
Abstract. This paper deals with the application of hybridized mixed methods for discretizing the Helmholtz problem. Starting from a mixed formulation, where the flux is considered a separate unknown, we use Raviart-Thomas finite elements to approximate the solution. We present two ways of hybridizing the problem, which means breaking the normal continuity of the fluxes and then imposing continuity weakly via functions supported on the element faces or edges. The first method is the Ultra-Weak Variational Formulation, first introduced by Cessenat and Després [7]; the second one uses Lagrange multipliers on element interfaces. We compare the two methods, and give numerical results. We observe that the iterative solvers applied to the two methods behave well for large wave numbers. 1.
Generalized DG-Methods for Highly Indefinite Helmholtz Problems based on the Ultra-Weak Variational Formulation
"... We develop a stability and convergence theory for the Ultra Weak Variational Formulation (UWVF) of a highly indefinite Helmholtz problem in R d, d ∈ {1, 2, 3}. The theory covers conforming as well as nonconforming generalized finite element methods. In contrast to conventional Galerkin methods where ..."
Abstract
-
Cited by 4 (3 self)
- Add to MetaCart
(Show Context)
We develop a stability and convergence theory for the Ultra Weak Variational Formulation (UWVF) of a highly indefinite Helmholtz problem in R d, d ∈ {1, 2, 3}. The theory covers conforming as well as nonconforming generalized finite element methods. In contrast to conventional Galerkin methods where a minimal resolution condition is necessary to guarantee the unique solvability, it is proved that the UWVF admits a unique solution under much weaker conditions. As an application we present the error analysis for the hp-version of the finite element method explicitly in terms of the mesh width h, polynomial degree p and wave number k. It is shown that the optimal convergence order estimate is obtained under the conditions that kh / √ p is sufficiently small and the polynomial degree p is at least O(log k). AMS Subject Classifications: 35J05, 65N12, 65N30 Key words: Helmholtz equation at high wavenumber, stability, convergence, discontinuous
Error analysis of Trefftz-discontinuous Galerkin methods for the time-harmonic Maxwell equations
- Math. Comp
, 2013
"... Abstract. In this paper, we extend to the time-harmonic Maxwell equations the p–version analysis technique developed in [R. Hiptmair, A. Moiola and I. Perugia, Plane wave discontinuous Galerkin methods for the 2D Helmholtz equation: analysis of the p-version, SIAM J. Numer. Anal., 49 (2011), 264-284 ..."
Abstract
-
Cited by 4 (0 self)
- Add to MetaCart
(Show Context)
Abstract. In this paper, we extend to the time-harmonic Maxwell equations the p–version analysis technique developed in [R. Hiptmair, A. Moiola and I. Perugia, Plane wave discontinuous Galerkin methods for the 2D Helmholtz equation: analysis of the p-version, SIAM J. Numer. Anal., 49 (2011), 264-284] for Trefftz-discontinuous Galerkin approximations of the Helmholtz problem. While error estimates in a mesh-skeleton norm are derived parallel to the Helmholtz case, the derivation of estimates in a mesh-independent norm re-quires new twists in the duality argument. The particular case where the local Trefftz approximation spaces are built of vector-valued plane wave functions is considered, and convergence rates are derived. 1.
On the negative-order norm accuracy of a local-structure-preserving LDG method
- J. Sci. Comput
"... Abstract The accuracy in negative-order norms is examined for a local-structure-preserving local discontinuous Galerkin method for the Laplace equation [Li and Shu, Methods and Applications of Analysis, v13 (2006), pp.215-233]. With its distinctive feature in using harmonic polynomials as local app ..."
Abstract
-
Cited by 3 (0 self)
- Add to MetaCart
(Show Context)
Abstract The accuracy in negative-order norms is examined for a local-structure-preserving local discontinuous Galerkin method for the Laplace equation [Li and Shu, Methods and Applications of Analysis, v13 (2006), pp.215-233]. With its distinctive feature in using harmonic polynomials as local approximating functions, this method has lower computational complexity than the standard local discontinuous Galerkin method while keeping the same order of accuracy in both the energy and the L 2 norms. In this note, numerical experiments are presented to demonstrate some accuracy loss of the method in negative-order norms.