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.

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.

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.

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

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

This paper presents how the FEMLab package can be used to perform the magnetic field analysis in the Active Magnetic Bearing (AMB). The AMB is an integral part of the industrial rotational machine laboratory model. The electromagnetic field distribution and density analysis allow verifying the designed AMB and the influence of the shaft and coil current changes on the bearing parameters.

Low order edge elements are widely used for electromagnetic field problems. Higher order edge approximations are receiving increasing interest but their definition become rather complex. In this paper we propose a simple definition for Whitney edge elements of polynomial degree higher than one. We give a geometrical localization of all degrees of freedom over particular edges and provide a basis for these elements on simplicial meshes. As for Whitney edge elements of degree one, the basis is...

The computation of nonlinear quasistationary two-dimensional magnetic fields leads to a nonlinear second order parabolic-elliptic initial-boundary value problem. Such a problem with a nonhomogeneous Dirichlet boundary condition on a part $\Gamma {}_1$ of the boundary is studied in this paper. The problem is discretized in space by the finite element method with linear functions on triangular elements and in time by the implicit-explicit method (the left-hand side by the implicit Euler method and the right-hand...

The present work aims at proposing a rigorous analysis of the mathematical and numerical modelling of ultrasonic piezoelectric sensors. This includes the well-posedness of the final model, the rigorous justification of the underlying approximation and the design and analysis of numerical methods. More precisely, we first justify mathematically the classical quasi-static approximation that reduces the electric unknowns to a scalar electric potential. We next justify the reduction of the computation...

The present work aims at proposing a rigorous analysis of the mathematical and numerical modelling of ultrasonic piezoelectric sensors. This includes the well-posedness of the final model, the rigorous justification of the underlying approximation and the design and analysis of numerical methods. More precisely, we first justify mathematically the classical quasi-static approximation that reduces the electric unknowns to a scalar electric potential. We next justify the reduction of the computation...

The present work aims at proposing a rigorous analysis of the mathematical and numerical modelling of ultrasonic piezoelectric sensors. This includes the well-posedness of the final model, the rigorous justification of the underlying approximation and the design and analysis of numerical methods. More precisely, we first justify mathematically the classical quasi-static approximation that reduces the electric unknowns to a scalar electric potential. We next justify the reduction of the computation...

We present and analyze an interior penalty method for the numerical discretization of the indefinite time-harmonic Maxwell equations in mixed form. The method is based on the mixed discretization of the curl-curl operator developed in [Houston et al., J. Sci. Comp. 22 (2005) 325–356] and can be understood as a non-stabilized variant of the approach proposed in [Perugia et al., Comput. Methods Appl. Mech. Engrg. 191 (2002) 4675–4697]. We show the well-posedness of this approach and derive optimal...

We present and analyze an interior penalty method for the numerical discretization of the indefinite time-harmonic Maxwell equations in mixed form. The method is based on the mixed discretization of the curl-curl operator developed in [Houston et al., J. Sci. Comp.22 (2005) 325–356] and can be understood as a non-stabilized variant of the approach proposed in [Perugia et al., Comput. Methods Appl. Mech. Engrg.191 (2002) 4675–4697]. We show the well-posedness of this approach and derive optimal...

We compute numerically the minimizers of the Dirichlet energy$$E\left(u\right)=\frac{1}{2}{\int }_{B^2}{\left|\nabla u\right|}^2\mathrm{d}x$$among maps $u:B^2\to S^2$ from the unit disc into the unit sphere that satisfy a boundary condition and a degree condition. We use a Sobolev gradient algorithm for the minimization and we prove that its continuous version preserves the degree. For the discretization of the problem we use continuous $P_1$ finite elements. We propose an original mesh-refining strategy needed to preserve the degree with the discrete version of the algorithm (which is a preconditioned...