DAE-Index and Convergence Analysis of Lumped Electric Circuits Refined by 3-D Magnetoquasistatic Conductor Models
|
|
- Randall Barton
- 5 years ago
- Views:
Transcription
1 Bergische Universität Wuppertal Fachbereich Mathematik un Naturwissenschaften Lehrstuhl für Angewante Mathematik un Numerische Mathematik Preprint BUW-AMNA 8/6 Sebastian Schöps DAE-Inex an Convergence Analysis of Lumpe Electric Circuits Refine by 3-D Magnetoquasistatic Conuctor Moels October 28
2 DAE-Inex an Convergence Analysis of Lumpe Electric Circuits Refine by 3-D Magnetoquasistatic Conuctor Moels Sebastian Schöps 1, Anreas Bartel 1, Herbert De Gersem 2, an Michael Günther 1 Abstract In this paper the fiel/circuit coupling is reconsiere for (non-linear) lumpe electric circuits refine by 3-D magnetoquasistatic conuctor moels, where the circuit is escribe by moifie noal analysis an the fiel is iscretize in the terms of the finite integration technique. The coupling of the two systems of ifferential-algebraic equations is iscusse an two numerical approaches are propose, the weak (co-simulation) an strong coupling (monolithic). The DAE-inex of the subproblems an of the full problem are analyze an the convergence properties of the co-simulation are stuie by inspecting the epenencies of algebraic constraints on solutions of previous iterations. Finally computational results of a simple half rectifier circuit are exemplarily given to prove the concepts. 1 Introuction Many basic elements in circuit analysis are escribe by (non-)linear relations, isregaring fiel effects. Sometimes, they are replace by more complex but more realistic companion moels to meet reality. These give, however, only a partial insight into fiel effects. In contrast refine moels irectly rely upon Maxwell s equations. We restrict our analysis to the refinement of two conuctor types, which are embee into electric circuits an exhibit proximity an skin effects with ey currents. The couple problem is a system of ifferential-algebraic equations (DAEs) originating from Kirchhof s Laws an the iscrete Maxwell Equations. It can be irectly aresse by solving one monolithic system using a fiel- or circuit-oriente approach. In the fiel approach, commonly the circuit is escribe using loop/branch techniques an is solve within the fiel simulator. This approach is quite successful an well unerstoo [1], but it is neither efficient for coupling with very large circuits nor usable within moern circuit simulators that are base on moifie noal 1 Bergische Universität Wuppertal, {schoeps,bartel,guenther}@math.uni-wuppertal.e 2 Katholieke Universitat Leuven, Herbert.DeGersem@kuleuven-kortrijk.be 1
3 2 Sebastian Schöps, Anreas Bartel, Herbert De Gersem, an Michael Günther analysis (MNA). The circuit-oriente approach relies on MNA. Although intensive research has been carrie out [2], companion moels are still wiesprea. Obviously the strongly couple approaches above will not have the avantages of both simulators. Therefore co-simulation becomes beneficial [3]. It allows the use of ifferent simulators for each subproblem an it provies a natural support for iversifying integration methos an time-stepping rates (multirate) towars the subproblems. The coupling is given mathematically by a ynamic iteration scheme. The paper is organize as follows: In Sections 2 an 3 we recall the circuit an fiel settings, in Section 4 we analyze the DAE-inex of the fiel-system, in Section 5 we introuce the monolithic coupling an its co-simulation with inex an convergence analysis, in Sections 6 an 7 we give an example an final conclusions. 2 Lumpe Electric Circuit Electric circuits are given by elements connecte by Kirchhoff s Laws an their escription is given here in terms of MNA, which yiels DAEs since it is base on reunant coorinates. In the charge-flux oriente formulation, [4], the system reas A C t q+a R r(a T R e,t)+a Li L + A V i V + A I i(t)+a λ i λ (A T λ e,t) =, t Φ AT L e =, AT Ve v(t) =, q q C (A T Ce,t) =, Φ Φ L (i L,t) =, (1) with incience matrices A, noe potentials e, inepenent an controlle current an voltage sources i, i λ an v, currents through voltage an flux controlle elements i V an i L, charges an fluxes q an Φ, functions of charges, fluxes an resistances q C, Φ L an r (with positive efinite erivatives), respectively. The numerical properties of (1) are well known, e.g. the DAE-inex has been iscusse by ecomposing the unknown (e,i V,i L,q,Φ) into algebraic an ifferential parts using a projector Q C onto the kernel of A T C, i.e., Q C kera T C = kera T C an A T CQ C =, an its complementary projection P C = I Q C. We assume in the terms above C1 No loops of capacitors an voltage sources, i.e., ker(a C,A R,A V ) T = {}. C2 C3 No cutsets of inuctors an/or current sources, i.e., kerqc T A V = {}. Voltage controlle current sources parallel to capacitors, i.e., QC T A λ i λ =. This splits the unknown into y := (P C e, j L ) T an z := (Q C e, j V,q,φ) T, such that t y = f 1(y,z,i λ ), = g 1 (y,z), (2) is an inex-1 escription of (1) since the erivative z g 1 can be shown to be nonsingular assuming C1-C3. This motivates the following important property [5]:
4 Analysis of Lumpe Electric Circuits Refine by 3-D MQS Conuctor Moels 3 Theorem 1. Let us consier a lumpe electric circuit in form (1) that respects C3, then the flux/charge oriente MNA leas to an inex-1 DAE iff C1-C2 hol, it leas otherwise to an inex-2 DAE. 3 Electromagnetic Fiel The electromagnetic fiel is escribe by Maxwell s equations. We assume that they are spatially iscretize by either the finite integration technique or the finite elements metho using lowest-orer Whitney elements an staggere gris [6,7]. This gives their iscrete counterpart, also calle the Maxwell s Gri Equations (MGE) C e = t b, C h = t + j, S = q, Sb =, (3) with iscrete curl operators C an C, ivergence operators S an S, electric an magnetic fiel strength e an h, current ensity, iscrete magnetic an electric flux j, b an, respectively. For linear materials, MGE are accomplishe by b = Mµ h, = Mε e, j = M σ e, (4) with symmetric positive (semi-)efinite matrices M µ, M ε an M σ for the permeabilities, permittivities an conuctivities, respectively; only M σ is generally singular ue to non-conucting regions [8]. We obtain from (3,4) the curl-curl equation M ε 2 t 2 a + M σ t a + K ν a = j src, (5) with a enoting the magnetic vector potential (MVP) an K ν := CM ν C. 4 Fiel Moels as Refine Network Elements i sol i str Conuctor moels for connecting fiel an circuit parts are well-known. Most common are soli an strane conuctors (Fig. 1). We use the given symbol for a (mulv sol v str (a) Soli (b) Strane (c) Symbol Fig. 1: Conuctor moels (a), (b) an evice symbol (c) that embes both into the circuit.
5 4 Sebastian Schöps, Anreas Bartel, Herbert De Gersem, an Michael Günther tiport) evice that consists of (multiple) conuctors of both types which are tightly couple by the fiel. The latter is escribe by the curl-curl equation an excite by the source current ensity j src ue to the connecte circuit [9]. Typically voltages of solis (v sol ) an currents of strane conuctors (i str ) are consiere to be given an thus the excitation reas j src = M σ Q sol v sol + Q str i str. (6) Here Q = [Q sol,q str ] enotes the coupling matrix. Each column correspons to a conuctor moel an imposes currents/voltages onto eges of the gri. The unknown currents i sol an voltages v str are obtaine by the aitional equations i sol = G sol v sol Q T sol M σ t a, v str = R str i str + Q T str t a, (7) with the iagonal matrices of soli conuctances G sol an strane resistances R str = G 1 str. Let the moel above fulfill the following assumptions F1 The fiel is magnetoquasistatic, i.e., max t max j src. F2 The solution is unique, i.e., [M σ,k ν ] := et(λm σ + K ν ) for a λ. F3 The moels are spatially isjunct, i.e., Q T (i) Q ( j) =, for all i j. F4 The excitation is consistent, i.e., ker(cq sol ) = {}, ker(cm σ,aniso + Q str) = {}. The first assumption F1 implies the neglect of M ε 2 t 2 a in the curl-curl equation (5), such that (5) becomes a first orer DAE, which is generally not uniquely solvable ue to the non-trivial nullspaces of M σ an C. Thus the integration requires the selection of one solution within the equivalent class efine by b = C a, therefore we nee the gauging assumption F2, [1]. From F3 we obtain the equivalent form M σ,fillin t a + K ν a = M σ Q sol v sol + Q str G str v str := j src, Q T sol K ν a = i sol, G str Q T strm + σ,aniso K ν a = i str, (8a) (8b) (8c) where M σ,aniso + is the pseuoinverse of the anisotropic conuctivity matrix for strane conuctors an M σ,fillin := M σ + Q str G str Q T str is a combine (ense) matrix for both conuctor types. Lemma 1. Let the fiel problem consist of soli an strane conuctors which fulfill F1-F3, then the curl-curl equation (8a) is (algebraic) inex-1 for given voltages an has only non-trivial components in the ifferential part. Proof. The matrix pencil of (8a) is regular ue to F2, so the symmetric positive semi-efiniteness of M σ,fillin implies that (8a) is inex-1 an there is the Kronecker form t a 1 (t)+u 1 K ν V 1 a 1 (t) = U 1 j src, a 2 (t) = U 2 j src, (9a) (9b)
6 Analysis of Lumpe Electric Circuits Refine by 3-D MQS Conuctor Moels 5 that splits the MVP a = V 1 a 1 +V 2 a 2 into ifferential an algebraic components by using the regular matrices U T = ( U1 T,U 2 T ) an V = (V1, V 2 ). From U 2 M σ,fillin = U 2 ( Mσ + Q str G str Q T str) = follows that both U 2 M σ an U 2 Q str G str vanish because the images of M σ an Q str are istinct, since F3 is assume. Hence we finally conclue that a 2 = U 2 j =. Let us now stuy the full system (8) in the abstract form t a = f 2a ( a,v λ ), = f 2b ( a,v λ ), = g 2 ( a,i λ ), (1) where the voltages v λ = (v sol,v str ) T an the currents i λ = (i sol,i str ) T are combine in vectors. The algebraic evaluation f 2b is trivial in our case because of Lemma 1 an the algebraic function g 2 consists of = g sol ( a,i sol ), = g str ( a,i str ). System (1) establishes a relation between currents (i sol, i str ) an voltages (v sol, v str ) an we can choose which quantity is treate as unknown for each conuctor type in the fiel system, since then the other quantity is efine by the couple electric circuit. Therefore we will istinguish between the possible sets in the following. Theorem 2. Let the fiel problem consist of soli an strane conuctors which fulfill F1-F4, then the system (8) has (ifferential) inex-2 an only for the case of given voltages v sol an v str it has inex-1. Proof. In the case of given voltages the currents i λ are obtaine by evaluations of the algebraic equation g 2. Thus one ifferentiation with respect to time yiels an ODEs, hence we have (ifferential) inex-1. In all other cases the arguments are analogue to the case of given i sol an v str. Now the function f 2a in (1) epens on the unknown v sol an one time erivative yiels the aitional hien constraint: = t g sol( a,i sol ) = a g sol f 2 ( a,v sol )+ t i sol =: h sol ( a,v sol, t i sol), an since the conuctivity matrices M σ an M σ,aniso reflect F3 (M σ,aniso Q sol = ), another ifferentiation of this constraint gives v sol h sol = Q T sol K νm + σ,fillin M σ Q sol = Q T sol K νq sol = Q T sol CT M ν CQ sol, which is non-singular ue to F4; thus it is inex-2. We conclue, if some voltage is consiere unknown, then (8) is an inex-2 Hessenberg system (with inex-1 evaluations), which is less ill-conitione than the general case. Aitionally the inex-2 variables enter the computation linearly an without time-epenence. Therefore the ifferential components are not affecte by the erivatives of perturbations, [11], an thus the numerical ifficulties correspon to inex-1.
7 6 Sebastian Schöps, Anreas Bartel, Herbert De Gersem, an Michael Günther 5 Coupling We assign the voltages v λ to ifferences of noe potentials e an map the controlle current sources i λ to the currents through the conuctors in (1) v λ = A T λ e, i λ = (i sol,i str ) T. (11) We finally obtain the monolithic system consisting of (1), (8) an (11). Theorem 3. Let us consier an electric circuit in the form (1) with C1-C2, which is monolithically couple via (11) to a fiel moel (8) of soli an strane conuctors fulfilling F1-F4, then the full system is (ifferential) inex-1. Proof. The algebraic part of the MVP is insignificant for soli an strane conuctors accoring to Lemma 1. Hence after embeing the fiel into the circuit system the separate unknowns of the full system rea y := (P C e, j L, a 1 ) T, z := (Q C e, j V,q,φ,i λ ) T. (12) The critical partial erivative of the algebraic equation z g consisting of g 1 an g 2 is non-singular, since the first is regular ue to C1-C2 an the secon is just an evaluation of a ifferential variable ( a 1 ). Thus we have inex-1. The assumption C3 is not require in the monolithic coupling because the algebraic part of the MVP was shown to vanish for the excitement of soli an strane conuctors. Alternatively the subproblems coul be treate separately by a waveform relaxation scheme (of Jacobi or Gauß-Seiel type). When applying these schemes to DAEs one has to pay attention to algebraic constraints to avoi numerical instabilities, [12]. We suggest the following Gauß-Seiel scheme t a (1) = f 2 ( a (1),v () ), v () := A T λ (y() + z () ), t y (1) = f 1 (y (1),z (1),i (1) λ ), = g 2 ( a (1),i (1) λ ), = g 1(y (1),z (1),i (1) λ ), where the fiel is couple via the current sources i λ to the circuit. The scheme above shows irectly that there is no epenence in algebraic constraints (g 1, g 2 ) on ol algebraic iterates (i λ, z) an hence the convergence is guarantee [13] an we obtain: Lemma 2. Let us consier an electric circuit in the form (1) that respects C1-C2 is couple via (11) to a fiel moel (8) of soli an strane conuctors fulfilling F1-F4, then both subproblems are inex-1 an the waveform-relaxation above will converge. The aitional assumption C3 can eliminate the i λ -epenence of the algebraic equation g 1 an allows us to exchanges the computational orer of the subproblems an thus we may compute the circuit first without losing the convergence guarantee.
8 Analysis of Lumpe Electric Circuits Refine by 3-D MQS Conuctor Moels v(t) C R loa -1 (a) Half rectifier: v eff = 25V, f = 5Hz, R loa = 1Ω an Shockley ioe I s =1µA (b) Voltage in noes 1 an 3, obtaine by mono, H = 5µs Fig. 2: Refine half rectifier circuit an its input an compute output voltages. 6 Numerical Experiments The numerical experiments have been obtaine with coe that is implemente within the COMSON DP ( using fiel moels constructe by EM Stuio from CST ( It is capable of both, the monolithic an the co-simulation of non-linear circuits refine by (linear) conuctor moels. The iscusse schemes have been supplie to the DP, where the time integration is kept simple by using Backwar Euler. The example of Fig. 2 is a refine half rectifier with a transformer consisting of two strane conuctors an a soli core. The co-simulation without iterations of time frames (cosim1) is slightly faster than the monolithic simulation (mono) using the same time steps H an it yiels comparable results if the accuracy requirement is quite low. For ecreasing step sizes cosim1 oes not linearly improve its accuracy as mono an cosim3 (3 iterations) o (Fig. 3), but the latter has an increase computational effort ue to the aitional iterations. Aaptive time-integrators in the co-simulation apply the same step size to the circuit an the fiel, as long as they o not have multirate potential itself. This is in line with the fact that the fiel reflects the ynamics of the couple circuit noes. 7 Conclusions The fiel problem is essentially an inex-1 DAE, the monolithic couple system is still inex-1 an the convergence of the propose co-simulation is guarantee, as ocumente by the computation of a refine rectifier circuit. The co-simulation uses problem-specific software packages an exploits multirate potentials if available in the circuit, but oes not have all possible benefits. Applying a step size an iteration control will improve efficiency, an the use of more complex equivalent circuits (e.g. aitional inuctivities) might require fewer fiel upates [14, 15]. Acknowlegements This work was partially supporte by the European Commission within the framework of the COMSON RTN project.
9 8 Sebastian Schöps, Anreas Bartel, Herbert De Gersem, an Michael Günther (a) mono, H = 1µs (b) cosim3, H = 1µs (c) cosim1, H = 1µs () mono, H = 1µs (e) cosim3, H = 1µs (f) cosim1, H = 1µs Fig. 3: Errors in the voltages compare to the results of mono, H = from Fig. 2b References 1. Benerskaya, G. et al.: Transient Fiel-Circuit Couple Formulation Base on the FIT an a Mixe Circuit Formulation. COMPEL 23(4), (24) 2. Dreher, T., Meunier, G.: 3D Line Current Moel of Coils an External Circuits. IEEE Trans. Mag. 31(3), (1995) 3. Berosian, G.: A New Metho for Coupling Finite Element Fiel Solutions with External Circuits an Kinematics. IEEE Trans. Mag. 29(2), (1993) 4. Günther, M., Rentrop, P.: Numerical Simulation of Elec. Circuits. GAMM 1/2, (2) 5. Estévez Schwarz, D., Tischenorf, C.: Structural Analysis of Electric Circuits an Consequences for MNA. Int. J. Circ. Theor. Appl. 28(2), (2) 6. Bossavit, A.: Computational Electromagnetism. Acaemic Press, San Diego (1998) 7. Clemens, M.: Large Systems of Equations in a Discrete Electromagnetism: Formulations an Numerical Algorithms. IEE Proc Sci Meas Tech 152(2), 5 72 (25) 8. Nicolet, A., Delincé, F.: Implicit Runge-Kutta Methos for Transient Magnetic Fiel Computation. IEEE Trans. Mag. 32(3), (1996) 9. De Gersem, H. et al.: Fiel-circuit Couple Moels in Electromagnetic Simulation. JCAM 168(1-2), (24) 1. Kettunen, L. et al.: Gauging in Whitney Spaces. IEEE Trans. Mag. 35(3), (1999) 11. Arnol, M. et al.: Errors in the Numerical Solution of Nonlinear Differential-Algebraic Systems of Inex 2. Martin-Luther-University Halle (1995) 12. Arnol, M., Günther, M.: Preconitione Dynamic Iteration for Couple Differential- Algebraic Systems. BIT 41(1), 1 25 (21) 13. Bartel, A.: Partial Differential-Algebraic Moels in Chip Design - Thermal an Semiconuctor Problems. Ph.D. thesis, TU Karlsruhe, VDI Verlag 14. Lange, E. et al.: A Circuit Coupling Metho Base on a Temporary Linearization of the Energy Balance of the Finite Element Moel. IEEE Trans. Mag. 44(6), (28) 15. Zhou, P. et al.: A General Co-Simulation Approach for Couple Fiel-Circuit Problems. IEEE Trans. Mag. 42(4), (26)
19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationA Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation
A Novel ecouple Iterative Metho for eep-submicron MOSFET RF Circuit Simulation CHUAN-SHENG WANG an YIMING LI epartment of Mathematics, National Tsing Hua University, National Nano evice Laboratories, an
More informationConservation laws a simple application to the telegraph equation
J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness
More informationAn inductance lookup table application for analysis of reluctance stepper motor model
ARCHIVES OF ELECTRICAL ENGINEERING VOL. 60(), pp. 5- (0) DOI 0.478/ v07-0-000-y An inuctance lookup table application for analysis of reluctance stepper motor moel JAKUB BERNAT, JAKUB KOŁOTA, SŁAWOMIR
More informationAdjoint Transient Sensitivity Analysis in Circuit Simulation
Ajoint Transient Sensitivity Analysis in Circuit Simulation Z. Ilievski 1, H. Xu 1, A. Verhoeven 1, E.J.W. ter Maten 1,2, W.H.A. Schilers 1,2 an R.M.M. Mattheij 1 1 Technische Universiteit Einhoven; e-mail:
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationAPPLICATION OF DIFFERENTIAL FORMS IN THE FINITE ELEMENT FORMULATION OF ELECTROMAGNETIC PROBLEMS
Technical article APPLICATION OF DIFFERENTIAL FORMS IN THE FINITE ELEMENT FORMULATION OF ELECTROMAGNETIC PROBLEMS INTRODUCTION In many physical problems, we have to stuy the integral of a quantity over
More informationRobust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k
A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine
More informationA Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges
Plasma Science an Technology, Vol.16, No.1, Oct. 214 A Simple Moel for the Calculation of Plasma Impeance in Atmospheric Raio Frequency Discharges GE Lei ( ) an ZHANG Yuantao ( ) Shanong Provincial Key
More informationConservation Laws. Chapter Conservation of Energy
20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action
More informationStable and compact finite difference schemes
Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long
More informationECE 422 Power System Operations & Planning 7 Transient Stability
ECE 4 Power System Operations & Planning 7 Transient Stability Spring 5 Instructor: Kai Sun References Saaat s Chapter.5 ~. EPRI Tutorial s Chapter 7 Kunur s Chapter 3 Transient Stability The ability of
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationA simple model for the small-strain behaviour of soils
A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationNoether s theorem applied to classical electrodynamics
Noether s theorem applie to classical electroynamics Thomas B. Mieling Faculty of Physics, University of ienna Boltzmanngasse 5, 090 ienna, Austria (Date: November 8, 207) The consequences of gauge invariance
More informationELEC3114 Control Systems 1
ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More informationA Comparison between a Conventional Power System Stabilizer (PSS) and Novel PSS Based on Feedback Linearization Technique
J. Basic. Appl. Sci. Res., ()9-99,, TextRoa Publication ISSN 9-434 Journal of Basic an Applie Scientific Research www.textroa.com A Comparison between a Conventional Power System Stabilizer (PSS) an Novel
More informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More informationApplication of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate
Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite
More informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
More informationTractability results for weighted Banach spaces of smooth functions
Tractability results for weighte Banach spaces of smooth functions Markus Weimar Mathematisches Institut, Universität Jena Ernst-Abbe-Platz 2, 07740 Jena, Germany email: markus.weimar@uni-jena.e March
More informationAll s Well That Ends Well: Supplementary Proofs
All s Well That Ens Well: Guarantee Resolution of Simultaneous Rigi Boy Impact 1:1 All s Well That Ens Well: Supplementary Proofs This ocument complements the paper All s Well That Ens Well: Guarantee
More informationExponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity
Preprints of the 9th Worl Congress The International Feeration of Automatic Control Cape Town, South Africa. August -9, Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity Zhengqiang
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationDesigning Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations
EECS 6B Designing Information Devices an Systems II Fall 07 Note 3 Secon Orer Differential Equations Secon orer ifferential equations appear everywhere in the real worl. In this note, we will walk through
More informationLagrangian and Hamiltonian Dynamics
Lagrangian an Hamiltonian Dynamics Volker Perlick (Lancaster University) Lecture 1 The Passage from Newtonian to Lagrangian Dynamics (Cockcroft Institute, 22 February 2010) Subjects covere Lecture 2: Discussion
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationIPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy
IPA Derivatives for Make-to-Stock Prouction-Inventory Systems With Backorers Uner the (Rr) Policy Yihong Fan a Benamin Melame b Yao Zhao c Yorai Wari Abstract This paper aresses Infinitesimal Perturbation
More informationInvestigation of local load effect on damping characteristics of synchronous generator using transfer-function block-diagram model
ORIGINAL ARTICLE Investigation of local loa effect on amping characteristics of synchronous generator using transfer-function block-iagram moel Pichai Aree Abstract of synchronous generator using transfer-function
More informationOptimal LQR Control of Structures using Linear Modal Model
Optimal LQR Control of Structures using Linear Moal Moel I. Halperin,2, G. Agranovich an Y. Ribakov 2 Department of Electrical an Electronics Engineering 2 Department of Civil Engineering Faculty of Engineering,
More informationModeling time-varying storage components in PSpice
Moeling time-varying storage components in PSpice Dalibor Biolek, Zenek Kolka, Viera Biolkova Dept. of EE, FMT, University of Defence Brno, Czech Republic Dept. of Microelectronics/Raioelectronics, FEEC,
More information4. CONTROL OF ZERO-SEQUENCE CURRENT IN PARALLEL THREE-PHASE CURRENT-BIDIRECTIONAL CONVERTERS
4. CONRO OF ZERO-SEQUENCE CURREN IN PARAE HREE-PHASE CURREN-BIDIRECIONA CONVERERS 4. A NOVE ZERO-SEQUENCE CURREN CONRO 4.. Zero-Sequence Dynamics he parallel boost rectifier moel in Figure.4 an the parallel
More informationPolynomial Inclusion Functions
Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl
More informationHarmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method
1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse
More informationSituation awareness of power system based on static voltage security region
The 6th International Conference on Renewable Power Generation (RPG) 19 20 October 2017 Situation awareness of power system base on static voltage security region Fei Xiao, Zi-Qing Jiang, Qian Ai, Ran
More informationHyperbolic Systems of Equations Posed on Erroneous Curved Domains
Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (
More informationRevisiting the Charge Concept in HBT/BJT Models
evisiting the Charge Concept in HBT/BJT Moels Zoltan Huska an Ehrenfrie Seebacher austriamicrosystems AG 23r Bipolar Arbeitkeis BipAK Meeting at STM Crolles, France, 5 October 2 Outline recalling the junction
More informationSliding mode approach to congestion control in connection-oriented communication networks
JOURNAL OF APPLIED COMPUTER SCIENCE Vol. xx. No xx (200x), pp. xx-xx Sliing moe approach to congestion control in connection-oriente communication networks Anrzej Bartoszewicz, Justyna Żuk Technical University
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More informationOn state representations of time-varying nonlinear systems
On state representations of time-varying nonlinear systems Paulo Sérgio Pereira a Silva a, Simone Batista a, a University of São Paulo, Escola Politécnicca PTC Av. Luciano Gualberto trav. 03, 158, 05508-900
More informationTHE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS
THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS
More informationPRACTICE 4. CHARGING AND DISCHARGING A CAPACITOR
PRACTICE 4. CHARGING AND DISCHARGING A CAPACITOR. THE PARALLEL-PLATE CAPACITOR. The Parallel plate capacitor is a evice mae up by two conuctor parallel plates with total influence between them (the surface
More informationChaos, Solitons and Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena
Chaos, Solitons an Fractals (7 64 73 Contents lists available at ScienceDirect Chaos, Solitons an Fractals onlinear Science, an onequilibrium an Complex Phenomena journal homepage: www.elsevier.com/locate/chaos
More informationIntroduction to variational calculus: Lecture notes 1
October 10, 2006 Introuction to variational calculus: Lecture notes 1 Ewin Langmann Mathematical Physics, KTH Physics, AlbaNova, SE-106 91 Stockholm, Sween Abstract I give an informal summary of variational
More informationLower Bounds for the Smoothed Number of Pareto optimal Solutions
Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.
More informationWe G Model Reduction Approaches for Solution of Wave Equations for Multiple Frequencies
We G15 5 Moel Reuction Approaches for Solution of Wave Equations for Multiple Frequencies M.Y. Zaslavsky (Schlumberger-Doll Research Center), R.F. Remis* (Delft University) & V.L. Druskin (Schlumberger-Doll
More informationThe canonical controllers and regular interconnection
Systems & Control Letters ( www.elsevier.com/locate/sysconle The canonical controllers an regular interconnection A.A. Julius a,, J.C. Willems b, M.N. Belur c, H.L. Trentelman a Department of Applie Mathematics,
More informationGoal of this chapter is to learn what is Capacitance, its role in electronic circuit, and the role of dielectrics.
PHYS 220, Engineering Physics, Chapter 24 Capacitance an Dielectrics Instructor: TeYu Chien Department of Physics an stronomy University of Wyoming Goal of this chapter is to learn what is Capacitance,
More informationDYNAMIC PERFORMANCE OF RELUCTANCE SYNCHRONOUS MACHINES
Annals of the University of Craiova, Electrical Engineering series, No 33, 9; ISSN 184-485 7 TH INTERNATIONAL CONFERENCE ON ELECTROMECHANICAL AN POWER SYSTEMS October 8-9, 9 - Iaşi, Romania YNAMIC PERFORMANCE
More informationTwo Dimensional Numerical Simulator for Modeling NDC Region in SNDC Devices
Journal of Physics: Conference Series PAPER OPEN ACCESS Two Dimensional Numerical Simulator for Moeling NDC Region in SNDC Devices To cite this article: Dheeraj Kumar Sinha et al 2016 J. Phys.: Conf. Ser.
More informationA new proof of the sharpness of the phase transition for Bernoulli percolation on Z d
A new proof of the sharpness of the phase transition for Bernoulli percolation on Z Hugo Duminil-Copin an Vincent Tassion October 8, 205 Abstract We provie a new proof of the sharpness of the phase transition
More informationOn the number of isolated eigenvalues of a pair of particles in a quantum wire
On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationSensors & Transducers 2015 by IFSA Publishing, S. L.
Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity
More informationThe effect of nonvertical shear on turbulence in a stably stratified medium
The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the
More informationOptimization of Geometries by Energy Minimization
Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.
More informationCONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS
BIT 0006-3835/00/4004-0001 $15.00 200?, Vol.??, No.??, pp.?????? c Swets & Zeitlinger CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS APPLIE TO THE SHALLOW WATER EQUATIONS JASON FRANK 1 an SEBASTIAN
More informationarxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003
Mass reistribution in variable mass systems Célia A. e Sousa an Vítor H. Rorigues Departamento e Física a Universiae e Coimbra, P-3004-516 Coimbra, Portugal arxiv:physics/0211075v2 [physics.e-ph] 23 Sep
More informationMartin Luther Universität Halle Wittenberg Institut für Mathematik
Martin Luther Universität alle Wittenberg Institut für Mathematik Weak solutions of abstract evolutionary integro-ifferential equations in ilbert spaces Rico Zacher Report No. 1 28 Eitors: Professors of
More informationLecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
More informationSYMMETRIC KRONECKER PRODUCTS AND SEMICLASSICAL WAVE PACKETS
SYMMETRIC KRONECKER PRODUCTS AND SEMICLASSICAL WAVE PACKETS GEORGE A HAGEDORN AND CAROLINE LASSER Abstract We investigate the iterate Kronecker prouct of a square matrix with itself an prove an invariance
More informationThe continuity equation
Chapter 6 The continuity equation 61 The equation of continuity It is evient that in a certain region of space the matter entering it must be equal to the matter leaving it Let us consier an infinitesimal
More information. Using a multinomial model gives us the following equation for P d. , with respect to same length term sequences.
S 63 Lecture 8 2/2/26 Lecturer Lillian Lee Scribes Peter Babinski, Davi Lin Basic Language Moeling Approach I. Special ase of LM-base Approach a. Recap of Formulas an Terms b. Fixing θ? c. About that Multinomial
More informationThe total derivative. Chapter Lagrangian and Eulerian approaches
Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function
More informationTIME-DELAY ESTIMATION USING FARROW-BASED FRACTIONAL-DELAY FIR FILTERS: FILTER APPROXIMATION VS. ESTIMATION ERRORS
TIME-DEAY ESTIMATION USING FARROW-BASED FRACTIONA-DEAY FIR FITERS: FITER APPROXIMATION VS. ESTIMATION ERRORS Mattias Olsson, Håkan Johansson, an Per öwenborg Div. of Electronic Systems, Dept. of Electrical
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationHomework 7 Due 18 November at 6:00 pm
Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine
More informationComputing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions
Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5
More informationarxiv: v1 [cond-mat.stat-mech] 9 Jan 2012
arxiv:1201.1836v1 [con-mat.stat-mech] 9 Jan 2012 Externally riven one-imensional Ising moel Amir Aghamohammai a 1, Cina Aghamohammai b 2, & Mohamma Khorrami a 3 a Department of Physics, Alzahra University,
More informationDynamics of the synchronous machine
ELEC0047 - Power system ynamics, control an stability Dynamics of the synchronous machine Thierry Van Cutsem t.vancutsem@ulg.ac.be www.montefiore.ulg.ac.be/~vct These slies follow those presente in course
More informationProblem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs
Problem Sheet 2: Eigenvalues an eigenvectors an their use in solving linear ODEs If you fin any typos/errors in this problem sheet please email jk28@icacuk The material in this problem sheet is not examinable
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationTRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS
TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS Francesco Bullo Richar M. Murray Control an Dynamical Systems California Institute of Technology Pasaena, CA 91125 Fax : + 1-818-796-8914 email
More informationThe effect of dissipation on solutions of the complex KdV equation
Mathematics an Computers in Simulation 69 (25) 589 599 The effect of issipation on solutions of the complex KV equation Jiahong Wu a,, Juan-Ming Yuan a,b a Department of Mathematics, Oklahoma State University,
More informationTowards One-Step Multirate Methods in Chip Design
Bergische Universität Wuppertal Fachbereich Mathematik und Naturwissenschaften Lehrstuhl für Angewandte Mathematik und Numerische Mathematik Preprint BUW-AMNA 04/09 M. Striebel, M.Günther Towards One-Step
More informationPure Further Mathematics 1. Revision Notes
Pure Further Mathematics Revision Notes June 20 2 FP JUNE 20 SDB Further Pure Complex Numbers... 3 Definitions an arithmetical operations... 3 Complex conjugate... 3 Properties... 3 Complex number plane,
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationDesign and Application of Fault Current Limiter in Iran Power System Utility
Australian Journal of Basic an Applie Sciences, 7(): 76-8, 13 ISSN 1991-8178 Design an Application of Fault Current Limiter in Iran Power System Utility M. Najafi, M. Hoseynpoor Department of Electrical
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More informationSlide10 Haykin Chapter 14: Neurodynamics (3rd Ed. Chapter 13)
Slie10 Haykin Chapter 14: Neuroynamics (3r E. Chapter 13) CPSC 636-600 Instructor: Yoonsuck Choe Spring 2012 Neural Networks with Temporal Behavior Inclusion of feeback gives temporal characteristics to
More informationNested Saturation with Guaranteed Real Poles 1
Neste Saturation with Guarantee Real Poles Eric N Johnson 2 an Suresh K Kannan 3 School of Aerospace Engineering Georgia Institute of Technology, Atlanta, GA 3332 Abstract The global stabilization of asymptotically
More informationNonlinear Dielectric Response of Periodic Composite Materials
onlinear Dielectric Response of Perioic Composite aterials A.G. KOLPAKOV 3, Bl.95, 9 th ovember str., ovosibirsk, 639 Russia the corresponing author e-mail: agk@neic.nsk.su, algk@ngs.ru A. K.TAGATSEV Ceramics
More informationGeneralization of the persistent random walk to dimensions greater than 1
PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,
More informationEVALUATING HIGHER DERIVATIVE TENSORS BY FORWARD PROPAGATION OF UNIVARIATE TAYLOR SERIES
MATHEMATICS OF COMPUTATION Volume 69, Number 231, Pages 1117 1130 S 0025-5718(00)01120-0 Article electronically publishe on February 17, 2000 EVALUATING HIGHER DERIVATIVE TENSORS BY FORWARD PROPAGATION
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More informationEnergy behaviour of the Boris method for charged-particle dynamics
Version of 25 April 218 Energy behaviour of the Boris metho for charge-particle ynamics Ernst Hairer 1, Christian Lubich 2 Abstract The Boris algorithm is a wiely use numerical integrator for the motion
More informationANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS
ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics
More informationLecture XVI: Symmetrical spacetimes
Lecture XVI: Symmetrical spacetimes Christopher M. Hirata Caltech M/C 350-17, Pasaena CA 91125, USA (Date: January 4, 2012) I. OVERVIEW Our principal concern this term will be symmetrical solutions of
More informationHow the potentials in different gauges yield the same retarded electric and magnetic fields
How the potentials in ifferent gauges yiel the same retare electric an magnetic fiels José A. Heras a Departamento e Física, E. S. F. M., Instituto Politécnico Nacional, México D. F. México an Department
More informationLeast-Squares Regression on Sparse Spaces
Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction
More informationarxiv: v1 [physics.class-ph] 20 Dec 2017
arxiv:1712.07328v1 [physics.class-ph] 20 Dec 2017 Demystifying the constancy of the Ermakov-Lewis invariant for a time epenent oscillator T. Pamanabhan IUCAA, Post Bag 4, Ganeshkhin, Pune - 411 007, Inia.
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationTransmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency
Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com
More information