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...

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...

As an example of a simple constrained geometric non-linear wave equation, we study a numerical approximation of the Maxwell Klein Gordon equation. We consider an existing constraint preserving semi-discrete scheme based on finite elements and prove its convergence in space dimension 2 for initial data of finite energy.

As an example of a simple constrained geometric non-linear wave equation, we study a numerical approximation of the Maxwell Klein Gordon equation. We consider an existing constraint preserving semi-discrete scheme based on finite elements and prove its convergence in space dimension 2 for initial data of finite energy.

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.

In this paper we prove the discrete compactness property for a discontinuous Galerkin approximation of Maxwell's system on quite general tetrahedral meshes. As a consequence, a discrete Friedrichs inequality is obtained and the convergence of the discrete eigenvalues to the continuous ones is deduced using the theory of collectively compact operators. Some numerical experiments confirm the theoretical predictions.

Two new time-dependent versions of div-curl results in a bounded domain $\Omega \subset {ℝ}^3$ are presented. We study a limit of the product $v_k{w}_k$, where the sequences $v_k$ and $w_k$ belong to ${Ł}_2\left(\Omega \right)$. In Theorem 2.1 we assume that $\nabla \times v_k$ is bounded in the $L_p$-norm and $\nabla ·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 ·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...

Three non-overlapping domain decomposition methods are proposed for the numerical approximation of time-harmonic Maxwell equations with damping (i.e., in a conductor). For each method convergence is proved and, for the discrete problem, the rate of convergence of the iterative algorithm is shown to be independent of the number of degrees of freedom.

Three non-overlapping domain decomposition methods are proposed for the numerical approximation of time-harmonic Maxwell equations with damping (i.e., in a conductor). For each method convergence is proved and, for the discrete problem, the rate of convergence of the iterative algorithm is shown to be independent of the number of degrees of freedom.

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...

In this paper we consider the Maxwell resolvent operator and its finite element approximation. In this framework it is natural the use of the edge element spaces and to impose the divergence constraint in a weak sense with the introduction of a Lagrange multiplier, following an idea by Kikuchi [14]. We shall review some of the known properties for edge element approximations and prove some new result. In particular we shall prove a uniform convergence in the $L^2$ norm for the sequence of discrete operators....

In this paper we consider the Maxwell resolvent operator and its finite element approximation. In this framework it is natural the use of the edge element spaces and to impose the divergence constraint in a weak sense with the introduction of a Lagrange multiplier, following an idea by Kikuchi [14]. We shall review some of the known properties for edge element approximations and prove some new result. In particular we shall prove a uniform convergence in the L2 norm for the sequence of discrete...

Since matrix compression has paved the way for discretizing the boundary integral equation formulations of electromagnetics scattering on very fine meshes, preconditioners for the resulting linear systems have become key to efficient simulations. Operator preconditioning based on Calderón identities has proved to be a powerful device for devising preconditioners. However, this is not possible for the usual first-kind boundary formulations for electromagnetic...

Since matrix compression has paved the way for discretizing the boundary integral equation formulations of electromagnetics scattering on very fine meshes, preconditioners for the resulting linear systems have become key to efficient simulations. Operator preconditioning based on Calderón identities has proved to be a powerful device for devising preconditioners. However, this is not possible for the usual first-kind boundary formulations for electromagnetic...

The electrowetting process is commonly used to handle very small amounts of liquid on a solid surface. This process can be modelled mathematically with the help of the shape optimization theory. However, solving numerically the resulting shape optimization problem is a very complex issue, even for reduced models that occur in simplified geometries. Recently, the second author obtained convincing results in the 2D axisymmetric case. In this paper, we propose and analyze a method that is suitable...

We present a numerical algorithm to solve the micromagnetic equations based on tangential-plane minimization for the magnetization update and a homothethic-layer decomposition of outer space for the computation of the demagnetization field. As a first application, detailed results on the flower-vortex transition in the cube of Micromagnetic Standard Problem number 3 are obtained, which confirm, with a different method, those already present in the literature, and validate our method and code. We...

We present a numerical algorithm to solve the micromagnetic equations based on tangential-plane minimization for the magnetization update and a homothethic-layer decomposition of outer space for the computation of the demagnetization field. As a first application, detailed results on the flower-vortex transition in the cube of Micromagnetic Standard Problem number 3 are obtained, which confirm, with a different method, those already present in the literature, and validate our method and...

We rigorously derive energy estimates for the second order vector wave equation with gauge condition for the electric field with non-constant electric permittivity function. This equation is used in the stabilized Domain Decomposition Finite Element/Finite Difference approach for time-dependent Maxwell’s system. Our numerical experiments illustrate efficiency of the modified hybrid scheme in two and three space dimensions when the method is applied for generation of backscattering data in the reconstruction...