Harmonic allocation to MV customers in rural distribution systems
|
|
- Iris Baldwin
- 5 years ago
- Views:
Transcription
1 University of Wollongong Researc Online Faculty of Engineering - Papers (Arcive) Faculty of Engineering and Information Sciences 2007 Harmonic allocation to MV customers in rural distribution systems Victor Gosbell University of Wollongong, vgosbell@uow.edu.au ttp://ro.uow.edu.au/engpapers/5380 Publication Details V. J. Gosbell, "Harmonic allocation to MV customers in rural distribution systems," in 2007 Australasian Universities Power Engineering Conference: AUPEC 2007, 2007, pp Researc Online is te open access institutional repository for te University of Wollongong. For furter information contact te UOW Library: researc-pubs@uow.edu.au
2 1 Harmonic Allocation to MV Customers in Rural Distribution Systems V. J. Gosbell, Member, IEEE Abstract-- Te armonic management of power systems in Australia is governed by AS/NZS Tis document is limited to a discussion of principles wit little advice on detailed calculation. A Standards Australia andbook sows ow te standard can be applied to find customer armonic allocations for city and urban power systems. Tis paper furter extends tis work to rural power systems aving te features of isolated lumped loads and radials wit spurs. A general armonic allocation equation is developed containing an allocation coefficient wic is required to be found for eac MV system. It is sown tat te computation to find tis coefficient can be simplified so tat a small spreadseet can be used. Index Terms armonics, standards, distribution systems, IEC. Symbol E Ii G MV I i k L MV L US R SBase S eq S i V MV x av x eq x i α I. NOMENCLATURE Meaning Harmonic current emission allocation to load "i" Harmonic voltage available to local MV loads Harmonic order Harmonic current from load "i" Allocation constant Local MV Planning Level Upstream Planning Level Ratio of fault levels at supply and downstream ends of feeder Base MVA Equivalent lumped max demand for uniformly distributed load Load "i" maximum demand Voltage at end of feeder due to MV loads Average impedance at pcc for uniformly distributed load Equivalent lumped impedance at pcc for uniformly distributed load Harmonic impedance at pcc of load "i" Summation exponent II. INTRODUCTION NTERNATIONAL armonic surveys sow tat THD I voltages are generally increasing at about 1% per decade and tere are many sites in Europe were armonic levels already exceed limits, also known as Planning Levels, for some armonics [1]. It is important tat utilities ave effective policies for ensuring tat te armonic loads being connected V. J. Gosbell is Honorary Professorial Fellow wit te University of Wollongong and Tecnical Adviser to te Integral Energy Power Quality & Reliability Centre ( v.gosbell@uow.edu.au). to teir system do not exceed teir system s capacity to absorb distortion. At LV, equipment armonic emission is controlled by equipment standards suc as [2], wile at MV, armonic current is allocated to eac customer s installation guided by te principles of AS/NZS [3]. Unfortunately, tis document is unclear in places, and is particular difficult for occasional users. Standards Australia commissioned a andbook [4] to give a guide to applying te standard. Tis guide was limited in application to MV distribution systems wit uniformly distributed loads and a simple radial topology, based on a correction factor as developed in [5]. Wit a growt in mining loads remote from towns, tere is now a need to extend te procedure to rural systems aving te following additional features of lumped loads and spur lines. Tis paper proposes an approac wic can be applied to rural systems and is illustrated wit a particular test system wic will be used to confirm te usefulness of te suggested approac. It allows te determination of te armonic current to be allocated to an installation under Stage 2 of te standard. It is assumed tat pfc capacitors are fitted wit series detuning reactors so tat te armonic impedance can be determined from simple fault level considerations. LV installations are ignored in tis paper since teir armonic effects on MV systems are not usually of great significance [4]. Te treatment of diversity given in Section III.A is based on a combination of experience and teory and is necessarily approximate. Hence little is to be gained by an attempt at an exact approac, especially as tis leads to te need for an enormous amount of data. Engineering judgment need to be exercised in most armonic allocations were te accuracy required is about 10%. III. STAGE 2 OVERVIEW Tis section gives te main aspects of Stage 2 as given in [3, 4]. It identifies te basic principles tat will be used in te paper for te development of analytical tecniques wic can address rural systems. A. Representation of time-varying quantities Under normal conditions, power system armonics are due to te interaction of many time-varying distorting loads. Te following principles are recommended for representing voltages and currents in tis situation [3]:
3 2 1. Time-varying voltages and currents are represented by teir 95% probability values. 2. Te 95% values of voltages and currents can be found from te 95% value of te independent contributions using te Summation Law to account for diversity. Tis law is Vtot α = V α 1 + V α α V n (1) were α varies wit te armonic order. A value of α = 1.4 is used for 5 10 and 2 for > 10. B. Allocation principles Te Standard adopts tese principles for armonic allocation under Stage 2: 1. Te armonic allowance to an installation increases wit maximum demand. 2. Te armonic allowance sould be suc tat te igest armonic voltage in te local power system just reaces te Planning Level wen te system is fully loaded and all installations are taking teir maximum allocation. Reference [3] discusses te coice of quantity to be allocated. In distribution systems wit long feeders, as is usually te case in Australia, armonic VA is sown to be a more suitable quantity tan voltage or current since it gives a useful allocation to customers at points of low fault level. Tis allocation quantity will be adopted ere. C. Interaction wit oter parts of te distribution system 33kV 11kV Fig kv power system wit study installation aving maximum demand S i In te local power system sown in Fig. 1, te armonic voltages are influenced by loads connected elsewere wit armonic effects propagating via te 33kV system. Similarly, te local power system influences armonic voltages upstream. Assumptions need to be made wic allow te user to avoid making any detailed study of upstream effects due to te 33kV system and oter more remote interconnected parts. It is assumed tat tat te wole power system is armonically fully loaded so tat te 33kV bus reaces its Planning Level. Te 33kV bus is represented as a constant voltage source equal to te Planning Level for eac armonic being studied. Tis is an approximation because some part of te 33kV bus armonic voltage is due to downstream effects of te local power system and it is incorrect to use a voltage source representation for tis contribution. However te assumption is pessimistic and gives some safety margin for te allocation calculation. If L MV is te MV Planning Level and L US te Planning Level at te next voltage level S i upstream, te armonic voltage available to local MV loads is, making use of (1), α α α G = L L (2) MV MV IV. HARMONIC ALLOCATION IN RURAL DISTRIBUTION SYSTEMS A. Overview Te system of Fig. 2 will be used to discuss te details of te allocation process. It as been cosen because it is simple but as all te features wic need to be treated to apply to rural system armonic allocations: Feeders 1 and 2 wit distributed loads. Feeders 2 and 3 wit lumped loads. Feeder 3 wit spurs 3A and 3B. V US S t, x t V MV impedance x s S D, x D S A (x sa to x da ) S B (x sb to x db ) Fig. 2. Example system US Feeder 1 Feeder 2 Feeder 3A S E, x E S C, x C Feeder 3B S F, x F One as to consider future as well as present loads in te allocation process (Section III.B, point 2), ence not all te loads sown in Fig. 2 need be connected at present. In some cases, te position of future loads is uncertain and engineering judgment needs to be made. Te allocation principles lead to an allocation law given in Section IV.B. Te armonic current drawn by MV loads is determined as sown in Sections IV.C, IV.D. Harmonic voltage can ten be calculated as sown in Section IV.F, owever tis calculation need only be done for te weakest feeder as defined in Section IV.E. Te value of k is ten determined by simple proportion as given in Section IV.G. B. Te allocation equation A load wit maximum demand S i is connected to te system at a point were te armonic impedance is x i. It is sown in [6] tat te principles of armonic allocation are followed if te load S i is given te following current allocation: 1/α k Si E Ii = (3) x i k is te allocation constant for all MV loads in te local power system and varies wit te armonic order.
4 3 C. Lumped MV loads If resistance is small at armonic frequencies, te armonic impedance at te point of connection of load S i is given closely by x i = x i1 (4) were x i1 is te fundamental reactance determined from fault level considerations. Te armonic current for installation i is taken as its full allocated emission from (3) 1/α k Si I i = (5) x D. Uniformly distributed MV loads Te metodology requires tat all loads are represented as lumped quantities. We ave sown in Appendix A tat uniformly distributed MV loads drawing equal armonic VA can be replaced wit adequate accuracy, for armonic calculations, by a lumped loads aving a somewat increased value of total MVA and situated at te appropriate pcc using (A.9, A.10). Te loads S A, S B are replaced by: Feeder 1: S Aeq is concentrated at x Aeq. Feeder 2: S Beq is concentrated at x Beq. Te system of Fig. 2 becomes transformed to tat of Fig. 3, wic also sows node numbers for future reference. x s 1 2 S D, x D i S Aeq, x Aeq S Beq, x Beq Feeder Feeder 3A S E, x E Fig. 3. Lumped model equivalent of Fig. 2 Feeder 2 S C, x C S F, x F Feeder 3B 7 E. Determination of weakest feeder It is not necessary to calculate te igest voltage at te end of eac feeder. Harmonic loading is related to fundamental loading since armonic currents relate to maximum demand and armonic impedance is related to fundamental impedances as given in (4). It is usually only necessary to determine te value of armonic voltage at te end of te feeder wit te igest fundamental loading. Te fundamental loading can be estimated for eac feeder by computing te sum Σx i S i and finding te one aving te largest value see Table I for te calculations for te system in Fig. 3. Uniform loads are given an impedance value corresponding to te average connection point. Tis is not exactly te same as te impedance of te equivalent lumped load used for armonic calculations (A.10). TABLE I FUNDAMENTAL LOADING FOR EACH FEEDER IN FIG. 2 Feeder 1: x Aav S A Feeder 2: x Bav S B + x C S C Feeder 3A: x D S D + x D S F + x E S E Feeder 3B: x D S D + x D S E + x F S F Notice tat for Spur 3A, it is necessary to allow for te impact of te load on Spur 3B at te common impedance x D. A similar remark applies to te impact on Spur 3B of te load on Spur 3A. Were tere is uniform construction trougout te local system (eiter overead open-wire, ABC or underground), reactances can be replaced by line lengts. F. Determination of armonic voltages We determine te armonic voltage appearing at te end of te weakest feeder as identified in Section IV.E. Te following steps are necessary: (a). Harmonic currents are determined for MV loads (5) wit k initially set to unity. (b). Harmonic impedances are cosen for eac load at te pcc wit te weakest feeder. For te Fig. 2 example, if Feeder 3A is te weakest, impedances are cosen as Load S Aeq : te supply impedance x s Load S F : te impedance to D, tat is x D Load S E : te impedance to E, tat is x E G. Determining te allocation constant Te armonic voltage contribution of MV loads to te armonic voltage at te end of te weakest feeder, wit k equal to unity, can be found by summing te separate contributions using te Summation Law (Section III.A). For oter values of k, linearity applies and te MV contribution will be V MV.weakest = k V MV.weakest ' (6) were te das ' indicates a value determined for k = 1. Hence, G k = MV V ' (7) MV.weakest H. Extension to oter frequencies We note tat tat k in (7) is te ratio of two quantities. Te numerator, G MV will need to be determined for eac frequency. Te denominator quantity only needs to be determined in detail at two armonic frequencies as discussed below. Te weakest feeder as been defined in terms of fundamental quantities (Section IV.E) and will not cange wit frequency. Te armonic voltage is determined using te Summation Law as te sum of terms suc as (I x ) α wit k initially at unity. From te allocation law (3), tese terms are of te form Sx α/2 = S α/2 x α/2 1.
5 4 Wen te resultant V is ten found by taking te α root, we find V ' (8) Hence G k MV (9) For a given value of α, if k 5, as been determined, oter values (for 10) can be found by proportion G 5 k = MV k 5 G (10) At = 11, α is cosen as 2, and k 11 cannot be found from (10) it needs to be determined from a detailed calculation, as for k 5. Once determined, oter values (for 11) can be found from G 11 k = MV k11 G (11) 5MV 11MV V. EXAMPLE CALCULATION A. System Data We now consider te case of Fig. 2 were te data in Table II applies. Te aim is to determine te allocation constant k for te 5 t armonic. ceck te proposed metodology. Te system as been represented wit eac loads replaced by te armonic current in Appendix B.D scaled up by k from Appendix B.E. Harmonic voltages were found as given in Table III. Te largest voltage is found at te end of Feeder 3B (Node 7) and just equals G MV as determined in Appendix B.C sowing tat te approac is sound. TABLE III V H AND FUNDAMENTAL LOADING FOR EACH NODE. NODES AT THE END OF FEEDERS ARE SHOWN IN BOLD. Node ΣS lengt V V as been plotted against fundamental loading in Fig. 4 to ceck te metod for finding te weakest feeder. It can be seen tat tere is close relationsip between tese two parameters so tat loading is a good predictor of te weakest feeder. V Feeder 1 Feeder 2 Feeder 3A Feeder 3B TABLE II DATA FOR FIRST EXAMPLE V US V MV S t x t α L US5 L MV5 x/km 33kV 11kV 25MV A 15% Ω/ km Load A B C D E F Type Uniform Uniform Lumped Lumped Lumped Lumpe d S- MVA Posnkm B. Outline of calculations Te full calculation is given in Appendix B. A base of 1MVA is cosen for per unit calculations. System armonic impedances are found in Appendix B.A, including equivalent impedances for lumped equivalents to S A and S B. Te igest armonic voltage is seen to occur at te end of Feeder 3B. In Appendix B.C, for te Planning Levels given in Table II, G MV, te armonic voltage available for local MV loads is found as 3.1%. Wit k = 1, MV loads would cause a voltage rise of 158% (Appendix B.D). Te ratio of G MV to tis voltage gives a value for k of As an example of its use, for load S C, E IC is determined from (3) as /1.4 / or pu. Relative to its fundamental current of 2.5 pu, tis is an allocated 5 t armonic current of 3.1%. C. MATLAB verification An independent MATLAB calculation as been made to D. Sensitivity Loading Fig. 4. Plot of V vs ΣS lengt Te sensitivity of k on system data as been determined from te MATLAB program. For system parameters, te sensitivity is 0.3%-per-% or less. Te igest sensitivity is for load caracteristics at D and F on te weakest feeder. Tis compares wit about 2%-per-% sensitivity to parameters defined in te standards suc as α and L MV. Tis demonstrates tat te proposed metod is robust. VI. CONCLUSIONS An approac as been developed to implement te basic principles of armonic standards for rural systems. Tis leads to te concept of a armonic allocation constant for a particular local power system. Calculations can be simplified by (i) determining te voltage only at te end of one feeder wic can be easily identified, (ii) making full calculations only at two armonic frequencies. Uniformly distributed MV loads can be lumped into a single equivalent load. Te calculations are sufficiently simple tat tey can be implemented by spreadseet.
6 5 VII. APPENDIX A - UNIFORMLY DISTRIBUTED MV LOADS Tis teory used in tis section is based on results from [5]. Suppose we ave a feeder wit te following features: Total distributed load S Load distributed between points aving fundamental reactance x s at supply end and x d at downstream end Harmonic order wit Summation Law coefficient α Load allocated armonic current following te principle of equal armonic VA to equal loads A specific installation S i is allocated a current I i based on te armonic VA principle and following (3). Te armonic current drawn by tis load can be closely approximated by [5]. I ~ k S (1/α) R -0.3 / (x s ) (A.1) Te armonic voltage caused by tis load is igest at te far end and can be approximated by V ~ k ( x s )S (1/α) R 0.33 (A.2) We now find an equivalent lumped load S eq and position x eq (given in terms of fundamental sort-circuit reactance) to give te same results. We tus need to solve te following simultaneous equations k S (1/α) eq / (x eq ) = k S (1/α) R -0.3 / (x s ) k ( x eq )S (1/α) eq = k ( x s )S (1/α) R 0.33 (A.3) (A.4) Multiplying (A.3) by (A.4) S 2/α eq = S 2/α R 0.03 giving S eq = SR 0.03α/2 (A.5) Dividing (A.4) by (A.3) 1/x eq = 1/(R 0.63 x s ) Giving x eq = R 0.63 x s (A.6) Tere are approximations in tis derivation and it was cecked and furter improved by means of a spreadseet study. A model was set up were te uniformly distributed load was assumed to be accurately represented by a 10-stage discrete model. Values of armonic current and voltage were determined for combinations of R = 6, 10, 28 and α = 1, 1.4, 2. Based on te above teory, te lumped model was assumed to be of te form S eq = SR x1.α (A.7) x eq = x s R x2 (A.8) Values were found for te parameters x1, x2 to minimise te maximum of te percentage errors between te accurate and te lumped equivalent model. Te following modified equations gives a worst-case error of less tan 11% and an average error of 5% over te studied range of R and α. S eq = SR 0.044α (A.9) x eq = x s R 0.64 (A.10) Te original equations (A.7, A.8) gave a worst-case error of 21% and an average error of 7%. VIII. APPENDIX B - NUMERICAL EXAMPLE Tese calculations ave been set out so tat tey can easily adapted to a spreadseet. Coose a base S Base = 1 MVA. Tis coice leads to eac value of S in per unit being identical to te MVA value as given in Table II. A. System impedances Find te relevant fundamental and armonic reactances for = 5. x s = x t = 0.15 (New base/old base) = /25 = pu x s = x s = pu At 11kV, Z Base = V 2 /S Base = 11 2 /1 = 121Ω Determine armonic reactances corresponding to eac lumped load connection point. For uniform loads, determine te quantity x d, te armonic reactance to te downstream end of te distributed load. At A: x Ad = x s + len A x/z Base = /121 = pu Similarly, x Bd = pu, x C = pu, x D = pu, x E = pu, x F = pu Determine equivalent MVA and point of connection for te distributed loads A and B. S Aeq = S A R A 0.044α = = 2.79 pu S Beq = S B R B 0.044α = = 2.70 pu x Aeq = x s R A 0.64 = = pu x Beq = x s R B 0.64 = = pu B. Weakest feeder Since in tis case line lengts are proportional to impedances, we use tem in te loading calculation for simplicity. 1: S A len Aav = /2 = : S B len Bav + S C len C = 2.5 5/ = A: S D len D + S E len E + S F len D = = 60 3B: S D len D + S E len D + S F len F = = 75 Feeder 3B is te weakest as it as te largest S lengt product. C. Harmonic voltage available to MV loads Determine te voltage available for MV loads. α α α 1.4 G MV = L L = = pu MV US 031 D. Harmonic voltage at end of weakest feeder Determine te voltage at te end of te weakest feeder 3B wit k provisionally put to "1", using te symbol ' as a reminder. We need to determine te armonic currents injected by eac load from (5). 1/α 1/1.4 SA 2.79 Feeder A: IA ' = = = pu x Aeq
7 6 Similarly, I B ' = 7.92 pu, I C ' = 3.87 pu, I D ' = 3.20 pu, I E ' = 2.36 pu, I F ' = 1.96 pu. We also need to determine te appropriate impedance for eac load to give te impact of its armonic current on te armonic voltage at te end of weakest feeder 3B. Tese are x s, x s, x s, x D, x D and x F for loads A-F respectively. Hence te armonic voltage at te end of Feeder 3B is made up of te following components wic will need to be combined by te Summation Law: A: V A ' = I A 'x s = = 0.21 pu Similarly, V B ' = 0.24 pu, V C ' = 0.12 pu, V D ' = 0.56 pu, V E ' = 0.41 pu, V F ' = 0.91 pu. Combining tese voltages V ' = α α α α α α α V A + VB + VC + VD + VE + VF = = 1.58 E. Allocation constant k is determined as te ratio of G MV from B.3 and V ' from B.4 k = /1.58 = As a quick ceck, we can compute te relative armonic current for a couple of te loads, for example B (te closest to te supply) and F (te furtest). For B: I B = k I B ' = = 0.16 pu Relative to I B1, I B = 0.16/S B = 0.16/2.5 = 6.2% For F: I F = k I F ' = = pu Relative to I F1, I F = 0.039/S F = 0.039/1.5 = 2.6%. Tese values are in te range of typical MV 5 t armonic current allocations and decrease for loads furter away from te supply as would be expected from a armonic VA allocation policy. IX. REFERENCES [1] Euroelectric Report, "PQ in European Electricity Supply Networks", February 2002 [2] AS/NZS :1998, "Electromagnetic compatibility (EMC) Part 3.2: Limits-Limits for armonic current emissions (equipment input current less tan or equal to 16A per pase)" [3] AS/NZS , "Limits Assessment of emission limits for distorting loads in MV and HV power systems", Standards Australia [4] V.J. Gosbell, S. Perera, V. Smit, D. Robinson and G. Sanders, "Power Quality Recommendations for te application of AS/NZS and AS/NZS ", Standards Au`stralia, HB , August 2003, ISBN [5] V.J. Gosbell and D Robinson, "Allocating armonic emission to MV customers in long feeder systems", AUPEC03, Sept-Oct, 2003, Cristcurc [6] V.J. Gosbell, D.A. Robinson and B.S.P. Perera and A. Baitc, "Te application of IEC to MV systems in Australia", ERA Conference, Tame, Feb 2001, pp X. BIOGRAPHIES Vic Gosbell (M'1975) obtained is BSc, BE and PD degrees from te University of Sydney. He as eld academic positions at te University of Sydney and te University of Wollongong were e became te foundation Professor of Power Engineering. He is now an Honorary Professorial Fellow and Tecnical Advisor to te Integral Energy Power Quality and Reliability Centre. He is currently working on armonic management, power quality monitoring and standards. He is a member of Australian standards and CIGRE sub-committees and is a Fellow of te Institution of Engineers, Australia.
HARMONIC ALLOCATION TO MV CUSTOMERS IN RURAL DISTRIBUTION SYSTEMS
HARMONIC ALLOCATION TO MV CUSTOMERS IN RURAL DISTRIBUTION SYSTEMS V Gosbell University of Wollongong Department of Electrical, Computer & Telecommunications Engineering, Wollongong, NSW 2522, Australia
More informationA new approach to harmonic allocation for medium-voltage installations
University of Wollongong Researc Online Faculty of ngineering and Information Sciences - Papers: Part A Faculty of ngineering and Information Sciences 03 A new approac to armonic allocation for medium-voltage
More informationThe application of IEC to MV systems in Australia
University of Wollongong Research Online Faculty of Informatics - Papers (Archive) Faculty of Engineering and Information Sciences 00 The application of IEC 6000-3-6 to MV systems in Australia Victor J.
More informationAn improved methodology for determining MV to LV voltage unbalance transfer coefficient
University of Wollongong Research Online Faculty of Engineering - Papers (Archive) Faculty of Engineering and Information Sciences 2008 An improved methodology for determining MV to LV voltage unbalance
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More informationAN IMPROVED WEIGHTED TOTAL HARMONIC DISTORTION INDEX FOR INDUCTION MOTOR DRIVES
AN IMPROVED WEIGHTED TOTA HARMONIC DISTORTION INDEX FOR INDUCTION MOTOR DRIVES Tomas A. IPO University of Wisconsin, 45 Engineering Drive, Madison WI, USA P: -(608)-6-087, Fax: -(608)-6-5559, lipo@engr.wisc.edu
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationSimulation and verification of a plate heat exchanger with a built-in tap water accumulator
Simulation and verification of a plate eat excanger wit a built-in tap water accumulator Anders Eriksson Abstract In order to test and verify a compact brazed eat excanger (CBE wit a built-in accumulation
More informationThe Zeckendorf representation and the Golden Sequence
University of Wollongong Researc Online Faculty of Engineering and Information Sciences - Papers: Part A Faculty of Engineering and Information Sciences 1991 Te Zeckendorf representation and te Golden
More information3.1 Extreme Values of a Function
.1 Etreme Values of a Function Section.1 Notes Page 1 One application of te derivative is finding minimum and maimum values off a grap. In precalculus we were only able to do tis wit quadratics by find
More informationSolution. Solution. f (x) = (cos x)2 cos(2x) 2 sin(2x) 2 cos x ( sin x) (cos x) 4. f (π/4) = ( 2/2) ( 2/2) ( 2/2) ( 2/2) 4.
December 09, 20 Calculus PracticeTest s Name: (4 points) Find te absolute extrema of f(x) = x 3 0 on te interval [0, 4] Te derivative of f(x) is f (x) = 3x 2, wic is zero only at x = 0 Tus we only need
More informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
More informationMaterial for Difference Quotient
Material for Difference Quotient Prepared by Stepanie Quintal, graduate student and Marvin Stick, professor Dept. of Matematical Sciences, UMass Lowell Summer 05 Preface Te following difference quotient
More informationContinuity and Differentiability of the Trigonometric Functions
[Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be te
More informationREVIEW LAB ANSWER KEY
REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g
More informationDerivation Of The Schwarzschild Radius Without General Relativity
Derivation Of Te Scwarzscild Radius Witout General Relativity In tis paper I present an alternative metod of deriving te Scwarzscild radius of a black ole. Te metod uses tree of te Planck units formulas:
More informationHOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS
HOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS Po-Ceng Cang National Standard Time & Frequency Lab., TL, Taiwan 1, Lane 551, Min-Tsu Road, Sec. 5, Yang-Mei, Taoyuan, Taiwan 36 Tel: 886 3
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationA MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES
A MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES Ronald Ainswort Hart Scientific, American Fork UT, USA ABSTRACT Reports of calibration typically provide total combined uncertainties
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationModeling of nonlinear loads in high-voltage network by measured parameters
nternational Conference on Renewable Energies and Power Quality (CREPQ 17) Malaga (Spain), 4 t to 6 t April, 017 Renewable Energy and Power Quality Journal (RE&PQJ) SSN 17-038 X, No.15 April 017 Modeling
More information. If lim. x 2 x 1. f(x+h) f(x)
Review of Differential Calculus Wen te value of one variable y is uniquely determined by te value of anoter variable x, ten te relationsip between x and y is described by a function f tat assigns a value
More informationStepped-Impedance Low-Pass Filters
4/23/27 Stepped Impedance Low Pass Filters 1/14 Stepped-Impedance Low-Pass Filters Say we know te impedance matrix of a symmetric two-port device: 11 21 = 21 11 Regardless of te construction of tis two
More informationMicrostrip Antennas- Rectangular Patch
April 4, 7 rect_patc_tl.doc Page of 6 Microstrip Antennas- Rectangular Patc (Capter 4 in Antenna Teory, Analysis and Design (nd Edition) by Balanis) Sown in Figures 4. - 4.3 Easy to analyze using transmission
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationChapter Seven The Quantum Mechanical Simple Harmonic Oscillator
Capter Seven Te Quantum Mecanical Simple Harmonic Oscillator Introduction Te potential energy function for a classical, simple armonic oscillator is given by ZÐBÑ œ 5B were 5 is te spring constant. Suc
More informationThe derivative function
Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationGraviton Induced Nuclear Fission through Electromagnetic Wave Flux Phil Russell, * Jerry Montgomery
Graviton Induced Nuclear Fission troug Electromagnetic Wave Flux Pil Russell, * Jerry Montgomery Nort Carolina Central University, Duram, NC 27707 Willowstick Tecnologies LLC, Draper, UT 84020 (Dated:
More informationLecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.
Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative
More informationExercises for numerical differentiation. Øyvind Ryan
Exercises for numerical differentiation Øyvind Ryan February 25, 2013 1. Mark eac of te following statements as true or false. a. Wen we use te approximation f (a) (f (a +) f (a))/ on a computer, we can
More informationDe-Coupler Design for an Interacting Tanks System
IOSR Journal of Electrical and Electronics Engineering (IOSR-JEEE) e-issn: 2278-1676,p-ISSN: 2320-3331, Volume 7, Issue 3 (Sep. - Oct. 2013), PP 77-81 De-Coupler Design for an Interacting Tanks System
More informationDefinition of the Derivative
Te Limit Definition of te Derivative Tis Handout will: Define te limit grapically and algebraically Discuss, in detail, specific features of te definition of te derivative Provide a general strategy of
More informationAN ANALYSIS OF AMPLITUDE AND PERIOD OF ALTERNATING ICE LOADS ON CONICAL STRUCTURES
Ice in te Environment: Proceedings of te 1t IAHR International Symposium on Ice Dunedin, New Zealand, nd t December International Association of Hydraulic Engineering and Researc AN ANALYSIS OF AMPLITUDE
More informationSECTION 1.10: DIFFERENCE QUOTIENTS LEARNING OBJECTIVES
(Section.0: Difference Quotients).0. SECTION.0: DIFFERENCE QUOTIENTS LEARNING OBJECTIVES Define average rate of cange (and average velocity) algebraically and grapically. Be able to identify, construct,
More informationMAT244 - Ordinary Di erential Equations - Summer 2016 Assignment 2 Due: July 20, 2016
MAT244 - Ordinary Di erential Equations - Summer 206 Assignment 2 Due: July 20, 206 Full Name: Student #: Last First Indicate wic Tutorial Section you attend by filling in te appropriate circle: Tut 0
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationch (for some fixed positive number c) reaching c
GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan
More informationEfficient algorithms for for clone items detection
Efficient algoritms for for clone items detection Raoul Medina, Caroline Noyer, and Olivier Raynaud Raoul Medina, Caroline Noyer and Olivier Raynaud LIMOS - Université Blaise Pascal, Campus universitaire
More information1watt=1W=1kg m 2 /s 3
Appendix A Matematics Appendix A.1 Units To measure a pysical quantity, you need a standard. Eac pysical quantity as certain units. A unit is just a standard we use to compare, e.g. a ruler. In tis laboratory
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationLIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION
LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y
More informationContinuity and Differentiability Worksheet
Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;
More informationPreface. Here are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed.
Preface Here are my online notes for my course tat I teac ere at Lamar University. Despite te fact tat tese are my class notes, tey sould be accessible to anyone wanting to learn or needing a refreser
More informationThese errors are made from replacing an infinite process by finite one.
Introduction :- Tis course examines problems tat can be solved by metods of approximation, tecniques we call numerical metods. We begin by considering some of te matematical and computational topics tat
More informationTHE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS. L. Trautmann, R. Rabenstein
Worksop on Transforms and Filter Banks (WTFB),Brandenburg, Germany, Marc 999 THE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS L. Trautmann, R. Rabenstein Lerstul
More informationExam 1 Review Solutions
Exam Review Solutions Please also review te old quizzes, and be sure tat you understand te omework problems. General notes: () Always give an algebraic reason for your answer (graps are not sufficient),
More informationTe comparison of dierent models M i is based on teir relative probabilities, wic can be expressed, again using Bayes' teorem, in terms of prior probab
To appear in: Advances in Neural Information Processing Systems 9, eds. M. C. Mozer, M. I. Jordan and T. Petsce. MIT Press, 997 Bayesian Model Comparison by Monte Carlo Caining David Barber D.Barber@aston.ac.uk
More informationMath 31A Discussion Notes Week 4 October 20 and October 22, 2015
Mat 3A Discussion Notes Week 4 October 20 and October 22, 205 To prepare for te first midterm, we ll spend tis week working eamples resembling te various problems you ve seen so far tis term. In tese notes
More informationNONLINEAR SYSTEMS IDENTIFICATION USING THE VOLTERRA MODEL. Georgeta Budura
NONLINEAR SYSTEMS IDENTIFICATION USING THE VOLTERRA MODEL Georgeta Budura Politenica University of Timisoara, Faculty of Electronics and Telecommunications, Comm. Dep., georgeta.budura@etc.utt.ro Abstract:
More informationCS522 - Partial Di erential Equations
CS5 - Partial Di erential Equations Tibor Jánosi April 5, 5 Numerical Di erentiation In principle, di erentiation is a simple operation. Indeed, given a function speci ed as a closed-form formula, its
More informationEOQ and EPQ-Partial Backordering-Approximations
USING A ONSTANT RATE TO APPROXIMATE A LINEARLY HANGING RATE FOR THE EOQ AND EPQ WITH PARTIAL BAKORDERING David W. Pentico, Palumo-Donaue Scool of Business, Duquesne University, Pittsurg, PA 158-18, pentico@duq.edu,
More informationINTRODUCTION AND MATHEMATICAL CONCEPTS
Capter 1 INTRODUCTION ND MTHEMTICL CONCEPTS PREVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips
More informationQuasiperiodic phenomena in the Van der Pol - Mathieu equation
Quasiperiodic penomena in te Van der Pol - Matieu equation F. Veerman and F. Verulst Department of Matematics, Utrect University P.O. Box 80.010, 3508 TA Utrect Te Neterlands April 8, 009 Abstract Te Van
More informationAalborg Universitet. Published in: Proceedings of The International conference on future Power Systems. Publication date: 2005
Aalborg Universitet Harmonic Domain Modelling of a Distribution System using te DigSilent PowerFactory Software Wasilewski, J.; Wiecowski, Wojciec Tomasz; Bak, Claus Let Publised in: Proceedings of Te
More information5.1 We will begin this section with the definition of a rational expression. We
Basic Properties and Reducing to Lowest Terms 5.1 We will begin tis section wit te definition of a rational epression. We will ten state te two basic properties associated wit rational epressions and go
More informationAverage Rate of Change
Te Derivative Tis can be tougt of as an attempt to draw a parallel (pysically and metaporically) between a line and a curve, applying te concept of slope to someting tat isn't actually straigt. Te slope
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More informationChapter 2 Limits and Continuity
4 Section. Capter Limits and Continuity Section. Rates of Cange and Limits (pp. 6) Quick Review.. f () ( ) () 4 0. f () 4( ) 4. f () sin sin 0 4. f (). 4 4 4 6. c c c 7. 8. c d d c d d c d c 9. 8 ( )(
More informationDerivatives of Exponentials
mat 0 more on derivatives: day 0 Derivatives of Eponentials Recall tat DEFINITION... An eponential function as te form f () =a, were te base is a real number a > 0. Te domain of an eponential function
More information2.11 That s So Derivative
2.11 Tat s So Derivative Introduction to Differential Calculus Just as one defines instantaneous velocity in terms of average velocity, we now define te instantaneous rate of cange of a function at a point
More information4.2 - Richardson Extrapolation
. - Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Definition Let x n n converge to a number x. Suppose tat n n is a sequence
More informationCHAPTER (A) When x = 2, y = 6, so f( 2) = 6. (B) When y = 4, x can equal 6, 2, or 4.
SECTION 3-1 101 CHAPTER 3 Section 3-1 1. No. A correspondence between two sets is a function only if eactly one element of te second set corresponds to eac element of te first set. 3. Te domain of a function
More informationINTRODUCTION AND MATHEMATICAL CONCEPTS
INTODUCTION ND MTHEMTICL CONCEPTS PEVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips of sine,
More informationDifferential Calculus (The basics) Prepared by Mr. C. Hull
Differential Calculus Te basics) A : Limits In tis work on limits, we will deal only wit functions i.e. tose relationsips in wic an input variable ) defines a unique output variable y). Wen we work wit
More information5. (a) Find the slope of the tangent line to the parabola y = x + 2x
MATH 141 090 Homework Solutions Fall 00 Section.6: Pages 148 150 3. Consider te slope of te given curve at eac of te five points sown (see text for figure). List tese five slopes in decreasing order and
More informationIntroduction to Derivatives
Introduction to Derivatives 5-Minute Review: Instantaneous Rates and Tangent Slope Recall te analogy tat we developed earlier First we saw tat te secant slope of te line troug te two points (a, f (a))
More informationDistribution of reynolds shear stress in steady and unsteady flows
University of Wollongong Researc Online Faculty of Engineering and Information Sciences - Papers: Part A Faculty of Engineering and Information Sciences 13 Distribution of reynolds sear stress in steady
More informationPositive Sequence Tracking Phase Locked Loops: A Unified Graphical Explanation
Positive Sequence Tracking Pase Locked Loops: A Unified Grapical Explanation L. Matakas Junior *, **, W. Komatsu *, and F. O. Martinz * * Polytecnic Scool of te University of Sao Paulo, Electrical Energy
More informationDiscriminate Modelling of Peak and Off-Peak Motorway Capacity
International Journal of Integrated Engineering - Special Issue on ICONCEES Vol. 4 No. 3 (2012) p. 53-58 Discriminate Modelling of Peak and Off-Peak Motorway Capacity Hasim Moammed Alassan 1,*, Sundara
More informationMath 34A Practice Final Solutions Fall 2007
Mat 34A Practice Final Solutions Fall 007 Problem Find te derivatives of te following functions:. f(x) = 3x + e 3x. f(x) = x + x 3. f(x) = (x + a) 4. Is te function 3t 4t t 3 increasing or decreasing wen
More informationFinite Difference Methods Assignments
Finite Difference Metods Assignments Anders Söberg and Aay Saxena, Micael Tuné, and Maria Westermarck Revised: Jarmo Rantakokko June 6, 1999 Teknisk databeandling Assignment 1: A one-dimensional eat equation
More informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
More informationSin, Cos and All That
Sin, Cos and All Tat James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 9, 2017 Outline Sin, Cos and all tat! A New Power Rule Derivatives
More informationExercise 19 - OLD EXAM, FDTD
Exercise 19 - OLD EXAM, FDTD A 1D wave propagation may be considered by te coupled differential equations u x + a v t v x + b u t a) 2 points: Derive te decoupled differential equation and give c in terms
More informationA general articulation angle stability model for non-slewing articulated mobile cranes on slopes *
tecnical note 3 general articulation angle stability model for non-slewing articulated mobile cranes on slopes * J Wu, L uzzomi and M Hodkiewicz Scool of Mecanical and Cemical Engineering, University of
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More information= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)
Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''
More informationPre-Calculus Review Preemptive Strike
Pre-Calculus Review Preemptive Strike Attaced are some notes and one assignment wit tree parts. Tese are due on te day tat we start te pre-calculus review. I strongly suggest reading troug te notes torougly
More informationHOMEWORK HELP 2 FOR MATH 151
HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,
More informationDepartment of Statistics & Operations Research, Aligarh Muslim University, Aligarh, India
Open Journal of Optimization, 04, 3, 68-78 Publised Online December 04 in SciRes. ttp://www.scirp.org/ournal/oop ttp://dx.doi.org/0.436/oop.04.34007 Compromise Allocation for Combined Ratio Estimates of
More informationMATH1151 Calculus Test S1 v2a
MATH5 Calculus Test 8 S va January 8, 5 Tese solutions were written and typed up by Brendan Trin Please be etical wit tis resource It is for te use of MatSOC members, so do not repost it on oter forums
More informationDigital Filter Structures
Digital Filter Structures Te convolution sum description of an LTI discrete-time system can, in principle, be used to implement te system For an IIR finite-dimensional system tis approac is not practical
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More informationChapter 3 Thermoelectric Coolers
3- Capter 3 ermoelectric Coolers Contents Capter 3 ermoelectric Coolers... 3- Contents... 3-3. deal Equations... 3-3. Maximum Parameters... 3-7 3.3 Normalized Parameters... 3-8 Example 3. ermoelectric
More informationLab 6 Derivatives and Mutant Bacteria
Lab 6 Derivatives and Mutant Bacteria Date: September 27, 20 Assignment Due Date: October 4, 20 Goal: In tis lab you will furter explore te concept of a derivative using R. You will use your knowledge
More informationESTIMATION OF TIME-DOMAIN FREQUENCY STABILITY BASED ON PHASE NOISE MEASUREMENT
ESTIMATION OF TIME-DOMAIN FREQUENCY STABILITY BASED ON PHASE NOISE MEASUREMENT P. C. Cang, H. M. Peng, and S. Y. Lin National Standard Time & Frequenc Laborator, TL, Taiwan, Lane 55, Min-Tsu Road, Sec.
More informationProblem Solving. Problem Solving Process
Problem Solving One of te primary tasks for engineers is often solving problems. It is wat tey are, or sould be, good at. Solving engineering problems requires more tan just learning new terms, ideas and
More informationQuantization of electrical conductance
1 Introduction Quantization of electrical conductance Te resistance of a wire in te classical Drude model of metal conduction is given by RR = ρρρρ AA, were ρρ, AA and ll are te conductivity of te material,
More informationTime (hours) Morphine sulfate (mg)
Mat Xa Fall 2002 Review Notes Limits and Definition of Derivative Important Information: 1 According to te most recent information from te Registrar, te Xa final exam will be eld from 9:15 am to 12:15
More informationINTRODUCTION TO CALCULUS LIMITS
Calculus can be divided into two ke areas: INTRODUCTION TO CALCULUS Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and
More informationPractice Problem Solutions: Exam 1
Practice Problem Solutions: Exam 1 1. (a) Algebraic Solution: Te largest term in te numerator is 3x 2, wile te largest term in te denominator is 5x 2 3x 2 + 5. Tus lim x 5x 2 2x 3x 2 x 5x 2 = 3 5 Numerical
More informationLines, Conics, Tangents, Limits and the Derivative
Lines, Conics, Tangents, Limits and te Derivative Te Straigt Line An two points on te (,) plane wen joined form a line segment. If te line segment is etended beond te two points ten it is called a straigt
More informationDifferentiation. Area of study Unit 2 Calculus
Differentiation 8VCE VCEco Area of stud Unit Calculus coverage In tis ca 8A 8B 8C 8D 8E 8F capter Introduction to limits Limits of discontinuous, rational and brid functions Differentiation using first
More informationOn my honor as a student, I have neither given nor received unauthorized assistance on this exam.
HW2 (Overview of Transport) (Print name above) On my onor as a student, I ave neiter given nor received unautorized assistance on tis exam. (sign name above) 1 Figure 1: Band-diagram before and after application
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More information