# An Algorithm for Evaluation of Lightning Electromagnetic Fields at Different Distances with respect to Lightning Channel

Authors: Mahdi Izadi; Mohd Zainal Abidin Ab Kadir; Maryam Hajikhani
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

The analytical field expressions are proposed to estimate the electromagnetic fields associated with vertical lightning channel at an observation point directly in the time domain using current and geometrical parameters. The proposed expressions embrace a large number of channel base current functions and widely used engineering current models. Moreover, they can be joined with different coupling models for evaluation of lightning induced voltage when they provide different required field components directly in the time domain. The data from measured fields and simulated fields from previous works were employed to validate the proposed field expressions. The simulated results were compatible with both computed field values from previous studies and the measured fields.


---

## Body

## 1. Introduction

Lightning is a natural phenomenon that can give effect on the power lines in two ways, that is, direct and indirect. In the direct effects, lightning strikes to power lines or towers or substations directly. On the other hand, the indirect effect considers electromagnetic fields associated with lightning channel and induced voltage due to coupling between generated electromagnetic fields and power system when lightning strikes to the ground or any objective near to power lines or substations [1, 2]. This study considers the evaluation of electromagnetic fields due to lightning channel. Several methods consider the calculation of electromagnetic fields due to lightning channel while they are validated using measured fields at one or a number of observation points [1] while the most common calculation methods can be concluded by(i)
the MONOPOLE method [3, 4],(ii)
the DIPOLE method [1, 3],(iii)
the FDTD method [5],(iv)
the FDTD-HYBRID methods [6, 7].It should be mentioned that the electromagnetic fields associated with the lightning channel rely basically on some geometrical, channel, and ground parameters. The basic input parameters in almost all calculation methods in perfect ground conductivity case are(i)
the channel base current that is usually prepared with direct measurement [2] from the triggered lightning technique or obtained from inverse procedure algorithms from measured electromagnetic fields [8, 9],(ii)
the current model for prediction of return stroke current wave shape along lightning channel,(iii)
the radial distance from lightning channel,(iv)
the observation point height with respect to ground surface,(v)
the return stroke current velocity along lightning channel.Two first calculation methods are expressed in the cylindrical [1] and spherical [3, 4] domains using current and charge density [4] parameters along the lightning channel while they have complexity in the calculation of induction partial when some realistic channel base current functions are used (as Heidler function) and they do not have an analytical solution of integral to time [7, 10]. On the other hand, the FDTD method just is validated for closed distances from lightning channel [5]. Moreover, the FDTD_HYBRID methods need to estimate magnetic fields at six points around observation point for calculation of electric fields at an observation point [6, 7, 11]. In order to provide the general analytical electromagnetic field expressions in the time domain, the widely used channel base current functions are introduced and categorized. Therefore, after generalization of current functions based on proposed expressions, the general electromagnetic field expressions due to vertical lightning channel on the prefect ground are proposed. These equations can be applied to a large number of leading current functions and their combinations and can support different generalized engineering current models directly in the time domain. The basic assumptions of this study are the following.(i)
The ground conductivity is infinite.(ii)
The ground surface is flat.(iii)
The lightning channel is perpendicular to the ground and has no branches.(iv)
The corona effect is negligible.In the next section, the widely used current functions and current models are reviewed. Therefore, the proposed electromagnetic field expressions are presented. Finally, a justification of the proposed method and its validation by measuring fields and simulated fields from previous studies is given.


## 2. Channel Base Current Functions

Previous studies presented several channel base current functions to simulate return stroke current wave shape at channel base on the ground surface where the fixed coefficients of the current functions were determined by the data from measured current or measured fields using inverse procedure algorithms. The significant channel base current functions are tabulated in Table1 where A, A1, A2, B, B1, B2, Γ1, Γ2, Γ11, Γ12, Γ21, Γ22, η, n1, n2 are constant coefficients and Ip, Ip1, Ip2, I01, I02 are current peaks. Note that similar symbols in different current functions do not have a similar concept.Table 1
The widely used channel base current functions.
Current model
Channel base current function
Bruce and Golde [29]
I
p
[
exp
⁡
(
-
A
t
)
-
exp
⁡
(
-
B
t
)
]
Improvement of Uman and McLain on Bruce and Golde function [30]
I
p
η
[
exp
⁡
(
-
A
t
)
-
exp
⁡
(
-
B
t
)
]
Improvement of Jones on Bruce and Golde function [31]
I
p
η
[
exp
⁡
(
-
t
*
Γ
1
)
-
exp
⁡
(
-
t
*
Γ
2
)
]
t
*
=
Γ
2
2
Γ
1
+
t
Pierce and Ciones [32]
I
p
1
[
exp
⁡
(
-
A
1
t
)
-
exp
⁡
(
-
B
1
t
)
]
+
I
p
2
[
exp
⁡
(
-
A
2
t
)
-
exp
⁡
(
-
B
2
t
)
]
Heidler function [33]
I
p
η
1
(
t
/
Γ
1
)
n
1
1
+
(
t
/
Γ
1
)
n
1
exp
⁡
⁡
(
-
t
Γ
2
)
Sum of two Heidler functions [11, 34]
i
01
η
1
(
t
/
Γ
11
)
n
1
1
+
(
t
/
Γ
11
)
n
1
exp
⁡
⁡
(
-
t
Γ
12
)
+
i
02
η
2
(
t
/
Γ
21
)
n
2
1
+
(
t
/
Γ
21
)
n
2
exp
⁡
⁡
(
-
t
Γ
22
)
Improvement of Nucci on Heidler function [35]
i
01
η
1
(
t
/
Γ
11
)
n
1
1
+
(
t
/
Γ
11
)
n
1
exp
⁡
(
-
t
Γ
12
)
+
i
02
[
exp
⁡
(
-
Γ
21
t
)
-
exp
⁡
(
-
Γ
22
t
)
]Different channel base current functions in Table1 can be determined from (1) and (2) by changing the constant coefficients or the combination of (1) and (2). Hence, this study adapts these two basic equations and considers the widely used current functions using a combination of them:
(1)I0,t=I0η0c1e-αt-c2e-γt,(2)I0,t=I01η01t/t1n1+t/t1nexp⁡⁡-tt2,
where
(3)η01=exp⁡⁡-t1t2nt2t11/n.


## 3. Return Stroke Current Models

One of the effective factors on electromagnetic field calculations is the behavior of return stroke currents at different heights along the lightning channel that is modeled via current models. Generally, the return stroke current models can be classified into four groups as follows [12, 13]:(i)
the gas-dynamic or physical models,(ii)
the current distributed models,(iii)
the electromagnetic models,(iv)
the engineering models.To develop general expressions for electromagnetic fields due to lightning channel, the engineering models are applied in this study. In the engineering models, the current at different heights along the lightning channel can be expressed by a special function based on the channel base current and an attenuation height dependent factor. Most of the engineering models can be presented by a general equation such as (4) [12, 13]. The internal parameters of (4) for different engineering current models are listed in Table 2 [1, 12] where c is considered as the speed of light in free space, λ as a decay constant (1~2 km), and Hc as the cloud height with respect to the ground surface [1]:
(4)Iz′,t=I0,t-z′v×Pz′×ut-z′vf,
where z′ is the temporary charge height along lightning channel, I(z′,t) is current distribution along the lightning channel at any height z′ and any time t, I(0,t) is channel base current, P(z′) is the attenuation height dependent factor, v is the current-wave propagation velocity, vf is the upward propagating front velocity, and u is the Heaviside function defended as
(5)ut-z′vf=1fort≥z′vf0fort<z′vf.Table 2
The internal parameters of (4) for different engineering models  [12].
Model
P
(
z
′
)
v
Bruce and Golde (BG) [29]
1
∞
Traveling Current Source (TCS) [36]
1
-
c
Transmission Line (TL) [1]
1
v
f
Modified Transmission Line with Linear Decay (MTLL) [1]
(
1
-
z
′
H
c
)
v
f
Modified Transmission Line with Exponential Decay (MTLE) [1]
exp
⁡
(
-
z
′
λ
)
v
fBased on (4), return stroke current at the height of z′ along the lightning channel can be estimated by using channel base current and the height dependent attenuation factor. Note that, in recent years, several models based on (4) and using measured electromagnetic fields and inverse procedure algorithms were presented [14–16]. Similarly, different current models can be examined by the comparison between simulated electromagnetic fields using current model and measured fields at different observation points [1]. In this model the assumed return stroke velocity is set at a constant value that is typically between c/2 and 2c/3 [17].


## 4. Electromagnetic Fields Associated with Lightning Channel

In this section, electromagnetic field expressions due to two expressed functions in (1) and (2) at an observation point above the ground surface are estimated whereas the current model is set on (4). According to the geometry of the situation (see Figure 1), the potential vector (A→) at the observation point above ground surface can be obtained from (6) when Maxwell’s equation [3, 18] and Lorentz gauge [19] are applied:
(6)∇2A→-μ0ε0∂2A→∂t2=-μ0J→,
where J→ is the current density, μ0=4π×10-7 V·s/(A·m), and ε0=8.85×10-12 F/m.Figure 1
Geometry of the problem.Moreover, Figure1 demonstrates that the solution to (6) can be obtained from (7) whereas the infinitesimal current source is located along z-axis at space in position r2→ from the origin. Also, the observation point is located at space and r1→ from the origin point; therefore, R→ can be estimated by R→=r1→-r2→:
(7)dA→=μ0[I→]dL′4πr1→-r2→,
where I→ is line current density and dL′ is unit vector along lightning channel.Likewise, by replacing return stroke current anddz′ with I→ and unit vector, respectively, in (7) can convert it to (8)(8)dAz→=μ0i(z′,t-R/c)dz′4πR,
where
(9)r=x2+y2,R=x2+y2+z-z′2.Further, the magnetic flux density and electric fields can be obtained from potential vectorA→ as expressed in (10) and (11), respectively [3, 18]:
(10)B→=∇×A→,(11)∇×E→=-∂B→∂t=-∂∇×A→∂t,
where B→ is the magnetic flux density and E→ is the electric field vector.Therefore, based on the general current model in (4), the field components for first current function as presented in (1) can be obtained from (12), when (8), (10), (11), the trapezoid, and FDTD methods are applied [20–22]. Note that the internal parameters are collected from Table 3. Also, B→z_RS1x,y,z,tn=0:
(12)B→x_RS1x,y,z,tn=∑i=1n∑m=1k+1amF1x,y,z,t=tn,z′=hm,ihhhhhhhhhih-am′F1x,y,z,t=tn,z′=hm,i′,B→y_RS1x,y,z,,tn=∑i=1n∑m=1k+1-xyamF1x,y,z,t=tn,z′=hm,ihhhhhhhhhihhhhhih-am′F1x,y,z,t=tn,z′=hm,i′,E→x_RS1x,y,z,tn=E→x_RS1x,y,z,tn-1+Δt×∑i=1n∑m=1k+1amF2x,y,z,t=tn,z′=hm,ihhhhhhhhiihhh-am′F2x,y,z,t=tn,z′=hm,i′,E→y_RS1x,y,z,tn=E→y_RS1x,y,z,tn-1+Δt×∑i=1n∑m=1k+1(yx)amF2x,y,z,t=tn,z′=hm,ihhhhhhhhhhhihhhhih-am′F2x,y,z,t=tn,z′=hm,i′,E→z_RS1x,y,z,tn=E→z_RS1x,y,z,tn-1+Δt×∑i=1n∑m=1k+1amF3x,y,z,t=tn,z′=hm,ihhhhhhhhhihhh-am′F3x,y,z,t=tn,z′=hm,i′,
where B→x_RS1 is the magnetic flux density at x-axis due to return stroke current function from (1), B→y_RS1 is the magnetic flux density at y-axis due to return stroke current function from (1), B→z_RS1 is the magnetic flux density at z-axis due to return stroke current function from (1), E→x_RS1 is the electric field at x-axis due to return stroke current function from (1), E→y_RS1 is the electric field at y-axis due to return stroke current function from (1), and E→z_RS1 is the electric field at z-axis due to return stroke current function from (1). Δt is the time step, n is the number of time steps,
(13)tn=r2+z2c+n-1Δtn=1,2,…,nmax⁡am=Δhi2×kform=1,m=k+1Δhikforothersam′=Δhi′2×kform=1,m=k+1Δhi′kforothers,k is division factor (≥2),(14)Δhi=βχ2cti-cti-1-βcti-z2+rχ2+βcti-1-z2+rχ2βχ2-βz-cti-βcti-z2+rχ2fori=1Δhi′=βχ2cti-1-cti+βcti+z2+rχ2-βcti-1+z2+rχ2βχ2-βz+cti+βcti+z2+rχ2fori=1hm,i=m-1×Δhik+hm=k+1,i-1m-1×Δhikfori=1hm,i′=m-1×Δhi′k+hm=k+1,i-1′m-1×Δhi′kfori=1.Table 3
The internal parameters of proposed electromagnetic fields expressions.
Parameter
Expression
F
1
10
-
7
I
0
y
P
(
z
′
)
η
0
{
[
c
1
α
exp
⁡
(
-
α
A
1
)
-
c
2
γ
exp
⁡
(
-
γ
A
1
)
c
R
2
]
-
[
c
1
exp
⁡
(
-
α
A
1
)
-
c
2
exp
⁡
(
-
γ
A
1
)
R
3
]
}
F
2
I
0
P
(
z
′
)
4
π
ε
0
η
0
(
z
-
z
′
)
x
{
-
3
(
c
1
α
exp
⁡
(
-
α
A
1
)
-
c
2
γ
exp
⁡
(
-
γ
A
1
)
)
c
R
4
+
c
1
α
2
exp
⁡
(
-
α
A
1
)
-
c
2
γ
2
exp
⁡
(
-
γ
A
1
)
c
2
R
3
+
3
(
c
1
exp
⁡
(
-
α
A
1
)
-
c
2
exp
⁡
(
-
γ
A
1
)
)
R
5
}
F
3
I
0
P
(
z
′
)
4
π
ε
0
η
0
{
3
(
x
2
+
y
2
)
(
c
1
α
exp
⁡
(
-
α
A
1
)
-
c
2
γ
exp
⁡
(
-
γ
A
1
)
)
c
R
4
+
2
(
-
c
1
α
exp
⁡
(
-
α
A
1
)
+
c
2
γ
exp
⁡
(
-
γ
A
1
)
)
c
R
2
+
(
x
2
+
y
2
)
(
-
c
1
α
2
exp
⁡
(
-
α
A
1
)
+
c
2
γ
2
exp
⁡
(
-
γ
A
1
)
)
c
2
R
3
-
3
(
x
2
+
y
2
)
(
c
1
exp
⁡
(
-
α
A
1
)
-
c
2
exp
⁡
(
-
γ
A
1
)
)
R
5
+
2
(
c
1
exp
⁡
(
-
α
A
1
)
-
c
2
exp
⁡
(
-
γ
A
1
)
)
R
3
}
F
4
10
-
7
I
01
P
(
z
′
)
y
A
2
η
01
[
(
A
1
/
t
1
)
n
+
1
]
{
(
A
1
/
t
1
)
n
c
t
2
R
2
-
(
A
1
/
t
1
)
n
R
3
-
n
(
A
1
/
t
1
)
n
-
1
c
t
1
R
2
+
n
(
A
1
/
t
1
)
2
n
-
1
c
t
1
R
2
[
(
A
1
/
t
1
)
n
+
1
]
}
F
5
I
01
P
(
z
′
)
x
(
z
-
z
′
)
A
2
4
π
ε
0
η
01
[
(
A
1
/
t
1
)
n
+
1
]
{
(
A
1
/
t
1
)
n
(
3
c
2
t
2
2
-
3
c
t
2
R
+
R
2
)
c
2
t
2
2
R
5
+
2
n
2
(
A
1
/
t
1
)
3
n
-
2
+
n
(
n
-
1
)
(
A
1
/
t
1
)
n
-
2
[
(
A
1
/
t
1
)
n
+
1
]
2
+
(
-
3
n
2
+
n
)
(
A
1
/
t
1
)
2
n
-
2
[
(
A
1
/
t
1
)
n
+
1
]
c
2
t
1
2
R
3
[
(
A
1
/
t
1
)
n
+
1
]
2
+
3
n
(
[
(
A
1
/
t
1
)
n
+
1
]
(
A
1
/
t
1
)
n
-
1
-
(
A
1
/
t
1
)
2
n
-
1
)
c
t
1
R
4
[
(
A
1
/
t
1
)
n
+
1
]
+
2
n
(
(
A
1
/
t
1
)
2
n
-
1
-
[
(
A
1
/
t
1
)
n
+
1
]
(
A
1
/
t
1
)
n
-
1
)
c
2
t
1
t
2
R
3
[
(
A
1
/
t
1
)
n
+
1
]
}
F
6
I
01
P
(
z
′
)
A
2
4
π
ε
0
η
01
[
(
A
1
/
t
1
)
n
+
1
]
{
(
x
2
+
y
2
)
[
(
A
1
/
t
1
)
n
(
-
3
c
2
t
1
t
2
2
+
3
c
t
1
t
2
R
-
t
1
R
2
)
+
n
(
A
1
/
t
1
)
n
-
1
(
-
3
c
t
2
2
R
+
2
t
2
R
2
)
]
c
2
t
1
t
2
2
R
5
+
2
(
A
1
/
t
1
)
n
(
c
t
1
t
2
-
t
1
R
)
+
2
n
(
A
1
/
t
1
)
n
-
1
t
2
R
c
t
1
t
2
R
3
-
n
(
n
-
1
)
(
x
2
+
y
2
)
(
A
1
/
t
1
)
n
-
2
c
2
t
1
2
R
3
+
(
x
2
+
y
2
)
[
2
n
2
t
2
R
(
A
1
/
t
1
)
2
n
-
2
+
3
n
c
t
1
t
2
(
A
1
/
t
1
)
2
n
-
1
+
n
(
n
-
1
)
t
2
R
(
A
1
/
t
1
)
2
n
-
2
-
2
n
t
1
R
(
A
1
/
t
1
)
2
n
-
1
]
-
2
n
c
t
1
t
2
R
2
(
A
1
/
t
1
)
2
n
-
1
c
2
t
1
2
t
2
R
4
[
(
A
1
/
t
1
)
n
+
1
]
-
2
n
2
(
x
2
+
y
2
)
(
A
1
/
t
1
)
3
n
-
2
c
2
t
1
2
R
3
[
(
A
1
/
t
1
)
n
+
1
]
2
}
R
=
x
2
+
y
2
+
(
z
-
z
′
)
2
;
A
1
=
t
-
R
c
-
|
z
′
|
v
;
A
2
=
exp
⁡
⁡
(
(
R
/
c
)
-
t
+
(
|
z
′
|
/
v
)
t
2
)In addition, by using the same methodology that was used for the first current function, the components of electromagnetic field for the second current function as expressed in (2) are proposed by (15). Note that B→z_RS2x,y,z,tn=0:
(15)B→x_RS2x,y,z,tn=∑i=1n∑m=1k+1amF4x,y,z,t=tn,z′=hm,ihhhhhhhhhhhh-am′F4x,y,z,t=tn,z′=hm,i′,B→y_RS2x,y,z,,tn=∑i=1n∑m=1k+1-xyamF4x,y,z,t=tn,z′=hm,ihhhhhhhhhhhhhhh-am′F4x,y,z,t=tn,z′=hm,i′,E→x_RS2x,y,z,tn=E→x_RS2x,y,z,tn-1+Δt×∑i=1n∑m=1k+1amF5x,y,z,t=tn,z′=hm,ihihhhhhhhhhhh-am′F5x,y,z,t=tn,z′=hm,i′,E→y_RS2x,y,z,tn=E→y_RS2x,y,z,tn-1+Δt×∑i=1n∑m=1k+1yxamF5x,y,z,t=tn,z′=hm,ihihhhhhhhhhhhhhhh-am′F5x,y,z,t=tn,z′=hm,i′,E→z_RS2x,y,z,tn=E→z_RS2x,y,z,tn-1+Δt×∑i=1n∑m=1k+1amF6x,y,z,t=tn,z′=hm,ihhhhhhhhhhhh-am′F6x,y,z,t=tn,z′=hm,i′,
where B→x_RS2 is the magnetic flux density at x-axis due to return stroke current function from (2), B→y_RS2 is the magnetic flux density at y-axis due to return stroke current function from (2), B→z_RS2 is the magnetic flux density at z-axis due to return stroke current function from (2), E→x_RS2 is the electric field at x-axis due to return stroke current function from (2), E→y_RS2 is the electric field at y-axis due to return stroke current function from (2), and E→z_RS2 is the electric field at z-axis due to return stroke current function from (2).Consequently, due to the whole return stroke current functions from Table1, the electromagnetic fields can be easily estimated via the proposed field expressions whereas the total fields due to two part functions are equal to the sum of the electromagnetic fields due to each part independently. On top of that the first part current functions in Table 3 should be converted to either (1) or (2) format before getting into proposed expressions. In other words, the proposed algorithm can support the whole current models related to (4). In the same way, the components of the electromagnetic fields in the Cartesian coordinates can be converted to cylindrical domain using (16) and (17) for the magnetic flux density and the electric field components, respectively:
(16)B→rr,z,ϕ,tnB→φr,z,ϕ,tnB→zr,z,ϕ,tn=cos⁡⁡π2-ϕsin⁡π2-ϕ0-sin⁡π2-ϕcos⁡⁡π2-ϕ0001×B→xx,y,z,tnB→yx,y,z,tnB→zx,y,z,tn,(17)E→rr,z,ϕ,tnE→φr,z,ϕ,tnE→zr,z,ϕ,tn=cos⁡⁡π2-ϕsin⁡π2-ϕ0-sin⁡π2-ϕcos⁡⁡π2-ϕ0001×E→xx,y,z,tnE→yx,y,z,tnE→zx,y,z,tn,
where ϕ=arccos(y/x2+y2), B→rr,z,ϕ,t/E→rr,z,ϕ,t is the magnetic flux density/electric field at r-direction (horizontal), B→φr,z,ϕ,t/E→φr,z,ϕ,t is the magnetic flux density/electric field at ϕ-direction, and B→zr,z,ϕ,t/E→zr,z,ϕ,t is the magnetic flux density/electric field at z-direction (vertical).


## 5. Result and Discussion

This section deals with the validation of the electromagnetic fields due to lightning channel using some channel base current functions. According to the proposed method, the widely used current functions can be converted into sum of individual parts using (1) and (2). Therefore, the electromagnetic fields due to each part can be calculated by (12) and (15). Then, the values of total electromagnetic field are equal to the sum of individual components of electromagnetic field. When the MTLE model is applied as current model, the simulated magnetic flux density and the vertical electric field caused by a sum of two Heidler function as a channel base current function at an observation point with 4.6 km distance from lightning channel are depicted in Figures 2 and 3, respectively. The measured magnetic flux density and vertical electric field with similar primary conditions as the ones in Figures 2 and 3 are demonstrated in Figures 4 and 5, respectively. Therefore, comparisons between the measured and simulated fields for magnetic flux density and vertical eclectic field are tabulated in Tables 4 and 5, respectively. It should be noted that the current parameters are i01=17kA, i02=8kA, Γ11=0.4×10-6, Γ12=4×10-6, Γ21=4×10-6, Γ22=50×10-6, n1=2, n2=2, v=1×108, and λ=1500 [23].Table 4
Numerical comparison between the simulated magnetic flux density and the corresponding measured fields from Figures2 and 4 after subtracting initial decay time.
Time (μs)
2
5
8
10
15
B
φ (measured), Wb/m2
2.75
×
10
-
7
2
×
10
-
7
1.7
×
10
-
7
1.5
×
10
-
7
1.3
×
10
-
7
B
φ (simulated), Wb/m2
2.7
×
10
-
7
1.94
×
10
-
7
1.6
×
10
-
7
1.4
×
10
-
7
1.25
×
10
-
7
Percentage difference, %
2
3
6
6
4Table 5
Numerical comparison between the simulated vertical electric field and the corresponding measured fields from Figures3 and 5 after subtracting initial decay time.
Time (μs)
3
10
15
20
25
-
E
z (measured), v/m
81
53
53
58
62
-
E
z (simulated), v/m
79
51
50
54
58
Percentage difference, %
3
4
6
7
6Figure 2
The simulated magnetic flux density using proposed method (z=10 m, r=4.6 km).Figure 3
The simulated vertical electric field using proposed method (z=10 m, r=4.6 km).Figure 4
The measured magnetic flux density for an observation point (z=10 m, r=4.6 km) [23, 24].Figure 5
The measured vertical electric field for an observation point (z=10 m, r=4.6 km) [23, 24].Additionally, the sum of two Heidler functions as a channel base current function can be expressed by the sum of two general current functions based on (2). Therefore, the total electromagnetic fields are estimated by using the proposed method and considering each current part separately. The comparison between simulated and measured fields shows the overall good agreement with experimental measurement for proposed field expressions. Similarly, the simulated fields are in accord with other previous method of calculation in [23, 24].The proposed method is also validated by Nucci current function in Table1 which is defined by the combination of (1) and (2) when the MTLL model is applied as a current model. The simulated magnetic flux density and vertical electric field are formed as in Figures 6 and 7, respectively. By the same token, Figure 8 depicts the measured vertical electric field.Figure 6
The simulated magnetic flux density using Nucci current function (z=0 m, r=50 m).Figure 7
The simulated vertical electric field using Nucci current function (z=0 m, r=15 m).Figure 8
The measured vertical electric field at an observation point (z=0 m, r=15 m) [5].In addition, the numerical comparison between the measured fields and simulated vertical fields is listed in Table6 suggesting that the simulated field is in agreement with the corresponding measured field. Moreover, magnetic flux densities simulated using the proposed field expressions are compatible with the simulated field calculated using other methods given in the reference [5]. Note that the current parameters are i01=3.25kA, i02=8.95kA, Γ11=0.072×10-6, Γ12=16.67×10-6, Γ21=100×10-6, Γ22=0.5×10-6, n1=2, and v=1.5×108 m/s [5]. The proposed expressions can directly be applied to widely used current functions in the time domain just by substituting the current and geometrical parameters in the proposed field expressions without needing to apply any extra algorithms and extra conversions (such as Fourier transform) when the realistic current functions are applied [1]. Besides, the proposed algorithm expressions can support different engineering current models, provided that they follow (4).Table 6
Numerical comparison between the simulated magnetic flux density and the corresponding measured fields from Figures7 and 8 after subtracting initial decay time.
Time (μs)
5
10
15
20
25
30
-
E
z (measured), kv/m
92
93
94.5
95.5
96
97.5
-
E
z (simulated), kv/m
90.5
91
92
93.5
94
95.5
Percentage difference, %
1.6
2.1
2.6
2
2
2The proposed method can corroborate most of the coupling models that compute lightning induced overvoltage (LIOV), directly in the time domain, because not only are all of the electromagnetic field components accessible but also the fundamental equations are in the Cartesian coordinates. This adds a new dimension to the development of LIOV coupling models which are restricted, in the present context, to apply a fast and direct computational block in the time domain for estimation of electromagnetic fields generated by lightning channel without needing to use extra computational stage as Fourier transform for a large number of channel base current functions and engineering current models. Likewise, the proposed field expressions can be used in the inverse procedure algorithms to evaluate lightning return stroke current using measured electromagnetic fields directly in the time domain [8, 9]. Moreover, the proposed field expressions can be developed for the case of nonperfect ground conductivity using Cooray-Rubinstein equation [25].Consequently, the evaluation of electromagnetic fields at different points along the power line can easily be computed with different values ofx or y [25–28]. The accuracy of the results can be increased by changing the values of k and Δt. However, the processing time will increase with such increasing of k. Further, the proposed field expressions can be used for close and medium distances from the lightning channel; however, some of the calculation methods in the time domain are merely confined to close distances from the lightning channel.


## 6. Conclusion

In this study, the general electromagnetic field expressions due to vertical lightning channel at an observation point are proposed directly in the time domain just by substituting of current and geometrical parameters in the proposed field expressions without needing to apply any algorithms and extra conversions. The proposed field expressions can support the widely used current functions and current models. Besides, the proposed method was validated with two return stroke current samples. The results showed that there is a good agreement between the electromagnetic fields computed by proposed method and both the measured fields and those calculated by previous simulation methods. In addition, the proposed method is applicable for both close and intermediate distances from the lightning channel and it is completely in accordance with many coupling models in evaluating the induced lightning voltage and also inverse procedure algorithms in the time domain. This would be very useful especially for the utility in assessing the lightning performance of an overhead line, particularly the distribution line.



---

*Source: 925463-2014-11-12.xml*