DOCSLIB.ORG

TOPICS in CELESTIAL MECHANICS 1. the Newtonian N-Body Problerm

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
TOPICS in CELESTIAL MECHANICS 1. the Newtonian N-Body Problerm TOPICS IN CELESTIAL MECHANICS RICHARD MOECKEL 1. The Newtonian n-body Problerm Celestial mechanics can be defined as the study of the solution of Newton's differ- ential equations formulated by Isaac Newton in 1686 in his Philosophiae Naturalis Principia Mathematica. The setting for celestial mechanics is three-dimensional space: 3 R = fq = (x; y; z): x; y; z 2 Rg with the Euclidean norm: p jqj = x2 + y2 + z2: A point particle is characterized by a position q 2 R3 and a mass m 2 R+. A motion of such a particle is described by a curve q(t) where t runs over some interval in R; the mass is assumed to be constant. Some remarks will be made below about why it is reasonable to model a celestial body by a point particle. For every motion of a point particle one can define: velocity: v(t) =q _(t) momentum: p(t) = mv(t): Newton formulated the following laws of motion: Lex.I. Corpus omne perservare in statu suo quiescendi vel movendi uniformiter in directum, nisi quatenus a viribus impressis cogitur statum illum mutare 1 Lex.II. Mutationem motus proportionem esse vi motrici impressae et fieri secundem lineam qua vis illa imprimitur. 2 Lex.III Actioni contrarium semper et aequalem esse reactionem: sive corporum duorum actiones in se mutuo semper esse aequales et in partes contrarias dirigi. 3 The first law is statement of the principle of inertia. The second law asserts the existence of a force function F : R4 ! R3 such that: p_ = F (q; t) or mq¨ = F (q; t): In celestial mechanics, the dependence of F (q; t) on t is usually indirect; the force on one body depends on the positions of the other massive bodies which in turn depend on t. The third law postulates the symmetry of the mutual interaction of two bodies which will apply, in particular, to the gravitational interaction. Date: January 21, 2020. 1Every body continues in its quiescent state or moves uniformly in direction, unless it is compelled by impressed forces to change its state. 2The change of momentum is proportional to the motive force impressed and takes place along the line where this force is impressed. 3To every action there is always an equal and opposite reaction: the actions of two bodies on one another are always equal and aimed in opposite directions. 1 2 RICHARD MOECKEL The n-body problem is about the motion of n point particles under the influence of their mutual gravitational attraction. Each particle has a mass mi > 0 and 3 position, velocity and momentum vectors qi; vi; pi 2 R . The whole system can be described using the vectors q; v; p 2 R3n where q = (q1; q2; : : : ; qn) v = (v1; v2; : : : ; vn) p = (m1v1; m2v2; : : : ; mnvn): According to Newton, the gravitational force acting on particle i due to the presence of particle j is Gmimj(qj − qi) Fij = 3 jqi − qjj where G is a constant. Note that Fij acts along the line containing the masses. It's proportional to the product of the two masses and inversely proportional to the distance between them (see Figure 1). The force produced on mj by mi is Fji = −Fij by Newton's third law. By choosing the units of mass, one can arrange that G = 1 and this will be assumed from now on (see Exercise 1.1). 1.0 m2 F21 0.5 -2.0 -1.5 -1.0 -0.5 0.5 1.0 m1 F12 -0.5 Figure 1. Newtonian gravitational forces. The force on the i-th mass due to the other n − 1 masses is: X X mimj(qj − qi) Fi = Fij = 3 : jqj − qij j6=i j6=i This can be written: Fi(q) = riU(q) where X mimj (1) U(q) = jq − q j (i;j) i j i<j and ri is the partial gradient operator: @U @U @U 3 riU = ; ; 2 R : @xi @yi @zi The function U(q) will be called the Newtonian gravitational potential function. V (q) = −U(q) is the gravitational potential energy. Newton's second law becomes p_i = miq¨i = riU(q) i = 1,. ,n: or, more concisely (2)p _ = rU(q) or Mq¨ = rU(q) TOPICS IN CELESTIAL MECHANICS 3 where r is the gradient operator in R3n and M is the 3n × 3n mass matrix M = diag(m1; m1; m1; : : : ; mn; mn; mn): It is worth digressing at this point to note two important, special features of the Newtonian interparticle potential: m m i j : jqi − qjj First of all, the presence of the factor mimj has the effect that the equation for the acceleration of the i-th mass, 1 q¨i = Fi mi is independent of mi. This corresponds to the observation, notably by Galileo, that the trajectory of a falling body is independent of its mass. Figure 2. Masses falling from a tower. Second, the fact that the potential is inversely proportional to the distance be- tween the particles provides some justification for the modeling of celestial bodies by point particles. While such bodies are not even approximately pointlike, they are approximately spherically symmetric. It turns out that with the Newtonian poten- tial, spherically symmetric bodies behave as if their total mass were concentrated at their centers. To see this, consider a more general massive body, specified by giving a bounded subset B ⊂ R3 together with a continuous mass density functions ρ. The gravita- tional force exerted by such a mass distribution on a point mass m at position q is F = rU(q) where Z m U(q) = dm B jq − pj where the triple integral is over p = (x; y; z) 2 B and dm = ρ(x; y; z)dxdydz. 4 RICHARD MOECKEL 3 Proposition 1.1. Suppose B is a ball of radius R centered at q0 2 R and ρ is a spherically symmetric density function. Then the mutual Newtonian potential of the ball and a mass m at any point q with jqj > R is m m U(q) = 0 jq0 − qj where m0 is the total mass in the ball. Proof. Using the symmetry of the Euclidean distance under rotations and trans- lations, it is no loss of generality to assume q0 = (0; 0; 0) and q = (0; 0; z); z > R. Using spherical coordinates (x; y; z) = r(cos θ sin φ, sin θ sin φ, cos φ), the spherical symmetry means that the density function depends only on r and then Z R Z π Z 2π m ρ(r)r2 sin φ dθ dφ dr U(q) = i : p 2 2 0 0 0 r + z − 2rz cos φ It is an exercise to carry out the first two integrals to show Z R m0m 2 U(q) = m0 = 4πr ρ(r) dr: z 0 QED There is an alternative proof of this result, based on the fact that f(x; y; z) = 1=jq − pj is a harmonic function of p = (x; y; z), that is, fxx + fyy + fzz = 0. The well-known mean value theorem for harmonic functions states that the average value of a harmonic function over a sphere is equal to the value at the center of the sphere, which can be proved using the divergence theorem (see exercise 1.3). Fixing a value of r and applying this to the function f(x; y; z) = mρ(r)=jq − pj and the sphere Sr = fjq − q0j = rg gives Z 4πr2ρ(r) m f = : Sr jq0 − qj Then integration over 0 ≤ r ≤ R completes the proof. Although R3 is the natural home of celestial mechanics, it is useful and interesting to consider the point-mass n-body problem in Rd for any positive integer d. In this d case the position vectors are q1; : : : ; qn 2 R and the vectors q; v; p are elements of Rdn. Newton's equations (2) form a system of real-analytic, second order differential equations on the configuration space, X = Rdn n ∆, where ∆ = fq : qi = qj for some i 6= jg is the collision set. It can be transformed in the usual way into a first-order system in the phase space: dn dn TX = X × R = f(q; v): q 2 X and v 2 R g namely: q_ = v (3) v_ = M −1rU(q): The notation TX takes note of the fact that the phase space is the tangent bundle of X. The Newtonian n-body problem is to study the solutions of equations (3). A solution to (3) is a differentiable curve (q(t); v(t)) where the time t lies in some interval I. Since the differential equation is given by real-analytic functions on phase space, the solutions will be real-analytic functions of time and of their TOPICS IN CELESTIAL MECHANICS 5 initial conditions, that is, they are given locally by convergent power series. Since the phase space is not compact, it may not be possible to extend solutions for all time t 2 R. In general the maximal interval of existence will be of the form I = (a; b) with −∞ ≤ a < b ≤ 1. By the general theory of ordinary differential equations, if b < 1 then as t ! b−,(q(t); v(t)) must leave every compact subset of Rdn n ∆ × Rdn and similarly for a > −∞. For example, this can happen if q(t) converges to a collision configurationq ¯ 2 ∆ (see exercise 1.5). Exercise 1.1. Using units of kilograms for mass, meters for distance and seconds −11 m3 for time, the gravitational constant is G ' 6:674 × 10 kg·sec2 .
Load more

Recommended publications
  • Parametric Surfaces and 16.6 Their Areas Parametric Surfaces
    Parametric Surfaces and 16.6 Their Areas Parametric Surfaces 2 Parametric Surfaces Similarly to describing a space curve by a vector function r(t) of a single parameter t, a surface can be expressed by a vector function r(u, v) of two parameters u and v. Suppose r(u, v) = x(u, v)i + y(u, v)j + z(u, v)k is a vector-valued function defined on a region D in the uv-plane. So x, y, and, z, the component functions of r, are functions of the two variables u and v with domain D. The set of all points (x, y, z) in R3 s.t. x = x(u, v), y = y(u, v), z = z(u, v) and (u, v) varies throughout D, is called a parametric surface S. Typical surfaces: Cylinders, spheres, quadric surfaces, etc. 3 Example 3 – Important From Book The vector equation of a plane through (x0, y0, z0) and containing vectors is rather 4 Example – Point on Surface? Does the point (2, 3, 3) lie on the given surface? How about (1, 2, 1)? 5 Example – Identify the Surface Identify the surface with the given vector equation. 6 Example – Identify the Surface Identify the surface with the given vector equation. 7 Example – Find a Parametric Equation The part of the hyperboloid –x2 – y2 + z = 1 that lies below the rectangle [–1, 1] X [–3, 3]. 8 Example – Find a Parametric Equation The part of the cylinder x2 + z2 = 1 that lies between the planes y = 1 and y = 3. 9 Example – Find a Parametric Equation Part of the plane z = 5 that lies inside the cylinder x2 + y2 = 16.
    [Show full text]
  • The Celestial Mechanics of Newton
    GENERAL I ARTICLE The Celestial Mechanics of Newton Dipankar Bhattacharya Newton's law of universal gravitation laid the physical foundation of celestial mechanics. This article reviews the steps towards the law of gravi­ tation, and highlights some applications to celes­ tial mechanics found in Newton's Principia. 1. Introduction Newton's Principia consists of three books; the third Dipankar Bhattacharya is at the Astrophysics Group dealing with the The System of the World puts forth of the Raman Research Newton's views on celestial mechanics. This third book Institute. His research is indeed the heart of Newton's "natural philosophy" interests cover all types of which draws heavily on the mathematical results derived cosmic explosions and in the first two books. Here he systematises his math­ their remnants. ematical findings and confronts them against a variety of observed phenomena culminating in a powerful and compelling development of the universal law of gravita­ tion. Newton lived in an era of exciting developments in Nat­ ural Philosophy. Some three decades before his birth J 0- hannes Kepler had announced his first two laws of plan­ etary motion (AD 1609), to be followed by the third law after a decade (AD 1619). These were empirical laws derived from accurate astronomical observations, and stirred the imagination of philosophers regarding their underlying cause. Mechanics of terrestrial bodies was also being developed around this time. Galileo's experiments were conducted in the early 17th century leading to the discovery of the Keywords laws of free fall and projectile motion. Galileo's Dialogue Celestial mechanics, astronomy, about the system of the world was published in 1632.
    [Show full text]
  • Multivariable and Vector Calculus
    Multivariable and Vector Calculus Lecture Notes for MATH 0200 (Spring 2015) Frederick Tsz-Ho Fong Department of Mathematics Brown University Contents 1 Three-Dimensional Space ....................................5 1.1 Rectangular Coordinates in R3 5 1.2 Dot Product7 1.3 Cross Product9 1.4 Lines and Planes 11 1.5 Parametric Curves 13 2 Partial Differentiations ....................................... 19 2.1 Functions of Several Variables 19 2.2 Partial Derivatives 22 2.3 Chain Rule 26 2.4 Directional Derivatives 30 2.5 Tangent Planes 34 2.6 Local Extrema 36 2.7 Lagrange’s Multiplier 41 2.8 Optimizations 46 3 Multiple Integrations ........................................ 49 3.1 Double Integrals in Rectangular Coordinates 49 3.2 Fubini’s Theorem for General Regions 53 3.3 Double Integrals in Polar Coordinates 57 3.4 Triple Integrals in Rectangular Coordinates 62 3.5 Triple Integrals in Cylindrical Coordinates 67 3.6 Triple Integrals in Spherical Coordinates 70 4 Vector Calculus ............................................ 75 4.1 Vector Fields on R2 and R3 75 4.2 Line Integrals of Vector Fields 83 4.3 Conservative Vector Fields 88 4.4 Green’s Theorem 98 4.5 Parametric Surfaces 105 4.6 Stokes’ Theorem 120 4.7 Divergence Theorem 127 5 Topics in Physics and Engineering .......................... 133 5.1 Coulomb’s Law 133 5.2 Introduction to Maxwell’s Equations 137 5.3 Heat Diffusion 141 5.4 Dirac Delta Functions 144 1 — Three-Dimensional Space 1.1 Rectangular Coordinates in R3 Throughout the course, we will use an ordered triple (x, y, z) to represent a point in the three dimensional space.
    [Show full text]
  • Math 1300 Section 4.8: Parametric Equations A
    MATH 1300 SECTION 4.8: PARAMETRIC EQUATIONS A parametric equation is a collection of equations x = x(t) y = y(t) that gives the variables x and y as functions of a parameter t. Any real number t then corresponds to a point in the xy-plane given by the coordi- nates (x(t); y(t)). If we think about the parameter t as time, then we can interpret t = 0 as the point where we start, and as t increases, the parametric equation traces out a curve in the plane, which we will call a parametric curve. (x(t); y(t)) • The result is that we can \draw" curves just like an Etch-a-sketch! animated parametric circle from wolfram Date: 11/06. 1 4.8 2 Let's consider an example. Suppose we have the following parametric equation: x = cos(t) y = sin(t) We know from the pythagorean theorem that this parametric equation satisfies the relation x2 + y2 = 1 , so we see that as t varies over the real numbers we will trace out the unit circle! Lets look at the curve that is drawn for 0 ≤ t ≤ π. Just picking a few values we can observe that this parametric equation parametrizes the upper semi-circle in a counter clockwise direction. (0; 1) • t x(t) y(t) 0 1 0 p p π 2 2 •(1; 0) 4 2 2 π 2 0 1 π -1 0 Looking at the curve traced out over any interval of time longer that 2π will indeed trace out the entire circle.
    [Show full text]
  • Parametric Curves You Should Know
    Parametric Curves You Should Know Straight Lines Let a and c be constants which are not both zero. Then the parametric equations determining the straight line passing through (b; d) with slope c=a (i.e., the line y − d = c=a(x − b)) are: x(t) = at + b ; −∞ < t < 1: y(t) = ct + d Note that when c = 0, the line is horizontal of the form y = d, and when a = 0, the line is vertical of the form x = b. 15 10 x(t)= 2t+ 3, y(t)=t+4 5 x(t)=3-2t, y(t)= 2t+1 �������� x(t)=t+ 2, y(t)=3- 3t x(t)= 3, y(t)=2+ 3t -10 -5 5 10 15 x(t)=2+3t, y(t)=3 -5 -10 Figure 1 A collection of lines whose parametric equations are given as above. Circles Let r > 0 and let x0 and y0 be real numbers. Then the parametric equations determining the circle of radius r centered at (x0; y0) are: x(t) = r cos t + x 0 ; 0 < t < 2π: (1) y(t) = r sin t + y0 To get a circular arc instead of the full circle, restrict the t-values in (1) to t1 < t < t2. 1 1 1 2 3 4 5 2 3 4 5 (3,-1) -1 -1 (3,-1) -2 -2 -3 -3 (a) 0 ≤ t ≤ 2π (b) π=4 ≤ t ≤ 7π=6 Figure 2 The circle x(t) = 2 cos t + 3, y(t) = 2 sin t − 1 and a circular arc thereof. Note that the center is at (3; 1).
    [Show full text]
  • Polynomial Curves and Surfaces
    Polynomial Curves and Surfaces Chandrajit Bajaj and Andrew Gillette September 8, 2010 Contents 1 What is an Algebraic Curve or Surface? 2 1.1 Algebraic Curves . .3 1.2 Algebraic Surfaces . .3 2 Singularities and Extreme Points 4 2.1 Singularities and Genus . .4 2.2 Parameterizing with a Pencil of Lines . .6 2.3 Parameterizing with a Pencil of Curves . .7 2.4 Algebraic Space Curves . .8 2.5 Faithful Parameterizations . .9 3 Triangulation and Display 10 4 Polynomial and Power Basis 10 5 Power Series and Puiseux Expansions 11 5.1 Weierstrass Factorization . 11 5.2 Hensel Lifting . 11 6 Derivatives, Tangents, Curvatures 12 6.1 Curvature Computations . 12 6.1.1 Curvature Formulas . 12 6.1.2 Derivation . 13 7 Converting Between Implicit and Parametric Forms 20 7.1 Parameterization of Curves . 21 7.1.1 Parameterizing with lines . 24 7.1.2 Parameterizing with Higher Degree Curves . 26 7.1.3 Parameterization of conic, cubic plane curves . 30 7.2 Parameterization of Algebraic Space Curves . 30 7.3 Automatic Parametrization of Degree 2 Curves and Surfaces . 33 7.3.1 Conics . 34 7.3.2 Rational Fields . 36 7.4 Automatic Parametrization of Degree 3 Curves and Surfaces . 37 7.4.1 Cubics . 38 7.4.2 Cubicoids . 40 7.5 Parameterizations of Real Cubic Surfaces . 42 7.5.1 Real and Rational Points on Cubic Surfaces . 44 7.5.2 Algebraic Reduction . 45 1 7.5.3 Parameterizations without Real Skew Lines . 49 7.5.4 Classification and Straight Lines from Parametric Equations . 52 7.5.5 Parameterization of general algebraic plane curves by A-splines .
    [Show full text]
  • TOPICS in CELESTIAL MECHANICS 1. the Newtonian N-Body Problem
    TOPICS IN CELESTIAL MECHANICS RICHARD MOECKEL 1. The Newtonian n-body Problem Celestial mechanics can be defined as the study of the solution of Newton's differ- ential equations formulated by Isaac Newton in 1686 in his Philosophiae Naturalis Principia Mathematica. The setting for celestial mechanics is three-dimensional space: 3 R = fq = (x; y; z): x; y; z 2 Rg with the Euclidean norm: p jqj = x2 + y2 + z2: A point particle is characterized by a position q 2 R3 and a mass m 2 R+. A motion of such a particle is described by a curve q(t) where t runs over some interval in R; the mass is assumed to be constant. Some remarks will be made below about why it is reasonable to model a celestial body by a point particle. For every motion of a point particle one can define: velocity: v(t) =q _(t) momentum: p(t) = mv(t): Newton formulated the following laws of motion: Lex.I. Corpus omne perservare in statu suo quiescendi vel movendi uniformiter in directum, nisi quatenus a viribus impressis cogitur statum illum mutare 1 Lex.II. Mutationem motus proportionem esse vi motrici impressae et fieri secundem lineam qua vis illa imprimitur. 2 Lex.III Actioni contrarium semper et aequalem esse reactionem: sive corporum duorum actiones in se mutuo semper esse aequales et in partes contrarias dirigi. 3 The first law is statement of the principle of inertia. The second law asserts the existence of a force function F : R4 ! R3 such that: p_ = F (q; t) or mq¨ = F (q; t): In celestial mechanics, the dependence of F (q; t) on t is usually indirect; the force on one body depends on the positions of the other massive bodies which in turn depend on t.
    [Show full text]
  • Perturbation Theory in Celestial Mechanics
    Perturbation Theory in Celestial Mechanics Alessandra Celletti Dipartimento di Matematica Universit`adi Roma Tor Vergata Via della Ricerca Scientifica 1, I-00133 Roma (Italy) ([email protected]) December 8, 2007 Contents 1 Glossary 2 2 Definition 2 3 Introduction 2 4 Classical perturbation theory 4 4.1 The classical theory . 4 4.2 The precession of the perihelion of Mercury . 6 4.2.1 Delaunay action–angle variables . 6 4.2.2 The restricted, planar, circular, three–body problem . 7 4.2.3 Expansion of the perturbing function . 7 4.2.4 Computation of the precession of the perihelion . 8 5 Resonant perturbation theory 9 5.1 The resonant theory . 9 5.2 Three–body resonance . 10 5.3 Degenerate perturbation theory . 11 5.4 The precession of the equinoxes . 12 6 Invariant tori 14 6.1 Invariant KAM surfaces . 14 6.2 Rotational tori for the spin–orbit problem . 15 6.3 Librational tori for the spin–orbit problem . 16 6.4 Rotational tori for the restricted three–body problem . 17 6.5 Planetary problem . 18 7 Periodic orbits 18 7.1 Construction of periodic orbits . 18 7.2 The libration in longitude of the Moon . 20 1 8 Future directions 20 9 Bibliography 21 9.1 Books and Reviews . 21 9.2 Primary Literature . 22 1 Glossary KAM theory: it provides the persistence of quasi–periodic motions under a small perturbation of an integrable system. KAM theory can be applied under quite general assumptions, i.e. a non– degeneracy of the integrable system and a diophantine condition of the frequency of motion.
    [Show full text]
  • Lecture 1: Explicit, Implicit and Parametric Equations
    4.0 3.5 f(x+h) 3.0 Q 2.5 2.0 Lecture 1: 1.5 Explicit, Implicit and 1.0 Parametric Equations 0.5 f(x) P -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.5 1.0 1.5 2.0 2.5 3.0 3.5 x x+h -0.5 -1.0 Chapter 1 – Functions and Equations Chapter 1: Functions and Equations LECTURE TOPIC 0 GEOMETRY EXPRESSIONS™ WARM-UP 1 EXPLICIT, IMPLICIT AND PARAMETRIC EQUATIONS 2 A SHORT ATLAS OF CURVES 3 SYSTEMS OF EQUATIONS 4 INVERTIBILITY, UNIQUENESS AND CLOSURE Lecture 1 – Explicit, Implicit and Parametric Equations 2 Learning Calculus with Geometry Expressions™ Calculus Inspiration Louis Eric Wasserman used a novel approach for attacking the Clay Math Prize of P vs. NP. He examined problem complexity. Louis calculated the least number of gates needed to compute explicit functions using only AND and OR, the basic atoms of computation. He produced a characterization of P, a class of problems that can be solved in by computer in polynomial time. He also likes ultimate Frisbee™. 3 Chapter 1 – Functions and Equations EXPLICIT FUNCTIONS Mathematicians like Louis use the term “explicit function” to express the idea that we have one dependent variable on the left-hand side of an equation, and all the independent variables and constants on the right-hand side of the equation. For example, the equation of a line is: Where m is the slope and b is the y-intercept. Explicit functions GENERATE y values from x values.
    [Show full text]
  • Plotting Circles and Ellipses Parametrically Example 1: the Unit Circle
    Lesson 1 Plotting Circles and Ellipses Parametrically Example 1: The Unit Circle Let's compare traditional and parametric equations for the unit circle : Traditional : Parametric : x22 y 1 x(t) cos(t) y(t) sin(t) P(t) cos(t),sin(t) t is called a parameter 0 t 2 You can see that the parametric equation satisfies the traditional equation by substituting one into the other: x22 y 1 (cos(t))22 (sin(t)) 1 11 Created by Christopher Grattoni. All rights reserved. Example 2: Circle of Radius r Let's compare traditional and parametric equations for a circle of radius r centeredTraditional at the : origin : Parametric : 2 2 2 x y r x(t) r cos(t) y(t) r sin(t) P(t) r cos(t),sin(t) 0 t 2 Orientation: Counterclockwise Think of this as dilating the unit circle by a factor of r. Created by Christopher Grattoni. All rights reserved. Example 3: Recentering the Circle Let's compare traditional and parametric equations for a circle of radius r centeredTraditional at (h,k) : : Parametric : 2 2 2 (x h) (y k) r x(t) r cos(t) h y(t) r sin(t) k P(t) r cos(t),sin(t) (h,k) 0 t 2 Think of this as dilating the unit circle by a factor of r and translating by the point (h,k). Let's add the orientation: Created by Christopher Grattoni. All rights reserved. Example 4: Ellipses Let's compare traditional and parametric equations for an ellipse centeredTraditional at (h,k) : : Parametric : 2 2 xh yk x(t) acos(t) h 1 ab y(t) bsin(t) k P(t) acos(t),bsin(t) (h,k) 0 t 2 Think of this as dilating the unit circle by a factor of "a" in the x-direction, a factor of "b" in the y-direction, and translated by the Created by Christopherpoint Grattoni.
    [Show full text]
  • The Celestial Mechanics Approach: Theoretical Foundations
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by RERO DOC Digital Library J Geod (2010) 84:605–624 DOI 10.1007/s00190-010-0401-7 ORIGINAL ARTICLE The celestial mechanics approach: theoretical foundations Gerhard Beutler · Adrian Jäggi · Leoš Mervart · Ulrich Meyer Received: 31 October 2009 / Accepted: 29 July 2010 / Published online: 24 August 2010 © Springer-Verlag 2010 Abstract Gravity field determination using the measure- of the CMA, in particular to the GRACE mission, may be ments of Global Positioning receivers onboard low Earth found in Beutler et al. (2010) and Jäggi et al. (2010b). orbiters and inter-satellite measurements in a constellation of The determination of the Earth’s global gravity field using satellites is a generalized orbit determination problem involv- the data of space missions is nowadays either based on ing all satellites of the constellation. The celestial mechanics approach (CMA) is comprehensive in the sense that it encom- 1. the observations of spaceborne Global Positioning (GPS) passes many different methods currently in use, in particular receivers onboard low Earth orbiters (LEOs) (see so-called short-arc methods, reduced-dynamic methods, and Reigber et al. 2004), pure dynamic methods. The method is very flexible because 2. or precise inter-satellite distance monitoring of a close the actual solution type may be selected just prior to the com- satellite constellation using microwave links (combined bination of the satellite-, arc- and technique-specific normal with the measurements of the GPS receivers on all space- equation systems. It is thus possible to generate ensembles crafts involved) (see Tapley et al.
    [Show full text]
  • Parametric Curves in the Plane
    Parametric Curves in the Plane Purpose The purpose of this lab is to introduce you to curve computations using Maple for para- metric curves and vector-valued functions in the plane. Background By parametric curve in the plane, we mean a pair of equations x = f(t) and y = g(t) for t in some interval I. A vector-valued function in the plane is a function r(t) that associates a vector in the plane with each value of t in its domain. Such a vector valued function can always be written in component form as follows, r(t) = f(t)i + g(t)j where f and g are functions defined on some interval I. From our definition of a para- metric curve, it should be clear that you can always associate a parametric curve with a vector-valued function by just considering the curve traced out by the head of the vector. Defining parametric curves and vector valued functions simply in Maple The easiest way to define a vector function or a parametric curve is to use the Maple list notaion with square brackets[]. Strictly speaking, this does not define something that Maple recognizes as a vector, but it will work with all of the commands you need for this lab. >f:=t->[2*cos(t),2*sin(t)]; You can evaluate this function at any value of t in the usual way. >f(0); This is how to access a single component. You would use f(t)[2] to get the second component. >f(t)[1] Plotting and animating curves in the plane The ParamPlot command is in the CalcP package so you have to load it first.
    [Show full text]