Coupling of the ultra-weak variational formulation and an integral representation using a fast multipole method in electromagnetism
type: Item
instance of: scholarly article
title: Coupling of the ultra-weak variational formulation and an integral representation using a fast multipole method in electromagnetism
zbMATH Open document ID: 1136.78008
published in: Journal of Computational and Applied Mathematics
publication date: 2007-06-04T00:00:00Z
review text: The ultra-weak variational formulation (UWVF) is a volume-based numerical method for solving the time-harmonic Maxwell system on a bounded domain developed by Despres and Cessenat. It uses local plane wave solutions on a finite element mesh to approximate the field. By varying the number of plane wave basis functions from element to element the UWVF can discretize the electromagnetic field with a coarser volume mesh in comparison to more classical methods like low order finite elements or finite differences. However, to approximate scattering on an unbounded domain, the UWVF requires an artificial boundary \(\Gamma_{\text{ext}}\), sufficiently far from the obstacle. A simple absorbing boundary condition on \(\Gamma_{\text{ext}}\), (as used in the original UWVF) implies a large domain around the obstacle and so a large number of degrees of freedom. In this paper the authors suggest using an integral representation of the unknown field on \(\Gamma_{\text{extt}}\) due to Hazard and Lenoir. In particular, their idea is to use an integral representation of the unknown on the artificial boundary \(\Gamma_{\text{ext}}\) thanks to the unknown field values on a third boundary \(\Sigma\) taken closer to the boundary \(\Gamma_{\text{int}}\). This couples the degrees of freedom on \(\Sigma\) to those on \(\Gamma_{\text{ext}}\). The main constraint is that the domain between \(\Sigma\) and \(\Gamma_{\text{ext}}\) be homogeneous (i.e. the background medium). The artificial boundary \(\Gamma_{\text{ext}}\) can then be taken very close to the boundary of the obstacle. However, this method requires the evaluation of integral operators which are expensive by direct means. They show that the use of the integral representation does not greatly increase the cost of the UWVF if the integral calculation is performed using a fast multipole method (FMM). Indeed, for a given accuracy, the numerical complexity of the new algorithm is better than the complexity of the standard UWVF with a simp…
Mathematics Subject Classification ID: 78A45
Mathematics Subject Classification ID: 78M30
zbMATH DE Number: 5161213
zbMATH Keywords: ultra-weak variational formulation
zbMATH Keywords: absorbing boundary conditions (ABC)
zbMATH Keywords: integral representation
zbMATH Keywords: fast multipole method
author: E. Darrigrand
author: Peter Monk
MaRDI profile type: Publication
full work available at URL: j.cam.2006.03.030
OpenAlex ID: W2075552902
cites work: Q4341079
cites work: Using Plane Waves as Base Functions for Solving Time Harmonic Equations with the Ultra Weak Variational Formulation
cites work: Inverse acoustic and electromagnetic scattering theory.
cites work: Coupling of fast multipole method and microlocal discretization for the 3-D Helmholtz equation
cites work: The fast multipole method: Numerical implementation
cites work: On the Solution of Time-Harmonic Scattering Problems for Maxwell’s Equations
cites work: A novel hybridization of higher order finite element and boundary integral methods for electromagnetic scattering and radiation problems
DOI: 10.1016/J.CAM.2006.03.030
description: scientific article; zbMATH DE number 5161213