Next Article in Journal
Investigating West African Monsoon Features in Warm Years Using the Regional Climate Model RegCM4
Previous Article in Journal
Chemical Characterization of PM2.5 at Rural and Urban Sites around the Metropolitan Area of Huancayo (Central Andes of Peru)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Interpretation of the Electromagnetic Fields of Lightning Return Strokes

1
Department of Engineering Sciences, Uppsala University, 752 37 Uppsala, Sweden
2
Karolinska Institute, 171 77 stockholm, Sweden
*
Author to whom correspondence should be addressed.
Atmosphere 2019, 10(1), 22; https://doi.org/10.3390/atmos10010022
Submission received: 11 December 2018 / Revised: 24 December 2018 / Accepted: 27 December 2018 / Published: 9 January 2019
(This article belongs to the Section Meteorology)

Abstract

:
Electric and/or magnetic fields are generated by stationary charges, uniformly moving charges and accelerating charges. These field components are described in the literature as static fields, velocity fields (or generalized Coulomb field) and radiation fields (or acceleration fields), respectively. In the literature, the electromagnetic fields generated by lightning return strokes are presented using the field components associated with short dipoles, and in this description the one–to-one association of the electromagnetic field terms with the physical process that gives rise to them is lost. In this paper, we have derived expressions for the electromagnetic fields using field equations associated with accelerating (and moving) charges and separated the resulting fields into static, velocity and radiation fields. The results illustrate how the radiation fields emanating from the lightning channel give rise to field terms varying as 1 / r and 1 / r 2 , the velocity fields generating field terms varying as 1 / r 2 , and the static fields generating field components varying as 1 / r 2 and 1 / r 3 . These field components depend explicitly on the speed of propagation of the current pulse. However, the total field does not depend explicitly on the speed of propagation of the current pulse. It is shown that these field components can be combined to generate the field components pertinent to the dipole technique. However, in this conversion process the connection of the field components to the physical processes taking place at the source that generate these fields (i.e., static charges, uniformly moving charges and accelerating charges) is lost.

1. Introduction

Estimating the field strength and temporal features of electromagnetic fields from lightning return strokes at a given distance is of interest both in engineering studies, where the protection of electrical installations from induced voltages from lightning is concerned [1,2,3], and in physics studies, where the properties of lightning return strokes are extracted from the features of electromagnetic fields measured under conditions where the propagation effects are minimal [4,5]. The procedure used in these studies is to specify the spatial and temporal variation of the return stroke current I ( z , t ) by appealing to a return stroke model and from that calculate the electromagnetic fields [6,7]. Once I ( z , t ) is specified, there are four methods to estimate the electromagnetic fields. Three of these methods, namely, the dipole technique, monopole technique and the apparent charge density technique, are described in [6,7]. The fourth method, based on the field equations pertinent to the moving and accelerating charges is described in [8]. For a given I ( z , t ) , all these techniques generate the same total electromagnetic fields but the various components that constitute the total fields differ in different techniques.
The goals of the present paper are the following. First, it is standard practice today to describe the electromagnetic fields of lightning return strokes in terms of static (field terms decreasing with distance as 1 / r 3 , where r is the distance from the source to the point of observation), induction (field terms decreasing with distance as 1 / r 2 ) and radiation (field terms decreasing as 1 / r ) [9]. In the sections to follow we will call these field components dipole-static, dipole-induction and dipole-radiation. Except in the case of distant radiation fields, this division of the field components cannot be directly attached to the physical processes that generate the electromagnetic fields. In reality, there are only two types of electromagnetic fields. These are the Coulomb and radiation fields. Coulomb fields are produced by stationary and uniformly moving charges. The Coulomb field produced by stationary or static charges is called the static field. When the charges are moving, the Coulomb field has to be modified to take into account their motion and this modified field, to separate it from the static field, is called the velocity field (or generalized Coulomb field). The radiation fields (or acceleration fields) are produced by accelerating charges. The magnetic field is generated either by moving charges or by accelerating charges. Thus, the magnetic field consists of either the velocity fields or the radiation fields, or both. Unfortunately, when one divides the total field into dipole-static (i.e., the field components changing as 1 / r 3 ), dipole-induction (i.e., the field components changing as 1 / r 2 ) and dipole-radiation (i.e., the field components changing as 1 / r ), as is the standard practice, except in the case of distant radiation field, the direct association of the field components with the physical process that generate electromagnetic fields is lost. In this paper, we derive expressions for the electromagnetic fields of a return stroke where each field component is directly associated with the physical process that gives rise to it, i.e., stationary charges, moving charges and accelerating charges. Second, we will show analytically and illustrate by example, that the resulting field components can be combined together to produce the field components that are identified in the literature as dipole-static, dipole-induction and dipole-radiation.

2. Mathematical Analysis

Some of the mathematical formulations used in the present paper are identical to those already presented previously by Cooray and Cooray [10,11]. For example, the approximations listed in Equations (1)–(11) and the field Equations (12)–(23) can be extracted directly from references [10] and [11]. However, since the other equations given in this paper are constructed using these equations, for the sake of completeness and for easy reference, they are also given here.
The problem under consideration is the calculation of electromagnetic fields from a return stroke channel when I ( z , t ) is specified. The normal procedure for such a calculation is to divide the channel into infinitesimal channel sections of length d z and first estimate the electromagnetic fields from such a channel element located at height z , where z is the height of the channel element. Once the electric and magnetic fields produced by this channel element are known, the total field can be obtained by summing the contributions from all these elements. Let us consider the element d z located at height z . The first step is to estimate the electromagnetic fields generated by the channel element. For convenience we will treat the problem in frequency domain first and later convert it into time domain. Assume that the current flowing along the channel element is given by i ( z ) e j ω t . We consider the case where this current travels along the channel element with speed u and is absorbed at the other end of the channel element. In other words, the current that appears at the bottom of the channel element at any time t will appear at the top of the channel element after a time delay of d z / u . The geometry necessary for the calculation of the electromagnetic fields from this channel element using the technique pertinent to the moving and accelerating charges is shown in Figure 1.
Before we proceed with the analysis, let us consider some of the geometrical simplifications that can be used in the analysis. First of all, we assume that the distance to the point of observation, r , is such that r > > d z . When this condition is satisfied one can also make the following simplifications:
δ θ 2 = ( θ 2 θ ) = d z sin θ 2 r   ;   δ θ 1 = ( θ θ 1 ) = d z sin θ 2 r
cos δ θ 1 1   ;   cos δ θ 2 1
r 1 = r + d z cos θ 2 ;   r 2 = r d z cos θ 2
1 r 1 = 1 r { 1 d z cos θ 2 r } ;   1 r 2 = 1 r { 1 + d z cos θ 2 r }
1 r 1 2 = 1 r 2 { 1 d z cos θ r } ;   1 r 2 2 = 1 r 2 { 1 + d z cos θ r }
sin θ 1 = sin θ { 1 d z cos θ 2 r } ;   sin θ 2 = sin θ { 1 + d z cos θ 2 r }
cos θ 1 = cos θ + d z sin 2 θ 2 r ;   cos θ 2 = cos θ d z sin 2 θ 2 r
1 u cos θ 1 c = { 1 u c cos θ } { 1 u d z sin 2 θ 2 r c ( 1 u c cos θ ) }
1 u cos θ 2 c = { 1 u c cos θ } { 1 + u d z sin 2 θ 2 r c ( 1 u c cos θ ) }
1 1 u cos θ 1 c = 1 1 u cos θ c { 1 + u d z sin 2 θ 2 r c ( 1 u c cos θ ) }
1 1 u cos θ 2 c = 1 1 u cos θ c { 1 u d z sin 2 θ 2 r c ( 1 u c cos θ ) }
Now we are in a position to write down the expressions for the electromagnetic fields. The electromagnetic fields generated by the channel element can be divided into different components as follows: (a) The electric and magnetic radiation fields generated at the initiation and termination of the current at the end points of the channel element due to charge acceleration and deceleration, respectively; (b) the electric and magnetic velocity fields generated by the movement of charges along the channel element; (c) the static fields generated by the accumulation of charges at the two ends of the channel element. Let us consider these different field components separately. In writing down these field components, we will depend heavily on the results published previously by Cooray and Cooray [10,11]. The field expressions identified by Equations (12)–(23) can be constructed easily from the results presented in [10].

2.1. Radiation Field Generated by the Charge Acceleration and Deceleration at the Ends of the Channel Element

The electric radiation field generated by the initiation of current at the bottom of the channel element and by the termination of that current at the top of the channel element is given by
d e r a d = i ( z ) u 4 π ε o c 2 { e j ω ( t r 1 / c ) sin θ 1 r 1 [ 1 u cos θ 1 c ] a θ 1 e j ω ( t d z / u r 2 / c ) sin θ 2 r 2 [ 1 u cos θ 2 c ] a θ 2 }
The above expression is exact and does not contain any approximations. In order to extract the electric fields of an infinitesimal current element, we will write down the components of this electric field in the direction of a r and a θ using the geometrical approximations listed in Equations (1)–(11), which are valid when r > > d z . Moreover, we also assume that ω r / c < < 1 and ω d z / u < < 1 . Using these approximations and keeping only the first order terms with respect to d z , the components of this field in the directions of a r and a θ become (with t = t r / c ; see Appendix A for the derivation):
d e r a d , θ = i ( z ) e j ω t u sin θ 4 π ε o c 2 r ( 1 u c cos θ ) { j ω d z u j ω d z cos θ c 2 d z cos θ r + u d z sin 2 θ r c ( 1 u c cos θ ) } a θ
d e r a d , r = i ( z ) e j ω t 4 π ε o c 2 r 2 { u d z sin 2 θ ( 1 u c cos θ ) } a r
These two equations define the total radiation field produced by the channel element.

2.2. Velocity Field Generated by the Charges Moving from A to B

The velocity field generated as the current pulse propagates along the channel element can be written as
d e v e l = i ( z ) e j ω t d z 4 π ε o r 2 u [ 1 u c cos θ ] 2 [ 1 u 2 c 2 ] a r i ( z ) e j ω t d z 4 π ε o r 2 c [ 1 u c cos θ ] 2 [ 1 u 2 c 2 ] a z
The components of this field in the direction of a r and a θ after some mathematical manipulations are given by
d e v e l , r = i ( z ) e j ω t d z 4 π ε o r 2 u [ 1 u c cos θ ] [ 1 u 2 c 2 ] a r
d e v e l , θ = i ( z ) e j ω t sin θ d z 4 π ε o r 2 c [ 1 u c cos θ ] 2 [ 1 u 2 c 2 ]   a θ

2.3. Electrostatic Field Generated by the Accumulation of Charge at A and B

As the positive current leaves point A, negative charge accumulates at A and when the current is terminated at B, positive charge is accumulated there. The static Coulomb field produced by these stationary charges is given by
d e s t a t ( t ) = i ( z ) 4 π ε o j ω [ e j ω ( t r 1 / c ) r 1 2 a r 1 e j ω ( t d z / u r 2 / c ) r 2 2 a r 2 ]
After using the approximations given earlier and following a procedure similar to that given in Appendix A, the components of this field in the direction of a r and a θ can be written as
d e s t a t , r = i ( z ) e j ω t d z 4 π ε o r 2 j ω { j ω cos θ c + cos θ r j ω u } a r
d e s t a t , θ ( t ) = i ( z ) e j ω t d z 4 π ε o r 2 j ω { sin θ r } a θ

2.4. Magnetic Radiation Field Generated during the Initiation and Termination of the Current

The magnetic radiation field generated during the initiation and termination of the current at the ends of the channel element is given by
d b r a d , ϕ = i ( z ) u 4 π ε o c 3 { sin θ 1 e j ω ( t r 1 / c ) r 1 [ 1 u cos θ 1 c ] sin θ 2 e j ω ( t d z / u r 2 / c ) r 2 [ 1 u cos θ 2 c ] } a ϕ
Utilizing the geometrical approximations mentioned earlier and following procedure almost identical to that presented in Appendix A (and keeping only the second order terms in d z ), one obtains
d b r a d , ϕ = i ( z ) e j ω t u sin θ d z 4 π ε o c 3 r 1 [ 1 u cos θ c ] { j ω u j ω cos θ c 2 cos θ r + u sin 2 θ r c ( 1 u c cos θ ) } a φ

2.5. Magnetic Velocity Field Generated as the Current Pulse Propagates Along the Channel Element

The magnetic velocity field generated during the passage of the current along the channel element is given by
d b v e l , ϕ = i ( z ) e j ω t d z sin θ 4 π ε o r 2 c 2 [ 1 u c cos θ ] 2 [ 1 u 2 c 2 ] a φ

3. Electromagnetic Fields of a Channel Element in the Time Domain

From the frequency domain equations, the time domain equations for the electric and magnetic fields generated by the channel element can be written directly. The results are the following.

3.1. Radiation Fields

d E r a d , θ ( t ) = d z sin θ   4 π ε o c 2 r { I ( z , t ) t 2 u cos θ r ( 1 u c cos θ ) I ( z , t ) + u 2 sin 2 θ r c ( 1 u c cos θ ) 2 I ( z , t ) } a θ
d E r a d , r ( t ) = I ( z , t ) d z 4 π ε o c 2 r 2 { u sin 2 θ ( 1 u c cos θ ) } a r
d B r a d , ϕ ( t ) = d z sin θ 4 π ε o c 3 r { I ( z , t ) t 2 u cos θ r ( 1 u c cos θ ) I ( z , t ) + u 2 sin 2 θ r c ( 1 u c cos θ ) 2 I ( z , t ) } a φ

3.2. Velocity Fields

d E v e l , θ ( t ) = I ( z , t ) d z sin θ 4 π ε o r 2 c [ 1 u c cos θ ] 2 [ 1 u 2 c 2 ] a θ
d E v e l , r ( t ) = I ( z , t ) d z 4 π ε o r 2 u [ 1 u c cos θ ] [ 1 u 2 c 2 ] a r
d B v e l , ϕ ( t ) = I ( z , t ) d z sin θ 4 π ε o r 2 c 2 [ 1 u c cos θ ] 2 [ 1 u 2 c 2 ] a φ

3.3. Static Fields

d E s t a t , r ( t ) = d z 4 π ε o r 2 { cos θ c I ( z , t ) 1 u I ( z , t ) + 2 cos θ   r 0 t I ( z , τ ) d τ } a r
d E s t a t , θ ( t ) = d z 4 π ε o r 2 { sin θ r 0 t I ( z , τ ) d τ } a θ
In the above equations τ = τ r / c .

4. Comparison of the Fields of the Channel Element with Fields of a Short Dipole

The set of equations given in the previous section describes the radiation, velocity and static fields generated by the channel element. Each of these field terms are associated with the particular process that generates these fields. Note also that each term is associated with the speed of propagation of the current pulse. Now, let us sum up the fields in the directions of a r and a θ without any regard to the physical mechanism of the field generation. This generates the following field components
d E θ ( t ) = d z sin θ 4 π ε o { 1 c 2 r I ( z , t ) t + 1 r 3 0 t I ( z , τ ) d τ + I ( z , t )   u 2 sin 2 θ r 2 c 3 ( 1 u c cos θ ) 2 I ( z , t ) 2 u cos θ c 2 r 2 ( 1 u c cos θ ) + I ( z , t ) ( 1 u 2 / c 2 ) c r 2 ( 1 u c cos θ ) 2 } a θ
d E r ( t ) = d z 4 π ε o { 2 cos θ   r 3 0 t I ( z , τ ) d τ + I ( z , t ) cos θ r 2 c I ( z , t ) u r 2 + I ( z , t ) u sin 2 θ c 2 r 2 ( 1 u c cos θ ) + I ( z , t ) ( 1 u 2 / c 2 ) u r 2 ( 1 u c cos θ ) } a r
d B ϕ ( t ) = d z sin θ 4 π ε o { 1 c 3 r I ( z , t ) t I ( z , t ) 2 u cos θ c 3 r 2 ( 1 u c cos θ ) + I ( z , t ) u 2 sin 2 θ r 2 c 4 ( 1 u c cos θ ) 2 + I ( z , t ) ( 1 u 2 / c 2 ) c 2 r 2 ( 1 u c cos θ ) 2 } a ϕ
With some mathematical manipulations, one can show that these equations will reduce to (the mathematical steps necessary for this reduction are given in Appendix B)
d E θ ( t ) = d z sin θ 4 π ε o { 1 r 3 0 t I ( z , τ ) d τ + I ( z , t ) c r 2 + 1 c 2 r I ( z , t ) t } a θ
d E r ( t ) = d z 4 π ε o { 2 cos θ   r 3 0 t I ( z , τ ) d τ + 2 cos θ r 2 c I ( z , t ) t } a r
d B ϕ ( t ) = d z sin θ 4 π ε o { 1 c 3 r I ( z , t ) t + I ( z , t ) r 2 c 2 } a ϕ
The Equations (35)–(37) are identical to the fields of a short dipole [7,12]. This proves the equivalence of the results pertinent to dipole technique and the results obtained using equations of accelerating and moving charges. This equivalence of the two techniques was previously shown for frequency domain fields in reference [10].
Observe that the field components given by Equations (24)–(31) depend on the speed of propagation of the pulse along the channel element, i.e., u . These field components will change their amplitudes if the speed of propagation is changed. For example, note that the velocity fields go to zero when the speed of propagation becomes equal to the speed of light in free space. However, the field components given by Equations (35)–(37), which were obtained by summing all the field terms given by Equations (24)–(31), are independent of the speed of propagation of the pulse. They depend only on the current waveform exciting the current element. This is the reason why it is not necessary to specify the speed of propagation of the current when writing down the electromagnetic fields of a short dipole. This also means that the field terms belonging to the radiation, velocity and electrostatics, which depend on the speed of propagation of the pulse, cancel out during the summation, leaving behind a total field which is independent of the speed of propagation. Thus, the speed of propagation enters into the electromagnetic field expressions of the channel element only if we wish to separate the total field into its physical constituents.

5. The Time Domain Fields of the Lightning Channel

The electromagnetic fields generated by the lightning channel (without the contribution of the ground) can be obtained by summing the contributions of all the channel elements located between z = 0 and z = H . The results are the following.

5.1. Radiation Fields

E r a d , θ ( t ) = 0 H sin θ   d z 4 π ε o c 2 r { I ( z , t ) t 2 u cos θ r ( 1 u c cos θ ) I ( z , t ) + u 2 sin 2 θ r c ( 1 u c cos θ ) 2 I ( z , t ) } a θ
E r a d , r ( t ) = 0 H I ( z , t ) d z 4 π ε o c 2 r 2 { u sin 2 θ ( 1 u c cos θ ) } a r
B r a d , ϕ ( t ) = 0 H d z sin θ 4 π ε o c 3 r { I ( z , t ) t 2 u cos θ r ( 1 u c cos θ ) I ( z , t ) + u 2 sin 2 θ r c ( 1 u c cos θ ) 2 I ( z , t ) } a φ

5.2. Velocity Fields

E v e l , θ ( t ) = 0 H I ( z , t ) d z sin θ 4 π ε o r 2 c [ 1 u c cos θ ] 2 [ 1 u 2 c 2 ] a θ
E v e l , r ( t ) = 0 H I ( z , t ) d z 4 π ε o r 2 u [ 1 u c cos θ ] [ 1 u 2 c 2 ] a r
B v e l , ϕ ( t ) = 0 H I ( z , t ) d z sin θ 4 π ε o r 2 c 2 [ 1 u c cos θ ] 2 [ 1 u 2 c 2 ] a φ

5.3. Static Fields

E s t a t , r ( t ) = 0 H d z 4 π ε o r 2 { cos θ c I ( z , t ) 1 u I ( z , t ) + 2 cos θ   r 0 t I ( z , τ ) d τ } a r
E s t a t , θ ( t ) = 0 H d z 4 π ε o r 2 { sin θ r 0 t I ( z , τ ) d τ } a θ

6. Electromagnetic Fields of the Lightning Channel over Perfectly Conducting Ground

The electric field at any point over perfectly conducting ground can be obtained by replacing the ground plane by an image channel. The relevant geometry is shown in Figure 2. The electric field at point P located at the surface of the perfectly conducting ground has only a z component and the corresponding radiation, velocity and static field terms are given by
E r a d , z ( t ) = 0 H   d z 2 π ε o c 2 r { I ( z , t ) sin 2 θ t 2 u sin 2 cos θ r ( 1 u c cos θ ) I ( z , t ) + u 2 sin 4 θ r c ( 1 u c cos θ ) 2 I ( z , t ) u cos θ sin 2 θ ( 1 u c cos θ ) I ( z , t ) } a z
E v e l ( t ) = 0 H I ( z , t ) d z 2 π ε o r 2 c [ 1 u c cos θ ] 2 { sin 2 θ c cos θ u } [ 1 u 2 c 2 ] a z
E s t a t ( t ) = 0 H d z 2 π ε o r 2 { cos 2 θ c I ( z , t ) + cos θ u I ( z , t ) { 3 sin 2 θ 2 }   r 0 t I ( z , τ ) d τ } a z
B r a d ( t ) = 0 H d z sin θ 2 π ε o c 3 r { I ( z , t ) t 2 u cos θ r ( 1 u c cos θ ) I ( z , t ) + u 2 sin 2 θ r c ( 1 u c cos θ ) 2 I ( z , t ) } a φ
B v e l ( t ) = 0 H I ( z , t ) d z sin θ 4 π ε o r 2 c 2 [ 1 u c cos θ ] 2 [ 1 u 2 c 2 ] a φ
These field expressions can be further simplified as
E r a d ( t ) = 0 H   d z 2 π ε o c 2 { I ( z , t ) sin 2 θ r   t + I ( z , t ) r 2 [ u 2 sin 4 θ c ( 1 u c cos θ ) 2 + u cos θ sin 2 θ { cos θ 2 } ( 1 u c cos θ ) ] } a z
E v e l ( t ) = 0 H I ( z , t ) d z 2 π ε o r 2 c [ 1 u c cos θ ] 2 [ 1 u 2 c 2 ] { cos θ sin 2 θ } a z
E s t a t ( t ) = 0 H d z 2 π ε o { cos θ u r 2 { 1 u c cos θ } I ( z , t ) { 3 sin 2 θ 2 }   r 3 0 t I ( z , τ ) d τ } a z
B r a d ( t ) = 0 H d z sin θ 2 π ε o c 3 { I ( z , t ) r   t u I ( z , t ) r 2 ( 1 u c cos θ ) { u sin θ c ( 1 u c cos θ ) 2 cos θ } } a φ
B v e l ( t ) = 0 H I ( z , t ) d z sin θ 4 π ε o r 2 c 2 [ 1 u c cos θ ] 2 [ 1 u 2 c 2 ] a φ
Based on the analysis presented in Section 4, these field terms can be reduced to dipole fields, resulting in
E r a d ( t ) = 0 H   d z 2 π ε o c 2 { I ( z , t ) sin 2 θ r   t + c { 3 sin 2 θ 2 } I ( z , t ) r 2 + c 2 { 3 sin 2 θ 2 }   r 3 0 t I ( z , τ ) d τ } a z
B r a d ( t ) = 0 H d z sin θ 2 π ε o c 3 { I ( z , t ) r   t + c I ( z , t ) r 2 } a φ
In order to illustrate this point further, let us consider the electric fields generated by a lightning return stroke as simulated using the modified transmission line model with exponential current attenuation (MTLE) [13,14]. The spatial and temporal distribution of the current, I ( z , t ) , associated with this model are given by
I ( z , t ) = 0 ,   if   t < z / u
I ( z , t ) = e z / λ I ( 0 , t ) ,   if   t > z / u
In the above equations, z is the height of the point of observation along the vertical channel, λ is the current decay height constant, u is the speed of propagation of the return stroke front and I ( 0 , t ) is the current at the channel base. The three field terms and the total field associated with the accelerating charge technique (Equations (51)–(53)) and the dipole technique (Equation (56)) calculated using the above current distribution are shown in Figure 3. In the calculations, the distance to the point of observation located at ground from the strike point is taken to be 5 km, λ = 2 km and u = 1.5 × 10 8 m/s. The channel base current with a peak 12 kA is identical to the one suggested by Nucci et al. [15] for a subsequent return stroke.
Note that the total field generated by both techniques is identical, whereas the three terms that build up this total field is different in the two techniques. For example, observe that the dipole-radiation term at this distance is bipolar, whereas the radiation field generated by accelerating charges is unipolar. As one can observe the rest of the terms are also different in the two techniques.

7. Discussion

In converting the field equations pertinent to pure radiation (fields produced by accelerating charges) and pure Coulomb (electric fields produced by stationary and uniformly moving charges) to those pertinent to dipoles, the connection of the field components to the processes which generate them is lost. Dipole fields have three field components, namely, dipole-static, dipole-induction and dipole-radiation. None of these separate field components completely describe either the radiation (fields generated by accelerating charges), velocity or static field components (Coulomb fields). Of course, the distant dipole-radiation field describes the energy radiated to infinity by the accelerating charges. However, this is not the only field component that is being produced by accelerating charges. The dipole-static term describes only part of the static field, whereas the dipole-induction field is a combination of static, velocity and radiation terms, and the dipole-radiation field describes only part of the radiation field. Note that the pure radiation field has terms both varying as 1 / r and 1 / r 2 .
It is important to point out that not only the 1 / r field terms carry energy from one place to another, but the field terms varying as 1 / r 2 terms do as well. The difference is that while the 1 / r terms carry energy to infinity, the 1 / r 2 terms carry energy into space in the vicinity of the dipole. As the current flows along the channel, charges are transported along it and these displaced charges generate electric fields in space. These electric fields (and also magnetic fields) carry energy and it is the 1 / r 2 terms that supply this energy. However, if the displaced charges reunite by coming back to the original location, these field components will transport this energy—located in space—back to the source. That means the flow of energy caused by the 1 / r 2 terms can change direction when the current flowing along the channel reverses its direction of flow or polarity. Indeed, in the case of oscillating dipoles, these terms carry energy back and forth from the source to space and vice-versa.
Finally, note that the total field calculated for a channel element or a lightning return stroke is the same irrespective of whether the equations pertinent to accelerating charges or the dipole equations are used. Even though this point has been demonstrated numerically in previous publications [3], it is demonstrated analytically in the present paper.

8. Conclusions

In this paper we have presented analytical expressions for the electromagnetic fields generated by return strokes in terms of the field components generated by accelerating charges, moving charges and stationary charges. It is shown analytically and demonstrated using return stroke fields that these three field components can be reduced to those of the field expressions pertinent to the dipole technique, but in the reduction process the association of the individual field components to the physical processes that generate the electromagnetic fields is lost.

Author Contributions

V.C. conceived the idea; both authors contributed equally to the analysis and writing of the paper.

Funding

This work was supported by Swedish Research Council grant VR-2015-05026 and partly by the fund from the B. John F. and Svea Andersson donation at Uppsala University.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Mathematical Steps Necessary to Convert (Equation (12) to Equations (13) and (14))

Let us start with Equation (12).
d e r a d = i ( z ) u 4 π ε o c 2 { e j ω ( t r 1 / c ) sin θ 1 r 1 [ 1 u cos θ 1 c ] a θ 1 e j ω ( t d z / u r 2 / c ) sin θ 2 r 2 [ 1 u cos θ 2 c ] a θ 2 }
This can be written as
d e r a d = i ( z ) u 4 π ε o c 2 { F 1 F 2 }
Let us start with the expression for F 1 .
F 1 = e j ω ( t r 1 / c ) sin θ 1 r 1 [ 1 u cos θ 1 c ] a θ 1
Since d z is an elementary length, the exponential term after substituting for r 1 from Equation (3) can be expanded as follows:
e j ω ( t r 1 / c ) = e j ω ( t r / c ) { 1 j ω d z cos θ 2 c }
Using this expression and substituting into 1 / r 1 , sin θ 1 and 1 / [ 1 u cos θ 1 / c ] from Equations (4), (6) and (10) respectively, we obtain
F 1 = e j ω ( t r / c ) { 1 j ω d z cos θ 2 c } sin θ { 1 d z cos θ 2 r } 1 r { 1 d z cos θ 2 r } { 1 + u d z sin 2 θ 2 r c ( 1 u c cos θ ) } [ 1 u cos θ c ] a θ 1
Reorganizing the terms, we can write
F 1 = e j ω ( t r / c ) sin θ { 1 j ω d z cos θ 2 c } { 1 d z cos θ 2 r } { 1 d z cos θ 2 r } { 1 + u d z sin 2 θ 2 r c ( 1 u c cos θ ) } r [ 1 u cos θ c ] a θ 1
Keeping only the first order terms with respect to d z , we obtain
F 1 = e j ω ( t r / c ) sin θ { 1 j ω d z cos θ 2 c d z cos θ r + u d z sin 2 θ 2 r c ( 1 u c cos θ ) } r [ 1 u cos θ c ] a θ 1
This vector can be resolved into the directions of a θ and a r using Equations (1) and (2) (note that sin δ θ 1 δ θ 1 ) as follows:
F 11 = e j ω ( t r / c ) sin θ { 1 j ω d z cos θ 2 c d z cos θ r + u d z sin 2 θ 2 r c ( 1 u c cos θ ) } r [ 1 u cos θ c ] a θ
F 12 = e j ω ( t r / c ) sin θ { 1 j ω d z cos θ 2 c d z cos θ r + u d z sin 2 θ 2 r c ( 1 u c cos θ ) } ( d z sin θ 2 r ) r [ 1 u cos θ c ] a r
Again, keeping only the first order terms with respect to d z we obtain
F 12 = e j ω ( t r / c ) sin 2 θ 2 r [ 1 u cos θ c ] a r
Let us now consider the term F 2 .
F 2 = e j ω ( t d z / u r 2 / c ) sin θ 2 r 2 [ 1 u cos θ 2 c ] a θ 2
Let us consider the exponential term. Substituting for r 2 from Equation (3) and expanding the exponential term (since d z is an elementary length) we obtain
e j ω ( t r 2 / c d z / u ) = e j ω ( t r / c + d z cos θ c d z u ) = e j ω ( t r / c ) { 1 j ω d z u + j ω d z cos θ 2 c }
This can be written as
e j ω ( t r 2 / c d z / u ) = e j ω ( t r / c ) { 1 j ω d z u ( 1 u 2 c cos θ ) }
Using this expression and substituting to 1 / r 2 , sin θ 2 and 1 / [ 1 u cos θ 2 / c ] from Equations (4), (6) and (11) respectively, we obtain
F 2 = e j ω ( t r / c ) { 1 j ω d z u ( 1 u 2 c cos θ ) } sin θ { 1 + d z cos θ 2 r } { 1 u d z sin 2 θ 2 r c ( 1 u c cos θ ) } { 1 + d z cos θ 2 r } r { 1 u c cos θ } a θ 2
Keeping only the first order terms in d z we obtain
F 2 = e j ω ( t r / c ) sin θ { 1 j ω d z u ( 1 u 2 c cos θ ) + d z cos θ r u d z sin 2 θ 2 r c ( 1 u c cos θ ) } r [ 1 u cos θ c ] a θ 2
This vector can be resolved into the directions of a θ and a r using Equations (1) and (2) (note that sin δ θ 2 δ θ 2 ) as follows:
F 21 = e j ω ( t r / c ) sin θ { 1 j ω d z u ( 1 u 2 c cos θ ) + d z cos θ r u d z sin 2 θ 2 r c ( 1 u c cos θ ) } r [ 1 u cos θ c ] a θ
F 22 = e j ω ( t r / c ) sin θ { 1 j ω d z u ( 1 u 2 c cos θ ) + d z cos θ r u d z sin 2 θ 2 r c ( 1 u c cos θ ) } ( d z sin θ 2 r ) r [ 1 u cos θ c ] a θ
Keeping only the first order terms in d z we obtain
F 22 = e j ω ( t r / c ) sin 2 θ d z 2 r [ 1 u cos θ c ] a r
Now, the net field component in the a θ direction is given by
d e r a d , θ = i ( z ) u 4 π ε o c 2 { F 11 F 21 }
After substituting for F 11 and F 2 1 and replacing t r / c by t the expression reduces to Equation (13), i.e.,
d e r a d , θ = i ( z ) e j ω t u sin θ 4 π ε o c 2 r ( 1 u c cos θ ) { j ω d z u j ω d z cos θ c 2 d z cos θ r + u d z sin 2 θ r c ( 1 u c cos θ ) } a θ
The net component in the a r direction is given by
d e r a d , r = i ( z ) u 4 π ε o c 2 { F 12 F 22 }
After substituting for F 1 2 and F 22 the expression reduces to Equation (14), i.e.,
d e r a d , r = i ( z ) e j ω t 4 π ε o c 2 r 2 { u d z sin 2 θ ( 1 u c cos θ ) } a r

Appendix B. Mathematical Steps Necessary for the Reduction of Electromagnetic Fields Obtained from Accelerating and Moving Charges to the Dipole Fields

(a) Reduction of d E θ to the dipole fields:
The mathematical steps necessary for the reduction starting from Equation (32) are shown below.
d E θ ( t ) = d z sin θ 4 π ε o { 1 c 2 r i ( z , t ) t + 1 r 3 0 t I ( z , τ ) d τ + i ( z , t )   u 2 sin 2 θ r 2 c 3 ( 1 u c cos θ ) 2 i ( z , t ) 2 u cos θ c 2 r 2 ( 1 u c cos θ ) + i ( z , t ) ( 1 u 2 / c 2 ) c r 2 ( 1 u c cos θ ) 2 } a θ
Let us combine the third and fifth terms inside the bracket. The result is given by
d E θ ( t ) = d z sin θ 4 π ε o { 1 c 2 r i ( z , t ) t + 1 r 3 0 t I ( z , τ ) d τ + i ( z , t )   r 2 c ( 1 u c cos θ ) 2 [ u 2 sin 2 θ c 2 + 1 u 2 c 2 ] i ( z , t ) 2 u cos θ c 2 r 2 ( 1 u c cos θ ) } a θ
d E θ ( t ) = d z sin θ 4 π ε o { 1 c 2 r i ( z , t ) t + 1 r 3 0 t I ( z , τ ) d τ + i ( z , t )   r 2 c ( 1 u c cos θ ) 2 [ u 2 c 2 u 2 cos 2 θ c 2 + 1 u 2 c 2 ] i ( z , t ) 2 u cos θ c 2 r 2 ( 1 u c cos θ ) } a θ
d E θ ( t ) = d z sin θ 4 π ε o { 1 c 2 r i ( z , t ) t + 1 r 3 0 t I ( z , τ ) d τ + i ( z , t )   r 2 c ( 1 u c cos θ ) 2 [ 1 u 2 cos 2 θ c 2 ] i ( z , t ) 2 u cos θ c 2 r 2 ( 1 u c cos θ ) } a θ
d E θ ( t ) = d z sin θ 4 π ε o { 1 c 2 r i ( z , t ) t + 1 r 3 0 t I ( z , τ ) d τ + i ( z , t )   r 2 c ( 1 u c cos θ ) 2 [ ( 1 u c cos θ ) ( 1 + u c cos θ ) ] i ( z , t ) 2 u cos θ c 2 r 2 ( 1 u c cos θ ) } a θ
d E θ ( t ) = d z sin θ 4 π ε o { 1 c 2 r i ( z , t ) t + 1 r 3 0 t I ( z , τ ) d τ + i ( z , t ) ( 1 + u c cos θ )   r 2 c ( 1 u c cos θ ) i ( z , t ) 2 u cos θ c 2 r 2 ( 1 u c cos θ ) } a θ
d E θ ( t ) = d z sin θ 4 π ε o { 1 c 2 r i ( z , t ) t + 1 r 3 0 t I ( z , τ ) d τ + i ( z , t )   r 2 c ( 1 u c cos θ ) [ ( 1 + u c cos θ ) 2 u cos θ c ] } a θ
d E θ ( t ) = d z sin θ 4 π ε o { 1 c 2 r i ( z , t ) t + 1 r 3 0 t I ( z , τ ) d τ + i ( z , t )   r 2 c ( 1 u c cos θ ) [ ( 1 u c cos θ ) ] } a θ
d E θ ( t ) = d z sin θ 4 π ε o { 1 c 2 r i ( z , t ) t + 1 r 3 0 t I ( z , τ ) d τ + i ( z , t )   r 2 c } a θ
(b) Reduction of d E r to the dipole fields:
The mathematical steps necessary for the reduction starting from Equation (33) are shown below.
d E r ( t ) = d z 4 π ε o { 2 cos θ   r 3 0 t i ( z , τ ) d τ + i ( z , t ) cos θ r 2 c i ( z , t ) u r 2 + i ( z , t ) u sin 2 θ c 2 r 2 ( 1 u c cos θ ) + i ( z , t ) ( 1 u 2 / c 2 ) u r 2 ( 1 u c cos θ ) } a r
d E r ( t ) = d z 4 π ε o { 2 cos θ   r 3 0 t i ( z , τ ) d τ + i ( z , t ) u r 2 [ u cos θ c 1 + u 2 sin 2 θ c 2 ( 1 u c cos θ ) + ( 1 u 2 / c 2 ) ( 1 u c cos θ ) ] } a r
d E r ( t ) = d z 4 π ε o { 2 cos θ   r 3 0 t i ( z , τ ) d τ + i ( z , t ) u r 2 [ u cos θ c 1 + 1 c 2 ( 1 u c cos θ ) ( u 2 sin 2 θ + c 2 u 2 ) ] } a r
d E r ( t ) = d z 4 π ε o { 2 cos θ   r 3 0 t i ( z , τ ) d τ + i ( z , t ) u r 2 [ u cos θ c 1 + 1 c 2 ( 1 u c cos θ ) ( c 2 u 2 cos 2 θ ) ] } a r
d E r ( t ) = d z 4 π ε o { 2 cos θ   r 3 0 t i ( z , τ ) d τ + i ( z , t ) u r 2 [ u cos θ c 1 + 1 ( 1 u c cos θ ) ( 1 u 2 c 2 cos 2 θ ) ] } a r
d E r ( t ) = d z 4 π ε o { 2 cos θ   r 3 0 t i ( z , τ ) d τ + i ( z , t ) u r 2 [ u cos θ c 1 + ( 1 + u cos θ c ) ] } a r
d E r ( t ) = d z 4 π ε o { 2 cos θ   r 3 0 t i ( z , τ ) d τ + i ( z , t ) r 2 [ 2 cos θ c ] } a r
(c) Reduction of d B ϕ to the dipole fields:
The mathematical steps necessary for the reduction starting from Equation (34) are shown below.
d B ϕ ( t ) = d z sin θ 4 π ε o { 1 c 3 r i ( z , t ) t i ( z , t ) 2 u cos θ c 3 r 2 ( 1 u c cos θ ) + i ( z , t ) u 2 sin 2 θ r 2 c 4 ( 1 u c cos θ ) 2 + i ( z , t ) ( 1 u 2 / c 2 ) c 2 r 2 ( 1 u c cos θ ) 2 } a ϕ
d B ϕ ( t ) = d z sin θ 4 π ε o { 1 c 3 r i ( z , t ) t i ( z , t ) 2 u cos θ c 3 r 2 ( 1 u c cos θ ) + i ( z , t ) r 2 c 2 ( 1 u c cos θ ) 2 ( u 2 c 2 sin 2 θ + ( 1 u 2 c 2 ) ) } a ϕ
d B ϕ ( t ) = d z sin θ 4 π ε o { 1 c 3 r i ( z , t ) t i ( z , t ) 2 u cos θ c 3 r 2 ( 1 u c cos θ ) + i ( z , t ) r 2 c 2 ( 1 u c cos θ ) 2 ( 1 u 2 c 2 cos 2 θ ) } a ϕ
d B ϕ ( t ) = d z sin θ 4 π ε o { 1 c 3 r i ( z , t ) t i ( z , t ) 2 u cos θ c 3 r 2 ( 1 u c cos θ ) + i ( z , t ) r 2 c 2 ( 1 u c cos θ ) ( 1 + u c cos θ ) } a ϕ
d B ϕ ( t ) = d z sin θ 4 π ε o { 1 c 3 r i ( z , t ) t + i ( z , t ) c 2 r 2 ( 1 u c cos θ ) ( 1 + u c cos θ 2 u c cos θ ) } a ϕ
d B ϕ ( t ) = d z sin θ 4 π ε o { 1 c 3 r i ( z , t ) t + i ( z , t ) c 2 r 2 } a ϕ

References

  1. Rachidi, F. A review of field-to-transmission line coupling models with special emphasis to lightning-induced voltages. IEEE Trans. Electromagn. Compat. 2012, 54, 898–911. [Google Scholar] [CrossRef]
  2. Nucci, C.A.; Rachidi, F.; Rubinstein, M. Interaction of Lightning-Generated Electromagnetic Fields with Overhead and Underground Cables. In Lightning Electromagnetics; Cooray, V., Ed.; IET: London, UK, 2012; pp. 678–718. [Google Scholar]
  3. Baba, Y.; Rakov, V. Voltages Induced on an Overhead Wire by Lightning Strikes to a Nearby Tall Grounded Object. IEEE Trans. Electromagn. Compat. 2006, 48, 212–224. [Google Scholar] [CrossRef] [Green Version]
  4. Ye, M.; Cooray, V. Propagation effects caused by a rough ocean surface on the electromagnetic fields generated by lightning return strokes. Radio Sci. 1994, 29, 73–85. [Google Scholar]
  5. Cooray, V.; Ye, M. Propagation effects on the lightning generated electromagnetic fields for homogeneous and mixed sea land paths. J. Geophys. Res. 1994, 99, 10641–10652. [Google Scholar] [CrossRef]
  6. Thottappillil, R.; Rakov, V.A. On different approaches to calculating lightning electric fields. J. Geophys. Res. 2001, 106, 14191–14205. [Google Scholar] [CrossRef] [Green Version]
  7. Thottappillil, R. Computation of electromagnetic fields from lightning discharges. In The Lightning Flash, 2nd ed.; Cooray, V., Ed.; IET: London, UK, 2014. [Google Scholar]
  8. Cooray, V.; Cooray, G. The electromagnetic fields of an accelerating charge: Applications in lightning return stroke models. IEEE Trans. Electromagn. Compat. 2010, 52, 944–955. [Google Scholar] [CrossRef]
  9. McLain, D.K.; Uman, M.A. Exact expression and moment approximation for the electric field intensity of the lightning return stroke. J. Geophys. Res. 1971, 76, 2101–2105. [Google Scholar] [CrossRef]
  10. Cooray, G.; Cooray, V. Electromagnetic fields of a short electric dipole in free space—Revisited. Prog. Electromagn. Res. 2012, 131, 357–373. [Google Scholar] [CrossRef]
  11. Cooray, V. Application of electromagnetic fields of an accelerating charge to obtain the electromagnetic fields of a propagating current pulse. In The Lightning Electromagnetics; Cooray, V., Ed.; IET: London, UK, 2012. [Google Scholar]
  12. Pannofsky, W.K.H.; Phillips, M. Classical Electricity and Magnetism; Addison-Wesley: Reading, MA, USA, 1962. [Google Scholar]
  13. Nucci, C.A.; Mazzetti, C.; Rachidi, F.; Ianoz, M. On lightning return stroke models for LEMP calculations. In Proceedings of the 19th International Conference on Lightning Protection, Graz, Austria, 25–29 April 1988. [Google Scholar]
  14. Rachidi, F.; Nucci, C.A. On the Master, Lin, Uman, Standler and the Modified Transmission Line lightning return stroke current models. J. Geophys. Res. 1990, 95, 20389–20394. [Google Scholar] [CrossRef]
  15. Nucci, C.A.; Diendorfer, G.; Uman, M.A.; Rachidi, F.; Ianoz, M.; Mazetti, C. Lightning return stroke models with specified channel base current: A review and comparison. J. Geophys. Res. 1990, 95, 20395–20408. [Google Scholar] [CrossRef]
Figure 1. Geometry, angles and unit vectors pertinent to the evaluation of electromagnetic fields generated by a channel element. The unit vector in the direction of positive z-axis is denoted by a z . The unit vectors in the radial directions r , r 1 and r 2 are denoted by a r , a r 1 and a r 2 , respectively. The unit vectors a θ , a θ 1 and a θ 2 are defined as a r × ( a r × a z ) , a r 1 × ( a r 1 × a z ) and a r 2 × ( a r 2 × a z ) , respectively. The unit vector a ϕ is in the direction of the vector a r × a θ (i.e., into the page). Note that point P can be located anywhere in space.
Figure 1. Geometry, angles and unit vectors pertinent to the evaluation of electromagnetic fields generated by a channel element. The unit vector in the direction of positive z-axis is denoted by a z . The unit vectors in the radial directions r , r 1 and r 2 are denoted by a r , a r 1 and a r 2 , respectively. The unit vectors a θ , a θ 1 and a θ 2 are defined as a r × ( a r × a z ) , a r 1 × ( a r 1 × a z ) and a r 2 × ( a r 2 × a z ) , respectively. The unit vector a ϕ is in the direction of the vector a r × a θ (i.e., into the page). Note that point P can be located anywhere in space.
Atmosphere 10 00022 g001
Figure 2. Geometry relevant to the calculation of electromagnetic fields from a return stroke.
Figure 2. Geometry relevant to the calculation of electromagnetic fields from a return stroke.
Atmosphere 10 00022 g002
Figure 3. (a) The three field components and the total field associated with accelerating and moving charges, and (b) the three field components and the total field associated with dipole fields. The electric fields at 5 km distance from a lightning channel are obtained using the MTLE model with a 12 kA current, 1.5 × 108 m/s return stroke speed and 2 km current decay height constant. Note that in depicting the electric fields, the field components directed into the ground (i.e., directed along the negative z-axis) are considered positive.
Figure 3. (a) The three field components and the total field associated with accelerating and moving charges, and (b) the three field components and the total field associated with dipole fields. The electric fields at 5 km distance from a lightning channel are obtained using the MTLE model with a 12 kA current, 1.5 × 108 m/s return stroke speed and 2 km current decay height constant. Note that in depicting the electric fields, the field components directed into the ground (i.e., directed along the negative z-axis) are considered positive.
Atmosphere 10 00022 g003

Share and Cite

MDPI and ACS Style

Cooray, V.; Cooray, G. A Novel Interpretation of the Electromagnetic Fields of Lightning Return Strokes. Atmosphere 2019, 10, 22. https://doi.org/10.3390/atmos10010022

AMA Style

Cooray V, Cooray G. A Novel Interpretation of the Electromagnetic Fields of Lightning Return Strokes. Atmosphere. 2019; 10(1):22. https://doi.org/10.3390/atmos10010022

Chicago/Turabian Style

Cooray, Vernon, and Gerald Cooray. 2019. "A Novel Interpretation of the Electromagnetic Fields of Lightning Return Strokes" Atmosphere 10, no. 1: 22. https://doi.org/10.3390/atmos10010022

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop