Consider the scattering of a time-harmonic acoustic incident wave by a bounded, penetrable, and isotropic elastic solid, which is immersed in a homogeneous compressible air or fluid. The paper concerns the numerical solution for such an acoustic-elastic interaction problem in three dimensions. An exact transparent boundary condition (TBC) is developed to reduce the problem equivalently into a boundary value problem in a bounded domain. The perfectly matched layer (PML) technique is adopted to truncate the unbounded physical domain into a bounded computational domain. The well-posedness and exponential convergence of the solution are established for the truncated PML problem by using a PML equivalent TBC. An a posteriori error estimate based adaptive finite element method is developed to solve the scattering problem. Numerical experiments are included to demonstrate the competitive behavior of the proposed method.

This paper is devoted to time domain numerical solutions of two-dimensional (2D) material interface problems governed by the transverse magnetic (TM) and transverse electric (TE) Maxwell's equations with discontinuous electromagnetic solutions. Due to the discontinuity in wave solutions across the interface, the usual numerical methods will converge slowly or even fail to converge. This calls for the development of advanced interface treatments for popular Maxwell solvers. We will investigate such interface treatments by considering two typical Maxwell solvers – one based on collocation formulation and another based on Galerkin formulation. To restore the accuracy reduction of the collocation finite-difference time-domain (FDTD) algorithm near an interface, the physical jump conditions relating discontinuous wave solutions on both sides of the interface must be rigorously enforced. For this purpose, a novel matched interface and boundary (MIB) scheme is proposed in this work, in which new jump conditions are derived so that the discontinuous and staggered features of electric and magnetic field components can be accommodated. The resulting MIB time-domain (MIBTD) scheme satisfies the jump conditions locally and suppresses the staircase approximation errors completely over the Yee lattices. In the discontinuous Galerkin time-domain (DGTD) algorithm – a popular GalerkinMaxwell solver, a proper numerical flux can be designed to accurately capture the jumps in the electromagnetic waves across the interface and automatically preserves the discontinuity in the explicit time integration. The DGTD solution to Maxwell interface problems is explored in this work, by considering a nodal based high order discontinuous Galerkin method. In benchmark TM and TE tests with analytical solutions, both MIBTD and DGTD schemes achieve the second order of accuracy in solving circular interfaces. In comparison, the numerical convergence of the MIBTD method is slightly more uniform, while the DGTD method is more flexible and robust.

In this work, we investigate wave propagation through a zero index meta-material (ZIM) waveguide embedded with triangular dielectric defects. We provide a theoretical guidance on how to achieve total reflection and total transmission (i.e., cloaking) by adjusting the defect sizes and/or permittivities of the defects. Our work provides a systematical way in manipulating wave propagation through ZIM in addition to the widely studied dielectric defects with cylindrical and rectangular geometries.

It is well known that conventional edge elements in solving vector Maxwell's eigenvalue equations by the finite element method will lead to the presence of spurious zero eigenvalues. This problem has been addressed for the first order edge element by Kikuchi by the mixed element method. Inspired by this approach, this paper describes a higher order mixed spectral element method (mixed SEM) for the computation of two-dimensional vector eigenvalue problem of Maxwell's equations. It utilizes Gauss-Lobatto-Legendre (GLL) polynomials as the basis functions in the finite-element framework with a weak divergence condition. It is shown that this method can suppress all spurious zero and nonzero modes and has spectral accuracy. A rigorous analysis of the convergence of the mixed SEM is presented, based on the higher order edge element interpolation error estimates, which fully confirms the robustness of our method. Numerical results are given for homogeneous, inhomogeneous, L-shape, coaxial and dual-inner-conductor cavities to verify the merits of the proposed method.

In this paper we present an efficient algorithm for the calculation of photonic crystal band structures and band structures of photonic crystal waveguides. Our method relies on the fact that the dispersion curves of the band structure are smooth functions of the quasi-momentum in the one-dimensional Brillouin zone. We show the derivation and computation of the group velocity, the group velocity dispersion, and any higher derivative of the dispersion curves. These derivatives are then employed in a Taylor expansion of the dispersion curves. We control the error of the Taylor expansion with the help of a residual estimate and introduce an adaptive scheme for the selection of nodes in the one-dimensional Brillouin zone at which we solve the underlying eigenvalue problem and compute the derivatives of the dispersion curves. These derivatives are then employed in a Taylor expansion of the dispersion curves. We control the error of the Taylor expansion with the help of a residual estimate and introduce an adaptive scheme for the selection of nodes in the one-dimensional Brillouin zone at which we solve the underlying eigenvalue problem and compute the derivatives of the dispersion curves. The proposed algorithm is not only advantageous as it decreases the computational effort to compute the band structure but also because it allows for the identification of crossings and anti-crossings of dispersion curves, respectively. This identification is not possible with the standard approach of solving the underlying eigenvalue problem at a discrete set of values of the quasi-momentum without taking the mode parity into account.

In this paper, the time-dependent Maxwell’s equations used to modeling wave propagation in dispersive lossy bi-isotropic media are investigated. Existence and uniqueness of the modeling equations are proved. Two fully discrete finite element schemes are proposed, and their practical implementation and stability are discussed.

Electromagnetic modeling provides an interesting context to present a link between physical phenomena and homology and cohomology theories. Over the past twenty-five years, a considerable effort has been invested by the computational electromagnetics community to develop fast and general techniques for defining potentials. When magneto-quasi-static discrete formulations based on magnetic scalar potential are employed in problemswhich involve conductive regionswith holes, cuts are needed to make the boundary value problem well defined. While an intimate connection with homology theory has been quickly recognized, heuristic definitions of cuts are surprisingly still dominant in the literature.

The aim of this paper is first to survey several definitions of cuts together with their shortcomings. Then, cuts are defined as generators of the first cohomology group over integers of a finite CW-complex. This provably general definition has also the virtue of providing an automatic, general and efficient algorithm for the computation of cuts. Some counter-examples show that heuristic definitions of cuts should be abandoned. The use of cohomology theory is not an option but the invaluable tool expressly needed to solve this problem.

In this paper we present a rigorous derivation of the material parameters for both the cylinder and rectangle cloaking structures. Numerical results using these material parameters are presented to demonstrate the cloaking effect.

In this paper, we propose several solar cell designs based on metamaterials. Extensive numerical simulations of various designs with different materials are carried out. Our tests show that metamaterial solar cells are quite efficient, and over 80% and 90% absorption rates can be attained for solar spectrum and visible rays, respectively.

A computational study on the enhancement of the second harmonic generation (SHG) in one-dimensional (1D) photonic crystals is presented. The mathematical model is derived from a nonlinear system of Maxwell’s equations, which partly overcomes the shortcoming of some existing models based on the undepleted pump approximation. We designed an iterative scheme coupled with the finite element method which can be applied to simulate the SHG in one dimensional nonlinear photonic band gap structures in our previous work. For the case that the nonlinearity is strong which is desirable to enhance the conversion efficiency, a continuation method is introduced to ensure the convergence of the iterative procedure. The convergence of our method is fast. Numerical experiments also indicate the conversion efficiency of SHG can be significantly enhanced when the frequencies of the fundamental and the second harmonic wave are tuned at the photonic band edges. The maximum total conversion efficiency available reaches more than 50% in all the cases studied.

We develop and analyze an adaptive hybridized Interior Penalty Discontinuous Galerkin (IPDG-H) method for H(curl)-elliptic boundary value problems in 2D or 3D arising from a semi-discretization of the eddy currents equations. The method can be derived from a mixed formulation of the given boundary value problem and involves a Lagrange multiplier that is an approximation of the tangential traces of the primal variable on the interfaces of the underlying triangulation of the computational domain. It is shown that the IPDG-H technique can be equivalently formulated and thus implemented as a mortar method. The mesh adaptation is based on a residual-type a posteriori error estimator consisting of element and face residuals. Within a unified framework for adaptive finite element methods, we prove the reliability of the estimator up to a consistency error. The performance of the adaptive symmetric IPDG-H method is documented by numerical results for representative test examples in 2D.