Lightning electromagnetic pulse simulation using 3D-FDTD method (Comparison between PEC and UPLMBoundary Conditions)

In this paper we present an extensive comparison between the uses of two formulations related to the finite differences time domain method in three dimensions (3D-FDTD) applied to the lightning electromagnetic pulses (LEMP) analysis. Electromagnetic models are implemented for the representation of the lightning return stroke current. Thus, the first formulation is based on the implementation of the Yee algorithm using perfect electric conductor (PEC) boundary conditions. The second formulation consists on the integration of Taflove formulation in the 3D-FDTD method using uniaxial perfectly matched layers (UPML) boundary conditions. For that effect, two computational electromagnetic codes have been developed, in Matlab programming environment, in order to determine the lightning return stroke current distribution and associated electromagnetic field components. Finally, for validation needs, the obtained simulation results, especially the lightning vertical electric field and the magnetic flux density, are compared with measured results taken from specialized literature.

[1]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[2]  Tran Huu Thang,et al.  FDTD Simulations of Corona Effect on Lightning-Induced Voltages , 2014, IEEE Transactions on Electromagnetic Compatibility.

[3]  Dennis M. Sullivan,et al.  Electromagnetic Simulation Using the FDTD Method , 2000 .

[4]  Yoshihiro Baba,et al.  Power System Transients: Theory and Applications , 2013 .

[5]  S. Yokoyama,et al.  Thin Wire Representation in Finite Difference Time Domain Surge Simulation , 2002, IEEE Power Engineering Review.

[6]  Allen Taflove,et al.  Computational Electrodynamics the Finite-Difference Time-Domain Method , 1995 .

[7]  Vladimir A. Rakov,et al.  Applications of the FDTD Method to Lightning Electromagnetic Pulse and Surge Simulations , 2014, IEEE Transactions on Electromagnetic Compatibility.

[8]  N. Nagaoka,et al.  Parametric Study on Unit Step Responses of Impulse Voltage Measuring Systems Based on FDTD Simulations , 2013, IEEE Transactions on Power Delivery.

[9]  Vernon Cooray,et al.  Evaluation of Lightning Return Stroke Current Using Measured Electromagnetic Fields , 2012 .

[10]  Yoshihiro Baba,et al.  On the use of lumped sources in lightning return stroke models , 2005 .

[11]  R. Stephenson.  A and V , 1962, The British journal of ophthalmology.

[12]  M. Kuzyk,et al.  Three-dimensional optical pulse simulation using the FDTD method , 2000 .

[13]  Vladimir A. Rakov,et al.  A new lightning return stroke model based on antenna theory , 2000 .

[14]  Y. Baba,et al.  On the mechanism of attenuation of current waves propagating along a vertical perfectly conducting wire above ground: application to lightning , 2005, IEEE Transactions on Electromagnetic Compatibility.

[15]  Yoshihiro Baba,et al.  Characteristics of electromagnetic return-stroke models , 2003 .

[16]  Matthew N. O. Sadiku,et al.  Numerical Techniques in Electromagnetics , 2000 .

[17]  Yoshihiro Baba,et al.  Electromagnetic models of the lightning return stroke , 2007 .

[18]  Y. Baba,et al.  Applications of Electromagnetic Models of the Lightning Return Stroke , 2008, IEEE Transactions on Power Delivery.

[19]  Roger F. Harrington,et al.  Field computation by moment methods , 1968 .

[20]  N. Nagaoka,et al.  An Improved Thin Wire Representation for FDTD Computations , 2008, IEEE Transactions on Antennas and Propagation.

[21]  Y. Baba,et al.  Electric and Magnetic Fields Predicted by Different Electromagnetic Models of the Lightning Return Stroke Versus Measured Fields , 2009, IEEE Transactions on Electromagnetic Compatibility.

[22]  Vladimir A. Rakov,et al.  Review and evaluation of lightning return stroke models including some aspects of their application , 1998 .

[23]  Annett Baier,et al.  Electromagnetic Simulation Using The Fdtd Method , 2016 .

[24]  Yoshihiro Baba,et al.  An Improvement of a Thin Wire Representation for FDTD Electromagnetic and Surge Calculations , 2009 .

[25]  J. Bérenger.  Three-Dimensional Perfectly Matched Layer for the Absorption of Electromagnetic Waves , 1996 .

[26]  S.H.H. Sadeghi,et al.  Analysis of lightning-radiated electromagnetic fields in the vicinity of lossy ground , 2005, IEEE Transactions on Electromagnetic Compatibility.

[27]  Yoshihiro Baba,et al.  On the transmission line model for lightning return stroke representation , 2003 .

[28]  N. Nagaoka,et al.  Modification on a Thin-Wire Representation for FDTD Calculations in Nonsquare Grids , 2008, IEEE Transactions on Electromagnetic Compatibility.