An Efficient Finite-Element Time-Domain Method via Hierarchical Matrix Algorithm for Electromagnetic Simulation
Keywords:
An Efficient Finite-Element Time-Domain Method via Hierarchical Matrix Algorithm for Electromagnetic SimulationAbstract
An efficient finite-element timedomain (FETD) method based on the hierarchical (H-) matrix algorithm is presented. The FETD method is on the basis of the second-order vector wave equation, obtained by eliminating one of the field variables from Maxwell’s equations. The time-dependent formulation employs the Newmark-beta method which is known as an unconditional stable time-integration method. H- matrix algorithm is introduced for the direct solution of a large sparse linear system at each time step, which is a serious handicap in conventional FETD method. H-matrix algorithm provides a data-sparse way to approximate the LU triangular factors of the FETD system matrix. Using the H-matrix arithmetic, the computational complexity and memory requirement of H-LU decomposition can be significantly reduced to almost logarithmic-linear. Once the H-LU factors are obtained, the FETD method can be computed very efficiently at each time step by the H-matrix formatted forward and backward substitution (H- FBS). Numerical examples are provided to illustrate the accuracy and efficiency of the proposed FETD method for the simulation of three-dimension (3D) electromagnetic problems.
Downloads
References
F. L. Teixeira, “A Summary Review on 25
Years of Progress and Future Challenges in
FDTD and FETD Techniques,” Applied
Computational Electromagnetic Society
(ACES) Journal, vol. 25, no. 1, pp. 1-14, Jan.
J. F. Lee, R. Lee, and A. C. Cangellaris,
“Time-Domain Finite Element Methods,”
IEEE Transactions on Antennas and
Propagation., vol. 45, pp. 430-442, 1997.
M. F. Wongm, O. Picon, and V. F. Hanna, “A
Finite-Element Method Based on Whitney
Forms to Solve Maxwell Equations in the
Time-Domain,” IEEE Trans. Magn., vol. 31,
pp. 1618-1621, 1995.
B. Donderici and F. L. Teixeira, “Mixed
Finite-Element Time-Domain Method for
Transient Maxwell Equations in Doubly
Dispersive Media,” IEEE Trans. Microwave
Theory Tech., vol. 56, no. 1, pp. 113-120,
H. –P. Tsai, Y. Wang, and T. Itoh, “An
Unconditionally Stable Extended (USE)
Finite-Element Time- Domain Solution of
Active Nonlinear Microwave Circuits Using
Perfectly Matched Layers,” IEEE Tran.
Microwave Theory Tech., vol. 50, no. 10, pp.
-2232, Oct. 2002.
N. Marais and Davidson, “Numerical
Evaluation of High-Order Finite Element
Time Domain Formulations in
Electromagnetics,” IEEE Trans. Antennas
Propagat., vol. 56, vol. 12, pp. 3743-3751,
S. D. Gedney and U. Navsariwala, “An
Unconditionally Stable Finite Element TimeDomain Solution of the Vector Wave
Equation,” IEEE Tran. Microwave and
Guided Wave Letters, vol. 5, no. 10, pp. 332-
, Oct. 1995.
M. Feliziani and F. Maradei, “Hybrid Finite
Element Solution of Time Dependent Maxwell’s Curl Equations,” IEEE Trans. Magn.,
vol. 31, no. 3, pp. 1330-1335, May 1995.
R. N. Rieben, G. H. Rodrigue, and D. A.
White, “A High-Order Mixed Vector Finite
Element Method for Solving the Time
Dependent Maxwell Equations on
Unstructured Grids,” J. Comp. Phys., vol.
, pp. 490-519, 2005.
S. D. Gedney, C. Luo, J. A. Roden, R. D.
Crawford, B. Guernsey, J. A. Miller, T.
Kramer, and E. W. Lucas, “The
Discontinuous Galerkin Finite-Element TimeDomain Method Solution of Maxwell’s
Equation,” Applied Computational
Electromagnetic Society (ACES) Journal, vol.
, no. 2, pp. 129-142, April 2009.
N. V. Kantartzis and T. D. Tsiboukis, Modern
EMC Analysis Techniques - Volume I: TimeDomain Computational Schemes. San Rafael,
CA, USA: Morgan & Claypool Publishers,
T. V. Yioultsis, N. V. Kantartzis, C. S.
Antonopoulos, and T. D. Tsiboukis, “A Fully
Explicit Whitney-Element Time-Domain
Scheme with Higher Order Vector Finite
Elements for Three-Dimensional HighFrequency Problems,” IEEE Trans. Magn.,
vol. 34, no. 5, pp. 3288-3291, Sept. 1998.
B. He and F. L. Teixeira, “A Sparse and
Explicit FETD via Approximate Inverse
Hodge (Mass) Matrix,” IEEE Microw. Wireless Comp. Lett., vol. 16, no. 6, pp. 348-350,
Y. Saad, Iterative Methods for Sparse Linear
Systems. New York: PWS Publishing, 1996.
Z. Jia and B. Zhu, “A Power Sparse Approximate Inverse Preconditioning Procedure for
Large Sparse Linear Systems,” Numer.
Linear Algebra Appl., vol. 16, pp. 259-299,
A. George, “Nested Dissection of a Regular
Finite Element Mesh,” SIAM J. on Numerical
Analysis, 10(2):345-363, April 1973.
W. Hackbusch, “A Sparse Matrix Arithmetic
Based on -Matrices. I. Introduction to -
Matrices,” Computing, 62 (2):89-108, 1999.
W. Hackbusch and B. Khoromskij, “A Sparse
-Matrix Arithmetic. Part II: Application to
Multi-Dimensional Problems,” Computing, 6,
pp. 21-47, 2000.
S. Börm and L. Grasedyck, “Low-Rank
Approximation of Integral Operators by Interpolation,” Computing, vol. 72, pp. 325-332,
M. Bebendorf and S. Rjasanow, “Adaptive
Low-Rank Approximation of Collocation
Matrices,” Computing, 70, pp. 1-24, 2003.
M. Bebendorf and W. Hackbusch, “Existence
of -Matrix Approximants to the Inverse FE
Matrix of Elliptic Operators with L∞-Coefficients,” Numer. Math., 95 (2003), pp.1-28.
H. Liu and D. Jiao, “A Direct Finite-ElementBased Solver of Significantly Reduced
Complexity for Solving Large-Scale ElectroWAN, CHEN, SHE, DING, FAN: AN EFFICIENT FETD METHOD VIA HIERARCHICAL MATRIX ALGORITHM 593
magnetic Problems,” IMS 2009, pp. 177-180,
H. Liu and D. Jiao, “Existence of -Matrix
Representations of the Inverse Finite-Element
Matrix of Electrodynamic Problems and -
Based Fast Direct Finite-Element Solvers,”
IEEE Trans. on Microwave Theory and
Techniques, vol. 58, no. 12, pp. 3697-3709,
Dec. 2010.
S. Borm, L. Grasedyck, and W. Hackbusch,
“Induction to Hierarchical Matrices with
Applications,” Engineering Analysis with
Boundary Elements, no. 27, pp. 405-422,
L. Grasedyck and W. Hackbusch, “Construction and Arithmetics of -Matrices,”
Computing, vol. 70, no. 4, pp. 295-344,
August 2003.
M. Bebendorf, “Why Finite Element
Discretizations can be Factored by Triangular
Hierarchical Matrices,” SIAM J. Matrix Anal.
Appl., 45(4):1472-1494, 2007.
J. P. Berenger, “A Perfectly Matched Layer
for the Absorption of Electromagnetic
Waves,” J. Compru. Phys., vol. 114, pp. 185-
, Oct. 1994.
S. D. Gedney, “An Anisotropic Perfectly
Matched Layer-Absorbing Medium for the
Truncation of FDTD Lattices,” IEEE
Transactions on Antennas and Propagation,
vol. 44, no. 11, pp. 1630-1639, 1996.
A. Bossavit, “Whitney Forms: A Class of
Finite Elements for Three Dimensional
Computations in Electromagnetism,” IEE
Proc. Pt. A, vol. 135, no. 8, pp. 493-500,
Nov. 1988.
Lei Du, R. S. Chen, and Z. B. Ye, “Perfectly
Matched Layers Backed with the First Order
Impedance Boundary Condition for the TimeDomain Finite-Element Solution of
Waveguide Problems,” Microwave and
Optical Technology Letters, vol. 50, no. 3, pp.
-843, March 2008.
R. S. Chen, E. K. N. Yung, C. H. Chan, D. X.
Wang, and D. G. Fang, “Application of the
SSOR Preconditioned CG Algorithm to the
Vector FEM for 3-D Full-Wave Analysis of
Electromagnetic-Field Boundary-Value Problems,” IEEE Trans. Microwave Theory
Tech., vol. 50, no. 4, pp. 1165-1172, 2002.