Program19 2D FDTD Unitless Open Boundary Suboptimal Courant factor 1.0

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 timestep, in a unitless space, where all physical parameters of free space are normalized to 1. The update is done 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 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, 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 three alternating for-loop updates in groups of two (for H) and one (for E) spanning entire spatial grid inside a main for-loop for time update spanning the entire time-grid. The courant stability factor is given a suboptimal value of 0.9*(1/sqrt(2)).Written for Course :-Computational Electromagnetics, Fall 2011Department of Electrical EngineeringIndian Institute of Technology Madras (IITM)Chennai - 600036, IndiaAuthors :-Sathya Swaroop Ganta, B.Tech., M.Tech. Electrical Engg.Kayatri, M.S. Engineering DesignPankaj, 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 KrishnanAssistant ProfessorDepartment of Electrical EngineeringIndian Institute of Technology Madrashttp://www.ee.iitm.ac.in/~ananthk

Requirements:	No special requirements
Platforms:	Matlab
Keyword:	Algorithm,  Alternating,  Coordinate,  Forloop,  Hfields,  Leapfrog,  Shown,  Timestepping,  Xcoordinate
Users rating:	0/10