Modelling of Distorted Electrical Power and Its Practical Compensation in Industrial Plant
MODELLING OF DISTORTED ELECTRICAL POWER AND ITS PRACTICAL COMPENSATION IN INDUSTRIAL PLANT
By
JAN HARM CHRISTIAAN PRETORIUS
THESIS presented in partial fulfilment of the requirements for the degree
DOCTOR 1NGENERIAE ' (D. Ing)
in the
FACULTY OF ENGINEERING
of the
RAND AFRIKAANS UNIVERSITY
SUPERVISOR: PROF. J.D. van WYK CO-SUPERVISOR: DR. P.H. SWART
DECEMBER 1997 11
PREFACE
Alternating current systems employing single-frequency sinusoidal waveforms render optimal service when the currents in that system are also sinusoidal and have a fixed phase relationship to the voltages that drive them. Under unity- power factor conditions, the currents are in phase with the voltages and optimal net-energy transfer takes place under minimum loading conditions, i.e. with the lowest effective values of current and voltage in the system.
The above conditions were realised in the earlier years, because supply authorities generated 50 Hz sinusoidal voltages and consumers drew 50 Hz sinusoidal currents with fixed phase relationships to these voltages. Static and rotating electrical equipment like transformers, motors, heating and lighting equipment were equally compatible with this requirement and well-behaved AC networks were more the rule than the exception. The fact that three-phase systems conveyed the bulk of the power from one topographical location to the next did not constrain the utilisation of that concept at all, even though poly-phase transmission systems were necessary to increase the economy of transmission and to furnish non-pulsating power transfer. Also, additional theory had to be developed to handle unbalanced conditions in these multi-phase systems and to take care of complex network analysis and fault conditions.
Difficulties begin to manifest themselves when equipment not meeting these requirements is connected to the network and when the currents it draws are not sinusoidal. An increasing number of applications demand DC-voltage supplies from which DC-currents are to be drawn. Because power transmission is carried out by means of AC networks, the DC is furnished by converting or rectifying the AC-supply. Power-electronic circuits, of which the R 2P2 power supplies the AEC employs are no exception, employ line-commutated AC/DC converters in their front-ends, and fall into that category.
Although these line-commutated, phase-controlled AC/DC converters are capable of handling giga-watt power levels, line-frequency commutation causes the currents they draw on the AC-side to be distorted, even though still to be periodic. These non-sinusoidal currents, drawn from the source, along the transmission lines and through other distribution system immittances, also give rise to non-sinusoidal voltage drops between the source and the load, which results in distorted voltage waveforms at other nodes and at the load.
Harmonic penetration studies are essential to evaluate the performance of transmission systems in the presence of current distortion sources. These sources do not only bring about voltage distortion within the confines of their own borders, but extend their influence outside into those of other consumers as well. Supplyutilities are wary of the distortion introduced into their networks by consumers and initial recommendations have now given way to rigid standards for curbing harmonic pollution by consumers. 111
Because conventional steady-state alternating current circuit theory fails in the presence of distortion there are only two ways in which harmonic penetration studies can be carried out. Numerical integration methods are mandatory in the study of transient performance of electrical networks during switching and similar occurrences, but become cumbersome when the networks contain more than just a few nodes and are impossible to use when several tens or hundreds of nodes are encountered. Fortunately, harmonic penetration studies can be confined to steady-state operating conditions in a network in which voltages and currents are distorted but remain periodic and are therefore Fourier transformable.
When viewed in the frequency-domain, non-sinusoidal but periodic current and voltage waveforms can be represented by discrete frequency spectra. Frequency-domain analysis offers a number of advantages. From the frequency-domain point of view, distortion can be quantified in terms of complex phasor values of voltages and currents at discrete harmonic frequencies that individually lend themselves to conventional circuit theory, permitting calculations to be carried out in extensive networks. Solutions that apply to these individual harmonic frequencies can then be summated across the spectrum to furnish aggregate or joint parameters of currents, voltages and powers and can also be transformed back into the time-domain for the reconstruction of the relevant time-dependent waveforms.
Both the frequency and time-domain waveforms, of voltage and current, constructed in the above manner are concise and convey the same numerical information. When attempting, however, to quantify the circuit behaviour in terms of the classical definitions of active, reactive and apparent power, it is soon discovered that different definitions are possible. The different definitions, unfortunately, lead to divergent results and it is impossible to assess the utility of each different theory on a general basis. Only by applying the different theories in dedicated measurements, can their relative worth be established in terms of specific circumstances. That is the main theme of this dissertation. iv
SUMMARY
An increasing number of applications, ranging from the very small power supplies in computers to multi-megawatt traction systems and DC-arc-furnaces, require DC-voltage supplies in their front-ends. The resonant regulating pulse power (R2P2) units, employed by the Atomic Energy Corporation of SA (AEC) to drive their pulse lasers, are no exception. To compensate the harmonic distortion in the plant to acceptable levels and to facilitate the search for the most economical combination of power supplies, filters and conditioning equipment, different compensation schemes are studied. In order to do this, a steady-state harmonic penetration model, based on harmonic superposition through nodal bus-admittance matrix formulation, has been developed on a Mathcad package, to furnish the flexibility by means of which the different configurations can be modelled. In addition to the availability of all the necessary parameters such as THD-values, rms-voltage and current values and power levels, the model also furnishes time-domain waveforms of the voltages and currents against which actual measurements are experimentally contrasted. The accuracy of this model is verified after which it is used to study larger networks.
The main alternative schemes investigated are those of compensation by means of:
Individual L.V. passive filters, distributed through the plant and of local design;
PWM-controlled, L.V.-connected power-electronic dynamic compensators.
To evaluate the findings a set of power definitions is necessary. Three sets of power definitions have been chosen in this study, namely: the theories of Budeanu, the IEEE Working Group and the Czarnecki power definitions. These definitions will be analysed numerically in the admittance matrix model and on measurements to evaluate their utility and drawbacks in a practical environment. V
OPSOMMING
`n Toenemende aantal toepassings, vanaf kragbronne vir rekenaars tot megawatt aandryfstelsels benodig gelykspanning as inset. Die resonante laaikragbron wat deur die Atoomenergie- korporasie as laser kragbronne gebruik word, is ook 'n voorbeeld hiervan. Om die harmoniese distorsie na aanvaarbare vlakke te kompenseer en om die ondersoek te vergemaklik vir die mees ekonomiese kombinasie van kragbronne, passiewe filters en dinamiese kompenseerders, word verskillende kompensasietopologiee bestudeer. Om dit sinvol te kan doen is 'n bestendige toestand harmoniese penetrasiemodel gebaseer op harmoniese superposisie d.m.v. die nodebus- admittansiematriksformulering op Mathcad ontwikkel. Met die model kan verskillende parameters byvoorbeeld THD, wgk-waardes en drywings bereken word. Die verskillende tydgolfvorms kan ook gegenereer word. Die akkuraatheid van die model word bewys d.m.v. eksperimetele metings.
Kompensering d.m.v. laagspanning passiewe filters en PWM ingevoerde dinamiese kompenseerders word ondersoek. Die drywingsteoriee van Budeanu, die IEEE Werksgroep en die Czarnecki definiesies word bestudeer. Die praktiese toepassings van hierdie teoried word ondersoek. vi
BEDANKINGS
Die volgende persone en instansies het elk 'n unieke hydrae gelewer:
Piet Swart vir sy tegniese leiding, aanmoediging, volgehoue ondersteuning asook hulp — dit het my geinspireer en gemotiveer.
Prof. Daan van Wyk vir sy rigtinggewing; dit was 'n voorreg om die studie onder sy leiding to kon doen.
Personeel van die AEK wat my ondersteun het, in besonder Neill Truter vir sy hulp.
Anlia vir jou aanmoediging, ondersteuning en liefde, ook ons kinders Rachelle en Christiaan wat verstaan het.
The Atomic Energy Corporation of South Africa and COGEMA of France for making this work possible. vii
CONTENTS
THE MLIS PLANT ENVIRONMENT 1
1.1 INTRODUCTION 1 1.2 R2P2 PULSE POWER SUPPLY CONFIGURATION 3 1.3 POWER SUPPLY FRONT-END TOPOLOGY 6 1.4 PLANT LAYOUT 8 1.5 SUMMARY 10
FUNDAMENTAL DEFINITIONS 11
2.1 INTRODUCTION 11 2.2 FUNDAMENTAL DEFINITIONS UNDER SINUSOIDAL CONDITIONS 12 2.2.1 Fundamental considerations 12 2.2.2 Complex voltage and current as a function of time 12 2.2.3 RAIS value 14 2.2.4 Steady-state impedance 15 2.3 FUNDAMENTAL DEFINITIONS UNDER NON-SINUSOIDAL CONDITIONS 15 2.3.1 Complex voltage and current 15 2.3.2 Harmonic phasors 17 2.3.3 Space-vector manipulations 17 2.3.4 Effective value 18 2.3.5 Individual harmonic distortion and total harmonic distortion 19 2.4 SUMMARY 19
POWER DEFINITIONS IN THE PRESENCE OF DISTORTION 21
3.1 INTRODUCTION 21 3.2 POWER IN CIRCUITS WITH SINUSOIDAL VOLTAGE AND CURRENT WAVEFORMS 22 3.2.1 Power in single frequency circuits 22 3.2.2 Separate power components in a simple circuit 22 3.2.3 Real power as function of time 23 3.2.4 Average active power 25 3.2.5 Imaginary power as function of time 26 3.2.6 Average imaginary power 27 3.2.7 Complex power as function of time. 27 3.2.8 Conclusions 28 3.3 FUNDAMENTAL CONSIDERATIONS FOR THE DEFINITION OF POWER IN CIRCUITS WITH DISTORTED VOLTAGES AND CURRENTS 28 3.3.1 Fundamental considerations in the definition of power in multi frequency circuits 28 3.3.2 Real power as function of time 29 viii
3.3.3 Average power 30 3.3.4 Imaginary power as function of time 31 3.3.5 Average imaginary power 32 3.3.6 Joint complex power 33 3.4 THE APPROACH ACCORDING TO BUDEANU 34 3.5 THE IEEE WORKING GROUP DEFINITIONS 37 3.5.1 Single-phase Relations 37 3.5.2 Three-phase Relations 39 3.6 THE CZARNECKI DEFINITIONS 39 3.6.1 Single-phase Relations 40 3.6.2 Three-phase Relations 44 3.6.3 Decomposition of load current with source impedance not zero 45 3.7 SUMMARY 46
DISTORTION COMPENSATION 47
4.1 FUNDAMENTAL CONCEPTS 47 4.1.1 Introduction 47 4.1.2 Frequency-domain view 47 4.2 COMPENSATION PRINCIPLES 49 4.2.1 Local and ambient distortion 49 4.2.2 Parallel distortion compensation 51 4.2.3 Series distortion compensation 52 4.2.4 Combined parallel and series compensation 53 4.3 GENERIC COMPENSATOR TYPES 54 4.4 PASSIVE FILTERS 55 4.5 DYNAMIC COMPENSATION 56 4.5.1 Basic considerations 56 4.5.2 Voltage-fed and current-fed topologies 56 4.5.3 Dynamic compensator control 58 4.6 HYBRID COMPENSATOR TOPOLOGIES 66 4.6.1 Series hybrid compensator topology 66 4.6.2 Parallel hybrid compensator topology 67 4.6.3 Series-dynamic hybrid filter topology 68 4.6.4 Series-parallel hybrid topologies 69 4.6.5 Resonance suppression 69 4.7 SUMMARY 70
HARMONIC SUPERPOSITION MODELLING 71
5.1 INTRODUCTION 71 ix
5.2 THREE-PHASE MODELLING TECHNIQUES 71 5.2.1 Compound admittances 72 5.2.2 Nodal matrix 74 5.2.3 Primitive network 74 5.2.4 Connection matrix 75 5.3 SUPPLY AUTHORITY SOURCE 76 5.4 TRANSMISSION SYSTEM COMPONENTS 77 5.5 6-PULSE AC-DC CONVERTER LOADS 79 5.6 PASSIVE FILTER MODEL 80 5.7 DYNAMIC COMPENSATOR 82 5.8 SUMMARY 82
PRACTICAL VERIFICATION OF THE SUPERPOSITION MODEL 83
6.1 INTRODUCTION 83 6.2 MEASUREMENT 1 AT THE PCC WITHOUT COMPENSATION 84 6.3 MEASUREMENT 2 - AT THE PCC WITH PASSIVE FILTER AT THE LOAD 86 6.4 MEASUREMENT 3 - SERIES INDUCTOR ON SUPPLY-SIDE OF FILTER 88 6.5 MEASUREMENT 4 - DYNAMIC COMPENSATOR AT THE LOAD 90 6.6 SUMMARY 92
MODELLING FOR ALTERNATIVE COMPENSATION SCHEMES AND POWER THEORIES 93
7.1 INTRODUCTION 93 7.2 THE MLISX4 PLANT NETWORK TOPOLOGY 94 7.3 MODELLING RUNS 96 7.4 CIRCUIT CONFIGURATION 1 96 7.4.1 Configuration 1 - Modelling Run 1 - Measuring at load, node 14. 96 7.4.2 Configuration 1 - Modelling Run 2 - Measuring at node 3 100 7.4.3 Configuration 1 - Modelling Run 3 - Measuring at node 2 102 7.4.4 Configuration 1 - Modelling Run 4 - Measuring at node 16 106 7.4.5 Summary of the measurements of circuit configuration 1 109 7.5 CIRCUIT CONFIGURATION 2 110 7.5.1 Configuration 2 - Modelling Run 5 - Measuring at the primary side of the step-down transformer, node 2 110 7.5.2 Configuration 2 - Modelling Run 6 - Measuring at node 14 113 7.5.3 Summary for circuit configuration 2 measurements 115 7.6 CIRCUIT CONFIGURATION 3 116 7.61 Configuration 3 - Modelling Run 7 - Measuring at node 14. 116 7.62 Configuration 3 - Modelling Run 8 - Measuring at node 7 119 7.6.3 Configuration 3 - Modelling Run 9 - Measuring at node 3 121 7.6.4 Configuration 3 - Modelling Run 10 - Measuring at node 2 124 7.6.5 Summary of circuit configuration 3 126 7.7 CIRCUIT CONFIGURATION 4 127 7.7.1 Configuration 4 - Modelling Run 11 - Measuring at node 14 127 7.7.2 Configuration 4 - Modelling Run 12 - Measuring at node 3 130 7.7.3 Configuration 4 - Modelling Run 13 - Measuring at the primary side of the step-down transformer, node 2 .132 7.7.4 Summary for the measurements in circuit configuration 4 135 7.8 CIRCUIT CONFIGURATION 5 135 7.8.1 Configuration 5 - Modelling Run 14 - Measuring at load, node 14. 135 7.8.2 Configuration 5 - Modelling Run 15 - Measuring at node 3 138 7.8.3 Configuration 5 - Modelling Run 16 - Measuring at node 2 140 7.8.4 Summary for circuit configuration 5 142 7.9 CIRCUIT CONFIGURATION 6 142 7.9.1 Configuration 6 - Modelling Run 17 - Measuring at node 16. 142 7.9.2 Configuration 6 - Modelling Run 18 - Measuring at node 3 145 7.9.3 Summary for circuit configuration 6 147 7.10 CIRCUIT CONFIGURATION 7 147 7.10.1 Configuration 7 - Modelling Run 19 - Measuring at node 2 148 7.10.2 Configuration 7 - Modelling Run 20 - Measuring at load, node 16. 150 7.10.3 Summary for circuit configuration 7 153 7.11 CONCLUSIONS REGARDING MODELLING CONFIGURATIONS 153 7.12 ALTERNATIVE POWER THEORIES 155 7.12.1 Budeanu 156 7.12.2 Czarnecki 157 7.12.3 IEEE power definitions 158 7.13 MEASUREMENTS AT THE PCC WITH LASERS IN OPERATION 158 7.14 SUMMARY 162
8. CONCLUSION AND EVALUATION 163
8.1 INTRODUCTION 163 8.2 FIELD COVERED AND RESULTS OBTAINED 163 8.3 EVALUATION 165 8.4 RECOMMENDATIONS 166
ANNEXURE A HARMONIC SUPERPOSITION MODELLING OF THE AC/DC CONVERTER Al
A.1 INTRODUCTION Al xi
A.2 DEFINITIONS Al A.3 MODELLING THE PULSES IN THE TIME DOMAIN A2 A.3.1 Positive-going pulses A2 A.3.2 Negative-going pulses A3 A.4 FOURIER TRANSFORMATION A4 A.4.1 Positive- going pulses A4 A.4.2 Negative-going pulses A6 A.4.3 Synthesis of the complete line-line current waveform: A7 A.4.4 Three-phase representation A8
ANNEXURE B - HARMONIC SUPERPOSITION MODEL B1
INTRODUCTION B1 SPECIFICATION OF UNITS B2 PER-UNIT SPECIFICATIONS B2 DISPLAY CONTROL B4 SUPPLY AUTHORITY EQUIVALENT SOURCE B5 TRANSFORMER MODELS B8 CABLE MODELS B11 PASSIVE FILTER MODELS B15 VOLTAGE VECTOR B19 SINGLE LASER MODEL BUS ADMITTANCE MATRIX B20 ELEVEN LASER BUS ADMITTANCE MATRIX B21 CURENT VECTOR FOR SINGLE LASER SYSTEM B25 CURRENT VECTOR FOR ELEVEN LASER SYSTEM B25 CALCULATIONS FOR TIME WAVEFORMS B27 LASER CURRENT INPUT PARAMETERS B27 GENERAL DISPLAY SECTION B31 MEASUREMENT RESULTS ON EXPERIMENTAL SETUP B31 NODE-NODE CURRENT CALCULATIONS B32 POWER CALCULATIONS B34 Budeanu distortion power B34 Czarnecki power definitions B35 IEEE power calculations B41 CALCULATIONS FOR TIME WAVEFORMS B45
REFERENCES R1 xii
LIST OF SYMBOLS AND DEFINITIONS a Voltage phase angle /3 Current phase angle 'an) The nth order lock admittance matrix Susceptance Susceptance for the different harmonic components Bne Czarnecki equivalent three-phase susceptance at harmonic order-n Db Budeanu distortion power Dg Czarnecki generated power Ds Czarnecki scattered power f(t) Time-dependent rms complex space-vector Fe Effective value offn(t) Fn Harmonic phasor of order-n fn(t) Time-dependent rms complex space-vector of the nth harmonic component Conductance Ge Fryze equivalent conductance Gn Conductance for the different harmonic components Gne Czarnecki equivalent three-phase conductance at harmonic order-n h(n) nth order symmetric component operator I Current phasor I rms current i (t) * Conjugate time-dependent rms scaled complex current 1(n) Current vector of order -n and dimension k i(t) Time-dependent rms scaled complex current i(t) Time-dependent scalar current II Fundamental current phasor II rms fundamental current rms value of the positive-sequence fundamental current IQ rms active current IA Primary generated three phase current ia(t) Time-dependent active current 113 Load generated three phase current Ig Czarnecki rms generated current is(t) Time-dependent Czarnecki generated current Ih Czarnecki rms unbalanced current 111 rms harmonic current ih(t) Time-dependent Czarnecki unbalanced current IHD Individual harmonic distortion Inc Czarnecki equivalent three-phase current I„ nth order harmonic current phasor in(t) Time-dependent rms scaled complex current at order-n in(t) nth order simple-harmonic, time-dependent complex current /„* Conjugate current phasor of order-n inR nth harmonic phase R current Ins Source harmonic current phasor of order-n. Ins nth harmonic phase S current InT nth harmonic phase T current Ire Czarnecki rms reactive current Time-dependent Czarnecki reactive current Is Czarnecki rms scattered current Time-dependent Czarnecki scattered current k Upper range of the order of harmonic summation n Harmonic order Total non-active power NA Set of orders n of harmonics for which P 0 NB Set of orders n of harmonics for which Pn<0 NH Total harmonic non-active power Average active power p(t) Real power as function of time p(t) Time-dependent real power in the nth harmonic component pg(t) Time-dependent real power in conductance PH Total harmonic active power ph(t) Time-dependent real power in susceptance P, Scalar-joint value of the active power in the harmonic components P, Average active power at harmonic order-n P„ Harmonic active power of order-n Q Average imaginary (reactive) power q(t) Imaginary power as function of time (time-dependent imaginary power) qb(t) Imaginary time-dependent power supplied to the susceptance qg(t) Imaginary time-dependent power supplied to the conductance Scalar joint value of the reactive power in the harmonic components Qn Reactive power at harmonic order-n Qn Harmonic reactive power of order-n q„(0 Time-dependent imaginary power in the nth harmonic component Qrc Czarnecki reactive power Apparent power s(t) Time-dependent complex or apparent power Fundamental apparent power Si+ Positive-sequence fundamental apparent power Sel Three-phase fundamental apparent power SF Forced apparent power SH Total harmonic apparent power Vector joint complex power
Sjh Joint complex harmonic power SN Non-fundamental apparent power S, Harmonic apparent power of order-n Sm„ Harmonic phasor apparent power of unequal harmonic orders Unbalanced fundamental apparent power Time T Period of fundamental component THD Total harmonic distortion Voltage phasor rms voltage u(t) Time-dependent rms scaled complex voltage u(t) Scalar time-dependent voltage Ul Fundamental voltage phasor rms fundamental voltage rms value of the positive-sequence fundamental voltage UH rms harmonic voltage Czarnecki equivalent three-phase voltage U(n) Voltage vector of order-n and dimension k Un Harmonic voltage phasor of order-n Un(t) Tithe-dependent rms scaled complex voltage at order-n un(t) nth order simple-harmonic, time-dependent complex voltage x v
Un* Conjugate voltage phasor of order-n UnR nth harmonic phase R voltage Uns nth harmonic phase S voltage U;„, Source harmonic voltage phasor of order-n. UnT nth harmonic phase T voltage V11 Passive filter 11th harmonic admittance Y5 Passive filter 5th harmonic admittance Y7 Passive filter 7th harmonic admittance Complex admittance of order-n Yp Passive filter parallel admittance
Ypnm Primitive admittance matrix Vs Passive filter series admittance Zim(n) Complex impedance between nodes / and m at harmonic order-n On Phase-angle of the nth harmonic phasor Fundamental harmonic frequency XV
LIST OF FIGURES Figure 1.2.1 - Power diagram of R2P2 unit 4 Figure 1.2.2 - Block diagram of laser power supply 5 Figure 1.3.1 - Line-current waveform of 6-pulse converter at full load 6 Figure 1.3.2 - Normalised harmonic line current spectrum for the 6-pulse converter at full load 6 Figure 1.3.3 and Figure 1.3.4 - Respective relationships of THD and output DC-voltage against inductor Lr value 7 Figure 1.4.1 - Admittance diagram for the plant 8 Figure 2.2.1 - Complex s-plane representation of voltage u(t) and current i(t) 14 Figure 2.3.1 - Graphical instantaneous representation of the voltage u(t) 16 Figure 3.2.1 - Parallel circuit with conductance and susceptance 22 Figure 3.4.1 - Representation of active, reactive and distortion power in cartesian space 36 Figure 3.6.1 — Single-phase load current decomposition according to Czarnecki 40 Figure 4.1.1 - Simplified transmission system model 48 Figure 4.1.2 - Circuit for fundamental and order-n harmonics 49 Figure 4.2.1 - Principle of parallel distortion compensation 51 Figure 4.2.2 - Fundamental and harmonic parallel compensation 52 Figure 4.2.3 - Principle of series compensation 53 Figure 4.2.4 - Combined parallel and series compensation 54 Figure 4.4.1 - Single-phase example of a passive filter network 55 Figure 4.5.1 - Essential power electronic structure of a voltage-fed dynamic compensator 57 Figure 4.5.2 - Basic current-fed compensator 58 Figure 4.5.3 - Simple single-phase steady-state compensator 60 Figure 4.5.4 - Current and voltage signals for dynamic compensator 61 Figure 4.5.5 - Dynamic compensator reference signal generator 61 Figure 4.5.6 - Steady-state per-phase control strategy for a dynamic compensator 64 Figure 4.6.1 - Series hybrid compensator topology 67 Figure 4.6.2 - Parallel hybrid compensator topology 67 Figure 4.6.3 - Series-dynamic compensator topology 68 Figure 4.6.4 - Series-parallel hybrid topology 69 Figure 4.6.5 - Resonance suppression example 70 Figure 5.2.1 - Admittance of a three-phase series element 73 Figure 5.2.2 - Single compound admittance representation for three series elements 73 Figure 5.2.3 - Actual connected network 74 Figure 5.2.4 - Primitive three-phase network with compound admittances 75 Figure 5.3.1 - Equivalent circuit for supply authority model 76 Figure 5.4.1 - Primitive network for three-phase transformer 78 Figure 5.5.1 - 6-pulse AC-DC converter line current waveform 79 xvi
Figure 6.1.1 - Per-phase admittance-diagram for dummy-load measurement setup 84 Figure 6.2.1- Measurement 1 - PCC without compensation 84 Figure 6.2.2 - Measurement 1 - Comparison of modelled and measured results 85 Figure 6.3.1- Measurement 2 - PCC with passive filter at load 86 Figure 6.3.2 - Measurement 2 - Comparison of modelled and measured results 86 Figure 6.4.1- Measurement 3 - At PCC with series inductor installed on supply-side of Passive Filter 89 Figure 64.2 Measurement 3 - Comparison of modelled and measured results 89 Figure 6.5.1 - Measurement 4 - At PCC with dynamic filter connected as shown in admittance-diagram 91 Figure 65.2 - Measurement 4 - Comparison of modelled and measured results 91 Figure 7.2.1 - MLISX4 power supply network to lasers 95 Figure 7.4.1 - Measuring at load, node 14 97 Figure 7.4.2 - Line-line voltage and current at the R 2P2 terminals 97 Figure 7.4.3 - Normalised harmonic voltage spectrum for node 14 voltage 98 Figure 7.4.4 - Normalised harmonic current spectrum for node 7-14 current 98 Figure 7.4.5 — Measuring at the secondary side of the transformer 100 Figure 7.4.6 - Line-line voltage and current at the secondary side of the transformer 100 Figure 7.4.7 - Normalised harmonic voltage spectrum 101 Figure 7.4.8 - Normalised harmonic current spectrum 101 Figure 7.4.9 — Measuring on the primary side of step-down transformer 103 Figure 7.4.10 - Line-line voltage and current at the transformer primary 103 Figure 7.4.11 - Normalised harmonic voltage spectrum 104 Figure 7.4.12 - Normalised harmonic current spectrum 104 Figure 7.4.13 - Measuring on the input side of an R3P2 power supply unit 107 Figure 7.4.14 - Line-line voltage and current at the R 3P2 terminals 107 Figure 7.4.15 - Normalised harmonic voltage spectrum 108 Figure 7.4.16 - Normalised harmonic current spectrum at the R 3P2 terminals 108 Figure 7.5.1 - Measuring at the primary side of the transformers 111 Figure 7.5.2 - Line-line voltage and current primary of the transformer 111 Figure 7.5.3 - Normalised harmonic voltage spectrum 112 Figure 7.5.4 - Normalised harmonic current spectrum 112 Figure 7.5:5 - Measuring at the laser load 113 Figure 7.5.6 - Line-line voltage and current at the R2P2 terminals 114 Figure 7.5.7 - Normalised harmonic voltage spectrum 114 Figure 7.5.8 - Normalised harmonic current spectrum at R2P2 terminals 114 Figure 7.6.1 - Measuring at the laser load, node 14 117 Figure 7.6.2 - Line-line voltage and current at the R2P2 terminals 117 Figure 7.6.3 - Normalised harmonic voltage spectrum 118 xvii
Figure 7.6.4 - Normalised harmonic current spectrum at R2P2 terminals 118 Figure 7.6.5 - Measuring at the input side of the series inductor to the laser load 119 Figure 7.6.6 - Line-line voltage and current at the R2P2 terminals 120 Figure 7.6.7 - Normalised harmonic voltage spectrum 120 Figure 7.6.8 - Normalised harmonic current spectrum 120 Figure 7.6.9 - Measuring at the secondary side of step-down transformer 1 122 Figure 7.6.10 - Line-line voltage and current at the PCC (node 5) 122 Figure 7.6.11 - Normalised harmonic voltage spectrum 123 Figure 7.6.12 - Normalised harmonic current spectrum at R2P2 terminals 123 Figure 7.6.13 - Measuring at the primary side of step-down transformers 124 Figure 7.6.14 - Line-line voltage and current at the R2P2 terminals 124 Figure 7.6.15 - Normalised harmonic voltage spectrum 125 Figure 7.6.16 - Normalised harmonic current spectrum 125 Figure 7.7.1 - Measuring on the laser load 128 Figure 7.7.2 - Line-line voltage and current at the R2P2 terminals 128 Figure 7.7.3 - Normalised harmonic voltage spectrum 129 Figure 7.7.4 - Normalised harmonic current spectrum at R2P2 terminals 129 Figure 7.7.5 - Measuring at the secondary side of transformer 1 130 Figure 7.7.6 - Line-line voltage and current at the transformer secondary node 3 131 Figure 7.7.7 - Normalised harmonic voltage spectrum 131 Figure 7.7.8 - Normalised harmonic current spectrum at R2P2 terminals 131 Figure 7.7.9 - Measuring at the primary side of the step-down transformers 133 Figure 7.7.10 - Line-line voltage and current at the R2P2 terminals 133 Figure 7.7.11 - Normalised harmonic voltage spectrum 134 Figure 7.7.12 - Normalised harmonic current spectrum 134 Figure 7.8.1 - Measuring at the load 136 Figure 7.8.2 - Line-line voltage and current at the R2P2 terminals 136 Figure 7.8.3 - Normalised harmonic voltage spectrum 137 Figure 7.8.4 - Normalised harmonic current spectrum at R2P2 terminals 137 Figure 7.8.5 - Measuring at the secondary of step-down transformer I 138 Figure 7.8.6 - Line-line voltage and current at the R2P2 terminals 138 Figure 7.8.7 - Normalised harmonic voltage spectrum 139 Figure 7.8.8 - Normalised harmonic current spectrum 139 Figure 7.8.9 - Measuring at the primary side of the step-down transformers 140 Figure 7.8.10 - Line-line voltage and current at the R2P2 terminals 141 Figure 7.8.11 - Normalised harmonic voltage spectrum 141 Figure 7.8.12 - Normalised harmonic current spectrum 141 Figure 7.9.1 - Measuring at a double-laser load 143 xviii
Figure 7.9.2 - Line-line voltage and current at the R2P2 terminals 143 Figure 7.9.3 - Normalised harmonic voltage spectrum 144 Figure 7.9.4 - Normalised harmonic current spectrum at R2P2 terminals 144 Figure 7.9.5 - Measuring at the secondary side of step-down transformer 1 145 Figure 7.9.6 - Line-line voltage and current at the transformer secondary for modelling run 18 145 Figure 7.9.7 - Normalised harmonic voltage spectrum 146 Figure 7.9.8 - Normalised harmonic current spectrum 146 Figure 7.10.1 - Measuring at the primary side of the step-down transformers 148 Figure 7.10.2 - Line-line voltage and current at the transformer primary option 19 148 Figure 7.10.3 - Normalised harmonic voltage spectrum 149 Figure 7.10.4 - Normalised harmonic current spectrum 149 Figure 7.10.5 - Measuring at load 151 Figure 7.10.6 - Line-line voltage and current at the R2P2 terminals 151 Figure 7.10.7 - Normalised harmonic voltage spectrum 152 Figure 7.10.8 - Normalised harmonic current spectrum at R2P2 terminals 152 Figure 7.13.1- - Line-Line voltage and current at node 5 159 Figure 7.13.2 - Comparison of modelled and measured results 159 Figure 7.13.3 - Normalised harmonic voltage spectrum - modelled 160 Figure 7.13.4 - Normalised harmonic voltage spectrum - measured 160 Figure 7.13.5 - Normalised harmonic current spectrum - measured 160 Figure 7.13.6 - Normalised harmonic current spectrum - measured 160 Figure A 1 - Time domain representation of line-line A3 Figure A 2 - Time domain portion of negative portion A4 Figure A 3 —Positive portion of pulse plotted from frequency domain A6 Figure A 4 — Negative-going pulse portion synthesised from Fourier components A7 Figure A 5 — Effect of changing a on the line-line converter current A8 Figure A 6 — Typical line current profile together with line-neutral voltage for a 6-pulse conveter A9 Figure B 1— Twelve laser layout B3 xix
LIST OF TABLES Table 1.2.1 - Salient Specifications for one R 2P2 Module 5 Table 1.3.1 — Salient parameters of R 2P2 input 7 Table 5.6.1 - Characteristics offilter 81 Table 6.2.1 - Measurement 1 - PCC without compensation 85 Table 6.3.1 - Measurement 2 - PCC with passive filter at load 87 Table 6.4.1 - Measurement 3 - PCC with series inductor on supply-side ofpassive filter 90 Table 6.5.1 - Measurement 4 - Dynamic compensator at node 4 92 Table 7.1.1 - List of circuit configurations modelled 93 Table 7.4.1 - Modelling run 1 - Measurements at the R 2P2 terminals for Configuration 1 99 Table 7.4.1 - Modelling run 2 - Circuit Configuration 1 - Measuring at the secondary side of the step- down transformer 102 Table 7.4.1 - Modelling run 3 - Measurements at the primary side of the step-down transformer 104 Table 7.4.1 - Modelling run 4 - Measurements at the R 3P2 terminals 108 Table 7.5.1 - Modelling run 5 - Measurements at the primary side of the step-down transformer 112 Table 7.5.1 - Modelling run 6 - Measurements at the R 2P2 terminals for Configuration 2 115 Table 7.6.1 - Modelling run 7 - Measurements at the R 2P2 terminals for Configuration 3 118 Table 7.6.1 - Modelling run 8 - Measurements at the at the R2P2 terminals configuration 3 120 Table 7.6.1 - Modelling run 9 - Measurements at the FCC (node 5) for Configuration 3 123 Table 7.6.1 - Modelling run 10 - Measurements at the transformer primary for configuration 3 125 Table 7.7.1 - Modelling run 11 - Measurements at the R 2P2 terminals for Configuration 4 129 Table 7.7.1 - Modelling run 12 - Measurements at the transformer secondary for Configuration 4 132 Table 7.7.1 - Modelling run 13 - Measurements at primary of transformer for Configuration 4 134 Table 7.8.1 - Modelling run 14 - Measurements at the R 2P2 terminals for Configuration 5 137 Table 7.8.1 - Modelling run 15 - Measurements at transformer secondary for Configuration 5 139 Table 7.8.1 - Modelling run 16 - Measurements at the primary of the transformer for Configuration 5 141 Table 7.9.1 - Modelling run 17 - Measurements at the R 3P2 terminals for Configuration 6 144 Table 7.9.1 - Modelling run 18 - Measurements at the secondary of the transformer for Configuration 6 146 Table 7.10.1 - Modelling run 19 - Measurements at the transformer primary for configuration 7 149 Table 7.10.1 - Modelling run 20 - Measurements at the R 2P2 terminals for configuration 7 152 Table 7.11.1 — Economic comparison of the different compensation topologies 155 Table 7.12.1 - Reactive powers for the different modelling runs 156 Table 7.12.1 - Harmonic active power components 157 Table 7.13.1 - Modelling and Measurements at PCC 161 Table 7.13.2 - Modelling and measurements results for the three power theories 161 1
1. THE MLIS PLANT ENVIRONMENT
1.1 INTRODUCTION
The Atomic Energy Corporation of South Africa (AEC) decided in the late seventies that it would invest in the Molecular Laser Isotope Separation (MLIS) process for uranium enrichment 1 .2 Research and Development has progressed, at the time of writing, to the point at which laser excitation systems, on the electrical side, have individually clocked multiple hours during endurance testing and are commercially and technically ready to be exploited in a plant environment.
The next stage in the enrichment program will be the composite testing of all the equipment (excitation systems, control schemes, lasers and isotope processing components) in a pilot plant. This pilot plant is currently in its final stages of preparation and will be used to demonstrate the engineering feasibility of the project and to highlight critical parameters that govern the capital and operating cost structures of such a plant. The final stage, planned upon the successful testing and evaluation of equipment in the pilot plant, will be the construction of a full-scale MLIS plant.
A key element in the development of the MLIS process is the industrial CO 2 transversely excited atmospheric (TEA) pulsed lasers that furnish the excitation in the separation process. These lasers have average optical powers in the kilowatt range and peak pulsed optical power outputs in the megawatt range3. Research in the development of these lasers and their associated electrical pulse excitation systems began in 1985 and has led to the establishment of a new technology of lasers and laser excitation systems 4.
The electrical pulse excitation systems, also referred to as pulsers, consist of a number of different elements. The resonant regulating pulse power (R 2P2) supply units, takes power from the mains network and carries out primary pulse conversion. The output pulses from the R 2P2 units are then further processed by thyratron pulsers or multi-stage pulse compressors 6, before they are compatible with the ultra-high peak current and ultra-short pulse durations required by the pulsed lasers.
These R2P2 units used for primary pulse modulation, employ intermediate DC-stages, in like manner to that of other power-electronic equipment. In this case, the DC-voltage levels are required to be very accurately and continuously adjustable to furnish highly repeatable and adjustable output pulse voltage levels.
Several options exist for the conversion of the supply network three-phase AC-power to DC. Conventional line-commutated phase control is presently employed, but the so-called "line- 2 friendly" or force-commutated PWM-modulated converter is also under investigation as an alternative'.
The use of phase-controlled AC-DC converters in the R2P2 front-ends leads to distortion in the currents drawn by these units from the AC network. On a smallscale, this distortion will not be a problem and can even be tolerated in a number of instances. At the scale envisaged for the main plant, which will extend to several hundred megawatt, and even in the case of the pilot plant, the distortion is expected to introduce problems that will extend from malfunctioning of equipment to infringements of supply-authority standards' for consumer harmonics 9.
Aside from the alternative supply-friendly front-end options, that may come in a number of forms and employ alternative topologies and control schemes, a number of external compensation options also exist for the mitigation of the distortion by the phase-controlled converter. The conventional line-commutated 6-pulse converter features very prominently in the present investigation, because of its relatively low cost, simplicity of operation, and its proven track record.
Distortion compensation technology is still in its infancy. Although passive filter networks have been employed for some time now to reduce distortion in supply currents to non-linear plants, the problems they introduce are many. They are most suited for installation on the H.V. side of consumer transformers and, because they tend to react with each other, multiple installations in close electrical proximity are usually impractical. Because of the spread-out nature of a large number of identical, relatively small distortion loads in the envisaged pilot and main plants, compensation of one form or another at the pulse power supply locations themselves may have to be considered to maintain acceptable limits of distortion in the feeders that will be spread throughout the plant.
Dynamic compensation offers a technical solution, but at the large power levels considered for the main plant, an expensive one. The total compensation solution, if that option is exercised, obviously lies in a hybrid compensation topology in which passive and dynamic compensation is integrated to yield both the most cost-effective and best technical results. This solution, and a search for the optimal combination of alternatives, form the theme of this study, and to which the answers are sought in the following chapters.
A comprehensive study of the problem outlined above necessitates extensive modelling. In this modelling, distorted voltage and current waveforms will have to be treated throughout. Established and concise definitions that are used in well-behaved sinusoidal systems are not applicable and other power definitions, that cater for the conditions of distortion, must be resorted to. Because this study concerns a practical plant and practical measurements and assessments, it also affords the ideal opportunity for the parallel testing and evaluation of the most prominent 3 power theories presently under consideration, namely that of Budeanu 15, the IEEE Working Group definitions3' and the Czarnecki power definitions 17.
Because other literature, 10" 1 already gives a basic background to the envisaged plant layout and operation, further explanations of the process will be avoided here. Schematics of the plant electrical layout furnishing details of the number and characteristics of the topographically distributed pulse power supply loads on the electrical distribution system will suffice here for the purpose of modelling and measurement.
The chapter will begin with a description of the R 2P2 units and their front-ends, that will constitute the main power consumers in the plant. Next, the simplified distribution layout for the envisaged plant will be sketched and typical load characteristics, modelled and measured, will be given.
1.2 R2P2 PULSE POWER SUPPLY CONFIGURATION
The pulsers driving the CO 2 pulse lasers in the programme, are driven by R2P2 modules that take three-phase AC power from the supply network at a nominal voltage of 400 V. Because the aggregate output pulse power requirement per laser exceeds the capacity of the individual modules, more than one module is used in a power supply configuration to time-multiplex the output pulses into the input of each laser pulser. In a typical case in which the laser requires to operate at 1 kHz, two R 2P2 modules can be employed at 500 Hz each or four modules can be employed at 250 Hz each. In practice, more than the minimum necessary number of modules will be used for high reliability and faulty modules will be automatically switched out of the pulse power supply upon failure without interrupting the operation of the power supply as a whole.
The R2P2 units at present employ conventional 6-pulse line-commutated phase-controlled converters in their front-ends to furnish an intermediate, continuously controllable DC-source that is in turn used for the generation of power pulses through primary pulse conversion by successive resonant energy transfer stages. The basic, simplified power diagram of an R 2P2 unit is shown in simplified form in Figure 1.2.1.
As shown, the 3-phase mains supply of 400 V is fed to a 6-pulse phase-controlled bridge in the input end of the R2P2 unit. The rectified output of the bridge is smoothed by the DC-choke inductor Lr and reservoir capacitor C I and is delivered to the power circuit at a nominal DC- voltage that can be controlled by adjustment of the firing angle of the phase-control thyristors from very low values to the maximum value of about 530 V (no-load value). 4
R 3
Rb Dm1 1: 2t>11-1 CI b D m 2 ->F C o Dm3 -[>[ L a MMI Tb 1 Dm4 Lr Lm TX • Cr
I I • D De
Figure 1.2.1 - Power diagram of R2P2 unit
Operation of the R2P2 unit is best explained by designating of three successive but separate phases' 2 to its operation.
The three phases are respectively the charging phase, the transfer phase and the regulating (also called the de - queuing) phase. The charging phase is initiated by firing thyristor Tb and is employed to charge the primary capacitor C o from the rectified supply by resonant energy transfer in the loop formed by Cr, La, Tb Co and the primary winding of transformer TX.
During the transfer phase, energy is resonantly transferred (through the step-up pulse transformer TX) from the primary capacitor C o, through the leakage inductance of the pulse transformer TX and the multiplexing diode, to the secondary (high voltage) capacitor C, in Figure 1.2.1.
The regulating phase is responsible for interrupting the charging phase under controlled conditions by the triggering of thyristor Ta, to terminate the delivery of charge to the primary capacitor Co, thereby accurately controlling the pulse voltage (and the energy) that is delivered to the load with each pulse.
Triggering of thyristor Ta during the second half of the charging transfer cycle, when the current gradient through La is negative, will turn Tb off, thereby terminating the transfer of charge to C o. In the same way, triggering of Tb when Co is discharged, after the next transfer cycle, will initiate the next charging cycle and terminate the current flow in the loop formed by L a and Ta. 5 The relevant specifications for the present configuration of the R 2P2 modules, used in the present study, is given in table 1.2.1.
Table 1.2.1 - Salient Specifications for one R 2P2 Module
DESCRIPTION OF PARAMETER UNIT VALUE Peak output pulse voltage kV 24 Repeatability of output pulse voltage % 0.01 Peak output pulse energy J 20 Peak output pulse frequency Hz 500 Peak input power kW 15 Input three-phase voltage V 380 Input power factor 0.84 Input total harmonic distortion % 47 Overall efficiency % 71
The input to the pulse compressor utilises multiplexing diode D m' to Dm4 to connect to the four power supply modules. In practice three modules are used for an oscillator (first laser in a laser train, and delivering only in the region of 17 J of energy) and four for an amplifier. A high-power laser (80 kW) is also under development and this laser may employ as many as eight power modules each delivering 40 J at 250 Hz. The multiplexing configuration is shown in block- diagrammatic form in Figure 1.2.2.
Power supplies 1 3 4-8
Pulse circuit
Laser load
Figure 1.2.2 - Block diagram of laser power supply
The pilot plant laser train will initially employ 11 lasers, operating at 1 kHz. This train will make use of 41 R2P2 modules, of which 22 will be operating at full load, with each drawing 15 kW of active power from the supply, or drawing a total of 330 kW in total. When fully operational, the chain will employ all 41 modules at 2 kHz, requiring a total input active power of 615 kW. The layout may change in future to accommodate the best achievable process parameters. These
6 parameters will be refined during on-going development of the process and the establishment of the most optimal laser configurations, powers and other parameters. The electrical power requirements and plant layout will ultimately depend on the latest findings.
1.3 POWER SUPPLY FRONT - END TOPOLOGY
There is a space constraint in the maximum size of choke inductor L r that can be installed in the R2P2 units. The maximum size of choke that is now capable of accommodating the continuous load with air cooling is 3.2 mH. This inductance is only just capable of furnishing continuous output current under full-load conditions at 1 kHz. At full power, the line current amplitudes yield continuous output current as shown in Figure 1.3.1.
400 320 240 160 Line 80 current 0 -80 -160 -240 -320
-400 -0.03 -0.02 -0.01 0 0.01 0.02 Time in seconds
Figure 1.3.1 - Line-current waveform of 6-pulse converter at full load
The normalised harmonic current spectrum for the above case, assuming steady-state conditions, is shown in Figure 1.3.2.
40 36 32 Harmonic 28 current as a 24 percentage of 20 the fundamental 16 12 8 4 0 U 5 10 15 20 Harmonic order
Figure 1.3.2 - Normalised harmonic line current spectrum for the 6-pulse converter at full load 7 The most important per-phase parameters, assuming phase symmetry for the above conditions are given in table 1.3.1.
Table 1.3.1 — Salient parameters of R2P2 input'
Equivalent total node voltage 239.03 V Equivalent fundamental node voltage: 239.013 V RMS load current: 84.171 A RMS fundamental component of line current: 78.657 A RMS harmonic component of line current: 29.964 A Apparent power: 60.358 kVA Fundamental apparent power: 52.313+21.08 j kVA Current THD: 38.095 % Voltage THD: • 1.196 %
A PSpice simulation of the 6-pulse converter against different values of the DC-choke value furnished the line current THD and output DC voltages, shown respectively in Figure 1.3.3 and Figure 1.3.4.
Figure 1.3.3 and Figure 1.3.4 - Respective relationships of THD and output DC-voltage against inductor L,. value
As shown in Figure 1.3.3 and Figure 1.3.4 the distortion on the line current drops from 120% to 40% and the line voltage from 500 to 445 when L4 increases from 300 1.tH to 3.2 mH. The module operates in the discontinuous mode with Lr= 300 11H and changes to the continuous mode when L1 approaches 3 mH. Both THD and output DC voltage tend to their asymptotic values when increasing I., beyond 3 mH. The choice of inductance for L r of 3.2 mH is therefore a
i The fundamental and harmonic components are separated in accordance with the theory discussed in chapter 3.
8 practical one. The line-current THD and output DC voltage behaviour is also borne out in the literature 13.
1.4 PLANT LAYOUT
The simplified, per-unit, line diagram of the proposed electrical installation that will typically serve the plant is shown in figure 1.4.1.
r-NLTi t Load Single laser PASSIVE ,- at A-EBP0002 FILTER.
15 Load Double laser at PASSIVE‘ A-EBP00:13 FILTER
(-19Toir 16 Load PASSIVE t Double laser at A-EBP0004 FILTER
I-1 10 Load 000) _L 17 Single laser at A-EBP0005 PASSIVE FILTER
(-Nroon} 18 Load v I Single laser at PASSIVE A-EBP009 FILTER
Load (--,12 ino-o. Double laser at PASSIVE 1 A-EBP008 FILTER
Load r-, (130-u 20 Doublelaser at PASSIVE t A-EBP0007 FILTER
r121_(7)()0) 22 Load _I_ Single laser at PASSIVE A-EBP0006 FILTER
Figure 1.4.1 - Admittance diagram for the plant
This diagram has been derived from the electrical reactance diagrams for the installation and only equivalent compound admittances and current sources are employed. Admittances and node- current values are obtained by source transformation from their respective source, line, cable, 9 transformer and other transmission component values. Admittance values are taken from their fundamental frequency values and are scaled for other harmonic frequencies as explained later.
As the diagram shows, power is taken from a single source at 11 kV as in the single-load study above, which in fact consists of several substations in a ring network. The Norton equivalent single source represents this network of substations. A single radial feeder, consisting of a 150mm2 11 kV PEX cable of length 391 m, brings this power to the substation, where it is transformed down to 380 V by two separate 2 MVA 5,69% impedance 11 kV to 380 V AY transformers, each feeding a separate busbar. In the event of a transformer failure, a bus-coupler can link the busbars and the plant is then operated on one transformer only. This configuration will also have to be investigated separately, because of the higher levels of voltage distortion that it will bring about.
Each transformer secondary is coupled to a dynamic compensator as well as to a passive filter, at nodes 3 and 4 for transformers 1 and 2 respectively. This provision again enables these components to be enabled or inhibited as the study demands, merely by redefining the coupling admittances in the case of the passive filters, or by toggling the current that the dynamic compensators inject into the network.
Each substation busbar, located at nodes 5 and 6, supplies a set of lasers. The particular allocation of lasers is representative of the present cable installations in the building and corresponds to the information presently known with respect to the physical location of the lasers. Alterations will only introduce marginal differences in cable and can be accommodated easily when final layouts have been finalised. A number of load models cater for double lasers and other for single lasers as shown.
The R2P2 loads are supplied from the two busbar sections by twin 95 mm 2 PVCAPVC cables that are generally about 60 m long. Provision is also made in the model for accommodating a smaller passive filter at each laser, and to equip each of the filters with a series inductor. The passive filters can be optionally inserted at nodes 14 to 20 and 22 in the model to substitute the parallel admittances as shown. The series admittances between the previous nodes and the laser nodes, such as Yp7_14, permit the installation of series decoupling inductors between the passive filters and the supply cables at each laser load.
The series inductors have important functions and certain configurations may not be possible or feasible without them. Firstly, their use is mandatory when multiple filters are installed to reduce parallel resonance between the multiple filters. They may also assist in reducing harmonic current distortion in the cables, provided that the ancillaries at each laser can be powered from a point before them and provided that voltage distortion will be permitted on the inputs of the AC/DC converters themselves. 1 0
1.5 SUMMARY
The basic plant layout and electrical requirements for the laser plant was discussed in this chapter. The layout and the modelling make provision for the installation of passive and dynamic filters in different locations in the circuit. The essential loads consist of 6-pulse line-commutated phase-controlled converters operating under steady load conditions. The ensuing study, commencing in the next chapters, will comprise finding the optimal configuration and ratings of compensating equipment for this plant. This study will also ideally lend itself to an analysis of different power theories, by means of which the modelling and measurement results can be analysed. 11 2. FUNDAMENTAL DEFINITIONS
2.1 INTRODUCTION
Steady-state linear AC circuit theory is derived from the forced response of circuits when they are excited by sinusoidal forcing functions'''. The definitions of complex phasor quantities and impedances, that form the basis of this theory, however, are confined to the requirement that the forcing functions be sinusoidal and that the circuits be linear.
When the forcing functions in a linear circuit are no longer sinusoidal, or if the circuit is not linear or is time-variant, non-sinusoidal responses take place and either the currents or the voltages or both may be non-sinusoidal. This deviation from the sinusoidal form is referred to as distortion and the voltages and currents are said to be distorted. When these distorted, but periodic responses are analysed in the frequency domain, they are found to possess more than one discrete frequency, and such responses may also be referred to as multi frequency responses.
The distorted or multi-frequency behaviour of AC circuits prohibits the use of complex analysis and the defined concepts of phasors, impedances and power are no longer applicable. When the distorted responses of voltage and current are periodic and when the circuit has a steady-state response, Fourier analysis yields discrete responses in the frequency domain. Each of these responses, unique in phase and magnitude are, however, sinusoidal and analysis becomes amenable to all the first-named definitions for conventional AC circuits. This is where the similarity to conventional treatment ends, however, for each discrete frequency component. New difficulties arise when the collective or joint power components have to be calculated. To do that, new definitions are required.
Studies in the field of circuit-analysis in the presence of distortion is becoming more important every day and has developed into a major field of research, yielding a proliferation of publications in which new analytical theories are advanced. Notable work in this regard, but not necessarily in any order of importance, is that of Budeanu 15, Fryze16, Czarnecki",
Depenbrock' 8 Custers and Moore 19 , Enslin2° and Akagi and Nabae21 . Before beginning with a discussion of some of these theories, a convenient set of definitions will first be made to cover the underlying theory and to gradually lead the reader from the conventional definitions for the sinusoidal case through to multi-frequency topologies in a number of steps. 12
2.2 FUNDAMENTAL DEFINITIONS UNDER SINUSOIDAL CONDITIONS
2.2.1 Fundamental considerations
Steady-state conditions in an AC circuit here refers to periodic behaviour, in the time domain, of the time-dependent values of voltage, current and power. Transient conditions in an AC circuit refer to deviations from periodic behaviour that are relatively short-lived and disappear within a fraction of a period to a few periods. Steady-state theory, on which Fourier analysis is based, does not cater for transient behaviour and both the duration and magnitude of the transient occurrences are assumed to be negligible over the period of operation to which this study applies.
2.2.2 Complex voltage and current as a function of time
A circuit is fully characterised by the topology of its active and passive components. Active components are either voltage or current sources and passive components are those of resistance, capacitance and inductance which, collectively or individually, constitute the immittance of the circuit'''.
Under sinusoidal, steady-state conditions, the voltage and current as time functions, may be represented in a circuit respectively by:
u(t) = 42 U cos(cot + a) (2.2.1) and
i(t) = /coot + p) (2.2.2) in which U and I respectively represent the scalar rms-values of voltage and current, and a and p the phase angles.
The above, measurable, time-dependent quantities represent the values of the respective current and voltage that can be expressed in accordance with the given equations for u(t) and i(t) in the circuit at any instant t. These values of voltage and current will be referred to as the time- dependent voltage and the time-dependent current in the circuit, in preference to the designations of instantaneous voltage and instantaneous current.
Equations (2.2.1) and (2.2.2) may also be expressed in a more general form through Euler's relation as:
u(t) = Re{ U e j(6)t + a)} (2.2.3) and 13
i(t) = yL Re{ I e i(cot + 13)} (2.2.4)
The generalisation may be taken further, to represent the time-dependent quantities in equations (2.2.3) and (2.2.4) directly in their complex form as:
u(t) = U e J(cot + a)) (2.2.5) and
i(t) = / e .i(cot 13)) (2.2.6)
The quantities are now complex and do not represent the normal physically measurable values of voltage and current. They do, however, still represent time-dependent values of voltage and current and lend themselves to a number of operations that will be very convenient later in subsequent analysis. It is worthwhile to note that these complex quantities represent the more general form from which the scalar values of u(t) and i(t) can be constructed. They cannot be designated as instantaneous values. In this regard, u(t) and i(t) are only time-dependent values. Note also that u(t) and i(t) are scaled as rms-quantities. This formulation is expedient to avoid the inevitable .42 and its derivatives in power calculations. It now follows that the complex time-dependent values of equations (2.2.5) and (2.2.6) are related to the scalar values by:
u(t) = 42 Re { u(t)} (2.2.7) and
i(t) = 42 Re { i(t)} (2.2.8)
It also follows from the previous two equations (2.2.7) and (2.2.8) that
u(t) = u e ja e icot (2.2.9) and
i(t)=IeJP a Jwt (2.2.10) in which:
u(t) = U e jc" (2.2.11) and:
i(t) = / e Jwt (2.2.12) 14
The complex quantities U and I are immediately recognised to represent the respective complex rms -phasor values of voltage and current.
The sketches in Figure 2.2.1 show that u(t) and i(t) represent rms-scaled vectors of constant magnitude that rotate in an anti-clockwise direction in the complex s-plane with angular frequency co. Rotation will be centred around the origin in the absence of a DC-component, or be shifted by the appropriate value on the real axis in the presence of one. Note that the respective phase angles a and p of the voltage and the current represent the angles that these vectors make with the real positive axis at time t = 0. u t) imaginary axis imaginary axis i ( t)
t t
a real axis real axis
Figure 2.2.1 - Complex s-plane representation of voltage u(t) and current i(t)
Reflection on the above discussion will show that the measurable physical values of u(t) and i(t), as expressed by equations (2.2.7) and (2.2.8) are simply numerically equal to the projections of the rotating vectors on the real axes, multiplied by the rms-scaling factor 42 to convert the vector magnitudes to peak values.
The term complexor has been suggested for the quantities u(t) and i(t) in the literature22. The corresponding phasor voltage U and phasor current I are complex constants that do not change under steady-state conditions, but which may vary 'gradually' to take the circuit operation from one steady-state condition to the next. In accordance with the foregoing discussion, measurement is not defined during these transitions.
2.2.3 RMS value
It can be seen through inspection that the complex (RMS) value may be obtained by expressing u(t) or i(t) at t = 0:
u(0) = U or i(0) = / (2.2.13), (2.2.14)
Alternatively, it can easily be shown that the scalar (RMS) value of u(t) or i(t) may be directly calculated from the time-dependent complex value as: 15
U= qu(t)u(t) * or I =.\, (2.2.15), (2.2.16)
Because of their identical mathematical nature, both complexors and phasors conform to the definition of space vectors and identical mathematical operations are equally valid in both cases.
2.2.4 Steady-state impedance
The relationship of phasor voltage U to phasor current I in a circuit under sinusoidal conditions is a function of the steady-state or AC-impedance 14 :
Z = 11 (2.2.17)
Instead of expressing impedance through the relationship between U and I, it is analytically valid, for example, to define the AC - impedance directly in terms of the complex time- dependent voltage u(t) and current i(t). That this is valid can easily be shown by respectively substituting equations (2.2.9) and (2.2.10) for u(t) and i(t) in eq. (2.2.17) and cancelling out the e Jwt
u(t) (2.2.18) Z = i(t)
2.3 FUNDAMENTAL DEFINITIONS UNDER NON-SINUSOIDAL CONDITIONS
2.3.1 Complex voltage and current
Time-dependent, single frequency, simple harmonic basic theory was reviewed in the previous section for circuits in which the voltage and current are sinusoidal. Circuits operating under steady-state, but periodic conditions, also yield simple harmonic component waveforms under Fourier analysis. This approach is equally applicable to each separate harmonic frequency and furnishes a convenient vehicle by means of which existing power concepts and definitions can be extended to all frequencies to provide a composite distortion power theory in the frequency domain.
The analysis in the case of single-frequency circuits employ complex values. In steady-state distorted periodic circuit analysis, summations are introduced to represent all the relevant harmonic components in the equations. Consider an instantaneous complex function of a periodic steady-state voltage or current represented in general by f(t):
It can be shown by simple manipulation of the exponential Fourier series 23 that:
k fit) = fFn ekca (2.3.1) n=0 16 in which fit) and the Fn are again rms-scaled complex values, just as in the single-frequency case. Note that n ranges from 0 to k in eq. (2.3.1) in order to include a DC-value when it is present and for which F0 will be real. Because f(t) is scaled to its rms-value, it is necessary to multiply it by -42 to obtain the peak value. The value of k in the summation is usually chosen lower than oo for practical reasons, but at the sacrifice of accuracy of representation. The nth order harmonic phasor Fn is defined for the nth order harmonic, identically to that of the single-frequency case.
It follows from eq. (2.3.1) that the nth order complexor is given by:
= Fn ejncot (2.3.2)
In eq. (2.3.2), fn(t) represents a vector of constant magnitude that rotates in an anti-clockwise direction in the complex plane with angular velocity nco. The corresponding "measurable" part offn(t) is then expressed as:
fn( t) = Re{ fn(t)} (2.3.3)
Because fn(t) is sinusoidal and periodic, it follows that the complex AC-circuit theory developed in the previous section is applicable at each separate harmonic order-n.
The complex quantity fn(t) can, for example, be represented on the complex plane as a time- dependent vector that will be rotating about a point F0 on the real axis. F0 represents the average or DC-value of the waveform. F0 will always be real and will lie on the real axis. The complexor fit) is represented graphically in Figure 2.3.1 for k= 3. It is shown to be formed by the vector sum of the instantaneous values of the fundamental and two harmonic components.
Figure 2.3.1 - Graphical instantaneous representation of the voltage u(t) 17
Although each of the fn(t) components have constant magnitude and a constant angular velocity of rotation about the origin of nco, the complexor fn(t) has neither constant magnitude nor does it rotate about the origin with a constant angular velocity. The average angular velocity of _At) about its point of rotation is, however, still equal to co because it makes one revolution in one fundamental period T. When the real part off(t) is plotted against time, it will yield the periodic but distorted waveform that its Fourier series represents.
2.3.2 Harmonic phasors
The nth harmonic phasor Fn can be calculated directly from the complex Fourier series 23 and adapted for unilateral application as:
r Fn = —TEO) ejnwtdt (2.3.4)
It follows from eq. (2.3.4) that:
Fn = Fn ei(1)11 (2.3.5)
Equation (2.3.5) can conveniently be written in the conventional polar form as:
Fn = Fn EL° (2.3.6)
2.3.3 Space-vector manipulations
Because the harmonic quantities conform to sinusoidal requirements at each harmonic order, it is simple to show that the following operations are equivalent:
The magnitude of the complex quantity fn(t) can be calculated as:
ln(t)fn(0 * =1In(t)12 (2.3.7)
Iffn(t) is substituted in eq. (2.3.7) from eq. (2.3.2), the exponents cancel, reducing it to:
FnFn* = iFni2 = Fn2 (2.3.8)
in which fn(t) and Fn represent the rms-scaled magnitude of the two quantities, therefore:
4.1,2(tYn(0* = Fn (2.3.9) 18 2.3.4 Effective value
The effective value of a waveform is defined to be the root mean square value of the measurable time-dependent value over a representative interval of time. For the periodic but non-sinusoidal function fit) of period T, the scalar effective value can be obtained from:
= 11 —T1 SNERe{ f (t)}2 dt (2.3.10)
By substituting forfit) from equations (2.3.3) and (2.3.3) the following equation is obtained:
~/ = iJ Re{f„ dt (2.3.11)
Because fn(t) has constant magnitude and makes a full revolution on the complex plane in one period, it can be shown that eq. (2.3.11) is also equivalent to the following:
T fn (t) fn (t) * dt (2.3.12) n=0 T
The order of summation and integration may be exchanged in eq. (2.3.12) to give: