Lightning ReturnStroke Models 1 Four Classes of Lightning

  • Slides: 48
Download presentation
Lightning Return-Stroke Models 1

Lightning Return-Stroke Models 1

Four Classes of Lightning Return-Stroke Models 1. Gas-dynamic models 2. Electromagnetic models 3. Distributed-circuit

Four Classes of Lightning Return-Stroke Models 1. Gas-dynamic models 2. Electromagnetic models 3. Distributed-circuit models 4. Engineering models 2

Classes of Lightning Return-Stroke Models Gas-dynamic models: The gas dynamic or "physical" models are

Classes of Lightning Return-Stroke Models Gas-dynamic models: The gas dynamic or "physical" models are primarily concerned with the radial evolution of a short segment of the lightning channel and its associated shock wave. These models typically involve the solution of three gas-dynamic equations (sometimes called hydrodynamic equations) representing the conservation of mass, of momentum, and of energy, coupled to two equations of state. Principal model outputs include temperature, pressure, and mass density as a function of the radial coordinate and time. 3 3

Gas-Dynamic Model of Paxton et al. (1990) 74 ns 91 μs Paxton et al.

Gas-Dynamic Model of Paxton et al. (1990) 74 ns 91 μs Paxton et al. (1990) 4 4

Electromagnetic models : The electromagnetic models are usually based on a thin-wire antenna approximation

Electromagnetic models : The electromagnetic models are usually based on a thin-wire antenna approximation to the lightning channel. These models involve a numerical solution of Maxwell's equations to find the current distribution along the channel, from which the remote electric and magnetic fields can be computed. [For example, the Electric Field Integral Equation (EFIE), equating to zero the tangential component of E-field on the surface of good conductors, is solved using Mo. M or the Faraday’s Law and Ampere’s Law equations are solved using the FDTD method. ] 5

Electromagnetic (antenna-theory) model of Moini et al (2000) Moini et al. (2000) 6 6

Electromagnetic (antenna-theory) model of Moini et al (2000) Moini et al. (2000) 6 6

Role of corona in making the return-stroke speed lower than the speed of light

Role of corona in making the return-stroke speed lower than the speed of light 2 rcore L Channel core (no corona) rcore C 2 rcorona 2 rcore Channel core (~104 S/m) 7 Corona sheath (~10− 6– 10− 5 S/m) rcore L rcorona C

8

8

Classes of Lightning Return-Stroke Models Distributed-circuit models: The distributed-circuit models represent the lightning discharge

Classes of Lightning Return-Stroke Models Distributed-circuit models: The distributed-circuit models represent the lightning discharge as a transient process on a vertical transmission line characterized by resistance (R), inductance (L), and capacitance (C), all per unit length. The governing equations are telegrapher’s equations. The distributedcircuit models (also called R-L-C transmission line models) are used to determine the channel current versus time and height and can therefore also be used for the computation of remote electric and magnetic fields. 9 9

Distributed-Circuit Models (usually R=f(t), L=const, C=const, G=0) Telegrapher’s Equations: 10 10

Distributed-Circuit Models (usually R=f(t), L=const, C=const, G=0) Telegrapher’s Equations: 10 10

Distributed-Circuit Model of Gorin and Markin (1975) t=1. 8 μs, V 0=50 MV. Instantaneously

Distributed-Circuit Model of Gorin and Markin (1975) t=1. 8 μs, V 0=50 MV. Instantaneously discharged corona sheath t=1. 8 μs, V 0=10 MV. No corona sheath Gorin and Markin (1975) 11 11

Engineering models: The engineering models specify a spatial and temporal distribution of the channel

Engineering models: The engineering models specify a spatial and temporal distribution of the channel current based on such observed lightning return-stroke characteristics as current at the channel base, the speed of the upward-propagating front, and the channel luminosity profile. In these models, the physics of the lightning return stroke is deliberately downplayed, and the emphasis is placed on achieving agreement between the model-predicted electromagnetic fields and those observed at distances from tens of meters to hundreds of kilometers. A characteristic feature of the engineering models is the small number of adjustable parameters, usually one or two besides the specified channel-base current. An engineering model is just a simple equation relating I(z’, t) and I(0, t). 12

Engineering Models Transmission Line (TL) Model vw I(z’, t) = I(z’=0, t - z’/vw)

Engineering Models Transmission Line (TL) Model vw I(z’, t) = I(z’=0, t - z’/vw) vw = z’/(z’/vw) = tanα, where α is the angle between the slanted broken line and the time axis. In this model, vf = vw = v, where vf is the upwardpropagating return-stroke front speed and vw is the current-wave propagation speed Subscript “w” is often dropped, so that: I(z’, t) = I(z’=0, t - z’/ v) 13 13 z’ 2/vw

Engineering Models Lightning return stroke models are used for calculating electromagnetic fields in various

Engineering Models Lightning return stroke models are used for calculating electromagnetic fields in various studies including: - lightning electromagnetic coupling with power and communication lines (e. g. , Zeddan and Degauque, 1990; Rachidi et al. 1996) - estimation of lightning properties from measured electric and magnetic fields (e. g. , Baker et al. 1990; Krider et al. 1996; Uman et al. 2002) - the production of transient optical emissions (elves) in the lower ionosphere (e. g. , Krider, 1994; Rakov and Tuni, 2003). An engineering return-stroke model is defined here as an equation relating the longitudinal channel current I(z', t) at any height z' and any time t to the current I(0, t) at the channel origin, z' = 0. Several simplest engineering models can be expressed by the following generalized current equation: I ( z’, t ) = u ( t - z’/vf ) P ( z’ ) I ( 0, t – z’/v ) where: u is the Heaviside function equal to unity for t ≥ z'/vf and zero otherwise P(z') is the height-dependent current attenuation factor vf is the upward-propagating return-stroke front speed v is the current-wave propagation speed 1 (1)

3. Engineering Models I ( z’, t ) = u ( t - z’/vf

3. Engineering Models I ( z’, t ) = u ( t - z’/vf ) P ( z’ ) I ( 0, t – z’/v ) Model Reference Transmission Line (TL) Model Uman and Mc. Lain (1969) Modified Transmission Line Model with Linear Current Decay with Height (MTLL) Rakov and Dulzon (1987) Modified Transmission Line Model with Exponential Current Decay with Height (MTLE) Nucci et al. (1988 a) Bruce-Golde (BG) Model Traveling Current Source (TCS) Model Bruce and Golde (1941) Heidler (1985) P (z’) v 1 vf 1 - z'/H vf exp(-z'/λ) vf 1 ∞ 1 -c H is the total channel height, c is the speed of light, λ is the current decay height constant (assumed by Nucci et al. (1988 a) to be 2000 m) 1

3. Engineering Models TCS Model BG Model TL Model Waveforms for current versus time

3. Engineering Models TCS Model BG Model TL Model Waveforms for current versus time t at ground (z′ = 0), and at two heights z 1′ and z 2′ above ground, for the TCS, BG, and TL return-stroke models. The slanting lines labeled vf represent the upward speed of the return-stroke front, and the lines labeled v represent the speed of the return-stroke current wave. The dark portions of the waveforms indicate when the current actually flows through a given channel section. Note that the current waveform at z′ = 0 and the front speed vf are the same for all three models. The Heaviside function u(t - z′/vf ) equals zero for t < z′/vf and unity for t ≥ z′/vf. 1

3. Engineering Models TL Model TCS Model Current versus height z′ above ground at

3. Engineering Models TL Model TCS Model Current versus height z′ above ground at an arbitrary instant of time t = t 1 for the TL and TCS models. Note that the current at z′ = 0 and vf are the same for both models. In both cases negative charge is transferred to the ground. 1

3. Engineering Models The most used engineering models can be grouped in two categories:

3. Engineering Models The most used engineering models can be grouped in two categories: - Transmission-line-type (lumped-source) models - Traveling-current-source-type (distributed-source) models Transmission-line-type models for t ≥ z'/vf TL Model (Uman and Mc. Lain, 1969) MTLL Model (Rakov and Dulzon, 1987) MTLE Model (Nucci et al. 1988 a) v = vf = const 1 H = const λ = const

3. Engineering Models The transmission-line-type models can be viewed as incorporating a lumped current

3. Engineering Models The transmission-line-type models can be viewed as incorporating a lumped current source at the channel base which injects a specified current wave into the channel. This wave propagates upward in the TL model without either distortion or attenuation and in the MTLL and MTLE models without distortion but with specified attenuation. I(z’, t) z’ I(0, t) 1

3. Engineering Models Traveling-current-source-type models for t BG Model (Bruce and Golde, 1941) TCS

3. Engineering Models Traveling-current-source-type models for t BG Model (Bruce and Golde, 1941) TCS Model (Heidler, 1985) DU Model (Diendorfer and Uman, 1990) v* = vf / (1+vf / c) 2 vf = const τD = const ≥ z'/vf

3. Engineering Models In the traveling-currentsource-type models, the return -stroke current may be viewed

3. Engineering Models In the traveling-currentsource-type models, the return -stroke current may be viewed as generated at the upwardmoving return-stroke front and propagating downward. These models can be also viewed as involving current sources distributed along the lightning channel that are progressively activated by the upward-moving return-stroke front. Distributed-source representation of the lightning channel in TCS-type models for the case of no strike object and no reflections at ground. Adapted from Rachidi et al. (2002). 21

Corona discharge from barbed wire Barbed wire subjected to a positive dc voltage of

Corona discharge from barbed wire Barbed wire subjected to a positive dc voltage of 400 -k. V magnitude in laboratory. Adapted from Gary (1994).

Channel core Corona sheath Leakage current DIRECTION OF CURRENT (positive charge transfer) (charge deposited)

Channel core Corona sheath Leakage current DIRECTION OF CURRENT (positive charge transfer) (charge deposited) DIRECTION OF NEGATIVE CHARGE TRANSFER Longitudinal current (charge transferred) MTLL and MTLE models 23

Equivalency Between the Lumped-Source and Distributed. Source Representations (Maslowski and Rakov, 2009). 24 24

Equivalency Between the Lumped-Source and Distributed. Source Representations (Maslowski and Rakov, 2009). 24 24

Model Description The distributions of current along the tower and along the lightning channel:

Model Description The distributions of current along the tower and along the lightning channel: Along the tower: Along the lightning channel: Transmission line representation of lightning strike to a tall object. 25

E-field due to an infinitesimal vertical dipole Idz’ R z’ P (z=0) Perfect conducting

E-field due to an infinitesimal vertical dipole Idz’ R z’ P (z=0) Perfect conducting ground z’ R Idz’ r Electrostatic Induction Radiation 26

Results Typical subsequent RS pulse Gaussian pulse 27

Results Typical subsequent RS pulse Gaussian pulse 27

Tower c Tower bottom Lightning Channel v=c/2 Tower top Height above ground (m) Time

Tower c Tower bottom Lightning Channel v=c/2 Tower top Height above ground (m) Time (µs) 28

Tower bottom Tower top Height above ground (m) Time (µs) 29

Tower bottom Tower top Height above ground (m) Time (µs) 29

Testing of Model Validity 30

Testing of Model Validity 30

Validation of Return-Stroke Models Using Measured Electric and Magnetic Fields “Typical-event” approach involves the

Validation of Return-Stroke Models Using Measured Electric and Magnetic Fields “Typical-event” approach involves the use of a typical channel-base current waveform i(0, t) and a typical front propagation speed vf as inputs to the model, and a comparison of the model-predicted fields with typically observed fields. “Individual-event” approach: In this approach, i(0, t) and vf, both measured for the same individual event are used to compute fields that are compared to the measured fields for that same event. 31 31

Testing Model Validity – Typical Event Approach Typical measured vertical electric field intensity (left

Testing Model Validity – Typical Event Approach Typical measured vertical electric field intensity (left column) and azimuthal magnetic flux density (right column) waveforms for first (solid line) and subsequent (dashed line) return strokes at distances of 1, 2, 5, 10, 15, 50, and 200 km. Adapted from Lin et al. (1979). 32

Testing Model Validity – Typical Event Approach 5 km Calculated vertical electric field (left-hand

Testing Model Validity – Typical Event Approach 5 km Calculated vertical electric field (left-hand scale, solid lines) and horizontal (azimuthal) magnetic field (right-hand scale, broken lines) for four return-stroke models at a distance r = 5 km displayed on (a) 100 µs and (b) 5 µs time scales. Adapted from Nucci et al. (1990). 100 km Calculated vertical electric field (left-hand scale) and horizontal (azimuthal) magnetic field (right-hand scale) for four return-stroke models at a distance r = 100 km displayed on (a) 100 µs and (b) 5 µs time scales. Adapted from Nucci et al. (1990). 33

Testing Model Validity – Typical Event Approach Leader Return Stroke Zero Crossing 15 m

Testing Model Validity – Typical Event Approach Leader Return Stroke Zero Crossing 15 m Leader 15 m 30 m DEL DERS DTL 30 m Electric field and electric field derivative (d. E/dt) waveforms for stroke 2 in rocket-triggered flash S 9918 measured at 15 and 30 m from the lightning channel at Camp Blanding, Florida. 34

Testing Model Validity – Typical Event Approach 50 m (a) The current at ground

Testing Model Validity – Typical Event Approach 50 m (a) The current at ground level and (b) the corresponding current derivative used by Nucci et al. (1990), Rakov and Dulzon (1991), and Thottappilil et al. (1997) for testing the validity of return-stroke models by the “typical-return-stroke” approach. The peak current is about 11 k. A, and peak current rate of rise is about 105 k. A/µs. Adapted from Nucci et al. (1990). 35 Calculated vertical electric fields for six returnstroke models at a distance r = 50 m. Adapted from Thottappillil et al. (1997).

Testing Model Validity – Specific Event Approach TL MTLE TCS DU 5. 15 km

Testing Model Validity – Specific Event Approach TL MTLE TCS DU 5. 15 km The calculated vertical electric fields (dotted lines) from the TL, MTLE, TCS, and DU models shown together with the measured field (solid lines) at 5. 16 km (right panel) for return stroke 8715_10. The measured current at the channel base and the measured return stroke speed are given in the left panel. Adapted from Thottappillil and Uman (1993). 36

Summary The overall results of the testing of the validity of the engineering models

Summary The overall results of the testing of the validity of the engineering models can be summarized as follows. • The relation between the initial field peak and the initial current peak is reasonably well predicted by the TL, MTLE, and DU models. • Electric fields at tens of meters from the channel after the first 10 -15 µs are reasonably reproduced by the MTLL, BG, TCS and DU model, but not by the TL and MTLE models. • From the standpoint of the overall field waveforms at 5 km all the models tested by Thottappillil and Uman (1993) should be considered less than adequate. Based on the entirety of the testing results and mathematical simplicity, Rakov and Uman (1998) ranked the engineering models in the following descending order: MTLL, DU, MTLE, TCS, BG, and TL. However, the TL model is recommended for the estimation of the initial field peak from the current peak or conversely the current peak from the field peak, since it is the mathematically simplest model with a predicted peak field/peak current relation that is equally or more accurate than that of the more mathematically complex models. 37

3. Engineering Models The TL model predicts that, as long as: (1) the height

3. Engineering Models The TL model predicts that, as long as: (1) the height above ground of the upwardmoving return-stroke front is much smaller than the distance r between the observation point on ground and the channel base, (2) the return-stroke front propagates at a constant speed, (3) the return-stroke front has not reached the top of the channel (the first 25 -50 us or so), and (4) the ground conductivity is high enough so that propagation effects are negligible, the vertical component Ezrad of the electric radiation field and the azimuthal component of the magnetic radiation field are proportional to the channel-base current I. The equation for the electric radiation field Ezrad is as follows, (2) where: 0 is the permittivity of free space, v is the upward propagation speed of the current wave, which is the same as the front speed vf in the TL as well as in the MTLL and MTLE models, and c is the speed of light. 38

3. Engineering Models Taking the derivative of with respect to time, one obtains (3)

3. Engineering Models Taking the derivative of with respect to time, one obtains (3) Equations 2 and 3 are commonly used, particularly the first one and its magnetic radiation field counterpart, found from |BΦrad| = |Ezrad|/c, for estimation of the peak values of return-stroke current and its time derivative, subject to the assumptions listed prior to Eq. 2. Equations 2 and 3 have been used for the estimation of v from measured Ep/Ip and (d. E/dt)p/(d. I/dt)p , respectively, where the subscript "z" and superscript "rad" are dropped, and the subscript “p” refers to peak values. Expressions relating channel base current and electric radiation field far from the channel for the BG, TCS, and MTLE models are given by Nucci et al. (1990). 39

Field-to-Current Conversion Equation

Field-to-Current Conversion Equation

Equation Based on the Transmission Line Model Absolute errors (ΔITL% = 100 × ΔITL

Equation Based on the Transmission Line Model Absolute errors (ΔITL% = 100 × ΔITL /ICB, where ΔITL = |ITL − ICB|) in estimating peak currents using the field-to-current conversion equation based on the TL model. v = c/3 v = c/2 Peak current (ITL) obtained using the TL model (for different return-stroke speeds) vs. directlymeasured channel-base peak current (ICB). The blue, red, green, and black lines are the best (least squares) fits to the data for return-stroke speed of c/3, c/2, 2 c/3, and c, respectively. Adapted from Mallick et al. [2014]. v = 2 c/3 v=c Mean 71% 20% 16% 43% Median 74% 19% 13% 42% Min 10% 0% 1. 6% 20% Max 129% 57% 71% The field-to-current conversion equation based on the TL model gives the best match with directly measured peak currents for return-stroke speeds between c/2 and 2 c/3, where c is the speed of light.

2. Return-Stroke Speed 42

2. Return-Stroke Speed 42

Why is the return-stroke speed lower than the speed of light? There are two

Why is the return-stroke speed lower than the speed of light? There are two possible reasons why the lightning returnstroke speed is lower than the speed of light: (1)The ohmic losses in the channel core that are sometimes represented in lightning models by the distributed constant or current-varying series resistance of the channel. The expected resistance per unit length of the dart-leader channel which is traversed by a subsequent return stroke is about 3. 5 /m. The expected resistance per unit length behind the returnstroke front is about 0. 035 Ω/m. (2)The effect of corona surrounding the current-carrying channel core. The radius of the charge-containing corona sheath is considerably larger than the radius of the core carrying the longitudinal channel current, so that (LC)-1/2 ( 0 0)-1/2 = c 43

Role of corona in making the return-stroke speed lower than the speed of light

Role of corona in making the return-stroke speed lower than the speed of light 2 rcore L Channel core (no corona) rcore C 2 rcorona 2 rcore Channel core (~104 S/m) 44 Corona sheath (~10− 6– 10− 5 S/m) rcore L rcorona C

The corona effect explanation is based on the following assumptions: (a) The longitudinal channel

The corona effect explanation is based on the following assumptions: (a) The longitudinal channel current flows only in the channel core, because the core conductivity, of the order of 104 S/m, is much higher than the corona sheath conductivity, of the order of 10− 6– 10− 5 S/m (Maslowski and Rakov, 2006). The longitudinal resistance of channel core is expected to be about 3. 5 Ω/m (Rakov, 1998), while that of a 2 -m radius corona sheath should be of the order of kiloohms to tens of kiloohms per meter. The corona current is radial (transverse) and hence cannot influence the inductance of the channel. 45

(b) The radial voltage drop across the corona sheath is negligible compared to the

(b) The radial voltage drop across the corona sheath is negligible compared to the potential of the lightning channel. According to Gorin (1985), the average radial electric field within the corona sheath is about 0. 5– 1. 0 MV/m, which results in a radial voltage drop of 1– 2 MV across a 2 -m radius corona sheath (expected for subsequent return strokes). The typical channel potential (relative to reference ground) is about 10– 15 MV for subsequent strokes (Rakov, 1998). For first strokes, both the corona sheath radius and channel potential are expected to be larger, so that about an order of magnitude difference between the corona sheath voltage drop and channel potential found for subsequent strokes should hold also for first strokes. 46

(c) The magnetic field due to the longitudinal current in channel core is not

(c) The magnetic field due to the longitudinal current in channel core is not significantly influenced by the corona sheath. For corona sheath conductivity of 10− 6– 10− 5 S/m and frequency of 1 MHz, the field penetration depth is 160 to 500 m (and more for lower frequencies), which is much larger than expected radii of corona sheath of a few meters. In summary, the corona sheath conductivity is low enough to neglect both the longitudinal current through the sheath and shielding effect of the sheath, but high enough to disregard the radial voltage drop across the sheath. 47

Lightning Return-Stroke Models 1. Introduction and Classification of Models 2. Return-Stroke Speed 3. Engineering

Lightning Return-Stroke Models 1. Introduction and Classification of Models 2. Return-Stroke Speed 3. Engineering Models 4. Equivalency Between the Lumped-Source and Distributed-Source Representations 5. Extension of Models to Include a Tall Strike Object 6. Testing of Model Validity 6. 1. Typical-Return-Stroke Approach 6. 1. 1. Distant (1 to 200 km) fields 6. 1. 2. Close (tens to hundreds of meters) fields 6. 2. Specific-Return-Stroke Approach 7. Summary 4