the solution of ()-(), an ecient numerical treatment requires variable steps. An alternative approach is to apply a time transformation of the form t
|
|
- Claude Thomas
- 5 years ago
- Views:
Transcription
1 Asymptotic Error Analysis of the Aaptive Verlet Metho Stephane Cirilli, Ernst Hairer Beneict Leimkuhler y May 3, 999 Abstract The Aaptive Verlet metho [7] an variants [6] are time-reversible schemes for treating Hamiltonian systems subject to a Sunman time transformation. These methos have been observe in computer experiments to exhibit superior numerical stability when implemente in a counterintuitive \reciprocal" formulation. Here we give a theoretical explanation of this behavior by examining the leaing terms of the moie equation (backwar error analysis) an those of the asymptotic error expansion. With this insight we are able to improve the algorithm by simply correcting the starting stepsize. keywors: Aaptive Verlet metho, time-reversible variable stepsizes, Hamiltonian systems, Sunman time-transformations, backwar error analysis, asymptotic expansions Introuction Consier a mechanical Hamiltonian system of the form t q = M p; q(t 0 ) = q 0 ; () t p = F (q); p(t 0) = p 0 ; () where q; p R N, M R N N is a positive enite matrix, an F = r q V for some smooth, real-value potential function V. If ierent time-scales are present in Section e mathematiques, Universite e Geneve, -4 rue e Lievre, CH- Geneve 4, Switzerlan. Stephane.Cirilli@math.unige.ch, Ernst.Hairer@math.unige.ch y Dept. of Mathematics, 405 Snow Hall, University of Kansas, Lawrence, KS 66045, USA. leimkuhl@math.ukans.eu. Visiting DAMTP, Cambrige University, AY Research supporte by NSF grant No an DAMTP.
2 the solution of ()-(), an ecient numerical treatment requires variable steps. An alternative approach is to apply a time transformation of the form t = G(q; p); with G a smooth, positive function (also calle Sunman transformation) so that, by viewing q an p as functions of the new variable, the equations of motion become q = gm p; (3) p = gf (q); (4) where g = G(q; p). Applying a numerical metho with constant stepsize h to (3)- (4) is essentially equivalent to applying the metho with stepsize hg n to ()-(). If G(q; p) = G(q; p), the change of the time variable preserves the reversing symmetry associate to ()-(), a property that has been foun to be important for recovering realistic qualitative behavior of the original system in numerical simulations [0, 8]. In this article we consier a variant of Verlet's metho written in the form q n+= = q n + h g nm p n+= ; (5) p n+= = p n + h g nf (q n ); (6) q n+ = q n+= + h g n+m p n+= ; (7) p n+ = p n+= + h g n+f (q n+ ): (8) Here h represents the ctive timestep (i.e. stepsize in ), an the inices inicate the timestep. The most obvious choice of g n is g n = G(q n ; p n ); (9) which, however, results in an implicit algorithm. In orer to avoi this implicitness, it was suggeste in [7] to put g 0 = G(q 0 ; p 0 ), an to upate the time-scaling factor g from a two-term recurrence relation (g n ; g n+ ; q n ; q n+= ; q n+ ; p n ; p n+= ; p n+ ) = 0 with symmetry preserving. Obviously, the latter equation has to be consistent with the equation g( n ) = G(q( n ); p( n )) in the limit of small h. A natural choice of is = g n+ + g n G(q n+= ; p n+= ); (0) but in experiments this has been foun to behave unreliably, particularly in the presence of large forces. A more robust approach is to use the corresponing equation for the reciprocal of g [6], = + g n+ g n G(q n+= ; p n+= ) : ()
3 Various semi-explicit an implicit variants of Aaptive Verlet propose in [7] also use the reciprocal relation. Until now, the reason for this counterintuitive moication of the recurrence relation has not been explaine. In Sect., we stuy for the various methos the leaing terms of the moie equation (backwar error analysis) an those of the asymptotic error expansions. Each of these has a (nonphysical) oscillating part etermine from a certain ierential or algebraic equation; a simple analysis of this leaing oscillatory term explains the observe ierences in stability of the methos. As an outcome of this analysis, we show in Sect. 3 how the oscillations can be attenuate by correcting the starting stepsize. Another extension of the Verlet scheme to variable stepsizes has recently been propose inepenently in [3] an in [9]. In aition to being symmetric, it is a symplectic metho, but it has the isavantage of being implicit in g n+. A comparison of the ierent extensions is given in []. Backwar Error Analysis an Asymptotic Expansions For the choice g n = G(q n ; p n ), the metho is a symmetric one-step metho applie to an orinary ierential equation. It is a classical result that the numerical solution can then be formally written as the exact solution of a moie ierential equation, an that the global error possesses an expansion in even powers of h (see [5, pp ] an [4, Sect. II.8]). We now turn to the more interesting case where g n is given by a two-term recursion, which allows the algorithm to be completely explicit. The following theorem shows that in this case the numerical solution is a superposition of a smooth function with oscillatory terms. Such a phenomenon is similar to the weak instability of two-step methos (see [4, Sect. III.9]). Theorem The numerical solution of the Aaptive Verlet metho (5)-(8) can formally be written as q n = eq( n )+( ) n bq( n ); p n = ep( n )+( ) n bp( n ); g n = eg( n )+( ) n bg( n ); () where n = nh an eq 0 = egm ep + h e Q () + : : : ; bq = h b Q () + : : : ; (3) ep 0 = egf (eq ) + h e P () + : : : ; bp = h b P () + : : : ; (4) eg = G(eq; ep ) + h e G () + : : : ; bg 0 = b G() + h b G () + : : : ; (5) with uniquely etermine initial values satisfying q 0 = eq(0) + bq(0); p 0 = ep(0) + bp(0); g 0 = eg(0) + bg(0): (6) All of the above expansions are formal an in even powers of h. The functions in (3)-(5) epen only on eq; ep an bg. The functions G; b G b ; : : : an also Q b ; P b ; : : : contain bg as factor, an bg(0) = O(h ). 3
4 Remark. The right han sie of the ierential equation for bg (5) contains bg as a factor, an bg(0) = O(h ); these facts imply that bg() = O(h ) for on compact intervals. Since, the functions Q b an P b also contain bg as factor, we have in aition bq() = O(h 4 ) an bp() = O(h 4 ). Proof. The iea is to insert () into the metho (5)-(8), to expan, an to compare like powers of h. This computation is signicantly simplie if we exploit the symmetric structure of the metho an if we expan the resulting expressions aroun := n +h=. Neglecting terms of orer O(h 4 ) an those of the form O(h bq ), O(h bp ), O(bq ), O(bq bp ), O( bp )) an omitting the obvious argument, we get (p n+ + p n ) = ep + h 8 ep 00 ( ) n h bp 0 + : : : h 4 (g n g n+ ) = h 4 eg h 0 + ( ) n h3 00 bg + bg + : : : : 6 Using these an similar relations for other variables, we obtain for the intermeiate approximations q n+= an p n+= the expansions p n+= = (p n+ + p n ) + h 4 g n F (q n ) g n+ F (q n+ ) = ep + h ep 00 (eg F (eq )) 0 (7) 8 ( ) h n bp 0 bgf (eq ) egf q (eq )bq + ( ) n h3 6 (bgf (eq ))00 + : : : ; q n+= = (q n+ + q n ) + h 4 (g n g n+ )M p n+= = eq + h eq 00 eg 0 M ep + bg M F (eq ) ( ) h n bq 0 bgm ep 8 +( ) n h3 6 bgm ep 00 + bg 00 M ep bg(eg M F (eq )) 0 bg eg 0 M F (eq ) + : : : ; (8) where F q (q) enotes the erivative of F with repect to q, an 0 enotes the erivative with respect to. We now write the main formulas of the metho (with the choice (0) for g n ) as q n+ q n = h (g n + g n+ )M p n+= ; (9) p n+ p n = h g n F (q n ) + g n+ F (q n+ ) ; (0) (g n+ + g n ) = G(q n+= ; p n+= ): () Inserting the relations () an (7)-(8), we get 4
5 heq 0 + h3 4 eq 000 ( ) n bq + : : : () = h eg + h 8 eg 00 ( ) h n M bg 0 ep + h 8 ( ep 00 (eg F (eq )) 0 ) +( ) n h bgf (eq ) + : : : ; h ep 0 + h3 4 ep 000 ( ) n bp + : : : (3) = h egf (eq ) + bgf q (eq )bq + h3 8 (egf (eq ))00 ( ) n h (bgf (eq ))0 + : : : ; eg + h 8 eg 00 ( ) n h bg 0 ( ) n h3 48 bg : : : (4) = G(eq; ep ) + G q (eq; ep )(q n+= eq ) + G p (eq; ep )(p n+= ep ) + : : : ; where G q an G p enote partial erivatives. Comparing the non-oscillating an oscillating parts in ()-(4) yiels equations for eq 0 ; ep 0 ; eg an bq; bp; bg 0, respectively. The higher erivatives appearing in the right-han sie have to be eliminate iteratively. This gives the moie equations (3)-(5). We obtain for example bg(eq; ep; bg ) = bg G q (eq; ep )M ep + G p (eq; ep )F (eq ) : (5) The algebraic relations for bq; bp; eg in (3)-(5) together with (6) constitute a set of six equations for the unknowns eq(0); bq(0); ep(0); bp(0); eg(0); bg(0). By the implicit function theorem they have a unique solution close to (q 0 ; 0; p 0 ; 0; g 0 ; 0),which can be written as a formal series in powers of h. Remark. For the choice () only the equations () an (4) have to be aapte, an one gets bg(eq; ep; bg ) = +bg G q (eq; ep )M ep + G p (eq; ep )F (eq ) ; (6) instea of (5). For the more general situation L(g n+ ) + L(g n ) = L G(q n+= ; p n+= ) ; which inclues (0) an () as special cases, we get bg(eq; ep; bg ) = bg + L gg(eg )eg G q (eq; ep )M ep + G p (eq; ep )F (eq ) : (7) L g (eg ) Corollary (Asymptotic Expansions) The functions eq(), bq(), ep(), bp(), eg(), bg() of Theorem all have asymptotic expansions in even powers of h. In particular, eg() = g() + h g () + O(h 4 ); bg() = h bg () + h 4 bg 4 () + O(h 6 ); (8) where g() = G(q(); p()) with (q(); p()) the solution of (3)-(4). 5
6 For the choice (0) of g n with G(q) only epening on q, bg () = G(q 0) G(q 0 )G qq (q 0 )(M p 0 ; M p 0 ) + (G q (q 0 )M p 0 ) ; (9) 8g() an for the choice () we have: bg () = g() 8 G(q 0)G qq (q 0 )(M p 0 ; M p 0 ): (30) Proof. We insert the formulas (8) an similar equations for eq(), bq(), ep(), bp() into the moie ierential-algebraic system (3)-(5), then compare like powers of h. This yiels ierential equations for eq i (), ep i (), bg i (), an algebraic relations for bq i (), bp i (), eg i (). A straightforwar computation gives (for (0)): g () = 8 g00 () + G q (q(); p()) 8 q00 () + q () 4 g0 ()M p() +G p (q(); p()) p () 8 p00 () ; (3) bg 0 () = bg ()g 0 ()=g(): (3) The ierential equation (3) can be solve for bg (), an we obtain bg ()g() = bg (0)g(0). Since q (0) = p (0) = 0, we get g (0) from (3) an then bg (0) from (6). This yiels the formula (9) for bg (). Formula (30) is obtaine in the same way. Example For an illustration of the oscillatory terms in the numerical solution an in the stepsizes we consier the problem q 0 = p; p 0 = =q ; G(q) = q ; with initial values q 0 = ; p 0 =. We apply the aaptive Verlet metho with ctive stepsize h = 0:08. The picture on the left in Fig. shows the values of g n as a function of the time t n, (i) for the choice (0) (inicate by small circles), an (ii) for the choice () (small squares). The secon an thir pictures of Fig. show the values g n g(t n ) together with the smooth curves h bg (t). For the choice (0) we observe increasing oscillations, an for t 0: the numerical solution becomes meaningless an the stepsizes hg n even become negative. We have joine consecutive points of fg n g by a polygon in orer to better illustrate this phenomenon. For the choice () the oscillations are very small (they can be observe only in the scale of the thir picture) an they are ecreasing. This can be explaine by the fact that bg (t) is inversely proportional to g(t) for the choice (0), whereas it is proportional to g(t) for the choice (), an g(t) = q(t) is approaching zero. The secon an thir pictures show excellent agreement of g n with the expansion g(t n )+h g (t n )+( ) n h bg (t n )+O(h 4 ), even for the rather large value h = 0:08 (observe that in the last two pictures the term h g (t) is not inclue in the smooth function). 6
7 Figure : The values of g n (left picture) an g n g(t n ) (secon an thir pictures) as a function of time t n ; circles inicate the results for the choice (0), an squares for the choice (). 3 Elimination of the Dominant Oscillatory Terms In the proof of Corollary, we have seen that the function bg () is the solution of a linear autonomous ierential equation. If we are able to achieve bg (0) = 0, the function bg () will remain zero for all an the ominant term in bg() will be eliminate, implying bg() = O(h 4 ). Since bg is a factor in the algebraic relations (3)-(4) for bq an bp, this will also imply that bq() = O(h 6 ) an bp() = O(h 6 ). The iea is to use g 0 = G(q 0 ; p 0 ) + h + 4 h 4 + : : : (33) for the initial stepsize instea of g 0 = G(q 0 ; p 0 ). Such a choice neither aects the moie ierential-algebraic system (3)-(5) nor the relations (6) for its initial values. If (33) is an expansion in even powers of h, the symmetry of the metho will not be estroye. If g 0 is compute from (33), the value of bg (0) is given by = g (0) + bg (0): Therefore, the choice = g (0) with g (0) given by (3) implies bg (0) = 0. Example For the problem of Example it is not icult to compute g (0) from (3). We have g (0) = 5 for the choice (0) an g (0) = for (). With the same ctive stepsize h = 0:08 an with the same initial values as before we apply the aaptive Verlet scheme, but we use g 0 = G(q 0 ) + h g (0) instea of g 0 = G(q 0 ). The result is shown in Fig.. The values of g n (rst picture) an g n g(t n ) for the choice (0) (secon picture) still show an oscillating behaviour, but the magnitue of the oscillations is smaller an proportional to h 4. For the choice () (small squares) the oscillations are completely eliminate. For general problems, the computation of the value g (0) with the help of a formula like (3) can be cumbersome, in particular if G also epens on p. It is of course possible to use automatic ierentiation []. Another possibility is to apply the aaptive Verlet scheme with g 0 = G(q 0 ; p 0 ) two steps with positive ctive stepsize 7
8 Figure : The same interpretation as in Fig., but with g 0 = G(q 0 ) + h g (0) instea of g 0 = G(q 0 ). an two steps with (for something like eps =4, where eps is the machine precision, the rouno an truncation errors are of the same orer). From the values g ; g ; g 0 ; g ; g, obtaine in this way, we then compute the fourth central ierence 4 = g 4g + 6g 0 4g + g ; which eliminates the smooth terms up to orer four an yiels 4 = 6 bg (0)+O( 4 ). With the starting stepsize g 0 = G(q 0 ; p 0 ) h 6 4 (34) we therefore eliminate the ominant oscillatory terms. We repeate the numerical experiment of Example with g 0 from (34) an observe the same eect as with g 0 = G(q 0 ) + h g (0). Acknowlegement. We are grateful to Gerhar Wanner for his useful comments on an earlier version of this paper. References [] M.P. Calvo, M.A. Lopez-Marcos an J.M. Sanz-Serna, Variable Step Implementation of Geometric Integrators, Appl. Numer. Math., to appear. [] G.F. Corliss an A. Griewank (es), Automatic Dierentiation of Algorithms. Theory, Implementation, an Application, SIAM, Philaelphia, PA, 99. [3] E. Hairer, Variable Timestep Integration with Symplectic Methos, Appl. Numer. Math. 5, 9-7, 997. [4] E. Hairer, S.P. Nrsett an G. Wanner, Solving Orinary Dierential Equations I. Nonsti Problems, Springer Series in Comput. Math., vol. 8, n eition, Springer-Verlag,
9 [5] E. Hairer an G. Wanner, Solving Orinary Dierential Equations II. Sti an Dierential-Algebraic Problems, Springer Series in Comput. Math., vol. 4, n eition, Springer-Verlag, 996. [6] Th. Holer, B. Leimkuhler an S. Reich, Explicit, Variable Stepsize, an Timereversible Integration, preprint, submitte, 998. [7] W. Huang an B. Leimkuhler, The Aaptive Verlet Metho, SIAM J. Sci. Comput. 8, 39-56, 997. [8] B. Leimkuhler, Reversible Aaptive Regularization I: Perturbe Kepler Motion an Classical Atomic Trajectories, Report 997/NA-08, DAMTP, University of Cambrige, 997. [9] S. Reich, Backwar Error Analysis for Numerical Integrators, SIAM J. Numer. Anal., to appear. [0] D. Stoer, Variable Steps for Reversible Integration Methos, Computing 55,-,
Energy 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 informationKonrad-Zuse-Zentrum für Informationstechnik Berlin
Konra-use-entrum für Informationstechnik Berlin Christof Schutte Smoothe Molecular Dynamics For Thermally Embee Systems SC 95-4 (May 995) Smoothe Molecular Dynamics For Thermally Embee Systems Christof
More information1 Introuction In the past few years there has been renewe interest in the nerson impurity moel. This moel was originally propose by nerson [2], for a
Theory of the nerson impurity moel: The Schrieer{Wol transformation re{examine Stefan K. Kehrein 1 an nreas Mielke 2 Institut fur Theoretische Physik, uprecht{karls{universitat, D{69120 Heielberg, F..
More informationSYMMETRIC PROJECTION METHODS FOR DIFFERENTIAL EQUATIONS ON MANIFOLDS
BIT 0006-3835/00/4004-0726 $15.00 2000, Vol. 40, No. 4, pp. 726 734 c Swets & Zeitlinger SYMMETRIC PROJECTION METHODS FOR DIFFERENTIAL EQUATIONS ON MANIFOLDS E. HAIRER Section de mathématiques, Université
More informationReducing round-off errors in symmetric multistep methods
Reducing round-off errors in symmetric multistep methods Paola Console a, Ernst Hairer a a Section de Mathématiques, Université de Genève, 2-4 rue du Lièvre, CH-1211 Genève 4, Switzerland. (Paola.Console@unige.ch,
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More informationGeometric Numerical Integration
Geometric Numerical Integration (Ernst Hairer, TU München, winter 2009/10) Development of numerical ordinary differential equations Nonstiff differential equations (since about 1850), see [4, 2, 1] Adams
More informationAdvanced Partial Differential Equations with Applications
MIT OpenCourseWare http://ocw.mit.eu 18.306 Avance Partial Differential Equations with Applications Fall 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.eu/terms.
More informationBackward error analysis
Backward error analysis Brynjulf Owren July 28, 2015 Introduction. The main source for these notes is the monograph by Hairer, Lubich and Wanner [2]. Consider ODEs in R d of the form ẏ = f(y), y(0) = y
More informationSOME PROPERTIES OF SYMPLECTIC RUNGE-KUTTA METHODS
SOME PROPERTIES OF SYMPLECTIC RUNGE-KUTTA METHODS ERNST HAIRER AND PIERRE LEONE Abstract. We prove that to every rational function R(z) satisfying R( z)r(z) = 1, there exists a symplectic Runge-Kutta method
More informationAccelerate Implementation of Forwaring Control Laws using Composition Methos Yves Moreau an Roolphe Sepulchre June 1997 Abstract We use a metho of int
Katholieke Universiteit Leuven Departement Elektrotechniek ESAT-SISTA/TR 1997-11 Accelerate Implementation of Forwaring Control Laws using Composition Methos 1 Yves Moreau, Roolphe Sepulchre, Joos Vanewalle
More informationEfficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators
MATEMATIKA, 8, Volume 3, Number, c Penerbit UTM Press. All rights reserved Efficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators Annie Gorgey and Nor Azian Aini Mat Department of Mathematics,
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 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 informationContents. 2.1 Motivation: Rates and Tangent Lines. Calculus I (part 2): Introduction to Dierentiation (by Evan Dummit, 2016, v. 2.
Calculus I (part 2): Introuction to Dierentiation (by Evan Dummit, 2016, v 250) Contents 2 Introuction to Dierentiation 1 21 Motivation: Rates an Tangent Lines 1 22 Formal Denition of the Derivative 3
More informationRank, Trace, Determinant, Transpose an Inverse of a Matrix Let A be an n n square matrix: A = a11 a1 a1n a1 a an a n1 a n a nn nn where is the jth col
Review of Linear Algebra { E18 Hanout Vectors an Their Inner Proucts Let X an Y be two vectors: an Their inner prouct is ene as X =[x1; ;x n ] T Y =[y1; ;y n ] T (X; Y ) = X T Y = x k y k k=1 where T an
More informationTHE Hamilton equations of motion constitute a system
Proceedings of the World Congress on Engineering 0 Vol I WCE 0, July 4-6, 0, London, U.K. Systematic Improvement of Splitting Methods for the Hamilton Equations Asif Mushtaq, Anne Kværnø, and Kåre Olaussen
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
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 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 information1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.
Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency
More informationd dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1
Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of
More informationThe Exact Form and General Integrating Factors
7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily
More informationinitial configuration initial configuration end
Design of trajectory stabilizing feeback for riftless at systems M. FLIESS y J. L EVINE z P. MARTIN x P. ROUCHON { ECC95 Abstract A esign metho for robust stabilization of riftless at systems aroun trajectories
More informationNOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,
NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which
More informationReducing round-off errors in rigid body dynamics
Reducing round-off errors in rigid body dynamics Gilles Vilmart INRIA Rennes, ENS Cachan Bretagne, avenue Robert Schuman, 357 Bruz, France Université de Genève, Section de mathématiques, -4 rue du Lièvre,
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 informationGLOBAL SOLUTIONS FOR 2D COUPLED BURGERS-COMPLEX-GINZBURG-LANDAU EQUATIONS
Electronic Journal of Differential Equations, Vol. 015 015), No. 99, pp. 1 14. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu ftp eje.math.txstate.eu GLOBAL SOLUTIONS FOR D COUPLED
More informationTMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments
Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary
More informationProblems Governed by PDE. Shlomo Ta'asan. Carnegie Mellon University. and. Abstract
Pseuo-Time Methos for Constraine Optimization Problems Governe by PDE Shlomo Ta'asan Carnegie Mellon University an Institute for Computer Applications in Science an Engineering Abstract In this paper we
More information19 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 informationREVERSIBLE LONG-TERM INTEGRATION WITH VARIABLE STEPSIZES. 1. Introduction. We consider the numerical treatment of systems
SIAM J. SCI. COMPUT. c 1997 Society for Industrial and Applied Mathematics Vol. 18, No. 1, pp. 257 269, January 1997 014 REVERSIBLE LONG-TERM INTEGRATION WITH VARIABLE STEPSIZES ERNST HAIRER AND DANIEL
More informationA Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential
Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix
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 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 information23 Implicit differentiation
23 Implicit ifferentiation 23.1 Statement The equation y = x 2 + 3x + 1 expresses a relationship between the quantities x an y. If a value of x is given, then a corresponing value of y is etermine. For
More informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
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 informationThe Sokhotski-Plemelj Formula
hysics 25 Winter 208 The Sokhotski-lemelj Formula. The Sokhotski-lemelj formula The Sokhotski-lemelj formula is a relation between the following generalize functions (also calle istributions), ±iǫ = iπ(),
More informationDAMTP 000/NA04 On the semi-norm of raial basis function interpolants H.-M. Gutmann Abstract: Raial basis function interpolation has attracte a lot of
UNIVERSITY OF CAMBRIDGE Numerical Analysis Reports On the semi-norm of raial basis function interpolants H.-M. Gutmann DAMTP 000/NA04 May, 000 Department of Applie Mathematics an Theoretical Physics Silver
More information11.7. Implicit Differentiation. Introduction. Prerequisites. Learning Outcomes
Implicit Differentiation 11.7 Introuction This Section introuces implicit ifferentiation which is use to ifferentiate functions expresse in implicit form (where the variables are foun together). Examples
More informationH. MOUNIER, J. RUDOLPH, M. FLIESS, AND P. ROUCHON. interior mass is considered. It is viewed as a linear delay system. A
ESAIM: Control, Optimisation an Calculus of Variations URL: http://www.emath.fr/cocv/ September 998, Vol. 3, 35{3 TRACKING CONTROL OF A VIBRATING STRING WITH AN INTERIOR MASS VIEWED AS DELAY SYSTEM H.
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 informationKramers Relation. Douglas H. Laurence. Department of Physical Sciences, Broward College, Davie, FL 33314
Kramers Relation Douglas H. Laurence Department of Physical Sciences, Browar College, Davie, FL 333 Introuction Kramers relation, name after the Dutch physicist Hans Kramers, is a relationship between
More information1 Heisenberg Representation
1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.
More informationSemianalytical method of lines for solving elliptic partial dierential equations
Chemical Engineering Science 59 (2004) 781 788 wwwelseviercom/locate/ces Semianalytical metho of lines for solving elliptic partial ierential equations Venkat R Subramanian a;, Ralph E White b a Department
More informationTHE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE
Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationthe electronic journal of combinatorics 4 (997), #R22 2 x. Introuction an the main results Asymptotic calculations are applie to stuy the egrees of ce
Asymptotics of oung Diagrams an Hook Numbers Amitai Regev Department of Theoretical Mathematics The Weizmann Institute of Science Rehovot 7600, Israel an Department of Mathematics The Pennsylvania State
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 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 informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationand from it produce the action integral whose variation we set to zero:
Lagrange Multipliers Monay, 6 September 01 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine
More informationThe Sokhotski-Plemelj Formula
hysics 24 Winter 207 The Sokhotski-lemelj Formula. The Sokhotski-lemelj formula The Sokhotski-lemelj formula is a relation between the following generalize functions (also calle istributions), ±iǫ = iπ(),
More informationDifferentiation ( , 9.5)
Chapter 2 Differentiation (8.1 8.3, 9.5) 2.1 Rate of Change (8.2.1 5) Recall that the equation of a straight line can be written as y = mx + c, where m is the slope or graient of the line, an c is the
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 informationPhysics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2
Physics 505 Electricity an Magnetism Fall 003 Prof. G. Raithel Problem Set 3 Problem.7 5 Points a): Green s function: Using cartesian coorinates x = (x, y, z), it is G(x, x ) = 1 (x x ) + (y y ) + (z z
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 informationAsymptotics of a Small Liquid Drop on a Cone and Plate Rheometer
Asymptotics of a Small Liqui Drop on a Cone an Plate Rheometer Vincent Cregan, Stephen B.G. O Brien, an Sean McKee Abstract A cone an a plate rheometer is a laboratory apparatus use to measure the viscosity
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson
JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises
More informationALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS
ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS MARK SCHACHNER Abstract. When consiere as an algebraic space, the set of arithmetic functions equippe with the operations of pointwise aition an
More informationSOLVING ODE s NUMERICALLY WHILE PRESERVING ALL FIRST INTEGRALS
SOLVING ODE s NUMERICALLY WHILE PRESERVING ALL FIRST INTEGRALS G.R.W. QUISPEL 1,2 and H.W. CAPEL 3 Abstract. Using Discrete Gradient Methods (Quispel & Turner, J. Phys. A29 (1996) L341-L349) we construct
More informationSome vector algebra and the generalized chain rule Ross Bannister Data Assimilation Research Centre, University of Reading, UK Last updated 10/06/10
Some vector algebra an the generalize chain rule Ross Bannister Data Assimilation Research Centre University of Reaing UK Last upate 10/06/10 1. Introuction an notation As we shall see in these notes the
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 informationExamining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing
Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June
More informationApplications of the Wronskian to ordinary linear differential equations
Physics 116C Fall 2011 Applications of the Wronskian to orinary linear ifferential equations Consier a of n continuous functions y i (x) [i = 1,2,3,...,n], each of which is ifferentiable at least n times.
More informationIERCU. Institute of Economic Research, Chuo University 50th Anniversary Special Issues. Discussion Paper No.210
IERCU Institute of Economic Research, Chuo University 50th Anniversary Special Issues Discussion Paper No.210 Discrete an Continuous Dynamics in Nonlinear Monopolies Akio Matsumoto Chuo University Ferenc
More information. ISSN (print), (online) International Journal of Nonlinear Science Vol.6(2008) No.3,pp
. ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.6(8) No.3,pp.195-1 A Bouneness Criterion for Fourth Orer Nonlinear Orinary Differential Equations with Delay
More informationQubit channels that achieve capacity with two states
Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March
More informationOne Dimensional Convection: Interpolation Models for CFD
One Dimensional Convection: Interpolation Moels for CFD ME 448/548 Notes Geral Recktenwal Portlan State University Department of Mechanical Engineering gerry@p.eu ME 448/548: D Convection-Di usion Equation
More informationThe Ehrenfest Theorems
The Ehrenfest Theorems Robert Gilmore Classical Preliminaries A classical system with n egrees of freeom is escribe by n secon orer orinary ifferential equations on the configuration space (n inepenent
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 informationWeb-Based Technical Appendix: Multi-Product Firms and Trade Liberalization
Web-Base Technical Appeni: Multi-Prouct Firms an Trae Liberalization Anrew B. Bernar Tuck School of Business at Dartmouth & NBER Stephen J. Reing LSE, Yale School of Management & CEPR Peter K. Schott Yale
More informationCalculus of Variations
Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,
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 informationParallel-in-time integrators for Hamiltonian systems
Parallel-in-time integrators for Hamiltonian systems Claude Le Bris ENPC and INRIA Visiting Professor, The University of Chicago joint work with X. Dai (Paris 6 and Chinese Academy of Sciences), F. Legoll
More informationAN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 14. Khyruddin Akbar Ansari, Ph.D., P.E.
AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD Mathca Release 14 Khyruin Akbar Ansari, Ph.D., P.E. Professor of Mechanical Engineering School of Engineering an Applie Science Gonzaga University SDC
More informationI. INTRODUCTION In this work we analize the slow roll approximation in cosmological moels of nonminimally couple theories of gravity (NMC) starting fr
Slow Roll Ination in Non-Minimally Couple Theories: Hyperextene Gravity Approach Diego F. Torres Departamento e Fsica, Universia Nacional e La Plata C.C. 67, 900, La Plata, Buenos Aires, Argentina Abstract
More informationPhysics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1
Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems
More informationGlobal Optimization for Algebraic Geometry Computing Runge Kutta Methods
Global Optimization for Algebraic Geometry Computing Runge Kutta Methos Ivan Martino 1 an Giuseppe Nicosia 2 1 Department of Mathematics, Stockholm University, Sween martino@math.su.se 2 Department of
More informationSparse Reconstruction of Systems of Ordinary Differential Equations
Sparse Reconstruction of Systems of Orinary Differential Equations Manuel Mai a, Mark D. Shattuck b,c, Corey S. O Hern c,a,,e, a Department of Physics, Yale University, New Haven, Connecticut 06520, USA
More informationMath 1B, lecture 8: Integration by parts
Math B, lecture 8: Integration by parts Nathan Pflueger 23 September 2 Introuction Integration by parts, similarly to integration by substitution, reverses a well-known technique of ifferentiation an explores
More informationu t v t v t c a u t b a v t u t v t b a
Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying
More informationepub WU Institutional Repository
epub WU Institutional Repository Michael A. Hauser Maximum Likelihoo Estimators for ARMA an ARFIMA Moels. A Monte Carlo Stuy. Paper Original Citation: Hauser, Michael A. (998) Maximum Likelihoo Estimators
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 informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
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 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 informationCentrum voor Wiskunde en Informatica
Centrum voor Wiskune en Informatica Moelling, Analysis an Simulation Moelling, Analysis an Simulation Conservation properties of smoothe particle hyroynamics applie to the shallow water equations J.E.
More informationA global Implicit Function Theorem without initial point and its applications to control of non-affine systems of high dimensions
J. Math. Anal. Appl. 313 (2006) 251 261 www.elsevier.com/locate/jmaa A global Implicit Function Theorem without initial point an its applications to control of non-affine systems of high imensions Weinian
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 informationMark J. Machina CARDINAL PROPERTIES OF "LOCAL UTILITY FUNCTIONS"
Mark J. Machina CARDINAL PROPERTIES OF "LOCAL UTILITY FUNCTIONS" This paper outlines the carinal properties of "local utility functions" of the type use by Allen [1985], Chew [1983], Chew an MacCrimmon
More informationConsider for simplicity a 3rd-order IIR filter with a transfer function. where
Basic IIR Digital Filter The causal IIR igital filters we are concerne with in this course are characterie by a real rational transfer function of or, equivalently by a constant coefficient ifference equation
More informationA Spectral Method for the Biharmonic Equation
A Spectral Metho for the Biharmonic Equation Kenall Atkinson, Davi Chien, an Olaf Hansen Abstract Let Ω be an open, simply connecte, an boune region in Ê,, with a smooth bounary Ω that is homeomorphic
More informationMarkov Chains in Continuous Time
Chapter 23 Markov Chains in Continuous Time Previously we looke at Markov chains, where the transitions betweenstatesoccurreatspecifietime- steps. That it, we mae time (a continuous variable) avance in
More informationBasic IIR Digital Filter Structures
Basic IIR Digital Filter Structures The causal IIR igital filters we are concerne with in this course are characterie by a real rational transfer function of or, equivalently by a constant coefficient
More informationOPTIMAL CONTROL OF A PRODUCTION SYSTEM WITH INVENTORY-LEVEL-DEPENDENT DEMAND
Applie Mathematics E-Notes, 5(005), 36-43 c ISSN 1607-510 Available free at mirror sites of http://www.math.nthu.eu.tw/ amen/ OPTIMAL CONTROL OF A PRODUCTION SYSTEM WITH INVENTORY-LEVEL-DEPENDENT DEMAND
More informationarxiv: v1 [math-ph] 5 May 2014
DIFFERENTIAL-ALGEBRAIC SOLUTIONS OF THE HEAT EQUATION VICTOR M. BUCHSTABER, ELENA YU. NETAY arxiv:1405.0926v1 [math-ph] 5 May 2014 Abstract. In this work we introuce the notion of ifferential-algebraic
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 information