### A stability analysis for finite volume schemes applied to the Maxwell system

Skip to main content (access key 's'),
Skip to navigation (access key 'n'),
Accessibility information (access key '0')

We present in this paper a stability study concerning finite volume schemes applied to the two-dimensional Maxwell system, using rectangular or triangular meshes. A stability condition is proved for the first-order upwind scheme on a rectangular mesh. Stability comparisons between the Yee scheme and the finite volume formulation are proposed. We also compare the stability domains obtained when considering the Maxwell system and the convection equation.

In this paper we discuss the approximate reconstruction of inhomogeneities of small volume. The data used for the reconstruction consist of boundary integrals of the (observed) electromagnetic fields. The numerical algorithms discussed are based on highly accurate asymptotic formulae for the electromagnetic fields in the presence of small volume inhomogeneities.

A new method for computation of the fundamental solution of electrodynamics for general anisotropic nondispersive materials is suggested. It consists of several steps: equations for each column of the fundamental matrix are reduced to a symmetric hyperbolic system; using the Fourier transform with respect to space variables and matrix transformations, formulae for Fourier images of the fundamental matrix columns are obtained; finally, the fundamental solution is computed by the inverse Fourier transform....

A Discontinuous Galerkin method is used for to the numerical solution of the time-domain Maxwell equations on unstructured meshes. The method relies on the choice of local basis functions, a centered mean approximation for the surface integrals and a second-order leap-frog scheme for advancing in time. The method is proved to be stable for cases with either metallic or absorbing boundary conditions, for a large class of basis functions. A discrete analog of the electromagnetic energy is conserved...

For solving the boundary-value problem for potential of a stationary magnetic field in two dimensions in ferromagnetics it is possible to use a linearization based on the succesive approximations. In this paper the convergence of this method is proved under some conditions.

Two new time-dependent versions of div-curl results in a bounded domain $\Omega \subset {\mathbb{R}}^{3}$ are presented. We study a limit of the product ${v}_{k}{w}_{k}$, where the sequences ${v}_{k}$ and ${w}_{k}$ belong to ${\u0141}_{2}\left(\Omega \right)$. In Theorem 2.1 we assume that $\nabla \times {v}_{k}$ is bounded in the ${L}_{p}$-norm and $\nabla \xb7{w}_{k}$ is controlled in the ${L}_{r}$-norm. In Theorem 2.2 we suppose that $\nabla \times {w}_{k}$ is bounded in the ${L}_{p}$-norm and $\nabla \xb7{w}_{k}$ is controlled in the ${L}_{r}$-norm. The time derivative of ${w}_{k}$ is bounded in both cases in the norm of $-1\left(\Omega \right)$. The convergence (in the sense of distributions) of ${v}_{k}{w}_{k}$ to the product $vw$ of weak limits...

During the development of a parallel solver for Maxwell equations by integral formulations and Fast Multipole Method (FMM), we needed to optimize a critical part including a lot of communications and computations. Generally, many parallel programs need to communicate, but choosing explicitly the way and the instant may decrease the efficiency of the overall program. So, the overlapping of computations and communications may be a way to reduce this drawback. We will see a implementation of this techniques...

During the development of a parallel solver for Maxwell equations by integral formulations and Fast Multipole Method (FMM), we needed to optimize a critical part including a lot of communications and computations. Generally, many parallel programs need to communicate, but choosing explicitly the way and the instant may decrease the efficiency of the overall program. So, the overlapping of computations and communications may be a way to reduce this drawback. We will see a implementation of this...

A simple and effective method based on Haar wavelets is proposed for the solution of Pocklington’s integral equation. The properties of Haar wavelets are first given. These wavelets are utilized to reduce the solution of Pocklington’s integral equation to the solution of algebraic equations. In order to save memory and computation time, we apply a threshold procedure to obtain sparse algebraic equations. Through numerical examples, performance of the present method is investigated concerning the...

In a posteriori error analysis of reduced basis approximations to affinely parametrized partial differential equations, the construction of lower bounds for the coercivity and inf-sup stability constants is essential. In [Huynh et al., C. R. Acad. Sci. Paris Ser. I Math.345 (2007) 473–478], the authors presented an efficient method, compatible with an off-line/on-line strategy, where the on-line computation is reduced to minimizing a linear functional under a few linear constraints. These constraints...

We investigate sufficient and possibly necessary conditions for the L2 stability of the upwind first order finite volume scheme for Maxwell equations, with metallic and absorbing boundary conditions. We yield a very general sufficient condition, valid for any finite volume partition in two and three space dimensions. We show this condition is necessary for a class of regular meshes in two space dimensions. However, numerical tests show it is not necessary in three space dimensions even on regular...

In order to get the fusion of small capsules containing a deuterium-tritium nuclear fuel, the MegaJoule laser (LMJ) will focus a large number of laser beams inside a cylinder (Hohlraum) which contains the fusion capsule. In order to control this process we have to know as well as possible the electromagnetic field created by the laser beams on both Hohlraum’s apertures. This article describes a numerical tool which computes this electromagnetic field...

For the Maxwell equations in time-dependent media only finite difference schemes with time-dependent conductivity are known. In this paper we present a numerical scheme based on the Magnus expansion and operator splitting that can handle time-dependent permeability and permittivity too. We demonstrate our results with numerical tests.

The topic of this paper is the numerical analysis of time periodic solution for electro-magnetic phenomena. The Limit Absorption Method (LAM) which forms the basis of our study is presented. Theoretical results have been proved in the linear finite dimensional case. This method is applied to scattering problems and transport of charged particles.