235 research outputs found
A simple preconditioned domain decomposition method for electromagnetic scattering problems
We present a domain decomposition method (DDM) devoted to the iterative
solution of time-harmonic electromagnetic scattering problems, involving large
and resonant cavities. This DDM uses the electric field integral equation
(EFIE) for the solution of Maxwell problems in both interior and exterior
subdomains, and we propose a simple preconditioner for the global method, based
on the single layer operator restricted to the fictitious interface between the
two subdomains.Comment: 23 page
A domain decomposition method for solving the three-dimensional time-harmonic Maxwell equations discretized by discontinuous Galerkin methods
We present here a domain decomposition method for solving the three-dimensional time-harmonic Maxwell equations discretized by a discontinuous Galerkin method. In order to allow the treatment of irregularly shaped geometries, the discontinuous Galerkin method is formulated on unstructured tetrahedral meshes. The domain decomposition strategy takes the form of a Schwarz-type algorithm where a continuity condition on the incoming characteristic variables is imposed at the interfaces between neighboring subdomains. A multifrontal sparse direct solver is used at the subdomain level. The resulting domain decomposition strategy can be viewed as a hybrid iterative/direct solution method for the large, sparse and complex coefficients algebraic system resulting from the discretization of the time-harmonic Maxwell equations by a discontinuous Galerkin method
Effective transmission conditions for domain decomposition methods applied to the time-harmonic curl-curl Maxwell's equations
The time-harmonic Maxwell equations describe the propagation of electromagnetic waves and are therefore fundamental for the simulation of many modern devices we have become used to in everyday life. The numerical solution of these equations is hampered by two fundamental problems: first, in the high frequency regime, very fine meshes need to be used in order to avoid the pollution effect well known for the Helmholtz equation, and second the large scale systems obtained from the vector valued equations in three spatial dimensions need to be solved by iterative methods, since direct factorizations are not feasible any more at that scale. As for the Helmholtz equation, classical iterative methods applied to discretized Maxwell equations have severe convergence problems.We explain in this paper a family of domain decomposition methods based on well chosen transmission conditions. We show that all transmission conditions proposed so far in the literature, both for the first and second order formulation of Maxwell's equations, can be written and optimized in the common framework of optimized Schwarz methods, independently of the first or second order formulation one uses, and the performance of the corresponding algorithms is identical. We use a decomposition into transverse electric and transverse magnetic fields to describe these algorithms, which greatly simplifies the convergence analysis of the methods. We illustrate the performance of our algorithms with large scale numerical simulations
An overlapping splitting double sweep method for the Helmholtz equation
We consider the domain decomposition method approach to solve the Helmholtz
equation. Double sweep based approaches for overlapping decompositions are
presented. In particular, we introduce an overlapping splitting double sweep
(OSDS) method valid for any type of interface boundary conditions. Despite the
fact that first order interface boundary conditions are used, the OSDS method
demonstrates good stability properties with respect to the number of subdomains
and the frequency even for heterogeneous media. In this context, convergence is
improved when compared to the double sweep methods in Nataf et al. (1997) and
Vion et al. (2014, 2016} for all of our test cases: waveguide, open cavity and
wedge problems
Discontinuous Galerkin discretizations of optimized Schwarz methods for solving the time-harmonic Maxwell equations
We show in this paper how to properly discretize optimized Schwarz methods for the time-harmonic Maxwell equations using a discontinuous Galerkin (DG) method. Due to the multiple traces between elements in the DG formulation, it is not clear a priori how the more sophisticated transmission conditions in optimized Schwarz methods should be discretized, and the most natural approach does not lead at convergence of the Schwarz method to the mono-domain DG discretization, which implies that for such discretizations, the DG error estimates do not hold when the Schwarz method has converged. We present an alternative discretization of the transmission conditions in the framework of a DG weak formulation, and prove that for this discretization the multidomain and mono-domain solutions for the Maxwell's equations are the same. We illustrate our results with several numerical experiments of propagation problems in homogeneous and heterogeneous media
Schwarz methods for second order Maxwell equations in 3D with coefficient jumps
We study non-overlapping Schwarz Methods for solving second order time-harmonic 3D Maxwell equations in heterogeneous media. Choosing the interfaces between the subdomains to be aligned with the discontinuities in the coefficients, we show for a model problem that while the classical Schwarz method is not convergent, optimized transmission conditions dependent on the discontinuities of the coefficients lead to convergent methods. We prove asymptotically that the resulting methods converge in certain cases independently of the mesh parameter, and convergence can even become better as the coefficient jumps increase
An Additive Schwarz Method Type Theory for Lions's Algorithm and a Symmetrized Optimized Restricted Additive Schwarz Method
International audienceOptimized Schwarz methods (OSM) are very popular methods which were introduced by P.L. Lions in [27] for elliptic problems and by B. Després in [8] for propagative wave phenomena. We give here a theory for Lions' algorithm that is the genuine counterpart of the theory developed over the years for the Schwarz algorithm. The first step is to introduce a symmetric variant of the ORAS (Optimized Restricted Additive Schwarz) algorithm [37] that is suitable for the analysis of a two-level method. Then we build a coarse space for which the convergence rate of the two-level method is guaranteed regardless of the regularity of the coefficients. We show scalability results for thousands of cores for nearly incompressible elasticity and the Stokes systems with a continuous discretization of the pressure
Recommended from our members
Schnelle Löser fĂŒr Partielle Differentialgleichungen
The workshop Schnelle Löser fĂŒr partielle Differentialgleichungen, organised by Randolph E. Bank (La Jolla), Wolfgang Hackbusch (Leipzig), and Gabriel Wittum (Frankfurt am Main), was held May 22ndâMay 28th, 2011. This meeting was well attended by 54 participants with broad geographic representation from 7 countries and 3 continents. This workshop was a nice blend of researchers with various backgrounds
- âŠ