Analysis of the energy‐conserved S‐FDTD scheme for variable coefficient Maxwell's equations in disk domains

In this paper, we analyze the energy‐conserved splitting finite‐difference time‐domain (FDTD) scheme for variable coefficient Maxwell's equations in two‐dimensional disk domains. The approach is energy‐conserved, unconditionally stable, and effective. We strictly prove that the EC‐S‐FDTD scheme for the variable coefficient Maxwell's equations in disk domains is of second order accuracy both in time and space. It is also strictly proved that the scheme is energy‐conserved, and the discrete divergence‐free is of second order convergence. Numerical experiments confirm the theoretical results, and practical test is simulated as well to demonstrate the efficiency of the proposed EC‐S‐FDTD scheme. Copyright © 2015 John Wiley & Sons, Ltd.     

Energy conservation and super convergence analysis of the EC‐S‐FDTD schemes for Maxwell equations with periodic boundaries

This paper is concerned with new energy analysis of the two dimensional Maxwell's equations and the symmetric energy‐conserved splitting finite difference time domain (EC‐S‐FDTD) method with the periodic boundary (PB) condition. New energy identities of the Maxwell's equations in terms of H¹ and H² norms are proposed and interpreted by considering the physical meanings of the H¹ and H² semi‐norms in the identities. It is found from these new identities that the first and second curls of the electromagnetic fields are conserved in terms their magnitudes. By the energy methods, the numerical energy identities of the symmetric EC‐S‐FDTD method are derived and shown to converge to the continuous energy identities of the Maxwell's equations. This proves that the symmetric EC‐S‐FDTD scheme is unconditionally stable and energy conserved in the discrete H¹ and H² norms. Also by the energy methods, it is proved that the symmetric EC‐S‐FDTD method with PB condition is of second order (super) convergence in the discrete H¹ and H² norms. Numerical experiments are carried out and confirm the analysis on energy conservation, stability and super convergence.     

A numerical study of electromagnetic waves in periodic waveguides

Here are considered time‐harmonic electromagnetic waves in a quadratic waveguide consisting of a periodic dielectric core enclosed by conducting walls. The permittivity function may be smooth or have jumps. The electromagnetic field is given by a magnetic vector potential in Lorenz gauge, and defined on a Floquet cell. The Helmholtz operator is approximated by a Chebyshev collocation, Fourier–Galerkin method. Laurent's rule and the inverse rule are employed for the representation of Fourier coefficients of products of functions. The computations yield, for known wavenumbers, values of the first few eigenfrequencies of the field. In general, the dispersion curves exhibit band gaps. Field patterns are identified as transverse electric, TE, transverse magnetic, TM, or hybrid modes. Maxwell's equations are fulfilled. A few trivial solutions appear when the permittivity varies in the guiding direction and across it. The results of the present method are consistent with exact results and with those obtained by a low‐order finite element software. The present method is more efficient than the low‐order finite element method. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 490–513, 2014     

On the Diffraction by Biperiodic Anisotropic Structures

This article studies the scattering of electromagnetic waves by a nonmagnetic biperiodic structure. The structure consists of anisotropic optical materials and separates two regions with constant dielectric coefficients. The time harmonic Maxwell equations are transformed to an equivalent strongly elliptic variational problem for the magnetic field in a bounded biperiodic cell with nonlocal boundary conditions. This guarantees the existence of quasiperiodic magnetic fields in H 1 and electric fields in H (curl) solving Maxwell's equations. The uniqueness is proved for all frequencies excluding possibly a discrete set. The analytic dependence of these solutions on frequency and incident angles is studied.     

Maxwell's equations in an unbounded structure

This paper is concerned with the mathematical analysis of the electromagnetic wave scattering by an unbounded dielectric medium, which is mounted on a perfectly conducting infinite plane. By introducing a transparent boundary condition on a plane surface confining the medium, the scattering problem is modeled as a boundary value problem of Maxwell's equations. Based on a variational formulation, the problem is shown to have a unique weak solution for a wide class of dielectric permittivity and magnetic permeability by using the generalized Lax–Milgram theorem. Copyright © 2016 John Wiley & Sons, Ltd.     

L p Estimates for Maxwell's Equations with Source Term

The goal of this paper is to establish interior and global L ^p -type estimates for the solutions of Maxwell's equations with source term in a domain filled with two different materials separated by a ² interface. The usual elliptic estimates cannot be applied directly, due to the singularity of the dielectric permittivity. A special curl-div decomposition is introduced for the electric field to reduce the problem to an elliptic equation in divergence form with jump coefficients. The potential analysis and the jump condition lead to the interior L ^p estimates which are superior to the straightforward Nash-Moser estimates. The reduction procedure is expected to be useful for future numerical simulation. Because of the natural physical requirements, the boundary condition is nonlocal and involves a first order pseudo-differential operator, the boundary estimate is established by delicate new maximum principles and Riesz convexity arguments. These estimates are then employed to solve a nonlinear optics problem that arises in the modeling of surface enhanced second-harmonic generation of nonlinear diffractive optics in periodic structures (gratings).     

An inverse source problem for Maxwell's equations in anisotropic media

In this article, we consider nonstationary Maxwell's equations in an anisotropic medium in the (x ₁,?x ₂,?x ₃)-space, where equations of the divergences of electric and magnetic flux densities are also unknown. Then we discuss an inverse problem of determining the x ₃-independent components of the electric current density from observations on the plane x ₃?=?0 over a time interval. Our main aim is, study conditional stability in the inverse problem provided the permittivity and the permeability are independent of x ₃. The main tool is a new Carleman estimate.     

Dynamic susceptibility computations for thin magnetic films

This paper deals with the dynamic modelling of thin ferromagnetic layers, based on the coupling of Maxwell's equations with the nonlinear Landau–Lifschitz–Gilbert law. A 2-D micromagnetic model is described which involves a FDTD code to determine equilibrium configurations and a finite element method to compute magnetostatic fields. Finally, after linearization, the susceptibility spectra of films supporting a weak-stripe-domain structure are computed and successfully compared to existing measurements without introducing any fitting parameter.     

Dissipative Maxwell's equations at low frequencies

A study of the low-frequency behaviour of solutions to dissipative Maxwell's equations with discontinuous coefficients is given.     

A hybrid technique for computing the power distribution generated in a lossy medium during microwave heating

Over the years researchers in the field of computational electromagnetics (CEM) have investigated and explored a number of different techniques to resolve electromagnetic fields inside waveguide and cavity structures. The equations that govern the fundamental behaviour of electromagnetic wave propagation in such structures are Maxwell's equations. In the literature, a number of different techniques have been employed to solve these equations and out of these methods, the classical finite-difference time-domain (FD-TD) scheme, which uses a staggered time and space discretisation, is the most well-known and widely used. However, this scheme is complicated to implement on an irregular computational domain using unstructured meshes.This research work builds upon previous work undertaken for a waveguide, where a coupled method was introduced for the solution of the governing electromagnetic equations. In that work, the free-space component of the solution was computed in the time-domain, whilst the power distribution in the load was resolved using the frequency dependent electric field Helmholtz equation. This methodology resulted in a time-frequency domain hybrid scheme. In this paper, the hybrid method has been tested further for both waveguide and cavity configurations that are loaded with a lossy dielectric material. Numerical tests highlight both the accuracy and computational efficiency of the proposed hybrid strategy for predicting the power distribution generated during microwave heating processes. The accuracy of the hybrid scheme is gauged by direct comparison with the FD-TD numerical solutions and previously published thermal images.     

A LOW-FREQUENCY ELECTROMAGNETIC NEAR-FIELD INVERSE PROBLEM FOR A SPHERICAL SCATTERER

The interior low-frequency electromagnetic dipole excitation of a dielectric sphere is uti- lized as a simplified but realistic model in various biomedical applications. Motivated by these considerations, in this paper, we investigate analytically a near-field inverse scatter- ing problem for the electromagnetic interior dipole excitation of a dielectric sphere. First, we obtain, under the low-frequency assumption, a closed-form approximation of the series of the secondary electric field at the dipole＇s location. Then, we utilize this derived approx- imation in the development of a simple inverse medium scattering algorithm determining the sphere＇s dielectric permittivity. Finally, we present numerical results for a human head model, which demonstrate the accurate determination of the complex permittivity by the developed algorithm.     

The a Posteriori Error Estimator of SDG Method for Variable Coefficients Time-Harmonic Maxwell's Equations

In this paper, we study the a posteriori error estimator of SDG method for variable coefficients time-harmonic Maxwell's equations. We propose two a posteriori error estimators, one is the recovery-type estimator, and the other is the residual-type estimator. We first propose the curl-recovery method for the staggered discontinuous Galerkin method (SDGM), and based on the super-convergence result of the postprocessed solution, an asymptotically exact error estimator is constructed. The residual-type a posteriori error estimator is also proposed, and it's reliability and effectiveness are proved for variable coefficients time-harmonic Maxwell's equations. The efficiency and robustness of the proposed estimators is demonstrated by the numerical experiments.     

Fully discrete T‐ψ finite element method to solve a nonlinear induction hardening problem

We study an induction hardening model described by Maxwell's equations coupled with a heat equation. The magnetic induction field is assumed a nonlinear constitutional relation and the electric conductivity is temperature‐dependent. The T ‐ψ method is to transform Maxwell's equations to the vector–scalar potential formulations and to solve the potentials by means of the finite element method. In this article, we present a fully discrete T ‐ψ finite element scheme for this nonlinear coupled problem and discuss its solvability. We prove that the discrete solution converges to a weak solution of the continuous problem. Finally, we conclude with two numerical experiments for the coupled system.     

Energy decay of solutions of maxwell's equations in bounded domains

In this paper the long time asymptotic behavior of solutions of Maxwell's equations in a bounded domain is investigated. A possibly nonlinear damping term arising from the electric conductivity is considered, whith vanishes on some part of the domain     

Analysis of stability and dispersion in a finite element method for Debye and Lorentz dispersive media

We study the stability properties of, and the phase error present in, a finite element scheme for Maxwell's equations coupled with a Debye or Lorentz polarization model. In one dimension we consider a second order formulation for the electric field with an ordinary differential equation for the electric polarization added as an auxiliary constraint. The finite element method uses linear finite elements in space for the electric field as well as the electric polarization, and a theta scheme for the time discretization. Numerical experiments suggest the method is unconditionally stable for both Debye and Lorentz models. We compare the stability and phase error properties of the method presented here with those of finite difference methods that have been analyzed in the literature. We also conduct numerical simulations that verify the stability and dispersion properties of the scheme. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Comparison of Taylor finite difference and window finite difference and their application in FDTD

The finite difference time domain (FDTD) method is an important tool in numerical electromagnetic simulation. There are many ways to construct a finite difference approximation such as the Taylor series expansion theorem, the filtering theory, etc. This paper aims to provide the comparison between the Taylor finite difference (TFD) scheme based on the Taylor series expansion theorem and the window finite difference (WFD) scheme based on the filtering theory. Their properties have been examined in detail, separately. In addition, the formula of the generalized finite difference (GFD) scheme is presented, which can include both the TFD scheme and the WFD scheme. Furthermore, their application in the numerical solution of Maxwell's equations is presented. The formulas for the stability criterion and the numerical dispersion relation are derived and analyzed. In order to evaluate their performance more accurately, a new definition of error is presented. Upon it, the effect of several factors including the grid resolution, the Courant number and the aspect ratio of the cell on the performance of the numerical dispersion is examined.     

Numerical solution of an inverse initial boundary-value problem for the full time-dependent Maxwell's equations in the presence of imperfections of small volume

We consider the numerical solution, in a three-dimensional bounded domain, of the inverse problem for identifying the location of small electromagnetic imperfections in a medium with homogeneous background. Our numerical algorithm is based on the coupling of a discontinuous Galerkin method for the time-dependent Maxwell's equations, on the exact controllability method and on a Fourier inversion. Several numerical results are given with one and two imperfections and the robustness and accuracy of the numerical method used for the dynamic detection problem are shown.     

Weak Galerkin finite element methods for electric interface model with nonhomogeneous jump conditions

In this paper, the weak Galerkin finite element method (WG-FEM) is applied to a pulsed electric model arising in biological tissue when a biological cell is exposed to an electric field. A fitted WG-FEM is proposed to approximate the voltage of the pulsed electric model across the physical media involving an electric interface (surface membrane), and heterogeneous permittivity and a heterogeneous conductivity. This method uses totally discontinuous functions in approximation space and allows the usage of finite element partitions consisting of general polygonal meshes. Optimal pointwise-in-time error estimates in L²-norm and H¹-norm are shown to hold for the semidiscrete scheme even if the regularity of the solution is low on the whole domain. Furthermore, a fully discrete approximation based on backward Euler scheme is analyzed and related optimal error estimates are derived.     

On the one- and one-half-dimensional relativistic Vlasov-Fokker-Planck-Maxwell system

In this paper, we study the relativistic Vlasov-Fokker-Planck-Maxwell system in one space variable and two momentum variables. This non-linear system of equations consists of a transport equation for the phase space distribution function combined with Maxwell's equations for the electric and magnetic fields. It is important in modelling distribution of charged particles in the kinetic theory of plasma. We prove the existence of a classical solution when the initial density decays fast enough with respect to the momentum variables. The solution which shares this same decay condition along with its first derivatives in the momentum variables is unique.