Program12 1D FDTD with units and conductivities

1D FDTD with units and conductivities

Youtube Links :-
http://www.youtube.com/watch?v=2ky6gkniVlA
http://www.youtube.com/watch?v=QSxIX0C_04c
http://www.youtube.com/watch?v=btK1bLA9Py4
http://www.youtube.com/watch?v=iNZx8TGcG4U

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.