2,600 research outputs found
Quadrature by Expansion: A New Method for the Evaluation of Layer Potentials
Integral equation methods for the solution of partial differential equations,
when coupled with suitable fast algorithms, yield geometrically flexible,
asymptotically optimal and well-conditioned schemes in either interior or
exterior domains. The practical application of these methods, however, requires
the accurate evaluation of boundary integrals with singular, weakly singular or
nearly singular kernels. Historically, these issues have been handled either by
low-order product integration rules (computed semi-analytically), by
singularity subtraction/cancellation, by kernel regularization and asymptotic
analysis, or by the construction of special purpose "generalized Gaussian
quadrature" rules. In this paper, we present a systematic, high-order approach
that works for any singularity (including hypersingular kernels), based only on
the assumption that the field induced by the integral operator is locally
smooth when restricted to either the interior or the exterior. Discontinuities
in the field across the boundary are permitted. The scheme, denoted QBX
(quadrature by expansion), is easy to implement and compatible with fast
hierarchical algorithms such as the fast multipole method. We include accuracy
tests for a variety of integral operators in two dimensions on smooth and
corner domains
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
- …