A full finite difference time domain methodology is developed for electromagnetic wave propagation in a plasma. The finite difference grid is consistent with central difference approximation of the curl, divergence, and gradient operators that appear in the joint equations of Euler and Maxwell, and the coupling effects between the fluid velocity and the electric field. To accomplish the time advancement, the central difference approximation is invoked for the time derivatives and leapfrog concepts are employed. The resulting difference equations converge to the exact equations, provided that the developed stability requirement is satisfied.

Finally, numerical results are provided and compared with the inverse fast Fourier transform results of closed-form, frequency domain solutions for the half space problem; the agreement between solutions is shown to be excellent. ¿ American Geophysical Union 1994