33,949 research outputs found
Rigorous analysis of internal resonances in 3-D hybrid FE-BIE formulations by means of the Poincaré-Steklov operator
3-D hybrid finite-element (FE) boundary integral equation (BIE) formulations are widely used because of their ability to simulate large inhomogeneous structures in both open and bounded simulation domains by applying each method where it is the most efficient. However, some formulations suffer from breakdown frequencies at which the solution is not uniquely defined and errors are introduced due to internal resonances. In this paper, we investigate the occurrence of spurious solutions resulting from these resonances by using the concept of the Poincare-Steklov or Dirichlet-to-Neumann operator, which provides a relation between the tangential electric field and the electric current on the boundary of a domain. By identifying this operator in both the FE and BIE method, several new properties of internal resonances in 3-D hybrid FE-BIE formulations are easily derived. Several conformal and nonconformal formulations are studied and the theory is then applied to a scattering problem
An introduction to Multitrace Formulations and Associated Domain Decomposition Solvers
Multitrace formulations (MTFs) are based on a decomposition of the problem
domain into subdomains, and thus domain decomposition solvers are of interest.
The fully rigorous mathematical MTF can however be daunting for the
non-specialist. We introduce in this paper MTFs on a simple model problem using
concepts familiar to researchers in domain decomposition. This allows us to get
a new understanding of MTFs and a natural block Jacobi iteration, for which we
determine optimal relaxation parameters. We then show how iterative multitrace
formulation solvers are related to a well known domain decomposition method
called optimal Schwarz method: a method which used Dirichlet to Neumann maps in
the transmission condition. We finally show that the insight gained from the
simple model problem leads to remarkable identities for Calderon projectors and
related operators, and the convergence results and optimal choice of the
relaxation parameter we obtained is independent of the geometry, the space
dimension of the problem{\color{black}, and the precise form of the spatial
elliptic operator, like for optimal Schwarz methods. We illustrate our analysis
with numerical experiments
Boundary element formulations for the numerical solution of two-dimensional diffusion problems with variable coefficients
This is the post-print version of the final paper published in Computers & Mathematics with Applications. The published article is available from the link below. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. Copyright @ 2012 Elsevier B.V.This paper presents new formulations of the radial integration boundary integral equation (RIBIE) and the radial integration boundary integro-differential equation (RIBIDE) methods for the numerical solution of two-dimensional diffusion problems with variable coefficients. The methods use either a specially constructed parametrix (Levi function) or the standard fundamental solution for the Laplace equation to reduce the boundary-value problem (BVP) to a boundary–domain integral equation (BDIE) or boundary–domain integro-differential equation (BDIDE). The radial integration method (RIM) is then employed to convert the domain integrals arising in both BDIE and BDIDE methods into equivalent boundary integrals. The resulting formulations lead to pure boundary integral and integro-differential equations with no domain integrals. Furthermore, a subdomain decomposition technique (SDBDIE) is proposed, which leads to a sparse system of linear equations, thus avoiding the need to calculate a large number of domain integrals. Numerical examples are presented for several simple problems, for which exact solutions are available, to demonstrate the efficiency of the proposed approaches
Fast Isogeometric Boundary Element Method based on Independent Field Approximation
An isogeometric boundary element method for problems in elasticity is
presented, which is based on an independent approximation for the geometry,
traction and displacement field. This enables a flexible choice of refinement
strategies, permits an efficient evaluation of geometry related information, a
mixed collocation scheme which deals with discontinuous tractions along
non-smooth boundaries and a significant reduction of the right hand side of the
system of equations for common boundary conditions. All these benefits are
achieved without any loss of accuracy compared to conventional isogeometric
formulations. The system matrices are approximated by means of hierarchical
matrices to reduce the computational complexity for large scale analysis. For
the required geometrical bisection of the domain, a strategy for the evaluation
of bounding boxes containing the supports of NURBS basis functions is
presented. The versatility and accuracy of the proposed methodology is
demonstrated by convergence studies showing optimal rates and real world
examples in two and three dimensions.Comment: 32 pages, 27 figure
- …