The aim of this paper is to present how to make a dedicaded computed language polymorphic and multi type, in C++to solve partial differential equations with the finite element method. The driving idea is to make the language as close as possible to the mathematical notation.
The aim of this paper is to present how to make a dedicaded computed language polymorphic and multi type, in C++ to solve partial differential equations with the finite element method. The driving idea is to make the language as close as possible to the mathematical notation.
This paper is devoted to the study of cloaking via anomalous localized resonance (CALR) in the two- and three-dimensional quasistatic regimes. CALR associated with negative index materials was discovered by Milton and Nicorovici [21] for constant plasmonic structures in the two-dimensional quasistatic regime. Two key features of this phenomenon are the localized resonance, i.e., the fields blow up in some regions and remain bounded in some others, and the connection between the localized resonance...
This paper is devoted to the practical computation of the magnetic potential induced by a distribution of magnetization in the theory of micromagnetics. The problem turns out to be a coupling of an interior and an exterior problem. The aim of this work is to describe a complete method that mixes the approaches of Ying [12] and Goldstein [6] which consists in constructing a mesh for the exterior domain composed of homothetic layers. It has the advantage of being well suited for catching the decay...
This paper is devoted to the practical computation of the magnetic potential induced by a distribution of magnetization in the theory of micromagnetics. The problem turns out to be a coupling of an interior and an exterior problem. The aim of this work is to describe a complete method that mixes the approaches of Ying [12] and Goldstein [6] which consists in constructing a mesh for the exterior domain composed of homothetic layers. It has the advantage of being well suited for catching the...
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....
We present a finite element method to compute guided modes in a stratified medium. The major difficulty to overcome is related to the unboundedness of the stratified medium. Our method is an alternative to the use of artificial boundary conditions and to the use of integral representation formulae. The domain is bounded in such a way we can write the solution on its lateral boundaries in terms of Fourier series. The series is then truncated for the computations over the bounded domain. The problem...
We present a finite element method to compute guided modes in a stratified medium. The major difficulty to overcome is related to the unboundedness of the stratified medium. Our method is an alternative to the use of artificial boundary conditions and to the use of integral representation formulae. The domain is bounded in such a way we can write the solution on its lateral boundaries in terms of Fourier series. The series is then truncated for the computations over the bounded domain. The problem...
The paper is devoted to an integral equation algorithm for studying the scattering of plane waves by multilayer diffraction gratings under oblique incidence. The scattering problem is described by a system of Helmholtz equations with piecewise constant coefficients in coupled by special transmission conditions at the interfaces between different layers. Boundary integral methods lead to a system of singular integral equations, containing at least two equations for each interface. To deal with...
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.