An Adaptive Model Order Reduction Approach for the Finite Element Method in Time Domain in Electromagnetics

TL;DR

This paper presents TA-ROMTD, an adaptive reduced order modeling technique for finite element time domain electromagnetic simulations that significantly cuts computational costs without sacrificing accuracy.

Contribution

It introduces a novel on-the-fly adaptive ROM method that switches between FEMTD and ROMTD based on error estimation, eliminating the need for prior problem knowledge.

Findings

Achieves substantial reduction in simulation time from hours to minutes.

Maintains high accuracy in complex electromagnetic structures.

Demonstrates effectiveness across various antenna and metamaterial examples.

Abstract

Time domain simulations are crucial for analyzing transient behavior and broadband responses in electromagnetic problems. However, conventional numerical methods such as finite element method in time domain (FEMTD) and finite difference time domain, can be computationally demanding due to their high-dimensional nature, making large-scale simulations impractical for design optimization and real-time analysis. This paper introduces TA-ROMTD, a time-adaptive reduced order model (ROM) for FEMTD simulations that significantly reduces computational costs while maintaining accuracy. The method alternates between FEMTD and a reduced order model in time domain (ROMTD), using an error estimator to detect when the ROMTD solution loses accuracy and switching back to FEMTD to update the ROM with new data. Thus, TA-ROMTD does not require prior knowledge of the problem, as the ROM is constructed on the fly using FEMTD data. A key feature of this approach is the use of a coarse time step during FEMTD time intervals, capturing essential system dynamics while minimizing computational overhead. By reducing the number of degrees of freedom, this method enables efficient electromagnetic simulations, making it a powerful tool for antenna and microwave circuit design. The efficiency of the TA-ROMTD strategy is demonstrated through numerical examples, including an antipodal Vivaldi antenna, a dielectric resonator antenna, a fully metallic dual-polarization frequency selective surface, and a mushroom-type electromagnetic bandgap structure. These cases show the capability of the proposed approach to achieve accurate solutions while significantly reducing the computation time from tens of hours to minutes.