Program13 1D FDTD with Perfect Magnetic Conductor boundary

Youtube Links :-
http://www.youtube.com/watch?v=sef5U-U9p4w
http://www.youtube.com/watch?v=Zw0fRDcaAIs
http://www.youtube.com/watch?v=2ac7apBUKss
http://www.youtube.com/watch?v=Ud3BPA09bDA

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 :-
The 1D TEM wave is x-directed z-polarized TEM wave containing the y-directed magnetic field Hy and z-directed electric field Ez. The time update in Yee Algorithm is done using Leapfrog time-stepping. Here, the H fields are updated every half time-step and E fileds are updated every full time-step. This is shown by two alternating vector updates spanning entire spatial grid inside a main for-loop for time update spanning the entire time-grid. Yes, there are No For Loops except for time steps. The vector updates span only a part of spatial grid where the wave, starting from source, has reached at that particular time instant (exploiting sparse vectors) avoiding field updates at all points in the grid which is unnecessary at that time instant. The spatial and temporal parameters are not unit less and are given real values. Mostly parameters of free space are used except that very less electrical and magnetic conductivities of 4 x 10^4 units are incorporated. Unlike in previous programs, which have open boundary condition, the boundary condition here is perfect magnetic boundary in the sense that the boundary grid points have zero magnetic field values irrespective of the influence of external fields.