Numerical modeling of Maxwells equations with dispersive materials

Numerical solver EMWSolver3D using the FDTD-method for 3D Maxwell equations is developed. Maxwell equations' approximation in the integral form on Yee lattice is used for the numerical solution. Problems in bounded and infinite domains are possible to be solved using UPML absorbing layer. Dispersive material model is included for the media with both negative dielectric permittivity and magnetic permeability research. EMWSolver3D involves modules for both single-processor systems and supercomputer IBM BlueGene /P that allow to perform calculations on grids of the order of 400 wave lengths in each direction. Absorbing layer parameters affecting the accuracy of solution, such that the boundary reflection is low and does not influence on the magnetic field distibution in calculation region, are obtained. In this work the illustration of plane wave impinging on metamaterial with $\varepsilon=\mu-1$ with different angles are shown. Surface plasmon waves on metamaterial interface can be observed.