We present a new method for evaluating the electromagnetic field created by a lightning discharge under the assumption of a finite conductivity ground. In this case the expressions for the vertical and radial components of the electric field and for the azimuthal component of the magnetic one involve the so-called Sommerfeld integrals, which require a huge computational effort due to their characteristics of slow convergence. From a mathematical standpoint the Sommerfeld’s integrals are improper integrals (i.e. over a semi-infinite domain). The presence in the integrand of the highly oscillating Bessel function and of an integrable singularity are the main causes of numerical problems. The proposed algorithm allows the extraction of the singularity and performs part of the integral analytically. Consequently, a standard quadrature formula can be applied to the remaining term. The method has been developed in the frequency domain and, for the evaluation of each field component, about 1000 Sommerfeld integrals are necessary, each of them requiring about six seconds of CPU time on a Pentium IV, 1.6 GHz PC. This method has been numerically tested in the simulation of a straight vertical lightning channel electromagnetic field, assuming typical values for the ground conductivity. In particular, the field has been evaluated at different distances from the channel in order to point out the distance ranges where conductivity effects cannot be neglected. All the simulations have been carried out assuming, for the channel base current, the well-known Heidler’s function and for the spatial-temporal distribution of the lightning current along the return stroke channel (i.e. the current return-stroke model) the classical MTLE expression.

Evaluation of lightning electromagnetic fields over a lossy ground

DELFINO, FEDERICO;PROCOPIO, RENATO;ROSSI, MANSUETO
2005-01-01

Abstract

We present a new method for evaluating the electromagnetic field created by a lightning discharge under the assumption of a finite conductivity ground. In this case the expressions for the vertical and radial components of the electric field and for the azimuthal component of the magnetic one involve the so-called Sommerfeld integrals, which require a huge computational effort due to their characteristics of slow convergence. From a mathematical standpoint the Sommerfeld’s integrals are improper integrals (i.e. over a semi-infinite domain). The presence in the integrand of the highly oscillating Bessel function and of an integrable singularity are the main causes of numerical problems. The proposed algorithm allows the extraction of the singularity and performs part of the integral analytically. Consequently, a standard quadrature formula can be applied to the remaining term. The method has been developed in the frequency domain and, for the evaluation of each field component, about 1000 Sommerfeld integrals are necessary, each of them requiring about six seconds of CPU time on a Pentium IV, 1.6 GHz PC. This method has been numerically tested in the simulation of a straight vertical lightning channel electromagnetic field, assuming typical values for the ground conductivity. In particular, the field has been evaluated at different distances from the channel in order to point out the distance ranges where conductivity effects cannot be neglected. All the simulations have been carried out assuming, for the channel base current, the well-known Heidler’s function and for the spatial-temporal distribution of the lightning current along the return stroke channel (i.e. the current return-stroke model) the classical MTLE expression.
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11567/269660
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus ND
  • ???jsp.display-item.citation.isi??? ND
social impact