Program 6A,6B,6C Unitless 1D FDTD Open Boundary with Different sources

1D unitless FDTD with open boundaries to compute propagation of fields

Youtube Links :- http://www.youtube.com/watch?v=2W5DOn2tNiM
http://www.youtube.com/watch?v=8Z0hkDHrFtM
http://www.youtube.com/watch?v=SpM7MEX9GFk

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 :-

1D FDTD with Open Boundaries with
a) Hard Always 'ON' source
b) Sinusoidal source
c) Gaussian Pulse source

Objective of the program is to solve for the Maxwell's equation for an x-directed z-polarized TEM wave containing the y-directed magnetic field Hy and z-directed electric field Ez. The fields are updated at every timestep, in a unitless space, where all physical parameters of free space are normalized to 1, using standard update equations obtained from the difference form of Maxwell's curl equations. The field points are defined in a grid described by Yee's algorithm. The H fields are defined at every half coordinate of a spacestep and E fields are defined at every coordinate point. Also, the time update 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 for-loop updates spanning entire spatial grid inside a main for-loop for time update spanning the entire time-grid.