Time-Domain Numerical Solutions of Maxwell Interface Problems with Discontinuous Electromagnetic Waves

ISCAS OpenIR

	Time-Domain Numerical Solutions of Maxwell Interface Problems with Discontinuous Electromagnetic Waves
	Zhang, Y; Nguyen, DD; Du, KW; Xu, J; Zhao, S
	2016
发表期刊	ADVANCES IN APPLIED MATHEMATICS AND MECHANICS
ISSN	2070-0733
卷号	8 期号:3 页码:353-385
摘要	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 Galerkin Maxwell 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.; 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 Galerkin Maxwell 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.
收录类别	SCI
关键词	Maxwell's Equations Finite-difference Time-domain (Fdtd) Discontinuous Galerkin Time-domain (Dgtd) Transverse Magnetic (Tm) And Transverse Electric (Te) sysTems High Order Interface Treatments Matched Interface And Boundary (Mib)
部门归属	Chinese Acad Sci, Inst Software, Beijing 100190, Peoples R China. Univ Alabama, Dept Math, Tuscaloosa, AL 35487 USA.
语种	英语
WOS记录号	WOS:000369433100001
引用统计
内容类型	期刊论文
URI标识	http://ir.iscas.ac.cn/handle/311060/17331
专题	中国科学院软件研究所
推荐引用方式 GB/T 7714	Zhang, Y,Nguyen, DD,Du, KW,et al. Time-Domain Numerical Solutions of Maxwell Interface Problems with Discontinuous Electromagnetic Waves[J]. ADVANCES IN APPLIED MATHEMATICS AND MECHANICS,2016,8(3):353-385.
APA	Zhang, Y,Nguyen, DD,Du, KW,Xu, J,&Zhao, S.(2016).Time-Domain Numerical Solutions of Maxwell Interface Problems with Discontinuous Electromagnetic Waves.ADVANCES IN APPLIED MATHEMATICS AND MECHANICS,8(3),353-385.
MLA	Zhang, Y,et al."Time-Domain Numerical Solutions of Maxwell Interface Problems with Discontinuous Electromagnetic Waves".ADVANCES IN APPLIED MATHEMATICS AND MECHANICS 8.3(2016):353-385.

条目包含的文件
文件名称/大小	文献类型	版本类型	开放类型	使用许可
div-class-title-time（3560KB）			开放获取	使用许可	请求全文