Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing
|
|
- Christopher Moody
- 5 years ago
- Views:
Transcription
1 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 1, 005 Abstract The precision of numerical integration greatly affects the ability to track an navigate a spacecraft throughout the solar system. In fact, precision of up to fifteen significant igits is require for accurate navigation of spacecraft missions heaing towar the outer planets. This paper explores the use of numerical integration techniques base symplectic an variational integrators to propagate orbit trajectories, plus an examination is mae of how to incorporate into these methos the non-conservative forces affecting the trajectories. These numerical techniques are of interest ue to their intrinsic ability to conserve quantites such as energy, momentum, the Hamiltonian, etc. from the inherent ynamics. Research Objective The ability to successfully navigate the solar system is irectly relate to how well a spacecraft state can be preicte. State etermination is accomplishe through both estimation from Earth-base observations, as well as through numeric propagation of governing equations of motion. Traitional techniques for numerical integration are RK (Runge-Kutta) an similar avance iscretization techniques that typically require small step sizes to ensure precision in the trajectory propagation of large solar istances. A benefit of these techniques is in their ability to incorporate non-conservative forces, along with starts an stops of the propagation in orer to inclue ifferent perturbations. The objective of this research is to implement symplectic an variational integrators for propagating orbit trajectories an to etermine the feasibility an methoology for incorporating non-conservative forces an perturbations into the evelope techniques. The approach to this research will be a step-wise increase in complexity of the orbit problem. The first step will be to formulate the Hamiltonian for a simple two-boy problem (involving a spacecraft in the gravity well of a single planet) with no external forcing other than the gravitational potential of the host planet (See simple schematic in figure 1). In aition, the integration scheme for a basic symplectic integrator will be evelope an formulate for this initial problem; the chosen symplectic integrator will be base on the mipoint rule. This first step will facilitate comparisons with classic RK methos (such as variable an fixe time step methos) since the two-boy problem has analytic solutions. Further, this step will enable an assessment of step size on numerical accuracy. The secon stage of this research will emonstrate that the chosen mipoint rule is not only a symplectic integrator but also a variational integrator. Through this extension, which comes 1
2 through the iscrete Euler-Lagrange equations (DEL) an the iscrete Legenre transform, a possible metho of incorporating non-conservative forces becomes apparant. The thir stage of the research will consier how non-conservative forcings enter variational integration schemes through the use of the Alembert Principle. Either the DEL-like approach just mentione will be consiere, or perhaps a return to the Hamiltonian-sie of mechanics will be necessary. Figure 1: Representation of Two-Boy Problem The secon step of the research will introuce impulsive forces at preset times in orer to etermine how symplectic versus RK propagation hanles the change in trajectory irection. Again, the analytic solutions for two-boy orbit mechanics are known, so meaningful comparisons can be mae. To implement the impulsive forces, both integrators will be stoppe, an a step change in velocity will be mae to current velocity values (thus imparting an impulsive force). The impulsive force is essentially bumping the energy (Hamiltonian) of the orbit to a new, constant level. The thir step of the research will be to etermine how to incorporate non-conservative forces into the equations of motion such that symplectic integration remains vali. The technique for oing this is not yet known, however, there are known classes of systems that conserve a Hamiltonian but not energy. Further, the notion of being locally Hamiltonian may play an important role in this implementation. This all remains to be etermine through the research. Equations of Motion The ynamics of an orbiting spacecraft uner the influence of only gravity obey a conservation of energy principal from which the equations of motion can be erive [1]. The energy (an Hamiltonian) for a spacecraft orbiting a planet is given by H(q, p) = 1 m p mµ q where q = q is the raius from the planet c.m. (center of mass), m is the spacecraft mass, µ is the gravitational parameter for a particular planets or boy being orbite, an p is the spacecraft momentum. The spacecraft mass is assume insignificant to the planet mass, so the spacecraft exerts no influence on the planet s orbit. Thus, the center of gravity of the two-boy system remains fixe at the planet c.m. The Hamiltonian in equation 1 governs the resulting equations of motion: q t = H p = 1 m p (1) ż = J H(z) () p t = H q = mµ q q 3 [ ] where z = (q, p) T 0 I an J =. Note also the secon-orer equation of motion that I 0 arises from the left han sie of equation : m q = mµ q q. 3
3 For spacecraft orbiting multi-boy systems (i.e.. systems with one spacecraft an multiple planets), a conservation of energy principal exists as well. However, the presence of nonconservative forces causes energy not to be conserve. Symplectic Integration For a n-imensional symplectic vector space (Z, Ω), let z = (q, p) Z enote canonical coorinates such that symplectic matrix J represents Ω. The evolution equations of a Hamiltonian system on phase space Z are given by [ where J = Proposition 1. 0 I I 0 ]. ż = J H(z), with z t=0 = z 0 The mapping φ : z n Z z n+1 Z base on the mipoint rule is symplectic, such that ż = z n+1 z n z = z n+1 + z n ( ) zn+1 + z n z n+1 = z n + J H (3) Proof. For this mapping to be symplectic, the Jacobian of the upate scheme must be symplectic: z n+1 = I + J ( ) zn+1 + z n H z n z n = I + 1 ( ) ( zn+1 + z n J H I + z ) n+1 (4) z n From equation 4, efine the following: where H is evaluate at A = 1 J H (5) ( ) zn+1+z n. Then, equation 4 can be re-written as (I A) z n+1 z n = I + A z n+1 z n = (I A) 1 (I + A) (6) where an assumption has been mae on the invertibility of I A. In other wors, the Hessian of H an the size of are such that I A is invertible. Assuming provie I A is invertible, the right-han sie of equation 6 is the Cayley transform, which is symplectic if an only if A is Hamiltonian [, p. 80]. Since A is a linear mapping, then A is Hamiltonian if an only if A is Ω-skew [, Prop..5.1, p. 77]. An, Ω- skewness of A is equivalent to the matrix JA being symmetric (i.e.. JA + A T J = 0) [, p. 78]. Thus, evaluate A from equation 5 as follows: 3
4 ( ) ( ) T 1 1 JA + A T J = J J H + J H J = 1 (JJ H + ( H ) ) T J T J = 1 ( I H + H I ) = 0 since H is symmetric. Thus, A is Ω-skew, implying A is Hamiltonian. Finally, by the Cayley transform, equation 6 is symplectic, inicating equation 3 efines a symplectic mapping. Implementing Symplectic Integration The ifficultly with implementing the symplectic integration scheme of equation 3 is that the metho is implicit. An initial metho for ealing with this problem is to utilize a Newton solve for fining the roots of the equation. Equation 3 can be viewe as ( ) zn+1 + z n f(z n+1 ) = z n+1 z n J H = 0 (7) where the roots of function f are the solutions z n+1 to equation 3. The Newton solve is initiate with a guess, such as guess = z n, such that f(z n+1 ) = 0 = f(guess + δ) f(guess) + f (guess)δ (8) with f (guess) being the Jacobian of f, as given in equation 6. The value of δ is etermine from equation 8: δ = (f (guess)) 1 f(guess) which is use to iterate on the initial guess (guess = guess + δ) until f(guess) < ɛ for some esire, small value of ɛ. Variational Algorithms an Variational Integration [3] The action integral on a Lagrangian L(q, q) is efine as S = b a L(q(t), q(t))t (9) where the Euler-Lagrange equations extremize S given fixe en points. The evelopment of a iscrete version of the Euler-Lagrange equations involves iscretizing the action integral S over a time step such that L (q k, q k+1 ) = tk+1 t k L(q(t), q(t))t L( q, q) (10) where L( q, q) is constant over the interval t (t k, t k+1 ) with q an q a iscrete approximation for q an q an t k+1 t k =. This efines a iscrete Lagrangian such that L : Q Q R. The associate iscrete analog to the action integral is then given by S = N 1 k=0 L (q k, q k+1 ), (11) 4
5 where S : Q N+1 R, an similar to the continuous case, extremizing S prouces the iscrete version of the Euler-Lagrange equations or DEL (Discrete Euler-Lagrange) equations: D 1 L (q k+1, q k+ ) + D L (q k, q k+1 ) = 0 (1) for k = 1,..., N 1. Note, equation 1 efines a variational algorithm. The Legenre transform has an analogous iscrete version, as efine in [3], FL : Q Q T Q; (q k, q k+1 ) (q k+1, D L (q k, q k+1 )) (13) which further implies p k+1 = D L (q k, q k+1 ). For the research conucte herein, the iscrete approximation to the continuous Lagrangian is chosen to be the symmetric form L sym,α (q k, q k+1 ) = L + L ( (1 α)q k + αq k+1, q ) k+1 q k ( αq k + (1 α)q k+1, q ) k+1 q k Theorem 1. [3] The mipoint rule in equation 3 is a variational algorithm through the iscrete Legenre tranformation with the iscrete Lagrangian L sym,α for α = 1. Proof. (in the reverse irection) Let q m = q k+1+q k an q = q k+1 q k, then for α = 1, ( L sym, 1 qk + q k+1 (q k, q k+1 ) = L, q ) k+1 q k = L (q m, q m ) (15) an (14) D 1 L sym, 1 (q k, q k+1 ) = D 1L (q m, q m ) D L (q m, q m ) (16) D L sym, 1 (q k, q k+1 ) = D 1L (q m, q m ) + D L (q m, q m ) (17) Through the DLT (iscrete Legenre transform) in equation 13 an the DEL equations (equation 1), p k DLT DEL = D L (q k 1, q k ) = D 1 L (q k, q k+1 ) (18) p k+1 = D L (q k, q k+1 ) (19) which efines an implicit algorithm such that (q k, p k ) (q k+1, p k+1 ). For L sym, 1, the variables p k an p k+1 are equivalent to equations 16 an 17, respectively. Aing an subtracting equations 18 an 19, plus incorporating the results of equations 16 an 17 for α = 1 provies the following algorithms: p k+1 p k p k+1 + p k Restrict the class of Lagrangian to those of canonical form = D l L (q m, q m ) (0) = D L (q m, q m ) (1) L(q, q) = 1 qt M q V (q) () 5
6 where M is a symmetric matrix, an V is a potential function. Note, this loses the generality in the proof, but this class of Lagrangian is vali for most mechanical systems, incluing those in the research within this paper. The continuous Legenre transform provies the following Hamiltonian for this class of canonical system: H(q, p) = p T q L(q, q) = 1 M 1 p T p + V (q) (3) with p = M q For the Lagrangian of equation, the first an secon slot erivatives are D 1 L(q, q) = V (q) an D L(q, q) = M q (4) which when combine with the algorithms in equations 0 an 1 result in the following equations p k+1 p k = V (q m ) = V ( q k+1 + q k ) p k+1 p k = H ( ) qk+1 + q k (5) q p k+1 + p k = M q m = M q k+1 q k q k+1 q k = H ( ) pk+1 + p k (6) q which are the exactly the form of equation 3 for z = (q, p) an the Hamiltonian of equation 3. Thus, the mipoint rule efines a variational algorithm with L sym,α when α = 1. Incorporating Non-Conservative Forces The results of Theorem 1 are quite powerful an provie a framework from which non-conservative forces might be incorporate into variational-like algorithms. No results will be shown from this brief section but it serves as a basis for how future work will continue. As iscusse in [3], the integral Lagrange- Alembert principal is given by δ L(q(t), q(t))t + F (q(t), q(t))t = 0 (7) which appears similar to the action integral in equation 9. In fact, a iscrete version of equation 7 exists from which force DEL equations can be erive: D 1 L (q k+1, q k+ ) + D L (q k, q k+1 ) + F (q k+1, q k+ ) + F + (q k, q k+1 ) = 0 (8) From these force DEL equations, an implicit algorithm is likely erivable in a manner analogous to that from the proof of Theorem 1. A key to the evelopment of this algorithm will require appropriate iscretization of the non-conservative forces, likely in a symmetric fashion as with L sym,α. Application of Symplectic/Variational Integration to Circular Orbits The first-orer symplectic integration scheme from equation 3 was implemente in Matlab for the Hamiltonian ynamics of equation 1 for a spacecraft orbiting a planet (two orbit types are analyze, one circular an one elliptic). The parameters chosen for the spacecraft an planet are given in table 1. Consequently, the initial conitions also set the value of the conserve Hamiltonian from equation 1 for the two orbit types; the Hamiltonian values for each orbit type are inclue in the table (they are negative, which correspons to circular an elliptic orbits). 6
7 Table 1: Simulation Parameter Data an Initial Conitions µ (km 3 /s ) m (kg) 1 q 0 (km) (0, , 0) T p 0 (circular) (km/s) ( , 0, 0) T H(q, p)(circular) (kg km /s ) p 0 (elliptic) (km/s) ( , 0, 0) T H(q, p)(elliptic) (kg km /s ) The Matlab coe also implements three aitional schemes for comparing numerical performance over roughly one ay (86400 secons) of orbit ata. The comparison schemes were the Matlab intrinsic RK45 scheme ODE45 (which has a variable time step) an two Mathworks fixe-step RK schemes calle ODE an ODE4 (RK an RK4, respectively), available from the company web site. Two orbit types were simulate with the numerical integrators: a circular orbit an an elliptical orbit. With both simulations, which are two-boy conic orbits, analytic solutions for the orbits are available, which facilitates comparison of numerical schemes an etermination of their accuracy. Circular Orbit Simulation Figure juxtaposes results from the symplectic, RK45, an RK4 integration schemes, along with an analytic solution for the two-boy circular orbit. Due to the choice of initial conitions, the orbit is planar. Figure : Comparison of Numerical Techniques Versus the Analytic Solution From figure, the variable-step RK45 scheme is seen to win towar the planet; the metho is essentially issipating energy from the orbit trajectory. The symplectic an RK4 (fixe-step) schemes track the exact orbit more closely an cannot be separate in figure. To istinguish these two schemes, figure 3 zooms in on a segment in the upper-left portion of figure. 7
8 Figure 3: Close-Up of Circular Orbit from Figure The close up of figure 3 reveals that both the RK4 an symplectic integrators track the analytic solution well, though both with an offset. Interestingly the symplectic integrator has a smaller offset than the RK4 metho, which is interesting given that RK4 is fourth orer an this particular symplectic scheme is secon orer. Utilizing an integrator that intrinsically accounts for the system geometry (conserving the Hamiltonian) oes have benefit over a more naive technique, at least with a simple circular orbit. The conservation of energy by the symplectic integrator is reveale in figure 4, which also clearly shows how RK4 introuces false issipation, something that woul cause a loss in accuracy over long-term integration. The issipation is in the 10 th igit only uring this time interval, so the plot scale is cryptic. Figure 4: Energy of Circular Orbit shoul be Conserve 8
9 Elliptical Orbit Simulation Figure 5 juxtaposes results from all four integration schemes (symplectic, RK, RK4, an RK45), along with an analytic solution for the two-boy elliptic orbit. The orbit is again planar ue to the choice of initial conitions; for the elliptic. Figure 5: Comparison of Numerical Techniques Versus the Analytic Solution From figure 5, the variable-step RK45 scheme is again seen to win towar the planet, artificially issipating energy as was the case with the circular orbit. The remaining schemes track the exact orbit more closely, which can be seen in the zoom of the orbit shown in figure 6. Figure 6: Close-Up of Ellipse Orbit from Figure 5 The close up of figure 6 shows that the RK4 integrator is significantly more accurate than the symplectic or RK technique. This is slightly iffering from the circular case but seems reasonable since the RK4 technique is 4th-orer, whereas the RK an chosen symplectic techniques are both secon orer. Outsie of the circular case, these plots woul appear to make the RK4 9
10 technique the integrator of choice. However, an examination of the energy of the system brings aitional unerstaning of the ability of the ifferent schemes to prouce longterm accuracy (See figures 7 an 8). Figure 7: Elliptic Orbit Energy Showing RK Decay an Symplectic Conservation Figure 8: Close-Up of Elliptic Orbit Energy Revealing Slow Decay in RK4 The results of figures 7 an 8 inicate that the symplectic integrator, though not constant in energy, oscillates about a mean value, offset from the actual energy. A correction factor is likely realizable for the numerical implementation of the symplectic integrator that woul make it oscillate about the exact value. The RK4 technique appears to conserve energy too in figure 7; however, a closer examination in figure 8 again shows that a very long-term integration woul exhibit energy ecay with RK4. Note, the number of significant igits isplaye in the MATLAB generate plots is cryptic as it oes not reflect the issipation occurring in the 10th igit for the RK4 energy, though the tren in the plot is unmistakable. 10
11 Future Steps One immeiate question to answer is how versatile is the symplectic/variational integrator approach? The RK methos an other similar techniques can incorporate iscretization of non-conservative forces an prouce certain known accuracies. However, can the force DEL equations prouce their own variational-like algorithms that preserve some element of the force ynamics, thus proviing an even higher accuracy than traitional techniques? This question is of immeiate interest to the research in this paper. The subsequent focus of this research will be to incorporate non-conservative forces, such as thruster firings an solar raiation pressure. From the aition of these forces, performance an precision comparisons can be mae with existing integration techniques. To assess precision an accuracy, simple two-boy orbits with impulsive thruster firings can be propagate to etermine how well ifferent numerical integration schemes track the known, analytic solutions. For an integration scheme to be useful in the real-worl spacecraft navigation, the scheme must have high precision an have timely performance. Future work can focus on the integration scheme itself; the mipoint rule is the lowest-orer symplectic integrator, so higher-orer schemes shoul be more accurate. Accuracy requires not only reproucing proper position an velocity solutions but also arriving at the solutions on a similar time scale; for instance, the perio of a conic orbit shoul be recoverable from the integration scheme. Regaring performance, explicit formulations, or goo approximations to implicit forms coul likely improve performance, so such ieas will be explore. References [1] R.R. Bate, D.D. Mueller, an J.E. White, Funamentals of Astroynamics, Dover Publications, [] J.E. Marsen an T.S. Ratiu, Introuction to Mechanics an Symmetry, Springer, n eition, 00. [3] C. Kane, J.E. Marsen, M. Ortiz, an M. West, Variational Integrators an the Newmark Algorithm for Conservative an Dissipative Mechanical Systems, Int. J. Num. Math. Eng., 49, pp ,
Lecture 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 informationDiscrete Hamilton Jacobi Theory and Discrete Optimal Control
49th IEEE Conference on Decision an Control December 15-17, 2010 Hilton Atlanta Hotel, Atlanta, GA, USA Discrete Hamilton Jacobi Theory an Discrete Optimal Control Tomoi Ohsawa, Anthony M. Bloch, an Melvin
More informationinflow outflow Part I. Regular tasks for MAE598/494 Task 1
MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the
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 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 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 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 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 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 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 informationA Second Time Dimension, Hidden in Plain Sight
A Secon Time Dimension, Hien in Plain Sight Brett A Collins. In this paper I postulate the existence of a secon time imension, making five imensions, three space imensions an two time imensions. I will
More informationVectors in two dimensions
Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication
More informationOptimization of a point-mass walking model using direct collocation and sequential quadratic programming
Optimization of a point-mass walking moel using irect collocation an sequential quaratic programming Chris Dembia June 5, 5 Telescoping actuator y Stance leg Point-mass boy m (x,y) Swing leg x Leg uring
More informationQuantum mechanical approaches to the virial
Quantum mechanical approaches to the virial S.LeBohec Department of Physics an Astronomy, University of Utah, Salt Lae City, UT 84112, USA Date: June 30 th 2015 In this note, we approach the virial from
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 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 informationApplications of First Order Equations
Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more
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 informationThe Principle of Least Action and Designing Fiber Optics
University of Southampton Department of Physics & Astronomy Year 2 Theory Labs The Principle of Least Action an Designing Fiber Optics 1 Purpose of this Moule We will be intereste in esigning fiber optic
More informationSurvey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013
Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing
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 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 informationVIRTUAL STRUCTURE BASED SPACECRAFT FORMATION CONTROL WITH FORMATION FEEDBACK
AIAA Guiance, Navigation, an Control Conference an Exhibit 5-8 August, Monterey, California AIAA -9 VIRTUAL STRUCTURE BASED SPACECRAT ORMATION CONTROL WITH ORMATION EEDBACK Wei Ren Ranal W. Bear Department
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 informationCalculus Class Notes for the Combined Calculus and Physics Course Semester I
Calculus Class Notes for the Combine Calculus an Physics Course Semester I Kelly Black December 14, 2001 Support provie by the National Science Founation - NSF-DUE-9752485 1 Section 0 2 Contents 1 Average
More informationPDE Notes, Lecture #11
PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =
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 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 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 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 informationDissipative numerical methods for the Hunter-Saxton equation
Dissipative numerical methos for the Hunter-Saton equation Yan Xu an Chi-Wang Shu Abstract In this paper, we present further evelopment of the local iscontinuous Galerkin (LDG) metho esigne in [] an a
More informationAngles-Only Orbit Determination Copyright 2006 Michel Santos Page 1
Angles-Only Orbit Determination Copyright 6 Michel Santos Page 1 Abstract This ocument presents a re-erivation of the Gauss an Laplace Angles-Only Methos for Initial Orbit Determination. It keeps close
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 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 informationNumerical Integrator. Graphics
1 Introuction CS229 Dynamics Hanout The question of the week is how owe write a ynamic simulator for particles, rigi boies, or an articulate character such as a human figure?" In their SIGGRPH course notes,
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 informationFree rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012
Free rotation of a rigi boy 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 1 Introuction In this section, we escribe the motion of a rigi boy that is free to rotate
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More informationDiophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations
Diophantine Approximations: Examining the Farey Process an its Metho on Proucing Best Approximations Kelly Bowen Introuction When a person hears the phrase irrational number, one oes not think of anything
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 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 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 informationVariational Integrators
Variational Integrators Thesis by Matthew West In Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy California Institute of Technology Pasaena, California 2004 (Submitte May
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
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 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 informationG4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const.
G4003 Avance Mechanics 1 The Noether theorem We alreay saw that if q is a cyclic variable, the associate conjugate momentum is conserve, q = 0 p q = const. (1) This is the simplest incarnation of Noether
More informationAN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 13. Khyruddin Akbar Ansari, Ph.D., P.E.
AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD Mathca Release 13 Khyruin Akbar Ansari, Ph.D., P.E. Professor of Mechanical Engineering School of Engineering Gonzaga University SDC PUBLICATIONS Schroff
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 information1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity
AP Physics Multiple Choice Practice Electrostatics 1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity. A soli conucting sphere is given a positive charge Q.
More informationensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y
Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay
More information6 Wave equation in spherical polar coordinates
6 Wave equation in spherical polar coorinates We now look at solving problems involving the Laplacian in spherical polar coorinates. The angular epenence of the solutions will be escribe by spherical harmonics.
More informationTime-of-Arrival Estimation in Non-Line-Of-Sight Environments
2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor
More informationOptimal Variable-Structure Control Tracking of Spacecraft Maneuvers
Optimal Variable-Structure Control racking of Spacecraft Maneuvers John L. Crassiis 1 Srinivas R. Vaali F. Lanis Markley 3 Introuction In recent years, much effort has been evote to the close-loop esign
More informationAN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A
AN INTRODUCTION TO AIRCRAFT WIN FLUTTER Revision A By Tom Irvine Email: tomirvine@aol.com January 8, 000 Introuction Certain aircraft wings have experience violent oscillations uring high spee flight.
More informationChapter 2 Governing Equations
Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement
More informationLinear First-Order Equations
5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)
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 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 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 information4. Important theorems in quantum mechanics
TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this
More informationSingle Arm, Centrifugal, Water Turbine for Low Head and Low Flow Application: Part 1- Theory and Design
Energy an Power 2018, 8(2): 51-55 DOI: 10.5923/j.ep.20180802.03 Single Arm, Centrifugal, Water Turbine for Low ea an Low Flow Application: Part 1- Theory an Design Kiplangat C. Kononen 1, Augustine B.
More informationLINEARIZATION METHODS FOR VARIATIONAL INTEGRATORS AND EULER-LAGRANGE EQUATIONS
LINEARIZATION METHODS FOR VARIATIONAL INTEGRATORS AND EULER-LAGRANGE EQUATIONS Tesis by AnilN.Hirani Avisors James R. Arvo, Jerrol E. Marsen In Partial Fulfillment of te Requirements for te Degree of Master
More informationG j dq i + G j. q i. = a jt. and
Lagrange Multipliers Wenesay, 8 September 011 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 informationThermal runaway during blocking
Thermal runaway uring blocking CES_stable CES ICES_stable ICES k 6.5 ma 13 6. 12 5.5 11 5. 1 4.5 9 4. 8 3.5 7 3. 6 2.5 5 2. 4 1.5 3 1. 2.5 1. 6 12 18 24 3 36 s Thermal runaway uring blocking Application
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 informationThe Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis
The Kepler Problem For the Newtonian 1/r force law, a miracle occurs all of the solutions are perioic instea of just quasiperioic. To put it another way, the two-imensional tori are further ecompose into
More informationON THE MEANING OF LORENTZ COVARIANCE
Founations of Physics Letters 17 (2004) pp. 479 496. ON THE MEANING OF LORENTZ COVARIANCE László E. Szabó Theoretical Physics Research Group of the Hungarian Acaemy of Sciences Department of History an
More informationarxiv:physics/ v4 [physics.class-ph] 9 Jul 1999
AIAA-99-2144 PROPULSION THROUGH ELECTROMAGNETIC SELF-SUSTAINED ACCELERATION arxiv:physics/9906059v4 [physics.class-ph] 9 Jul 1999 Abstract As is known the repulsion of the volume elements of an uniformly
More informationThe Press-Schechter mass function
The Press-Schechter mass function To state the obvious: It is important to relate our theories to what we can observe. We have looke at linear perturbation theory, an we have consiere a simple moel for
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 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 informationConstruction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems
Construction of the Electronic Raial Wave Functions an Probability Distributions of Hyrogen-like Systems Thomas S. Kuntzleman, Department of Chemistry Spring Arbor University, Spring Arbor MI 498 tkuntzle@arbor.eu
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 informationSYNCHRONOUS SEQUENTIAL CIRCUITS
CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents
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 informationThe Hamiltonian particle-mesh method for the spherical shallow water equations
ATMOSPHERIC SCIENCE LETTERS Atmos. Sci. Let. 5: 89 95 (004) Publishe online in Wiley InterScience (www.interscience.wiley.com). DOI: 10.100/asl.70 The Hamiltonian particle-mesh metho for the spherical
More informationPhysics 2212 GJ Quiz #4 Solutions Fall 2015
Physics 2212 GJ Quiz #4 Solutions Fall 215 I. (17 points) The magnetic fiel at point P ue to a current through the wire is 5. µt into the page. The curve portion of the wire is a semicircle of raius 2.
More informationCode_Aster. Detection of the singularities and computation of a card of size of elements
Titre : Détection es singularités et calcul une carte [...] Date : 0/0/0 Page : /6 Responsable : Josselin DLMAS Clé : R4.0.04 Révision : 9755 Detection of the singularities an computation of a car of size
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 informationTorque OBJECTIVE INTRODUCTION APPARATUS THEORY
Torque OBJECTIVE To verify the rotational an translational conitions for equilibrium. To etermine the center of ravity of a rii boy (meter stick). To apply the torque concept to the etermination of an
More information5.4 Fundamental Theorem of Calculus Calculus. Do you remember the Fundamental Theorem of Algebra? Just thought I'd ask
5.4 FUNDAMENTAL THEOREM OF CALCULUS Do you remember the Funamental Theorem of Algebra? Just thought I' ask The Funamental Theorem of Calculus has two parts. These two parts tie together the concept of
More informationExperiment 2, Physics 2BL
Experiment 2, Physics 2BL Deuction of Mass Distributions. Last Upate: 2009-05-03 Preparation Before this experiment, we recommen you review or familiarize yourself with the following: Chapters 4-6 in Taylor
More informationImproved Geoid Model for Syria Using the Local Gravimetric and GPS Measurements 1
Improve Geoi Moel for Syria Using the Local Gravimetric an GPS Measurements 1 Ria Al-Masri 2 Abstract "The objective of this paper is to iscuss recent research one by the author to evelop a new improve
More informationarxiv: v1 [math.na] 8 Feb 2017
Variational iscretization of te nonequilibrium termoynamics of simple systems arxiv:17.594v1 [mat.na] 8 Feb 17 François Gay-Balmaz Hiroaki Yosimura CNRS - LMD - IPSL Scool of Science an Engineering Ecole
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 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 informationSources and Sinks of Available Potential Energy in a Moist Atmosphere. Olivier Pauluis 1. Courant Institute of Mathematical Sciences
Sources an Sinks of Available Potential Energy in a Moist Atmosphere Olivier Pauluis 1 Courant Institute of Mathematical Sciences New York University Submitte to the Journal of the Atmospheric Sciences
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 information221A Lecture Notes Notes on Classica Mechanics I
1A Lecture Notes Notes on Classica Mechanics I 1 Precursor: Fermat s Principle in Geometric Optics In geometric optics, you talk about how light rays go. In homogeneous meiums, the light rays go straight.
More informationMATH , 06 Differential Equations Section 03: MWF 1:00pm-1:50pm McLaury 306 Section 06: MWF 3:00pm-3:50pm EEP 208
MATH 321-03, 06 Differential Equations Section 03: MWF 1:00pm-1:50pm McLaury 306 Section 06: MWF 3:00pm-3:50pm EEP 208 Instructor: Brent Deschamp Email: brent.eschamp@ssmt.eu Office: McLaury 316B Phone:
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 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 informationCode_Aster. Detection of the singularities and calculation of a map of size of elements
Titre : Détection es singularités et calcul une carte [...] Date : 0/0/0 Page : /6 Responsable : DLMAS Josselin Clé : R4.0.04 Révision : Detection of the singularities an calculation of a map of size of
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 informationChapter 4. Electrostatics of Macroscopic Media
Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1
More informationNON-ADIABATIC COMBUSTION WAVES FOR GENERAL LEWIS NUMBERS: WAVE SPEED AND EXTINCTION CONDITIONS
ANZIAM J. 46(2004), 1 16 NON-ADIABATIC COMBUSTION WAVES FOR GENERAL LEWIS NUMBERS: WAVE SPEED AND EXTINCTION CONDITIONS A. C. MCINTOSH 1,R.O.WEBER 2 ang.n.mercer 2 (Receive 14 August, 2002) Abstract This
More informationSolution to the exam in TFY4230 STATISTICAL PHYSICS Wednesday december 1, 2010
NTNU Page of 6 Institutt for fysikk Fakultet for fysikk, informatikk og matematikk This solution consists of 6 pages. Solution to the exam in TFY423 STATISTICAL PHYSICS Wenesay ecember, 2 Problem. Particles
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 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 information