### A combined method for computing the field of a point source in a surface waveguide.

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

In this paper, we study a model for the magnetization in thin ferromagnetic films. It comes as a variational problem for ${S}^{1}$-valued maps ${m}^{\text{'}}$ (the magnetization) of two variables ${x}^{\text{'}}$: ${E}_{\epsilon}\left({m}^{\text{'}}\right)=\epsilon \int {|{\nabla}^{\text{'}}\xb7{m}^{\text{'}}|}^{2}d{x}^{\text{'}}+\frac{1}{2}\int {\left||{\nabla}^{\text{'}}{|}^{-1/2}{\nabla}^{\text{'}}\xb7{m}^{\text{'}}\right|}^{2}d{x}^{\text{'}}$. We are interested in the behavior of minimizers as $\epsilon \to 0$. They are expected to be ${S}^{1}$-valued maps ${m}^{\text{'}}$ of vanishing distributional divergence ${\nabla}^{\text{'}}\xb7{m}^{\text{'}}=0$, so that appropriate boundary conditions enforce line discontinuities. For finite $\epsilon >0$, these line discontinuities are approximated by smooth transition layers, the so-called Néel walls. Néel...

A continuous finite element method to approximate Friedrichs' systems is proposed and analyzed. Stability is achieved by penalizing the jumps across mesh interfaces of the normal derivative of some components of the discrete solution. The convergence analysis leads to optimal convergence rates in the graph norm and suboptimal of order ½ convergence rates in the L2-norm. A variant of the method specialized to Friedrichs' systems associated with elliptic PDE's in mixed form and reducing the number...

In this paper, a new hybrid simulated annealing algorithm for constrained global optimization is proposed. We have developed a stochastic algorithm called ASAPSPSA that uses Adaptive Simulated Annealing algorithm (ASA). ASA is a series of modifications to the basic simulated annealing algorithm (SA) that gives the region containing the global solution of an objective function. In addition, Simultaneous Perturbation Stochastic Approximation (SPSA)...

We study in this paper the electromagnetic field generated in a conductor by an alternating current density. The resulting interface problem (see Bossavit (1993)) between the metal and the dielectric medium is treated by a mixed–FEM and BEM coupling method. We prove that our BEM-FEM formulation is well posed and that it leads to a convergent Galerkin method.

In this paper we construct a new H(div)-conforming projection-based p-interpolation operator that assumes only Hr(K) $\cap $$\tilde{\mathbf{H}}$-1/2(div, K)-regularity (r > 0) on the reference element (either triangle or square) K. We show that this operator is stable...

In this paper we construct a new H(div)-conforming projection-based p-interpolation operator that assumes only Hr(K) $\cap $$\tilde{\mathbf{H}}$-1/2(div, K)-regularity (r > 0) on the reference element (either triangle or square) K. We show that this operator is stable with...

We use the work of Milton, Seppecher, and Bouchitté on variational principles for waves in lossy media to formulate a finite element method for solving the complex Helmholtz equation that is based entirely on minimization. In particular, this method results in a finite element matrix that is symmetric positive-definite and therefore simple iterative descent methods and preconditioning can be used to solve the resulting system of equations. We also derive an error bound for the method and illustrate...

The paper deals with the application of a non-conforming domain decomposition method to the problem of the computation of induced currents in electric engines with moving conductors. The eddy currents model is considered as a quasi-static approximation of Maxwell equations and we study its two-dimensional formulation with either the modified magnetic vector potential or the magnetic field as primary variable. Two discretizations are proposed, the first one based on curved finite elements and the...

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.

Isogeometric analysis has been developed recently to use basis functions resulting from the CAO description of the computational domain for the finite element spaces. The goal of this study is to develop an axisymmetric Finite Element PIC code in which specific spline Finite Elements are used to solve the Maxwell equations and the same spline functions serve as shape function for the particles. The computational domain itself is defined using splines...

The aim of this paper is to analyze a formulation of the eddy current problem in terms of a time-primitive of the electric field in a bounded domain with input current intensities or voltage drops as source data. To this end, we introduce a Lagrange multiplier to impose the divergence-free condition in the dielectric domain. Thus, we obtain a time-dependent weak mixed formulation leading to a degenerate parabolic problem which we prove is well-posed. We propose a finite element method for space...

We present the current Reduced Basis framework for the efficient numerical approximation of parametrized steady Navier–Stokes equations. We have extended the existing setting developed in the last decade (see e.g. [S. Deparis, SIAM J. Numer. Anal. 46 (2008) 2039–2067; A. Quarteroni and G. Rozza, Numer. Methods Partial Differ. Equ. 23 (2007) 923–948; K. Veroy and A.T. Patera, Int. J. Numer. Methods Fluids 47 (2005) 773–788]) to more general affine and nonaffine parametrizations (such as volume-based...