2,397 research outputs found
Generalized Debye Sources Based EFIE Solver on Subdivision Surfaces
The electric field integral equation is a well known workhorse for obtaining
fields scattered by a perfect electric conducting (PEC) object. As a result,
the nuances and challenges of solving this equation have been examined for a
while. Two recent papers motivate the effort presented in this paper. Unlike
traditional work that uses equivalent currents defined on surfaces, recent
research proposes a technique that results in well conditioned systems by
employing generalized Debye sources (GDS) as unknowns. In a complementary
effort, some of us developed a method that exploits the same representation for
both the geometry (subdivision surface representations) and functions defined
on the geometry, also known as isogeometric analysis (IGA). The challenge in
generalizing GDS method to a discretized geometry is the complexity of the
intermediate operators. However, thanks to our earlier work on subdivision
surfaces, the additional smoothness of geometric representation permits
discretizing these intermediate operations. In this paper, we employ both ideas
to present a well conditioned GDS-EFIE. Here, the intermediate surface
Laplacian is well discretized by using subdivision basis. Likewise, using
subdivision basis to represent the sources, results in an efficient and
accurate IGA framework. Numerous results are presented to demonstrate the
efficacy of the approach
Numerical methods for computing Casimir interactions
We review several different approaches for computing Casimir forces and
related fluctuation-induced interactions between bodies of arbitrary shapes and
materials. The relationships between this problem and well known computational
techniques from classical electromagnetism are emphasized. We also review the
basic principles of standard computational methods, categorizing them according
to three criteria---choice of problem, basis, and solution technique---that can
be used to classify proposals for the Casimir problem as well. In this way,
mature classical methods can be exploited to model Casimir physics, with a few
important modifications.Comment: 46 pages, 142 references, 5 figures. To appear in upcoming Lecture
Notes in Physics book on Casimir Physic
High-frequency asymptotic compression of dense BEM matrices for general geometries without ray tracing
Wave propagation and scattering problems in acoustics are often solved with
boundary element methods. They lead to a discretization matrix that is
typically dense and large: its size and condition number grow with increasing
frequency. Yet, high frequency scattering problems are intrinsically local in
nature, which is well represented by highly localized rays bouncing around.
Asymptotic methods can be used to reduce the size of the linear system, even
making it frequency independent, by explicitly extracting the oscillatory
properties from the solution using ray tracing or analogous techniques.
However, ray tracing becomes expensive or even intractable in the presence of
(multiple) scattering obstacles with complicated geometries. In this paper, we
start from the same discretization that constructs the fully resolved large and
dense matrix, and achieve asymptotic compression by explicitly localizing the
Green's function instead. This results in a large but sparse matrix, with a
faster associated matrix-vector product and, as numerical experiments indicate,
a much improved condition number. Though an appropriate localisation of the
Green's function also depends on asymptotic information unavailable for general
geometries, we can construct it adaptively in a frequency sweep from small to
large frequencies in a way which automatically takes into account a general
incident wave. We show that the approach is robust with respect to non-convex,
multiple and even near-trapping domains, though the compression rate is clearly
lower in the latter case. Furthermore, in spite of its asymptotic nature, the
method is robust with respect to low-order discretizations such as piecewise
constants, linears or cubics, commonly used in applications. On the other hand,
we do not decrease the total number of degrees of freedom compared to a
conventional classical discretization. The combination of the ...Comment: 24 pages, 13 figure
Explicit Solution of the Time Domain Volume Integral Equation Using a Stable Predictor-Corrector Scheme
An explicit marching-on-in-time (MOT) scheme for solving the time domain volume integral equation is presented. The proposed method achieves its stability by employing, at each time step, a corrector scheme, which updates/corrects fields computed by the explicit predictor scheme. The proposedmethod is computationally more efficient when compared to the existing filtering techniques used for the stabilization of explicit MOT schemes. Numerical results presented in this paper demonstrate that the proposed method maintains its stability even when applied to the analysis of electromagnetic wave interactions with electrically large structures meshed using approximately half a million discretization elements
- …