A Schwarz Method for the Magnetotelluric Approximation of Maxwell's equations

The magnetotelluric approximation of the Maxwell's equations is used to model the propagation of low frequency electro-magnetic waves in the Earth's subsurface, with the purpose of reconstructing the presence of mineral or oil deposits. We propose a classical Schwarz method for solving this magnetotelluric approximation of the Maxwell equations, and prove its convergence using maximum principle techniques. This is not trivial, since solutions are complex valued, and we need a new result that the magnetotelluric approximations satisfy a maximum modulus principle for our proof. We illustrate our analysis with numerical experiments.Comment: 9 pages, 3 figure

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

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

Parallelizing the Kolmogorov-Fokker-Planck Equation

We design two parallel schemes, based on Schwarz Waveform Relaxation (SWR) procedures, for the numerical solution of the Kolmogorov equation. The latter is a simplified version of the Fokker-Planck equation describing the time evolution of the probability density of the velocity of a particle. SWR procedures decompose the spatio- temporal computational domain into subdomains and solve (in parallel) subproblems, that are coupled through suitable conditions at the interfaces to recover the solution of the global problem. We consider coupling conditions of both Dirichlet (Classical SWR) and Robin (Optimized SWR) types. We prove well-posedeness of the schemes subproblems and convergence for the proposed algorithms. We corroborate our findings with some numerical tests

Frozen Gaussian Approximation based domain decomposition methods for the linear Schrödinger equation beyond the semi-classical regime

International audienceThe paper is devoted to develop efficient domain decomposition methods for the linear Schrödinger equation beyond the semiclassical regime, which does not carry a small enough rescaled Planck constant for asymptotic methods (e.g. geometric optics) to produce a good accuracy, but which is too computationally expensive if direct methods (e.g. finite difference) are applied. This belongs to the category of computing middle-frequency wave propagation, where neither asymptotic nor direct methods can be directly used with both efficiency and accuracy. Motivated by recent works of the authors on absorbing boundary conditions [X. Antoine et al, J. Comput. Phys., 277 (2014), 268–304] and [X. Yang and J. Zhang, SIAM J. Numer. Anal., 52 (2014), 808–831], we introduce Semiclassical Schwarz Waveform Relaxation methods (SSWR), which are seamless integrations of semiclassical approximation to Schwarz Waveform Relaxation methods. Two versions are proposed respectively based on Herman-Kluk propagation and geometric optics, and we prove the convergence and provide numerical evidence of efficiency and accuracy of these methods