3D FDTD Magnetized Plasma Example

Example code for the scattering of electromagnetic radiation off a Gaussian plasma ball in the presence of a background magnetic field.

This code is based off the works of [1-3]. Make sure you have an NVidia GPU otherwise you will have to select the option of non-GPU use, which is very slow. When the code runs you should be able to visualize electromagnetic radiation propagating towards a Gaussian plasma ball, which represents a meteor. Futher analysis such as radar cross section or separation of scattered fields will need to be implemented by you.

Take note, while an ADE-PML (similiar to CPML) absorbing boundary condition is used, using this technique can introduce errors into the simulation [4]. If you find this to be an issue, you may need to implement methods as prescribed in [5] - however this particular method requires more computations and memory.

References

[1] S. Pokhrel, V. Shankar and J. J. Simpson, "3-D FDTD Modeling of Electromagnetic Wave Propagation in Magnetized Plasma Requiring Singular Updates to the Current Density Equation," in IEEE Transactions on Antennas and Propagation, vol. 66, no. 9, pp. 4772-4781, Sept. 2018, doi: 10.1109/TAP.2018.2847601.

[2] Marshall, Robert & Close, S.. (2015). An FDTD model of scattering from meteor head plasma. Journal of Geophysical Research: Space Physics. 120. 10.1002/2015JA021238.

[3] Elsherbeni, Atef & Demir, V.. (2016). The finite-difference time-domain method for electromagnetics with MATLAB® simulations: ACES series, 2nd edition. 10.1049/SBEW514E.

[4] Y. Yu and J. J. Simpson, "A Magnetic Field-Independent Absorbing Boundary Condition for Magnetized Cold Plasma," in IEEE Antennas and Wireless Propagation Letters, vol. 10, pp. 294-297, 2011, doi: 10.1109/LAWP.2011.2139191.

[5] Cho, Jeahoon & Park, Min-Seok & Jung, Kyung-Young. (2020). Perfectly Matched Layer for Accurate FDTD for Anisotropic Magnetized Plasma. Journal of Electromagnetic Engineering and Science. 20. 277-284. 10.26866/jees.2020.20.4.277.