Analysis of FDTD to UPML for Maxwell equations in polar coordinates

An FDTD system associated with uniaxial perfectly matched layer(UPML) for an electromagnetic scattering problem in two-dimensional space in polar coordinates is considered.Particularly the FDTD system of an initial-boundary value problems of the transverse magnetic(TM) mode to Maxwell’s equations is obtained by Yee’s algorithm,and the open domain of the scattering problem is truncated by a circle with a UPML.Besides,an artificial boundary condition is imposed on the outer boundary of the UPML.Afterwards,stability of the FDTD system on the truncated domain is established through energy estimates by the Gronwall inequality.Numerical experiments are designed to approve the theoretical analysis.     

An overfilled cavity problem for Maxwell's equations

This paper is concerned with the mathematical analysis of the scattering of a time‐harmonic electromagnetic plane wave by an open and overfilled cavity that is embedded in a perfect electrically conducting infinite ground plane, where the electromagnetic wave propagation is governed by the Maxwell equations. Above the flat ground surface and the open aperture of the cavity, the space is assumed to be filled with a homogeneous medium with a constant permittivity and permeability, whereas the interior of the cavity is filled with some inhomogeneous medium with a variable permittivity and permeability. The scattering problem is modeled as a boundary value problem over a bounded domain, with transparent boundary condition proposed on the hemisphere enclosing the inhomogeneity represented by the cavity. The existence and uniqueness of the weak solution for the model problem are established by using a variational approach. The perfectly matched layer (PML) method is investigated to truncate the unbounded electromagnetic cavity scattering problem. It is shown that the truncated PML problem attains a unique solution. An explicit error estimate is given between the solution of the original scattering problem and that of the truncated PML problem. The error estimate implies that the PML solution converges exponentially to the original cavity scattering problem by increasing either the PML medium parameter or the PML layer thickness. The convergence result is expected to be useful for determining the PML medium parameter in the computational electromagnetic scattering problem. Copyright © 2012 John Wiley & Sons, Ltd.     

Initial-boundary value problem and integrodifferential equations of electrodynamics for an inhomogeneous conductive body

The initial-boundary value problem of determining the electromagnetic field in an inhomogeneous conducting sample and the surrounding external medium is solved under certain assumptions on the sample, the external medium, and the external current density. The existence of a classical solution to this problem is proved. The electromagnetic field under small variations in the sample’s electric conductivity is computed by applying perturbation theory.     

An initial-boundary value problem for the Maxwell equations

In this paper, by using methods from complex analysis and quaternionic analysis, we investigate an initial-boundary value problem for the Maxwell equations and obtain the general solutions and solvable conditions of the problem respectively in different cases. In addition, by using a similar method, we also discuss an initial-boundary value problem for a hyperbolic complex system of first order equations in R³.     

A convergence analysis for a sweeping preconditioner for block tridiagonal systems of linear equations

We study sweeping preconditioners for symmetric and positive definite block tridiagonal systems of linear equations. The algorithm provides an approximate inverse that can be used directly or in a preconditioned iterative scheme. These algorithms are based on replacing the Schur complements appearing in a block Gaussian elimination direct solve by hierarchical matrix approximations with reduced off‐diagonal ranks. This involves developing low rank hierarchical approximations to inverses. We first provide a convergence analysis for the algorithm for reduced rank hierarchical inverse approximation. These results are then used to prove convergence and preconditioning estimates for the resulting sweeping preconditioner. Copyright © 2014 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 Stability Estimate for a Solution to a Three-Dimensional Inverse Problem for the Maxwell Equations

We consider the problem of determining dielectric permittivity and conductivity in the Maxwell equations. As additional information we prescribe the traces of the tangential components of the electromagnetic field on the lateral surface of a cylindric domain. We establish a stability estimate for a solution to the inverse problem and a uniqueness theorem.     

Analysis of a finite PML approximation for the three dimensional time-harmonic Maxwell and acoustic scattering problems

We consider the approximation of the frequency domain three-dimensional Maxwell scattering problem using a truncated domain perfectly matched layer (PML). We also treat the time-harmonic PML approximation to the acoustic scattering problem. Following work of Lassas and Somersalo in 1998, a transitional layer based on spherical geometry is defined, which results in a constant coefficient problem outside the transition. A truncated (computational) domain is then defined, which covers the transition region. The truncated domain need only have a minimally smooth outer boundary (e.g., Lipschitz continuous). We consider the truncated PML problem which results when a perfectly conducting boundary condition is imposed on the outer boundary of the truncated domain. The existence and uniqueness of solutions to the truncated PML problem will be shown provided that the truncated domain is sufficiently large, e.g., contains a sphere of radius . We also show exponential (in the parameter ) convergence of the truncated PML solution to the solution of the original scattering problem inside the transition layer.

Our results are important in that they are the first to show that the truncated PML problem can be posed on a domain with nonsmooth outer boundary. This allows the use of approximation based on polygonal meshes. In addition, even though the transition coefficients depend on spherical geometry, they can be made arbitrarily smooth and hence the resulting problems are amenable to numerical quadrature. Approximation schemes based on our analysis are the focus of future research.

     

Asymptotic Stability of Equilibrium State to the Mixed Initial-Boundary Value Problem for Quasilinear Hyperbolic Systems

Under the internal dissipative condition, the Cauchy problem for inhomogeneous quasilinear hyperbolic systems with small initial data admits a unique global C1 solution, which exponentially decays to zero as t → +∞, while if the coefficient matrixΘ of boundary conditions satisfies the boundary dissipative condition, the mixed initialboundary value problem with small initial data for quasilinear hyperbolic systems with nonlinear terms of at least second order admits a unique global C1 solution, which also exponentially decays to zero as t → +∞. In this paper, under more general conditions, the authors investigate the combined effect of the internal dissipative condition and the boundary dissipative condition, and prove the global existence and exponential decay of the C1 solution to the mixed initial-boundary value problem for quasilinear hyperbolic systems with small initial data. This stability result is applied to a kind of models, and an example is given to show the possible exponential instability if the corresponding conditions are not satisfied.     

Global existence and asymptotic behaviour for systems of nonlinear hyperbolic equations

We consider the initial-boundary value problem for a class of nonlinear hyperbolic equations system in a bounded domain. Using the potential well theory, the existence of global solutions is investigated. We also established the asymptotic behaviour of global solutions as t?→?+?∞ by applying the multiplier method.     

Propagation of TM waves in a Kerr nonlinear layer

TM electromagnetic waves propagating through a nonlinear homogeneous isotropic unmagnetized dielectric layer located between two homogeneous isotropic half-spaces are studied. The nonlinearity in the layer obeys the Kerr law. The problem is reduced to a system of nonlinear ordinary differential equations. A dispersion relation for the propagation constants is derived. The results are compared with those in the case of a linear layer.     

New energy identities and super convergence analysis of the energy conserved splitting FDTD methods for 3D Maxwell's equations

The energy‐conserved splitting finite‐difference time‐domain (EC‐S‐FDTD) method has recently been proposed to solve the Maxwell equations with second order accuracy while numerically keep the L² energy conservation laws of the equations. In this paper, the EC‐S‐FDTD scheme for the 3D Maxwell equations is proved to be energy‐conserved and unconditionally stable in the discrete H¹ norm. The EC‐S‐FDTD scheme is of second‐order accuracy both in time step and spatial steps, which suggests the super‐convergence of this scheme in the discrete H¹ norm. And the divergence of the electric field of the EC‐S‐FDTD scheme in the discrete L² norm is second‐order accurate. Numerical experiments confirm our theoretical analysis. Copyright © 2012 John Wiley & Sons, Ltd.     

Simple geometric constructions of quadratically and cubically convergent iterative functions to solve nonlinear equations

In this paper, we derive one-parameter families of Newton, Halley, Chebyshev, Chebyshev-Halley type methods, super-Halley, C-methods, osculating circle and ellipse methods respectively for finding simple zeros of nonlinear equations, permitting f ′ (x) = 0 at some points in the vicinity of the required root. Halley, Chebyshev, super-Halley methods and, as an exceptional case, Newton method are seen as the special cases of the family. All the methods of the family and various others are cubically convergent to simple roots except Newton’s or a family of Newton’s method.      

Energy Identities and Stability Analysis of the Yee Scheme for 3D Maxwell Equations

In this paper numerical energy identities of the Yee scheme on uniform grids for three dimensional Maxwell equations with periodic boundary conditions are proposed and expressed in terms of the $L^2$, $H^1$ and $H^2$ norms. The relations between the $H^1$ or $H^2$ semi-norms and the magnitudes of the curls or the second curls of the fields in the Yee scheme are derived. By the $L^2$ form of the identity it is shown that the solution fields of the Yee scheme is approximately energy conserved. By the $H^1$ or $H^2$ semi norm of the identities, it is proved that the curls or the second curls of the solution of the Yee scheme are approximately magnitude (or energy)-conserved. From these numerical energy identities, the Courant-Friedrichs-Lewy (CFL) stability condition is re-derived, and the stability of the Yee scheme in the $L^2$, $H^1$ and $H^2$ norms is then proved. Numerical experiments to compute the numerical energies and convergence orders in the $L^2$, $H^1$ and $H^2$ norms are carried out and the computational results confirm the analysis of the Yee scheme on energy conservation and stability analysis.     

On the solutions of a parametric family of cubic Thue equations

In this paper, we solve a family of Diophantine equations associated with families of number fields of degree 3. In fact, we use Baker’s method find all solutions to the Thue equation

Stability analysis of Runge-Kutta methods for nonlinear neutral delay integro-differential equations

The sufficient conditions for the stability and asymptotic stability of Runge-Kutta methods for nonlinear neutral delay integro-differential equations are derived. A numerical test that confirms the theoretical results is given in the end.     

Continuation of a solution to a homogeneous system of Maxwell equations

We consider the problem on the continuation of a solution to a system of Maxwell equations, using its values on a part of the domain boundary.     

Some classes of inverse evolution problems for parabolic equations

We consider the problem of simultaneously determining coefficients of a second order nonlinear parabolic equation and a solution to this equation. The unknown coefficients occur in the main part and in the nonlinear summand as well. The overdetermination conditions are conditions of the Dirichlet type on a family of planes of arbitrary dimension. It is demonstrated that the problem in question is solvable locally in time in Hölder spaces. When the unknown functions enter the right-hand side and the equation is linear, the theorem of global unique existence (in time) is established.     

On a class of nonlinear integral equations

In this paper we consider abstract equations of the typeK _ν ν +ν =w ₀, in a closed convex subset of a separable Hilbert spaceH. For eachv in the closed convex subset,K _v :H →H is a bounded linear map. As an application of our abstract result we obtain an existence result for nonlinear integral equations of the typeν(s)+ν(s)∫ ₀ ¹ k(s,t)ν(t)dt =W ₀(s) in the spaceL ² [0,1].