HARMONIC ALLOCATION TO MV CUSTOMERS IN RURAL DISTRIBUTION SYSTEMS
|
|
- Carol Curtis
- 5 years ago
- Views:
Transcription
1 HARMONIC ALLOCATION TO MV CUSTOMERS IN RURAL DISTRIBUTION SYSTEMS V Gosbell University of Wollongong Department of Electrical, Computer & Telecommunications Engineering, Wollongong, NSW 2522, Australia 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. 1. INTRODUCTION International armonic surveys sow tat THD 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 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 2.1 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 judgement need to be exercised in most armonic allocations were te accuracy required is about 10%. A list symbols is given after te References section. 2. 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. 2.1 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]: 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 V tot = V 1 + V V n were varies wit te armonic order. A value of = 1.4 is used for 5 10 and 2 for > 10. (1) 122
2 2.2 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. [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. 2.3 Interaction wit oter parts of te distribution system 33kV 11kV 3. HARMONIC ALLOCATION IN RURAL DISTRIBUTION SYSTEMS 3.1 Overview Te system of Figure 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 ) Feeder 1 Feeder 2 Feeder 3A S E, x E S C, x C Feeder 3B S F, x F Figure 1-11 kv power system wit study installation aving maximum demand S i In te local power system sown in Figure 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 upstream, te armonic voltage available to local MV loads is, making use of (1), G MV MV US = L L (2) S i Figure 2 - Example system One as to consider future as well as present loads in te allocation process (Section 2.2, point 2), ence not all te loads sown in Figure 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 3.2. Te armonic current drawn by MV loads is determined as sown in Sections 3.3, 3.4. Harmonic voltage can ten be calculated as sown in Section 3.6, owever tis calculation need only be done for te weakest feeder as defined in Section 3.5. Te value of k is ten determined by simple proportion as given in Section 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 k is te allocation constant for all MV loads in te local power system and varies wit te armonic order. i 123
3 3.3 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) I i = k 1/ Si 3.4 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 Figure 2 becomes transformed to tat of Figure 3, wic also sows node numbers for future reference. x s 1 2 x i S Aeq, x Aeq Feeder Feeder 2 (5) 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 Figure 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 eac feeder in Figure 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. 3.6 Determination of armonic voltages We determine te armonic voltage appearing at te end of te weakest feeder as identified in Section 3.5. 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 Figure 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 S Beq, x Beq S C, x C 3.7 Determining te allocation constant S D, x D 5 6 Feeder 3A S E, x E S F, x F Feeder 3B 7 Figure 3 - Lumped model equivalent of Figure 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 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 2.1). 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 3.8 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 124
4 to be determined in detail at two armonic frequencies as discussed below. Te weakest feeder as been defined in terms of fundamental quantities (Section 3.5) 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 1 /2. Wen te resultant V is ten found by taking te root, we find Hence V ' (8) 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) 5MV 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) 11MV 4. EXAMPLE CALCULATION 4.1. System Data We now consider te case of Figure 2 were te data in Table II applies. Te aim is to determine te allocation constant k for te 5 t armonic. Table II - Data for first example V US V MV S t x t L US5 L MV5 x/km 33kV 11kV 25MVA 15% Ω/ km Load A B C D E F Type Uniform Uniform Lumped Lumped Lumped Lumped S- MVA Posn km 4.2. 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 Section B.1, 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 Section B.3, 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% (Section B.4). 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% MATLAB verification An independent MATLAB calculation as been made to ceck te proposed metodology. Te system as been represented wit eac loads replaced by te armonic current in Section B.4 scaled up by k from Section B.5. 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 Section B.3 sowing tat te approac is sound. Table III - V and fundamental loading for eac node. Nodes at te end of feeders are sown in bold. Node ΣS lengt V V as been plotted against fundamental loading in Figure 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 Loading 4.4. Sensitivity Figure 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. 125
5 5. 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. REFERENCES [1] Euroelectric Report, February 2002 "PQ in European Electricity Supply Networks" [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 Australia, 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 LIST OF MAIN SYMBOLS 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 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 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 sortcircuit reactance) to give te same results. We tus need to solve te following simultaneous equations k S eq (1/) / (x eq ) = k S (1/) R -0.3 / (x s ) k ( x eq )S eq (1/) = k ( x s )S (1/) R 0.33 Multiplying (A.3) by (A.4) S eq 2/ = S 2/ R 0.03 (A.3) (A.4) 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. x eq = x s R x2 (A.7) (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. 126
6 S eq = SR (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%. 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. B.1 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 = = 2.79 pu S Beq = S B R B = = 2.70 pu x Aeq = x s R A 0.64 = = pu x Beq = x s R B 0.64 = = pu B.2 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. B.3 Harmonic voltage available to MV loads Determine te voltage available for MV loads. 1.4 G MV = L L = = pu MV US 031 B.4 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 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 + V + V + V + V + V A B C D E F + = 1.58 B.5 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. 127
Harmonic allocation to MV customers in rural distribution systems
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
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
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 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 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 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 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 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 informationPrecalculus Test 2 Practice Questions Page 1. Note: You can expect other types of questions on the test than the ones presented here!
Precalculus Test 2 Practice Questions Page Note: You can expect oter types of questions on te test tan te ones presented ere! Questions Example. Find te vertex of te quadratic f(x) = 4x 2 x. Example 2.
More informationIntegral Calculus, dealing with areas and volumes, and approximate areas under and between curves.
Calculus can be divided into two ke areas: Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and minima problems Integral
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 informationSection 2.1 The Definition of the Derivative. We are interested in finding the slope of the tangent line at a specific point.
Popper 6: Review of skills: Find tis difference quotient. f ( x ) f ( x) if f ( x) x Answer coices given in audio on te video. Section.1 Te Definition of te Derivative We are interested in finding te slope
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 informationSection 15.6 Directional Derivatives and the Gradient Vector
Section 15.6 Directional Derivatives and te Gradient Vector Finding rates of cange in different directions Recall tat wen we first started considering derivatives of functions of more tan one variable,
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 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 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 informationNumerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1
Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.
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 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 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 informationDerivatives. By: OpenStaxCollege
By: OpenStaxCollege Te average teen in te United States opens a refrigerator door an estimated 25 times per day. Supposedly, tis average is up from 10 years ago wen te average teenager opened a refrigerator
More informationFlavius Guiaş. X(t + h) = X(t) + F (X(s)) ds.
Numerical solvers for large systems of ordinary differential equations based on te stocastic direct simulation metod improved by te and Runge Kutta principles Flavius Guiaş Abstract We present a numerical
More informationDiffraction. S.M.Lea. Fall 1998
Diffraction.M.Lea Fall 1998 Diffraction occurs wen EM waves approac an aperture (or an obstacle) wit dimension d > λ. We sall refer to te region containing te source of te waves as region I and te region
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 informationThe Krewe of Caesar Problem. David Gurney. Southeastern Louisiana University. SLU 10541, 500 Western Avenue. Hammond, LA
Te Krewe of Caesar Problem David Gurney Souteastern Louisiana University SLU 10541, 500 Western Avenue Hammond, LA 7040 June 19, 00 Krewe of Caesar 1 ABSTRACT Tis paper provides an alternative to te usual
More informationPolynomials 3: Powers of x 0 + h
near small binomial Capter 17 Polynomials 3: Powers of + Wile it is easy to compute wit powers of a counting-numerator, it is a lot more difficult to compute wit powers of a decimal-numerator. EXAMPLE
More informationMath 1241 Calculus Test 1
February 4, 2004 Name Te first nine problems count 6 points eac and te final seven count as marked. Tere are 120 points available on tis test. Multiple coice section. Circle te correct coice(s). You do
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 information