Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userConvergence analysis of oversampled collocation boundary element methods in 2D
92
0
0.0
(
0
)
Added by
Georg Maierhofer
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Collocation boundary element methods for integral equations are easier to implement than Galerkin methods because the elements of the discretization matrix are given by lower-dimensional integrals. For that same reason, the matrix assembly also requires fewer computations. However, collocation methods typically yield slower convergence rates and less robustness, compared to Galerkin methods. We explore the extent to which oversampled collocation can improve both robustness and convergence rates. We show that in some cases convergence rates can actually be higher than the corresponding Galerkin method, although this requires oversampling at a faster than linear rate. In most cases of practical interest, oversampling at least lowers the error by a constant factor. This can still be a substantial improvement: we analyze an example where linear oversampling by a constant factor $J$ (leading to a rectangular system of size $JN times N$) improves the error at a cubic rate in the constant $J$. Furthermore, the oversampled collocation method is much less affected by a poor choice of collocation points, as we show how oversampling can lead to guaranteed convergence. Numerical experiments are included for the two-dimensional Helmholtz equation.
rate research
Read More
In this paper, we examine the effectiveness of classic multiscale finite element method (MsFEM) (Hou and Wu, 1997; Hou et al., 1999) for mixed Dirichlet-Neumann, Robin and hemivariational inequality boundary problems. Constructing so-called boundary correctors is a common technique in existing methods to prove the convergence rate of MsFEM, while we think not reflects the essence of those problems. Instead, we focus on the first-order expansion structure. Through recently developed estimations in homogenization theory, our convergence rate is provided with milder assumptions and in neat forms.
A Sinc-collocation method has been proposed by Stenger, and he also gave theoretical analysis of the method in the case of a `scalar equation. This paper extends the theoretical results to the case of a `system of equations. Furthermore, this paper proposes more efficient method by replacing the variable transformation employed in Stengers method. The efficiency is confirmed by both of theoretical analysis and numerical experiments. In addition to the existing and newly-proposed Sinc-collocation methods, this paper also gives similar theoretical results for Sinc-Nystr{o}m methods proposed by Nurmuhammad et al. From a viewpoint of the computational cost, it turns out that the newly-proposed Sinc-collocation method is the most efficient among those methods.
In this paper we discuss a hybridised method for FEM-BEM coupling. The coupling from both sides use a Nitsche type approach to couple to the trace variable. This leads to a formulation that is robust and flexible with respect to approximation spaces and can easily be combined as a building block with other hybridised methods. Energy error norm estimates and the convergence of Jacobi iterations are proved and the performance of the method is illustrated on some computational examples.
We design and analyze a coupling of a discontinuous Galerkin finite element method with a boundary element method to solve the Helmholtz equation with variable coefficients in three dimensions. The coupling is realized with a mortar variable that is related to an impedance trace on a smooth interface. The method obtained has a block structure with nonsingular subblocks. We prove quasi-optimality of the $h$- and $
A number of non-standard finite element methods have been proposed in recent years, each of which derives from a specific class of PDE-constrained norm minimization problems. The most notable examples are $mathcal{L}mathcal{L}^*$ methods. In this work, we argue that all high-order methods in this class should be expected to deliver substandard uniform h-refinement convergence rates. In fact, one may not even see rates proportional to the polynomial order $p > 1$ when the exact solution is a constant function. We show that the convergence rate is limited by the regularity of an extraneous Lagrange multiplier variable which naturally appears via a saddle-point analysis. In turn, limited convergence rates appear because the regularity of this Lagrange multiplier is determined, in part, by the geometry of the domain. Numerical experiments support our conclusions.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
