Οι πύργοι ελέγχου αεροδρομίων αποτελούν βασικό συστατικό για την ομαλή και συνεχή εξυπηρέτηση των πτήσεων πολιτικών και στρατιωτικών αεροσκαφών. Βρίσκονται εγκατεστημένοι σε κάθε αεροδρόμιο και αποτελούν το υψηλότερο κτίριο στο περιβάλλοντα χώρο, πράγμα που δικαιολογεί και τον μεγάλο αριθμό κεραυνικών πληγμάτων που μπορούν να δεχτούν. Σκοπός αυτής της διπλωματικής εργασίας είναι η μελέτη ενός συστήματος αντικεραυνικής προστασίας που βρίσκεται εγκατεστημένο σε έναν πύργο ελέγχου αεροδρομίου. Η απόδοση του συστήματος αντικεραυνικής προστασίας διαδραματίζει σημαντικό ρόλο στην αποτελεσματική λειτουργία του πύργου ελέγχου και στη γενικότερη αποτελεσματική λειτουργία του αερολιμένα. Με τη χρήση του προγράμματος εξομοίωσης αναλογικών και ψηφιακών κυκλωμάτων Orcad-Pspice, προσομοιώνεται το σύστημα αντικεραυνικής προστασίας, με βάση τον πύργο ελέγχου που του διεθνούς αερολιμένα Ελ. Βενιζέλος, μετράται το δυναμικό στη βάση του πύργου ελέγχου και τα ρεύματα στους αγωγούς καθόδου, όταν κεραυνός πλήττει το σύστημα συλλεκτήριων αγωγών. / The airport traffic control towers are a key component for the smooth and continuous service flights for both political and military aircraft. They are located at each airport and are the highest building in the surroundings, which justifies the number of lightning strikes that accepts. The purpose of this thesis is the study of a lightning protection system which is installed in an airport control tower. The efficiency of the lightning protection system constitute an important role in the effective operation of the control tower and the overall efficient operation of the airport. Using the analog and digital circuits simulation program Orcad-Pspice, the lightning protection system simulated , based on the control tower to the International airport Eleftherios Venizelos , the measured potential at the base of the control tower and the currents in the downconductors , when lightning strikes the collectors system

Lightning is a natural phenomenon involving transient high current discharge in the atmosphere. Cloud-to-ground lightning, wherein the discharge occurs between the cloud and the ground is quite hazardous to systems on the ground. Apart from threat to life, the devastating effects of lightning can be mainly of thermal, mechanical and electromagnetic origin. Many a times, thermal and electromagnetic effects are of main concern. A direct hit, wherein the system under consideration becomes a part of the lightning path, could be quite catastrophic to many vulnerable systems like oil rigs, chemical factories, missile/satellite launch pads. From the safety and operational point of view, lightning is of serious concern for electrical systems including transmission lines and substations, nuclear power stations, telecommunication station and data banks. Lightning cannot be avoided, however, by employing a suitable Lightning Protection System (LPS), adequate protection against a direct hit can be provided to ground based systems. A typical lightning protection system involves: 1) Air termination network, which is responsible for stroke interception, 2) Down conductor system, which provides to the stroke current a minimal impedance path to the ground and 3) Earth termination network, for safe dissipation of current into the ground. Similarly, for the indirect effects, which are basically of electromagnetic origin, suitable protection can be designed. The key factors in a protective action involve interception of the dangerous strokes, minimization of the consequential potential rise on down conductors, as well as, at earth termination and keeping the field in the protective volume within an acceptable level. The last aspect can be generally categorized into secondary level protection. For critical systems, the lightning protection system is generally isolated from it. In such designs, potential rise on LPS governs the physical isolation required between the protected and protection system. For a given level of bypass strokes, cost of the LPS increases with the amount of physical separation employed. All most all of the earlier works have concentrated on lightning surge response of power transmission line towers. Apart from their relatively moderate heights, the intention was to arrive at a model, which can be incorporated in circuit simulation software like EMTP. Consequently, they envisage or approximate the mode of propagation to be TEM. In reality, for down conductors of height greater than say 30 m, only TM mode prevails during the initial critical time period. Hence the earlier models cannot be extended to general lightning protection schemes and for down conductor of larger lengths. Only limited literature seems to be available on the characteristics of general down conductor configurations. The problem in hand is very important and some serious research efforts are very much essential. In view of the above, the present work aims to evaluate the rise in potential as well as current injected into the soil at the base for: (i) practical range of down conductor configurations involving single down conductor (with height exceeding 30 m) and (ii) pertinent values of stroke current parameters. The protection schemes considered are isolated vertical down conductor, isolated tower (both square and triangular cross-section) and, tower with insulated lightning mast carrying ground wires. The parameters under consideration are: (i) height and cross section for the down conductor, (ii) clearance between the down conductor and the protected system, (iii) channel geometry, wherein only inclination is to be considered, (iv) velocity of current along the channel and (v) wave shape and rise time for the stroke current. For the evaluation of lightning surge response of transmission line towers, many theoretical and experimental approaches are found in the literature. However, works considering the TM mode of current propagation is relatively limited. In that both experimental and theoretical approaches have been adopted. Theoretical approach invariably adopted numerical field computation in frequency domain using Numerical Electromagnetic code (NEC-2). Fourier Transform techniques are employed to extract the time domain quantities. This approach is very economical, free from experimental errors and least time consuming. Hence it is selected for the present work. However, there are certain limitations in this approach. In NEC simulation, there is a restriction on the size and the arrangement of individual elements. Therefore, although fairly complex tower structures can be simulated, some simplification in the geometry is unavoidable. Such an approximation has been reported to cause insignificant error. NEC is not accurate for calculations in low frequency regime. But in the present work, the initial time regime is of concern wherein the high frequency components dominate. Therefore the above said limitation is not of any serious concern. In order to validate the approach, potential rise is computed for 120 m tall cylindrical down conductor and tower. Results are compared favorably with earlier works, which are based on potential lead wire method. A careful re-look into the ’potential rise’ on the down conductors reveal several things. The electric field in the region between the protection system and protected system is the root cause for the breakdown/flashover. For a given geometry, the integral of the electric field along the shortest path between the two systems must be representing the overall stress on the air gap. Further, for the later time periods, this integral coincides with the well-known quasi-static potential. All the available data and models for breakdown of long air gaps are basically in terms of this quasi-static potential. In view of this, the above path integral is defined as ’equivalent potential rise’ (which will be hereafter termed as ’potential rise’), and taken as the index for surge response. Further, observation of the computed spatio-temporal radial electric field around the down conductor reveals some additional features, which are not common in the quasi-static regime. Electric field reverses its polarity in space, which is due to the opposite current flowing in the lightning channel. Therefore, ’potential rise’, which is taken as the representative for the dielectric stress on the air, should not be evaluated for larger distances. Considering this and noting that the protected system generally lies well within a distance of 50% of the H, height of the down conductor, potential rise is evaluated by integrating electric field within this distance (12.5%H, 25%H, 50%H). Three heights (100%H, 75%H, 50%H) are considered for the evaluation of the potential. The influences of various down conductor and lightning channel parameters are analyzed. Finally vertical channel with full velocity for current propagation is arrived for the investigations. Also, the influence of neighboring conducting objects is briefly studied. It is argued that it needs to be ignored for the general study. Analysis is carried out for a range of down conductor configurations of heights ranging from 45 m to 120 m. Cylindrical down conductor is selected for the detailed study on the overall characteristics and its dependency on pertinent parameters. The characteristics of potential rise are found to be significantly different from that given by the commonly employed uniform transmission line model. In the regime of very fast front currents, down conductor of comparable heights have comparable potential rise. For the larger time to crest, behavior tends more to wards that for quasi-static regime. The dependency of the potential rise on radius of the down conductor seems to be logarithmic in nature. Surge response of isolated towers of both square and triangular cross sections is studied for heights ranging from 45 m to 120 m. The overall characteristics are found to be similar to cylindrical down conductor. Dispersive propagation is found to exist on towers. As a result, the base currents are slightly lower and potential rise exhibits less oscillations. Data curves on potential rise at three different heights and for three different spatial extents are generated for the range of down conductor heights with rise time of the stroke current as the variable. Several interesting observations have been made. Next the investigation is taken up for the insulated mast scheme. The parameters of the study are taken as the number of ground wires, grounding location of ground wires and length of the insulation cylinder. Potential across the insulation, tower base currents, and ground wire end currents are deduced. The basic characteristics of the potential rise are shown to be quite similar to that for the transmission line. For fast front currents the temporal variation is bipolar with a smooth decay. In other words, oscillations are sustained for considerably longer duration. Voltage stress across the insulation surface for one ground wire design is found to be higher by 1.4 - 2.4 times than that for isolated tower. The highest amplification of the ground end current, which occurs for fast front currents, is about 1.8 times. Potential difference across the insulation for two-ground wire design is higher by a factor of 1.3 - 1.85 than that for isolated tower. For the design with four ground wires, potential across the insulation is comparable with that for the tower. However, the mechanical strength of the insulating support should also be considered in the selection of number of ground wires. There exists, especially for fast front strokes, significant induction to the supporting tower. The height of the insulation seems to possess no appreciable influence on the potential rise and base currents. Several issues need to be considered before selecting this design. The contribution made by the present work can be summarized as follows. It basically deals with lightning surge response of isolated down conductors of height in the range 45 - 120 m. The configurations considered are, cylindrical down conductor, tower with both square and triangular cross section and insulated mast scheme. It makes a careful study on the ’potential rise’ on down conductors and a suitable definition for the same is proposed. Basic characteristics of potential rise and ground end currents are studied for the above-mentioned designs. Their salient features are enumerated. For the towers, design data curves are provided for relevant range of stroke current rise time. The issues that need to be considered in the insulated mast scheme are discussed along with the data on potential rise and base currents. The findings of this work are believed to be very useful for the design of lightning protection scheme involving isolated down conductor. Further the results are useful in analyzing the consequential lightning generated threat of being close to tall towers.

Satellites have become absolute necessity in the growing modern space technology. At present, launch pads are the only means for launching of satellites or any other space vehicles. Due to the large magnitude of current and the associated rate of rise, a lightning strike to launch pads can be quite disastrous. Satellite launch complex forms typically the tallest object in that region. This makes them the more vulnerable to cloud-to-ground lightning. In addition, most of the launch pads are situated near the coastal area, where the isokeraunic levels are quite high. In view of these, almost all the satellite launch pads are provided with suitable Lightning Protection Systems (LPS). The LPS is basically intended for protecting against a direct lightning hit. The present work is related with the LPS to Indian satellite launch pads, Pad-I and Pad-II. The protection system for Pad-I consists of three 120 m tall towers placed approximately at the vertices of an equilateral triangle of 180 m. The same for Pad-II consist of 120 m tall towers placed at vertices of rectangle of size 90 m x 105 m. Towers are interconnected by 6 shield wires at the top. A mast of 10m length forms the top of the tower. Significant work on the analysis of interception efficacy of these protection systems has been reported in the literature. The lightning surge response of these systems have also been analysed and reported. The interception efficacy of these LPS in field can be ascertained by pertinent measurements. Measuring the lightning current on LPS seems to be one of the most suitable choices for this purpose. It would also greatly facilitate collection of local lightning current statistics, data on which is almost absent. Several considerations suggest that the tower bases form ideal place for such measurement. However, such lightning current records would involve mainly the current resulting from stroke interception, as well as, induced current due to strokes nearby. Literature on categorisation of measured currents to the type of stroke and correlation of measured currents to the incident stroke currents is rather limited. This is especially true for interconnected protection system of the type dealt in the present work. Considering these the present work is taken up and its scope is defined as: (i) Evolve a suitable model for study of current distribution in LPS due to Lightning and using the same deduce the current due to stroke interception and that due to stroke nearby. (ii) For the purpose of categorization identify the salient characteristics of current due to the intercepted strokes and that due to bypass/nearby strokes (iii) For the intercepted strokes, develop a processor for estimating the injected stroke current from the measured tower base currents. Lightning event, apart from other associated physical phenomena, is strongly governed by electromagnetic fields. Any method employed for the analysis, either theoretical or experimental, should satisfy the governing electromagnetic equations. As experimentation on actual system, as well as, their laboratory simulation is nearly impossible, theoretical modelling approach is selected. Modelling involves modelling of the channel along with its excitation, modelling of the LPS and modelling of the ground. Channel, following the literature, is represented as a loaded conductor with a lumped current source at the junction point. Such models have quite successfully predicted the electromagnetic fields and current in other places on the down conductor. For the LPS, some simplifications on the geometry are very essential. Tower lattice elements of dimensions much smaller than the wavelength of highest dominant frequency component of lightning current spectrum are neglected. Suitable modification is made for the tower top involving a plate and interconnection of several short members. For the close range within 200 – 400 m, even for the induced currents, the influence of ground in the literature has been reported to be small. Also, there is an extensive grounding network in these systems. In view of the same, a perfectly conducting ground along with suitable ground termination impedance is considered. Only the numerical solution of the problem is feasible and for the same, following the literature, NEC-2 is employed. All the guidelines of NEC are respected in the discretisation. Geometric mean radius is employed for modelling the complex tower elements. Fourier Transform Techniques are employed for time domain conversion of the computed frequency domain quantities. Occasionally, numerical inversion error of magnitude less than 5% is encountered. For the validation of the numerical modelling for both direct stroke and that nearby, time domain experimentation on electromagnetically reduced scale models (35:1) is employed. As the channel electrical and geometrical parameters are stochastic in nature, it is necessary to ensure that the deduction made using the model is practically relevant. For this, some parametric studies are conducted. The influence of channel length and inclination, stroke current velocity etc. has been shown to be insignificant for the case of intercepted strokes. Simulations are carried out for the stroke intercepted (i.e. direct strikes) by the LPS. The characteristics of the tower base currents are investigated. The base currents indicate a dispersive propagation along the towers and further a frequency dependent current division at the tower-shield wire junctions. Base currents contain superimposed oscillations, which basically originate from various junctions of the system. The magnitude of the oscillations is obviously dependent on the rise time of the incident currents. The tower base currents settle within about 10 -15 µs, which is shorter than that for isolated tower. Further, the full-frequency model could be limited to this time period. The corresponding current transfer functions are deduced. For the stroke interception by shield wires, based on the earlier work, only stroke to midspan is found to be relevant and hence it is considered. The nature of tower base currents for a stroke to midspan of the shield wires seem to be similar. However there are some distinct features, which are helpful in identifying the stroke location on the LPS. From the time correlated tower base currents, a suitable methodology for identifying the stroke interception location on LPS is developed. Next, simulations for induced current due to a bypass stroke, as well as, stroke to ground outside the LPS, however, within 1 km radius are taken up. In fact, it is estimated that latter is nearly 5 – 13 times higher than the strokes collected by LPS, indicating it as the most probable event. The objective here is characterization, rather than correlation. In this study, the influence of charge induced on the LPS by the descending leader is neglected and the upward leader activity is approximately considered. To the best of author’s knowledge, studies on such induced currents in down conductors are very scarce. Considering this and noting that the number of parameters is quite large, first the basic study is taken up on simple cylindrical down conductors. Many important and interesting deductions are made. The nature of the induced current is highly dependent on the rate of rise as well as the velocity of propagation of the stroke current. The magnitude and to some extent, the wave shape of the induced current is found to depend on the average as well as maximum di/dt of the stroke current. For a given wave shape, the magnitude of the induced current increases with rate of rise of the wave front; however, saturating trend will onset after some point. The height of the down conductor mainly governs the frequency of the oscillatory component of the induced current. The dependency of the induced current on the radius of the down conductor seems to be logarithmic (which is in accordance with the antenna theory). Based on these results, the parameters for the corresponding study on LPS under consideration, is chosen. The results of the investigation on the induced currents in LPS show that they have quite distinct waveform. They are basically bipolar and oscillatory in nature, with relatively short duration. These unique features facilitate clear distinction of the induced currents from that due to stroke interception. Basic characteristics are reasonably insensitive to the separation distance of the protection system and the channel, current propagation velocity along the channel, channel inclination and shape of the current front. The salient features of the induced current due to a bypass stroke are also enumerated. • The noise, if any, in the measured current can be addressed only after acquiring sufficient data. Based on the above, the following procedure is suggested for the stroke classification and estimation. • By employing the distinct features of the resulting tower base currents, analyze the measured tower base currents and classify the strokes into the intercepted stroke or stroke to ground. • For the latter case, using the salient features of the bypass strokes, further classify the strokes to bypass strokes and stroke to ground outside the protected volume. • For the intercepted strokes, using the relative strengths and wave shapes, identify the interception point to either tower top or the midspan of the shield wires. • Then by using the corresponding transfer functions and Fourier Transform techniques, compute the injected stroke current. • Using the above, other tower base currents are computed and compared with the measured currents. This gives quantification for the accuracy of the method. In summary the present work has made some original contribution to the classification and estimation of stroke currents measured on the interconnected LPS.

In India, the natural energy resources (thermal and hydro) are unevenly distributed and are mostly present in the remote areas and the load centers are distributed across various regions of the country. Therefore high voltage lines have become necessary for the devel-opment of large interconnected power networks and for the reliable and economic transfer of power. The increase in electric power demand due to the electric load growth has lead to the expansion of the transmission systems to ultra high voltage levels. Presently, Ultra High Voltage (UHV) power transmission lines are being built to transfer large electric power to distant load centers from the generating stations. Increasing the line voltage increases the surge impedance loading, stability and the thermal capacity of the line. Lightning is one of the major causes for the line outages and interruptions of UHV power lines. A lightning strike generates a very large voltage leading to insulator puncture, melting, burning and pitting of conductors and the supporting hardware. Lightning can lead to transient over-voltages thus leading to ash-over in the power transmission lines which are dangerous for the power equipments as well as for the human beings working in the vicinity. Ground wires are used for the protection of overhead power transmission lines against a lightning stroke. The overhead ground wires are installed such that the lightning attaches to it and shunts the lightning current to the ground through the tower, thus protecting the phase conductors. Shielding failure happens when the lightning strikes the phase conductor instead of the ground wires. Lightning shielding failure is a major con-cern in UHV lines due to their large height, very high operating voltage and wide exposure area of the phase conductors. The lightning over-voltages injected on the phase conductor (shielding failure) nally reaches the substation causing serious threat to the substation components and can lead to temporary or permanent outage of the power transmission system. There have been cases of very high shielding failure ash-overs of UHV lines and thus lightning attachment to power transmission lines need to be studied in detail to reduce the power system line outages. Several models such as electro-geometric model (EGM) and leader progression model (LPM) have been developed to study the shielding failure of power transmission lines. EGM has been extensively used to obtain lightning attachment to power transmission lines but in recent years it is seen that EGM is unable to accurately predict the lightning attach-ment to UHVAC lines. The shielding failure rates obtained by EGM does not match with the observed shielding failure rate for UHV lines. For this reason LPM is considered to obtain lightning attachment to UHV lines but LPM in its initial stage do not deal with the detailed physics of the upward leader inception, i.e., corona inception and unstable as well as stable upward leader inception from the object on the ground. In this thesis a model for the lightning attachment has been developed based on the present knowledge of the lightning physics. The thesis mainly focuses on the modelling of upward leader inception and propagation for lightning attachment to UHV power trans-mission lines. Upward leader inception is modeled based on the corona charge present near the conductor region and the upward leader propagation model is based on the correlation between the lightning induced voltage on the conductor and the voltage drop along the upward leader channel. The present model considers corona inception and modelling of unstable and stable upward leader inception from the ground object for the analysis of the lightning attachment process. The upward leader inception model developed is compared with the previous inception models and the results obtained using the present and previous models are found to be comparable. Lightning striking distances ( nal jump) for various lightning return stroke current were computed for di erent conductor heights using present lightning attachment model. It is seen that the striking distance increases with the increase in lightning re-turn stroke current and increases with increase in conductor heights. The striking distance computed using the present model matches with the value calculated using the equation proposed by the IEEE working group for the applicable conductor heights of up to 8 m. The in uence of the conductor operating voltage, cloud electric eld, lightning down-ward leader lateral distance, conductor length, transmission line tower and conductor sag on the upward lightning leader inception are analysed and reported in the thesis. It is found that the lightning attraction to power transmission line increases with increase in conductor positive operating voltage and decreases with increase in conductor negative op-erating voltage. The presence of transmission line tower reduces the lightning attachment to the conductor lines and the probability of lightning strike decreases with the increase in downward leader lateral distance from the conductor lines. The present lightning attachment model is applied to study the shielding failure of UHV power transmission lines rated for 1200 kV ac (delta and horizontal con guration) and for 800 kV dc (with and without a dedicated metallic return conductor) and thereby the lightning shielding failure ash-over rate is computed for the UHV power transmission lines. It is seen that the lightning shielding rate for UHV power transmission lines depend on the lateral distance of the downward leader channel, instantaneous 50 Hz voltage on the transmission line conductor, height of the transmission line conductor, induced voltages on the conductor and the lightning return stroke current.