Electromagnetic Computation Methods for Lightning Surge Protection Studies

Yoshihiro Baba, Vladimir A. Rakov

• Publisher: Wiley-IEEE Press
• Year: 2016
• ISBN:
  9781118275665
• Edition: 1
• Topic: Scienze fisiche
• Subtopic: Elettromagnetismo

Presents current research into electromagnetic computation theories with particular emphasis on Finite-Difference Time-Domain Method

This book is the first to consolidate current research and to examine the theories of electromagnetic computation methods in relation to lightning surge protection. The authors introduce and compare existing electromagnetic computation methods such as the method of moments (MOM), the partial element equivalent circuit (PEEC), the finite element method (FEM), the transmission-line modeling (TLM) method, and the finite-difference time-domain (FDTD) method. The application of FDTD method to lightning protection studies is a topic that has matured through many practical applications in the past decade, and the authors explain the derivation of Maxwell's equations required by the FDTD, and modeling of various electrical components needed in computing lightning electromagnetic fields and surges with the FDTD method. The book describes the application of FDTD method to current and emerging problems of lightning surge protection of continuously more complex installations, particularly in critical infrastructures of energy and information, such as overhead power lines, air-insulated sub-stations, wind turbine generator towers and telecommunication towers.

• Both authors are internationally recognized experts in the area of lightning study and this is the first book to present current research in lightning surge protection
• Examines in detail why lightning surges occur and what can be done to protect against them
• Includes theories of electromagnetic computation methods and many examples of their application
• Accompanied by a sample printed program based on the finite-difference time-domain (FDTD) method written in C++ program

1
Introduction

1.1 Historical Overview of Lightning Electromagnetic-Field and Surge Computations

Lightning return-stroke electromagnetic fields have been calculated using analytical expressions, derived for a vertical lightning channel (e.g., Uman et al. 1975). Effects of finite ground conductivity on lightning electromagnetic fields have also been studied using analytical expressions (e.g., Rachidi et al. 1996). These analytical expressions are still being used. Lightning-induced voltages on an overhead power distribution line or telecommunication line have been calculated using an engineering model of the lightning return stroke (e.g., Uman et al. 1975) and a field-to-wire coupling model (e.g., Rachidi 1993). Horizontal electric fields above a finitely conducting ground, which are needed for calculating lightning-induced voltages, have been evaluated using approximate expressions such as the Cooray–Rubinstein formula (Rubinstein 1996). Note that the Cooray–Rubinstein formula is given in the frequency domain, although its time-domain counterparts also exist. Lightning surges due to a direct lightning strike to an overhead power transmission or distribution line have been analyzed using distributed-circuit simulation methods such as the electromagnetic transients program (EMTP) (Dommel 1969). EMTP and other similar programs are still widely used in lightning surge simulations.

Around 1990, electromagnetic computation methods were first applied to lightning electromagnetic and surge simulations. One of the advantages of electromagnetic computation methods, in comparison with circuit simulation methods, is that they allow a self-consistent full-wave solution for both the transient current distribution in a 3D conductor system and resultant electromagnetic fields, although they are computationally expensive. Podgorski and Landt (1987) applied the method of moments (MoM) in the time domain (Van Baricum and Miller 1972; Miller et al. 1973) to analyze the lightning current along a tall object struck by lightning. Grcev and Dawalibi (1990) applied the MoM in the frequency domain (Harrington 1968) to analyze the surge characteristics of a grounding electrode. Since then, the MoM in the frequency domain has been frequently used in lightning surge simulations (e.g., Baba and Rakov 2008 and references therein).

Tanabe (2001) applied the finite-difference time domain (FDTD) method (Yee 1966), which is one of the electromagnetic computation methods, to studying the surge characteristics of a grounding electrode. Baba and Rakov (2003) used the FDTD method to compute lightning electromagnetic fields. More than 60 journal papers and a large number of conference papers, which use the FDTD method in lightning electromagnetic-field and surge simulations, have been published during the last 15 years (e.g., Baba and Rakov 2014 and references therein). Interest in using the FDTD method continues to grow. The FDTD method is presently the most widely used electromagnetic computation method in lightning electromagnetic-field and surge simulations.

Other electromagnetic computation methods such as the finite-element method (FEM) (e.g., Sadiku 1989), the partial-element equivalent-circuit (PEEC) method (Ruehli 1974), the hybrid electromagnetic model (HEM) (Visacro and Soares 2005), the transmission line matrix or modeling (TLM) method (Johns and Beurle 1971), and the constrained interpolation profile (CIP) method (Takewaki et al. 1985) have been recently applied to analyzing lightning electromagnetic fields and surge simulations (e.g., Yutthagowith et al. 2009 (PEEC); Silveira et al. 2010 (HEM); Smajic et al. 2011 (FEM); Tanaka et al. 2014a (TLM); Tanaka et al. 2014b (CIP)).

In the following section, we briefly introduce each of these electromagnetic computation methods.

1.2 Overview of Existing Electromagnetic Computation Methods

1.2.1 Method of Moments

The MoM in the time domain (Van Baricum and Miller 1972; Miller et al. 1973) has been used to analyze responses of thin-wire conducting structures to external transient electromagnetic fields. The entire conducting structure is modeled by a combination of cylindrical wire segments whose radii are much smaller than the wavelengths of interest. The so-called electric-field integral equation for a perfectly conducting thin wire in air (shown in Figure 1.1) is given below, assuming that current I and charge q are confined to the wire axis (thin-wire approximation) and that the boundary condition on the tangential electric field on the surface of the wire (this field must be equal to zero) is fulfilled:

Figure 1.1 Thin-wire cylindrical segment for method of moment (MoM)-based computation. Current is confined to the wire axis, and the tangential electric field on the surface of the wire is set to zero.

(1.1)

where C is an integration path along the wire axis; E_inc denotes the incident electric field that induces current I; r and t denote the observation location (a point on the wire surface) and time, respectively; r′ and t′ denote the source location (a point on the wire axis) and time, respectively; R = r − r′; s and s′ denote the distance along the wire surface at r and that along the wire axis at r′, respectively; ŝ and ŝ′ denote unit vectors tangential to path C in Eq. (1.1) at r and r′, respectively; μ₀ is the permeability of vacuum; and c is the speed of light. By numerically solving Eq. (1.1), which is based on Maxwell’s equations, the time-dependent current distribution along the wire structure excited by external field is obtained.

The thin-wire time domain (TWTD) code (Van Baricum and Miller 1972), developed at the Lawrence Livermore National Laboratory, is based on the MoM in the time domain. One of the advantages of the time-domain MoM is that it can incorporate nonlinear effects such as the ligh...

Table of contents

• Cover
• Title Page
• Table of Contents
• Preface
• 1 Introduction
• 2 Lightning
• 3 The Finite-Difference Time Domain Method for Solving Maxwell’s Equations
• 4 Applications to Lightning Surge Protection Studies
• Appendix: 3D-FDTD Program in C++
• Index
• End User License Agreement