Skip to main content Accessibility help
×
Home

A High-Order Time Domain Discontinuous Galerkin Method with Orthogonal Tetrahedral Basis for Electromagnetic Simulations in 3-D Heterogeneous Conductive Media

  • Jun Yang (a1) (a2), Wei Cai (a3) (a2) and Xiaoping Wu (a1)

Abstract

We present a high-order discontinuous Galerkin (DG) method for the time domain Maxwell's equations in three-dimensional heterogeneous media. New hierarchical orthonormal basis functions on unstructured tetrahedral meshes are used for spatial discretization while Runge-Kutta methods for time discretization. A uniaxial perfectly matched layer (UPML) is employed to terminate the computational domain. Exponential convergence with respect to the order of the basis functions is observed and large parallel speedup is obtained for a plane-wave scattering model. The rapid decay of the out-going wave in the UPML is shown in a dipole radiation simulation. Moreover, the low frequency electromagnetic fields excited by a horizontal electric dipole (HED) and a vertical magnetic dipole (VMD) over a layered conductive half-space and a high frequency ground penetrating radar (GPR) detection for an underground structure are investigated, showing the high accuracy and broadband simulation capability of the proposed method.

Copyright

Corresponding author

*Corresponding author. Email addresses: wcai@uncc.edu, wcai@csrc.ac.cn (W. Cai)

References

Hide All
[1] Bohren, C. F. and Huffman, D. R.. Absorption and Scattering of Light by Small Particles, John Wiley & Sons, 2008.
[2] Cai, W.. Computational Methods for Electromagnetic Phenomena: Electrostatics in Solvation, Scattering, and Electron Transport, Cambridge University Press, 2013.
[3] Cockburn, B. and Shu, C.-W.. TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws II: general framework, Math. Comp., 52 (1989), pp. 411435.
[4] Cockburn, B. and Shu, C.-W.. The Runge-Kutta local projection P1-Discontinuous Galerkin finite element method for scalar conservation laws, Mathematical Modelling and Numerical Analysis, 25 (1991), pp. 337361.
[5] Cockburn, B., Hou, S. and Shu, C.-W.. The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws IV: the multidimensional case, Math. Comp., 54 (1990), pp. 545581.
[6] Cockburn, B. and Shu, C.-W.. The Runge-Kutta discontinuous Galerkin method for conservation laws V: multidimensional systems, J. Comput. Phys., 141 (1998), pp. 199224.
[7] Dolean, V., Fahs, H., Fezoui, L., and Lanteri, S.. Locally implicit discontinuous Galerkin method for time domain electromagnetics, J. Comput. Phys., 229 (2010), pp. 512526.
[8] Gedney, S. D.. An anisotropic perfectly matched layer-absorbing medium for the truncation of FDTD lattices, IEEE Trans. Antennas Propagat., 44 (1996), pp. 16301639.
[9] Geuzaine, C. and Remacle, J. F.. Gmsh: a three-dimensional finite element mesh generator with built-in pre- and post-processing facilities, Int. J. Numr. Meth. Engng, 79 (2009), pp. 13091331.
[10] Ji, X., Lu, T., Cai, W. and Zhang, P.. Discontinuous Galerkin time domain (DGTD) methods for the study of 2-D waveguide-coupled microring resonators, Journal of lightwave technology, 23 (2005), pp. 38643874.
[11] Jin, J.-M.. The Finite Element Method in Electromagnetics, John Wiley & Sons, 2014.
[12] Kabakian, A. V., Shankar, V. and Hall, W. F.. Unstructured grid-based discontinuous Galerkin method for broadband electromagnetic simulations, J. Sci. Comput., 20 (2004), pp. 405431.
[13] Keast, P.. Moderate-degree tetrahedral quadrature formulas, Computer Methods in Applied Mechanics and Engineering, 55 (1986), pp. 339348.
[14] Kim, J. and Gopinath, A.. Simulation of a metamaterial containing cubic high dielectric resonators, Phys. Rev. B, 76 (2007), 115126
[15] Kopriva, D.A., Woodruff, S. L. and Hussaini, M. Y.. Computation of electromagnetic scattering with a non-onforming discontinuous spectral element method, Int. J. Numr. Meth. Engng, 53 (2002), pp. 105122.
[16] Li, J., Chen, Y. and Liu, Y.. Mathematical simulation of metamaterial solar cells, Advances in Applied Mathematics and Mechanics, 3 (2011), pp. 702715.
[17] Li, J., Waters, J.W. and Machorro, E. A.. An implicit leap-frog discontinuous Galerkin method for the time-domain Maxwell's equations in metamaterials, Computer Methods in Applied Mechanics and Engineering, 223 (2012), pp. 4354.
[18] Li, J. and Hesthaven, J. S.. Analysis and application of the nodal discontinuous Galerkin method for wave propagation in metamaterials, J. Comput. Phys., 258(2014), pp. 915930.
[19] Lu, T., Cai, W. and Zhang, P.. Discontinuous Galerkin Time-Domain Method for GPR Simulation in Dispersive Media, IEEE Trans Geoscience and Remote Sensing, 43 (2005), pp. 7280.
[20] Lu, T., Zhang, P. and Cai, W.. Discontinuous Galerkin methods for dispersive and lossy Maxwell's equations and PML boundary conditions, J. Comput. Phys., 200 (2004), pp. 549580.
[21] Nédélec, J. C.. Mixed finite elements in , Numerische Mathematik, 35 (1980), pp. 315341.
[22] Reed, W. H. and Hill, T. R.. Triangular mesh methods for the neutron transport equation, Technical Report, LA-UR-73-479, Los Alamos Scientific Laboratory, 1973.
[23] Russer, P.. Electromagnetics, Microwave Circuit and Antenna Design for Communications Engineering, Artech House, 2003.
[24] Sacks, Z. S., Kingsland, D. M., Lee, R. and Lee, J. F.. A perfectly matched anisotropic absorber for use as an absorbing boundary condition, IEEE trans. Antennas Propagat., 43 (1995), pp. 14601463.
[25] Simpson, J. J.. Current and future applications of 3-D global earth-ionosphere models based on the full-vector Maxwell's equations FDTD method, Surveys in Geophysics, 30 (2009), pp. 105130.
[26] Streich, R.. Controlled-source electromagnetic approaches for hydrocarbon exploration and monitoring on land, Surveys in Geophysics, 37 (2016), pp. 4780.
[27] Taflove, A. and Hagness, S. C.. Computational Electrodynamics: The Finite-Difference Time-Domain Method, Artech house publishers, 2005.
[28] Warburton, T.. Application of the discontinuous Galerkin method to Maxwell's equations using unstructured polymorphic hp-finite elements, In Discontinuous Galerkin Methods (pp. 451458), Springer Berlin Heidelberg, 2000.
[29] Ward, S.H. and Hohmann, G.W.. Electromagnetic theory for geophysical applications, in Electromagnetic Methods in Applied Geophysics, Vol. 1, chap. 4, pp. 131311, ed. Nabighian, M.N., Society of Exploration Geophysicists, 1991.
[30] Xie, Z., Wang, B. and Zhang, Z.. Space-Time discontinuous Galerkin method for Maxwell's equations, Commun. Comput. Phys., 14 (2013), pp. 916939.
[31] Xie, Z., Wang, J., Wang, B. and Chen, C.. Solving Maxwell's equation in meta-materials by a CG-DG method, Commun. Comput. Phys., 19 (2016), pp. 12421264.
[32] Xin, J. and Cai, W.. Well-conditioned orthonormal hierarchical bases on simplicial elements, J. Sci. Comput., 50 (2012), 446461.
[33] Xin, J., Guo, N. and Cai, W.. On the construction of well-conditioned hierarchical bases for tetrahedral H(curl)-conforming Nédélec elements, J. Comput. Math. 29 (2011), 526542.
[34] Yee, K. S.. Numerical solution of initial boundary value problems involving Maxwell's equations in isotropic media, IEEE Trans. Antennas Propagat., 14 (1966), pp. 302307.
[35] Youssef, N. N.. Radar cross section of complex targets, Proceedings of the IEEE, 77 (1989), pp. 722734.

Keywords

MSC classification

A High-Order Time Domain Discontinuous Galerkin Method with Orthogonal Tetrahedral Basis for Electromagnetic Simulations in 3-D Heterogeneous Conductive Media

  • Jun Yang (a1) (a2), Wei Cai (a3) (a2) and Xiaoping Wu (a1)

Metrics

Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed