Abstract
In this paper, we propose a new time-domain finite element method for solving the time dependent Maxwell’s equations coupled with the Lorentz metamaterial model. The Lorentz metamaterial Maxwell’s equations are much more complicated than the standard Maxwell’s equations in free space. Our fully discrete scheme uses edge elements to approximate the unknowns in space, and uses the leap-frog scheme in time discretization. Numerical stability and the optimal error estimate in the \(L^2\) norm are proved for our proposed scheme. Extensive numerical results are presented to confirm the theoretical analysis and applications of our scheme to model many interesting phenomena happened when wave propagates in the Lorentz metamaterials. Examples include the convergence effect happened in the concave lenses formed by the negative refraction index metamatrials, and total reflection and total transmission observed in the zero index metamaterials.
Similar content being viewed by others
References
Ainsworth, M., Coyle, J.: Hierarchic hp-edge element families for Maxwell’s equations on hybrid quadrilateral/triangular meshes. Comput. Methods Appl. Mech. Eng. 190, 6709–6733 (2001)
Banks, H.T., Bokil, V.A., Gibson, N.L.: Analysis of stability and dispersion in a finite element method for Debye and Lorentz media. Numer. Methods Partial Differ. Equ. 25, 885–917 (2009)
Bao, G., Li, P., Wu, H.: An adaptive edge element method with perfectly matched absorbing layers for wave scattering by biperiodic structures. Math. Comput. 79, 1–34 (2010)
Beck, R., Hiptmair, R., Hoppe, H.W., Wohlmuth, B.: Residual based a posteriori error estimators for eddy current computation. Math. Model. Numer. Anal. 34, 159–182 (2000)
Boffi, D., Fernandez, P., Perugia, I.: Computational models of electromagnetic resonators: analysis of edge element approximation. SIAM J. Numer. Anal. 36, 1264–1290 (1999)
Brenner, S.C., Gedicke, J., Sung, L.-Y.: Hodge decomposition for two-dimensional time-harmonic Maxwell’s equations: impedance boundary condition. Math. Methods Appl. Sci. (2015). doi:10.1002/mma.3398
Chen, Z., Du, Q., Zou, J.: Finite element methods with matching and nonmatching meshes for Maxwell equations with discontinuous coefficients. SIAM J. Numer. Anal. 37, 1542–1570 (2000)
Chung, E.T., Ciarlet Jr, P., Yu, T.F.: Convergence and superconvergence of staggered discontinuous Galerkin methods for the three-dimensional Maxwells equations on Cartesian grids. J. Comput. Phys. 235, 14–31 (2013)
Cockburn, B., Li, F., Shu, C.W.: Locally divergence-free discontinuous Galerkin methods for the Maxwell equations. J. Comput. Phys. 194, 588–610 (2004)
Creusé, E., Nicaise, S., Tang, Z., Le Menach, Y., Nemitz, N., Piriou, F.: Residual-based a posteriori estimators for the \(A-\varphi \) magnetodynamic harmonic formulation of the Maxwell system. Math. Models Methods Appl. Sci. 22, 1150028 (2012)
Cui, T.J., Smith, D.R., Liu, R. (eds.): Metamaterials: Theory, Design, and Applications. Springer, Berlin (2009)
Demkowicz, L., Kurtz, J., Pardo, D., Paszynski, M., Rachowicz, W., Zdunek, A.: Computing with hp finite elements. II. Frontiers: three-dimensional Elliptic and Maxwell problems with applications. Chapman & Hall/CRC, Boca Raton (2007)
Fernandes, P., Raffetto, M.: Well posedness and finite element approximability of time-harmonic electromagnetic boundary value problems involving bianisotropic materials and metamaterials. Math. Models Methods Appl. Sci. 19, 2299–2335 (2009)
Fu, Y., Xu, L., Hang, Z.H., Chen, H.: Unidirectional transmission using array of zero-refractive-index metamaterials. Appl. Phys. Lett. 104(19), 193509 (2014)
Gopalakrishnan, J., Pasciak, J.E., Demkowicz, L.F.: Analysis of a multigrid algorithm for time harmonic Maxwell equations. SIAM J. Numer. Anal. 42, 90–108 (2004)
Grote, M.J., Schneebeli, A., Schötzau, D.: Interior penalty discontinuous Galerkin method for Maxwell’s equations: energy norm error estimates. J. Comput. Appl. Math. 204, 375–386 (2007)
Hao, J., Yan, W., Qiu, M.: Super-reflection and cloaking based on zero index metamaterial. Appl. Phys. Lett. 96(10), 101109 (2010)
Hao, Y., Mittra, R.: FDTD Modeling of Metamaterials: Theory and Applications. Artech House, Boston (2008)
Hesthaven, J.S., Warburton, T.: Nodal high-order methods on unstructured grids: I. Time-domain solution of Maxwell’s equations. J. Comput. Phys. 181, 186–221 (2002)
Hiptmair, R.: Finite elements in computational electromagnetism. Acta Numer. 11, 237–339 (2002)
Houston, P., Perugia, I., Schötzau, D.: Mixed discontinuous Galerkin approximation of the Maxwell operator: non-stabilized formulation. J. Sci. Comput. 22, 315–346 (2005)
Huang, Y., Li, J., Yang, W.: Modeling backward wave propagation in metamaterials by the finite element time-domain method. SIAM J. Sci. Comput. 35, B248–B274 (2013)
Huang, Y., Li, J.: Total reflection and cloaking by triangular defects embedded in zero index metamaterials. Adv. Appl. Math. Mech. 7, 135–144 (2015)
Lanteri, S., Scheid, C.: Convergence of a discontinuous Galerkin scheme for the mixed time-domain Maxwell’s equations in dispersive media. IMA J. Numer. Anal. 33, 432–459 (2013)
Li, J.: Finite element study of the Lorentz model in metamaterials. Comput. Methods Appl. Mech. Eng. 200, 626–637 (2011)
Li, J.: Numerical convergence and physical fidelity analysis for Maxwell’s equations in metamaterials. Comput. Methods Appl. Mech. Eng. 198, 3161–3172 (2009)
Li, J., Huang, Y.: Time-Domain Finite Element Methods for Maxwel’s Equations in Metamaterials, vol. 43. Springer, New York (2013)
Li, J., Wood, A.: Finite element analysis for wave propagation in double negative metamaterials. J. Sci. Comput. 32, 263–286 (2007)
Lu, T., Zhang, P., Cai, W.: Discontinuous Galerkin methods for dispersive and lossy Maxwell’s equations and PML boundary conditions. J. Comput. Phys. 200, 549–580 (2004)
Markos, P., Soukoulis, C.M.: Wave Propagation: From Electrons to Photonic Crystals and Left-Handed Materials. Princeton University Press, Princeton (2008)
Monk, P.: Finite Element Methods for Maxwell’s Equations. Oxford University Press, Oxford (2003)
Monk, P., Süli, E.: A convergence analysis of Yee’s scheme on nonuniform grids. SIAM J. Numer. Anal. 32, 393–412 (1994)
Mu, L., Wang, J., Ye, X., Zhang, S.: A weak Galerkin finite element method for the Maxwell equations. J. Sci. Comput. 65, 363–386 (2015)
Nguyen, V.C., Chen, L., Halterman, K.: Total transmission and total reflection by zero index metamaterials with defects. Phys. Rev. Lett. 105(23), 233908 (2010)
Sármány, D., Botchev, M.A., van der Vegt, J.J.W.: Dispersion and dissipation error in high-order Runge–Kutta discontinuous Galerkin discretisations of the Maxwell equations. J. Sci. Comput. 33, 47–74 (2007)
Shi, D., Yao, C.: Nonconforming finite element approximation of time-dependent Maxwell’s equations in Debye medium. Numer. Methods Partial Differ. Equ. 30, 1654–1673 (2014)
Wang, B., Xie, Z., Zhang, Z.: Error analysis of a discontinuous Galerkin method for Maxwell equations in dispersive media. J. Comput. Phys. 229, 8552–8563 (2010)
Wu, Y., Li, J.: Total reflection and cloaking by zero index metamaterials loaded with rectangular dielectric defects. Appl. Phys. Lett. 102(18), 183105 (2013)
Yang, Z., Wang, L.-L.: Accurate simulation of circular and elliptic cylindrical invisibility cloaks. Commun. Comput. Phys. 17, 822–849 (2015)
Zhong, L., Chen, L., Shu, S., Wittum, G., Xu, J.: Convergence and optimality of adaptive edge finite element methods for time-harmonic Maxwell equations. Math. Comput. 81, 623–642 (2012)
Author information
Authors and Affiliations
Corresponding author
Additional information
This work was supported by NSFC Projects 11401506 and 11271310, NSFC Key Project 91430213, Hunan Education Department Project (15B236), and NSF Grant DMS-1416742.
Rights and permissions
About this article
Cite this article
Yang, W., Huang, Y. & Li, J. Developing a Time-Domain Finite Element Method for the Lorentz Metamaterial Model and Applications. J Sci Comput 68, 438–463 (2016). https://doi.org/10.1007/s10915-015-0144-y
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10915-015-0144-y
Keywords
- Maxwell’s equations
- Finite element method
- Lorentz metamaterial model
- Edge elements
- Total transmission
- Total reflection