10,237 research outputs found
Efficient Fast Multipole Accelerated Boundary Elements via Recursive Computation of Multipole Expansions of Integrals
In boundary element methods (BEM) in , matrix elements and
right hand sides are typically computed via analytical or numerical quadrature
of the layer potential multiplied by some function over line, triangle and
tetrahedral volume elements. When the problem size gets large, the resulting
linear systems are often solved iteratively via Krylov subspace methods, with
fast multipole methods (FMM) used to accelerate the matrix vector products
needed. When FMM acceleration is used, most entries of the matrix never need be
computed explicitly - {\em they are only needed in terms of their contribution
to the multipole expansion coefficients.} We propose a new fast method for the
analytical generation of the multipole expansion coefficients produced by the
integral expressions for single and double layers on surface triangles; charge
distributions over line segments and over tetrahedra in the volume; so that the
overall method is well integrated into the FMM, with controlled error. The
method is based on the per moment cost recursive computation of the
moments. The method is developed for boundary element methods involving the
Laplace Green's function in . The derived recursions are first
compared against classical quadrature algorithms, and then integrated into FMM
accelerated boundary element and vortex element methods. Numerical tests are
presented and discussed.Comment: 6 figures, preprin
Isogeometric Boundary Elements in Electromagnetism: Rigorous Analysis, Fast Methods, and Examples
We present a new approach to three-dimensional electromagnetic scattering
problems via fast isogeometric boundary element methods. Starting with an
investigation of the theoretical setting around the electric field integral
equation within the isogeometric framework, we show existence, uniqueness, and
quasi-optimality of the isogeometric approach. For a fast and efficient
computation, we then introduce and analyze an interpolation-based fast
multipole method tailored to the isogeometric setting, which admits competitive
algorithmic and complexity properties. This is followed by a series of
numerical examples of industrial scope, together with a detailed presentation
and interpretation of the results
Fast Multipole Method for the Symmetric Boundary Element Method in MEG/EEG
The accurate solution of the forward electrostatic problem is an essential first step before solving the inverse problem of magneto- and electro-encephalography (MEG/EEG). The symmetric Galerkin boundary element method is accurate but is difficule to use for very large problems because of its computational complexity and memory requirements. We describe a fast multipole-based acceleration for the symmetric BEM with complexity. It creates a hierarchical structure of the elements and approximates far interactions using spherical harmonics expansions. The accelerated method is shown to be as accurate as the direct method, yet for large problems it is both faster and more economical in terms of memory consumption
The development of a fast multipole boundary element method for coupled acoustic and elastic problems
This thesis presents a dual fast multipole boundary element method (FMBEM) for modelling 3D acoustic coupled fluid-structure interaction problems in the frequency domain. Boundary integral representations are used to represent both the exterior fluid and interior elastic solid domains and the fast multipole method is employed to accelerate the calculations in both domains. The dual FMBEM yields a similar solution accuracy to the conventional models, while its solution times and memory requirements are substantially reduced
An Adaptive Fast Multipole Boundary Element Method for Poisson−Boltzmann Electrostatics
The numerical solution of the Poisson−Boltzmann (PB) equation is a useful but a computationally demanding tool for studying electrostatic solvation effects in chemical and biomolecular systems. Recently, we have described a boundary integral equation-based PB solver accelerated by a new version of the fast multipole method (FMM). The overall algorithm shows an order N complexity in both the computational cost and memory usage. Here, we present an updated version of the solver by using an adaptive FMM for accelerating the convolution type matrix-vector multiplications. The adaptive algorithm, when compared to our previous nonadaptive one, not only significantly improves the performance of the overall memory usage but also remarkably speeds the calculation because of an improved load balancing between the local- and far-field calculations. We have also implemented a node-patch discretization scheme that leads to a reduction of unknowns by a factor of 2 relative to the constant element method without sacrificing accuracy. As a result of these improvements, the new solver makes the PB calculation truly feasible for large-scale biomolecular systems such as a 30S ribosome molecule even on a typical 2008 desktop computer
- …