2D FDTD of a region with Mur's absorbing boundary

Youtube links:-
http://www.youtube.com/watch?v=GMHwmWUMnbw&feature=related
http://www.youtube.com/watch?v=6ho03GQ6u04
http://www.youtube.com/watch?v=NqvVcuKR5Is&feature=related

Written for Course :-
Computational Electromagnetics, Fall 2011
Department of Electrical Engineering
Indian Institute of Technology Madras (IITM)
Chennai - 600036, India

Authors :-
Sathya Swaroop Ganta, B.Tech., M.Tech. Electrical Engg.
Kayatri, M.S. Engineering Design
Pankaj, M.S. Electrical Engg.
Sumantra Chaudhuri, M.S. Electrical Engg.
Projesh Basu, M.S. Electrical Engg.
Nikhil Kumar CS, M.S. Electrical Engg.

Instructor :-
Dr. Ananth Krishnan
Assistant Professor
Department of Electrical Engineering
Indian Institute of Technology Madras
http://www.ee.iitm.ac.in/~ananthk

Description :-

A 2D TM wave containing the xy-plane polarized magnetic field having components Hy and Hx and z-polarized electric field Ez. The fields are updated at every time step, in a space, where all physical parameters of free space are not normalized to 1 but are given real and known values. The update is done using standard update equations obtained from the difference form of Maxwell's curl equations with very low electric and magnetic conductivities of 4x10^(-4) units incorporated in them. The field points are defined in a grid described by Yee's algorithm. The H fields are defined at every half coordinate of spacesteps. More precisely, the Hx part is defined at every half y-coordinate and full x-coordinate and the Hy part is defined at every half x-coordinate and full y-coordinate and E fields i.e the Ez part is defined at every full x and full y-coordinate points.Also here, the space-step length is taken as 1 micron instead of 1 unit in unitless domain assumed in previous programs. Also, the time update is done using Leapfrog time-stepping. Here, H-fields i.e Hx and Hy are updated every half time-step and E fields i.e Ez are updated every full time-step. This is shown by two alternating vector updates spanning only a part of spatial grid where the wave, starting from source, has reached at that particular time instant avoiding field updates at all points in the grid which is unnecessary at that time instant. These spatial updates are inside the main for-loop for time update, spanning the entire time grid. Also, here, the matrices used as multiplication factors for update equations are initialized before the loop starts to avoid repeated calculation of the same in every loop iteration, a minor attempt at optimization. The
boundary condition here is Mur's Absorbing Boundary Condition (ABC) where the fields at the grid points have electric field values formulated using Engquist Majda one way wave equations where the boundaries give a sense of absorbing the total field incident on them and reflecting none back to the domain.