Program14 1D FDTD with Perfect Electric Conductor boundary

The PEC condition is introduced in 1D FDTD

Youtube Links :-
http://www.youtube.com/watch?v=Fa_0O2bMF7M
http://www.youtube.com/watch?v=fTUwNf1EDxc
http://www.youtube.com/watch?v=WWREgLHl8n4
http://www.youtube.com/watch?v=mpRoLNbXjQw

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 electric boundary in the sense that the boundary grid points have zero electric field values irrespective of the influence of external fields.