Show simple item record

dc.contributor.author
Casati, Daniele
dc.contributor.supervisor
Hiptmair, Ralf
dc.contributor.supervisor
Smajic, Jasmin
dc.contributor.supervisor
Tsukerman, Igor
dc.date.accessioned
2019-12-23T13:46:06Z
dc.date.available
2019-12-20T09:24:29Z
dc.date.available
2019-12-20T14:28:11Z
dc.date.available
2019-12-23T13:46:06Z
dc.date.issued
2020-03
dc.identifier.isbn
978-3-906916-80-4
en_US
dc.identifier.uri
http://hdl.handle.net/20.500.11850/386816
dc.identifier.doi
10.3929/ethz-b-000386816
dc.description.abstract
We propose a hybrid method to solve time-harmonic Maxwell's equations in R^3 through the Finite Element Method (FEM) in a bounded region encompassing parameter inhomogeneities, coupled with the Multiple Multipole Program (MMP) in the unbounded complement. FEM and MMP enjoy complementary capabilities. FEM requires a mesh of the computational domain of interest. This is expensive, but can treat inhomogeneous materials, shapes with corners, or other singularities. Moreover, FEM allows a purely local construction of the discrete system of equations. On the other hand, MMP belongs to the class of methods of auxiliary sources and of Trefftz methods, as it employs point sources that are exact solutions of the homogeneous equations as (global) basis functions: one only needs to evaluate them on hypersurfaces to obtain the discrete problem, and the resulting linear combination is valid in the whole domain where the equations hold, which can be unbounded. MMP performs well where the electromagnetic field is smooth, i.e. in the free space far from physical sources and material interfaces. Thus, a natural way to combine the strengths of these methods arises when one needs to simulate the electromagnetic field of components with inhomogeneous parameters or nonsmooth shapes surrounded by free space: use FEM on a mesh defined on the components and MMP in the unbounded domain outside. The boundary between the FEM and MMP domains can be nonphysical if one surrounds the components by a conforming mesh of an "airbox" also modeled by FEM. The interface conditions on the surface of the FEM domain are key to accurate coupled FEM--MMP solutions. Integrating by parts the variational form solved by FEM, surface integrals appear, through which one can impose interface conditions by substituting the ansatz of MMP. However, one interface condition (for example, the continuity of the tangential trace for Maxwell's equations) cannot be imposed in this way. To include this additional condition, we derive four methods from Lagrangian functionals that enforce both the variational form in the FEM domain and all the interface conditions between different discretizations: 1. Least-squares-based coupling using techniques from PDE-constrained optimization. This approach optimizes a functional for the additional condition, subject to a constraint expressed by the variational form of FEM. 2. Discontinuous Galerkin (DG) between the meshed FEM domain and the single-entity "MMP mesh", interpreting MMP as FEM with special basis functions. 3. Multi-field variational formulation in the spirit of mortar finite element methods, relying on the same interpretation of MMP as the DG-based coupling. 4. Coupling through the Dirichlet-to-Neumann operator: here we introduce a weak formulation of the additional condition where MMP basis functions are used as test functions. This approach is generalized to couple MMP with any numerical method based on volume meshes that fits co-chain calculus. Furthermore, to minimize the meshed region and for the sake of generalization, we assume that equation parameters are piecewise constant in the MMP domain (instead of constant everywhere). This induces a partition: we consider the case of one unbounded subdomain with other bounded, but possibly very large, subdomains, each requiring its own MMP discretization space (but without the need of meshes). Hence, on top of the interface conditions between the FEM and MMP domains, one should also impose interface conditions between the MMP subdomains themselves. We compare the FEM--MMP coupling approaches in a series of numerical experiments with different geometries and equation parameters, including examples that exhibit triple-point singularities. Convergence tests are performed for examples that cover the whole spectrum of computational electromagnetics, obtaining the expected results.
en_US
dc.format
application/pdf
en_US
dc.language.iso
en
en_US
dc.publisher
ETH Zurich
en_US
dc.rights.uri
http://rightsstatements.org/page/InC-NC/1.0/
dc.subject
Finite Element Method
en_US
dc.subject
Multiple Multipole Program
en_US
dc.subject
Method of Auxiliary Sources
en_US
dc.subject
Trefftz methods
en_US
dc.subject
computational electromagnetics
en_US
dc.title
Coupling Finite Elements and Auxiliary Sources
en_US
dc.type
Doctoral Thesis
dc.rights.license
In Copyright - Non-Commercial Use Permitted
dc.date.published
2019-12-23
ethz.size
151 p.
en_US
ethz.code.ddc
DDC - DDC::5 - Science::510 - Mathematics
en_US
ethz.grant
Coupling of Finite Elements and Multiple Multipoles for Computational Electromagnetism
en_US
ethz.identifier.diss
26522
en_US
ethz.publication.place
Zurich
en_US
ethz.publication.status
published
en_US
ethz.leitzahl
ETH Zürich::00002 - ETH Zürich::00012 - Lehre und Forschung::00007 - Departemente::02000 - Dep. Mathematik / Dep. of Mathematics::02501 - Seminar für Angewandte Mathematik / Seminar for Applied Mathematics
en_US
ethz.leitzahl
ETH Zürich::00002 - ETH Zürich::00012 - Lehre und Forschung::00007 - Departemente::02000 - Dep. Mathematik / Dep. of Mathematics::02501 - Seminar für Angewandte Mathematik / Seminar for Applied Mathematics::03632 - Hiptmair, Ralf / Hiptmair, Ralf
en_US
ethz.leitzahl.certified
ETH Zürich::00002 - ETH Zürich::00012 - Lehre und Forschung::00007 - Departemente::02000 - Dep. Mathematik / Dep. of Mathematics::02501 - Seminar für Angewandte Mathematik / Seminar for Applied Mathematics::03632 - Hiptmair, Ralf / Hiptmair, Ralf
en_US
ethz.grant.agreementno
165674
ethz.grant.fundername
SNF
ethz.grant.funderDoi
10.13039/501100001711
ethz.grant.program
ethz.date.deposited
2019-12-20T09:24:38Z
ethz.source
FORM
ethz.eth
yes
en_US
ethz.availability
Open access
en_US
ethz.rosetta.installDate
2019-12-23T13:46:17Z
ethz.rosetta.lastUpdated
2020-02-15T23:04:59Z
ethz.rosetta.exportRequired
true
ethz.rosetta.versionExported
true
ethz.COinS
ctx_ver=Z39.88-2004&rft_val_fmt=info:ofi/fmt:kev:mtx:journal&rft.atitle=Coupling%20Finite%20Elements%20and%20Auxiliary%20Sources&rft.date=2020-03&rft.au=Casati,%20Daniele&rft.isbn=978-3-906916-80-4&rft.genre=unknown&rft_id=info:doi/978-3-906916-80-4&rft.btitle=Coupling%20Finite%20Elements%20and%20Auxiliary%20Sources
 Search via SFX

Files in this item

Thumbnail

Publication type

Show simple item record