Time-Domain Analysis of Left-Handed Materials Based on a Dispersive Meshless Method with PML Absorbing Boundary Condition

TL;DR

This paper introduces a novel dispersive meshless method with PML boundary conditions for accurate time-domain electromagnetic analysis of left-handed materials, addressing a gap in simulation techniques for metamaterials.

Contribution

It presents the first dispersive meshless formulation incorporating PML boundaries for LH media, enhancing modeling flexibility and accuracy in wave propagation simulations.

Findings

The proposed method accurately simulates wave propagation in LH materials.

Inclusion of PML improves boundary absorption and reduces reflections.

The approach offers a flexible alternative to grid-based methods for complex geometries.

Abstract

In this paper, we have proposed a dispersive formulation of scalar-based meshless method for time-domain analysis of electromagnetic wave propagation through left-handed (LH) materials. Moreover, we have incorporated Berenger's perfectly matched layer (PML) absorbing boundary condition (ABC) into the dispersive formulation to truncate open-domain structures. In general, the LH medium as a kind of dispersive media can be described by frequency-dependent constitutive parameters. The most appropriate numerical techniques for analysis of LH media are dispersive formulations of conventional numerical methods. In comparison to the conventional grid-based numerical methods, it is proved that meshless methods not only are strong tools for accurate approximation of derivatives in Maxwell's equations but also can provide more flexibility in modeling the spatial domain of problems. However, we have not seen any reports on using dispersive forms of meshless methods for simulation of wave propagation in metamaterials and applying any PML ABCs to dispersive formulation of meshless method. The proposed formulation in this paper enables us to take advantage of meshless methods in analysis of LH media. For modeling the frequency behavior of the medium in the proposed dispersive formulation, we have used auxiliary differential equation (ADE) method based on the relations between electromagnetic fields intensities and current densities. Effectiveness of the proposed formulation is verified by a numerical example; also, some basic factors which affect the accuracy and computational cost of the simulations are studied.