Analytical & semi-analytical propagation of space orbits: The role of polar-nodal variables
|
|
- Aubrie Ford
- 6 years ago
- Views:
Transcription
1 Analytical & semi-analytical propagation of space orbits: The role of polar-nodal variables Martin Lara GRUCACI - Scientific Computation Group University of La Rioja, Spain Astronet II - Final meeting Tossa de Mar, June 15 19, 2015 IEEC/Marie-Curie Research Training Network on Astrodynamics support
2 Background Propagation of space orbits: analytical & S/A methods last 2 decades: eclipsed by numerical in practical apps. still may play a role in the propagation of space orbits Space debris mitigation guidelines: EoL disposal strategies guarantee some safe condition for hundreds of years S/A + global optimizers: suitable design strategy Armellin et al. Adv. Space Research 2015 Onboard orbit propagation: micro & nano satellite missions reduced power consumption is a critical requirement Efficiency in evaluating analytical expressions depends on variables used in the theory formulation Delaunay (action-angle of the Kepler problem), non-singular (Poincaré), polar-nodal (Whittaker/Hill) 1
3 Outline The two-body problem Newtonian formulation: vectorial reduction (G,e,t 0 ) Lagrangian, Hamiltonian formulation symmetries: Polar-nodal variables complete reduction by H-J: Delaunay variables Main problem of AST intermediaries Perturbed Kepler problems Recall perturbation methods Brouwer solution Normal forms, averaging Conclusions 2
4 The Kepler Problem Two bodies: Newton gravitational theory + 3rd law m 1 d 2 x 1 dt 2 = Gm 1 m 2 x 2 x 1 3(x 2 x 1 ) m 2 d 2 x 2 dt 2 = Gm 1m 2 x 2 x 1 3(x 2 x 1 ) 2nd order differential system, 6-DOF problem Add accelerations + quadratures: m 1 x 1 + m 2 x 2 = A t + B A, B, constant vectors, functions of the iicc preservation of the linear momentum of c.o.m. Relative motion: d2 x dt 2 = µ x { r = x = x2 x 1 r 3 µ = G(m 1 + m 2 ) singularity at the origin (collision orbits) 3-DOF, integrability needs 3 more independent integrals 3-1
5 By simple vector algebra (Danby s 1962 book, for instance) G = x (dx/dt) = const. (orbital plane) µe = (dx/dt) G µx/r = const. (shape, orientation) conic solution r = (G G/µ)/(1 + e cos f), f (e, x) Superintegrability: 5 independent integrals (G e = 0) corresponding to orbital elements: a, e, I, Ω, ω additive integration constant: time of periapsis passage Energy integral: E = T + V (kinetic plus potential energy) E = 1 2 dx dt dx dt µ r = µ2 (1 e e) 2G G bounded motion for E < 0 ( e < 1) 3-2
6 Lagrangian formulation: L = T V (scalar, state function) may use generalized coordinates L L(q, q), q = dq/dt eqs. of motion d ( ) L L = 0 (2nd order diff. system) dt q q Kepler problem: L = 1 dx 2 dt dx dt +µ r Hamiltonian formulation: H = Q q L, where Q = L/ q flow Hamilton equations: dq dt = H Q, dq dt = H q 1st order differential system Kepler problem: Cartesian coordinates Q X = dx/dt H = 1 2 X X µ dx r dt = X, dx dt = µ x r 3 the same as Newtonian equations, no advantage Good point: invariant under canonical transformations 3-3
7 Polar-nodal variables Polar-nodal variables (r, θ, ν, R, Θ, N), N = Θ cos I x and dx/dt define the (instantaneous) orbital plane ν not defined for I = 0 4-1
8 Match orbital (ˆx, Ĝ ˆx, Ĝ) and inertial frames: rotations x y z X Y Z = R 3 ( ν) R 1 ( I) R 3 ( θ) R 1, R 3, usual rotation matrices R 1 (φ) = cos φ sin φ 0 sin φ cos φ R 3 (φ) = r ṙ = R 0 r θ = Θ/r 0 0 cos φ sin φ 0 sin φ cos φ Canonical transformation (Whittaker s book 1904) H = 1 2 ( R 2 + Θ2 r 2 ) µ r N, Θ, ν constant decouples the differential system: 1-DOF in (r, R) dr dt = H R = R, dr dt = H r = Θ2 r 3 µ r 2 quadrature: dθ dt = H Θ = Θ r 2 θ = θ 0 + Θ 1 r(t) dt 4-2
9 Solution by H-J: Complete reduction Find canonical transform T : (x, X) (y, Y ) such that T : H(x, X) H(x(y, Y ), X(y, Y )) = Φ(, Y ) trivial integration: Y = Y 0, y = y 0 + ( Φ/ Y ) t H-J method: Generating function in mixed var. S(x, Y ) transformation: X = S/ x, y = S/ Y S solved from Jacobi eq.: H (x, S/ x) = Φ(, Y ) PDE; apparently this way complicates the problem Kepler problem: T : (r, θ, ν, R, Θ, N) (l, g, h, L, G, H) choose S S(r, θ, ν, L, G, H) = ν H + θ G + W(r, L, G) N = S ν = H, l = S L = W L, Θ = S θ = G, g = S G = θ + W G, R = S r = W r, h = S H = ν 5-1
10 Hamilton-Jacobi eq. 1 2 [ ( W r ) ] 2 + G2 r 2 µ r = Φ(L, G, ) W = 2Φ + 2µ/r G 2 /r 2 dr solved by quadrature Transformation in mixed variables (h = ν, N = H, Θ = G) R = W r = 2Φ + 2µ/r G 2 /r 2 l = W L = Φ L 1 R(r) dr g = θ + W G = θ + G 1 R(1/r) d ( ) 1 r + Φ G 1 R(r) dr Use anomalies: (1/r) = (1 + e cos f)/p, r = a(1 e cos u) p = G 2 /µ, a = µ/(2φ), e = 1 p/a l = Φ µ L ( 2Φ) 3/2(u e sin u), g = θ + f + Φ/ G Φ/ L l family of transformations parametrized by Φ 5-2
11 Standard Hamiltonian Φ Φ(L, G) = µ2 2G 2(1 e2 ) e e(l, G), nondimensional e e(η), η = G/L Φ L = µ2 ( G 2e e η ) Φ, η L G = µ2 G 3(1 e2 ) + µ2 G 2e e η family of transformations parametrized by e Action-angle variables: angle condition dl = 2π eη 2 de leads to: (1 e 2 ) 3/2 dη = 1 e2 = 1 η 2 = 1 G2 L 2 ( ) η Φ = µ2 g = θ f, l = u e sin u (Kepler equation), 2L2 Kepler probl. completely reduced Φ = Φ(L) degenerate! trivial solution g, h, L, G, H const., l = l 0 + (µ 2 /L 3 ) t l = M, g = ω, h = Ω, L = µa, G = Lη, H = G cos I Delaunay var.: singular for circular and equatorial orbits 5-3 G
12 Perturbed Kepler problems: AST Real problems: non-central potential V = µ r µ r α i i i 2 r i j=0 (C i,j cos jϑ + S i,j sin jϑ)p i,j (sin ϕ) (r, ϑ, ϕ) spherical coord., α equatorial radius P j,k associated Legendre func., C i,j, S i,j, harmonic coeff. Earth case J 2 = C 2,0 = O(10 3 ), C n,m, S n,m = O(J 2 ) 2 Main problem of AST ( H = 1 R 2 + Θ2 2 r 2 non-integrable ) [ ) 2 ( ) ] 3 µ ( α 1 J 2 r r 2 sin2 I sin 2 θ 1 2 Danby 1968, Irigoyen & Simó 93, Celleti & Negrini 95 still, a lot of info can be obtained 6-1
13 Axial symmetry: ν cyclic N constant decouple de motion of the orbital plane Invariant manifolds: polar and equatorial orbits Particular solutions: equatorial and rectilinear orbits Poincaré maps, periodic orbits,... Secular dynamics: constant term of H as a Fourier series slowly precessing ellipse regression of the perigee critical inclination: frozen-perigee orbits Short time spans: neglect second order effects intermediary solutions 6-2
14 Truncations of the main problem Expand the disturbing function: H = K + P Q + P R + P Z H = 1 2 ( µ 2r µ 2r µ 2r R 2 + Θ2 r 2 α 2 r 2 J 2 α 2 r 2 J 2 α 2 r 2 J 2 ) µ r ( 3 ) 2 sin2 I Kepler Equatorial main Cid intermediary 3 2 sin2 I cos 2θ Full problem Radial Hamiltonians: 1-DOF integrable term 1/r 3 : requires elliptic functions 6-3
15 Good intermediaries: secular effects O(J 2 ) Keplerian: K = 1 2 ( R 2 + Θ2 r 2 ) µ r, 1 2π Equatorial main problem E = K + P Q E = 1 2 ( R 2 + Θ2 r 2 ) µ r µ 2r α 2 r 2 J 2, 1 2π 2π 0 2π Cid & Lahulla intermediary C = K + P Q + P R C = 1 ( ) R 2 + Θ2 2 r 2 µ r µ α 2 ( 2r r 2 J sin2 I 2π 2π 1 1 (H C)dM = P Z dm= 0 2π 0 2π 0 paradigm of common intermediaries solution depends on elliptic functions 0 (H K)dM 0 (H E)dM 0 Radial, Θ = const: miss long-period effects in e and I ) 6-4
16 Zonal intermediaries Non-separability: related to 1/r 3 = (1 + e cos f)/(r 2 p) H = H 0 + H 1 (s sin I, c cos I, p = Θ 2 /µ) H 0 = 1 ) (R 2 + Θ2 µ 2 r 2 r µ ( α 2 2p r J s2 + 3 ) 2 s2 cos 2θ H 1 = µ ( ) ( r α 2 2r p 1 r J s2 + 3 ) 2 s2 cos 2θ H 0 2-DOF but integrable (Aksnes 1965, elliptic funct.) include long-period effects in e and I Good intermediary for low e: H 1 M = O(e 2 J 2 ) Other intermediaries Sterne 1957, Garfinkel 1958,..., Oberti 2005 Vinti 1959, Aksenov et al (2nd order effects of J 2 ) 6-5
17 Natural Intermediaries Common intermediaries: same vars. as the main problem Cid & Lahulla (1969): natural intermediary in new variables results from an infinitesimal contact transformation x = x + ɛa( x, X), X = X + ɛb( x, X) retains all secular and periodic effects up to O(J 2 ) Deprit (1981): common intermediaries can be naturalized DRI: H 0 = 1 ( ) R 2 + Θ2 2 r 2 µ r 1 Θ 2 α 2 2 r 2 p 2 J 2 (1 32 ) s2 integrable in trigonometric functions quasi-keplerian system, variable angular momentum useful for onboard orbit propagation (Gurfil & Lara 2014) If second order effects are required: perturbation theory 6-6
18 Perturbation theory Effects of non-centralities small, slightly distorted ellipses VOP method: osculating elem. x (a, e, i, Ω, ω), y = M ẋ = ɛ F (y, x, ɛ), ẏ = n(y, x) + ɛ Ψ(y, x, ɛ) Solution of VOP approached by Picard iterations recall that V = V (f), f f(m, e); Kaula expansions single iteration: main problem in closed form Hautesserres 2013 (elliptic integrals) Picard iterations may require expanding trig. functions mixed secular-periodic terms: M cos(k 1 M +... ) limit the solution to short time intervals Different perturbation methods avoid mixed terms issue books by: Giacaglia, Nayfeh, Ferraz-Mello,
19 H-J perspective: reduction by transformation of variables find (x, y) ( x, ỹ, ɛ) close to identity; ɛ small parameter x = x + j 1 ɛj j F j( x, ỹ), y = ỹ + j 1 ɛj j Ψ j( x, ỹ) such that ( x, ỹ) simpler than (ẋ, ẏ) if truncated to O(ɛ m ) Usual simplifications: reduce the number of DOF; v.g. x = 0<j<m ɛ j F j ( x) + ɛm j 0 ɛ j F j+m ( x, ỹ), ỹ = ñ( x) + 0<j<m ɛ j Ψ j ( x) + ɛm j 0 ɛ j Ψ j+m ( x, ỹ) After truncation, flow in original and tildes are different useful sol. constrained to some regions of phase space example: removing M keeps out tesseral resonances Hamiltonian perturbations: H = H 0 + ɛh 1, H 0 integrable select the form of the transformed Hamiltonian compute corresponding transformation equations 7-2
20 e Ω deg days Moon perturbations (GNSS-type orbit a = km) time history reveals a trend in evolution plus noise Averaging: filters the noise (elements const. on one period) simplifies solution of VOP by removing M days orbit evolution in mean variables evolution eqs. result from a transformation of variables Usual criterium for selecting the new H: Averaging 7-3
21 Brouwer solution to AST (AJ 1959) Model: disturbing effects of J 2, J 3, J 4 and J 5 axial symmetry, H integral, 2-DOF problem Analytical solution in Delaunay variables separates secular, long-, and short-period effects double averaging (based on von Zeipel s method) inclination resonances excluded because of g averaging trouble for low e: short-period corrections l, g Closed form solution using relation: a 2 η dm = r 2 df expansions of f = f(l, e) low convergence for high e Secular terms up to O(J 2 2 ), periodic terms up to O(J 2) Here, results for H = µ 2a + µ r m=2,3 α m r m J mp m (s sin θ) only 8-1
22 H l = µ 2a µ 1 2a 4 J 2 p η(4 2 6s2 ) + µ 3 2a 2 J 3 p (4 3 5s2 )s η e sin g [ ]. µ 3 2a 64 J 2 2 α 4 p η m 4 0,0 + m 0,1 η + m 0,2 η 2 + m 2,2 e 2 cos 2g α 2 m i,j m i,j (s) inclination polynomials s = 1 c 2, c = H/G, e = 1 η 2, η = G/L, a = L 2 /µ Hamiltonian in new (mean) variables; H and L constant Brouwer s 1st order corrections (Delaunay variables) = j k P (e)q(i) trig(jf + kg), trig sin, cos need to evaluate 20 different trigonometric terms! Brouwer s clever arrangement using r, a polar variable... still, evaluate 8 circular functions l, g: e in denominators, deteriorate for low e orbits α 3 8-2
23 H l,g = µ 2a µ 1 2a 2 J 2 p η(2 2 3s2 ) µ 2a J 2 2 α 2 ] α [m 4 p η 4 0,0 + m 0,1 η m 0,2 η 2 completely reduced Ham. in new (action-angle) variables H, L, G constant, secular terms: trivial integration Brouwer s 1st order corrections ( δg = α2 p J m 1 5c 2 0 (s) + m ) 2(s) 1 5c 2e2 sin 2g α p δl =..., δh =..., δg =..., δh = 0, δl = 0 J 3 1 ( s J 2 2 e e ) s c2 cos g trouble for low e (δl, δg), and also for low i orbits (δg, δh) Brouwer s theory does not apply to the critical inclination trouble for orbits very close to it: 1 5c 2 = O(J 2 ) essential singularity, contrary to virtual of e, s = 0 Limited in precision (desirable 2nd order corrections) 8-3
24 Non-singular variables: Lyddane s approach Poincaré canonical elements x 1 = L x 2 = 2 L G cos(g + h) x 3 = 2 G H cos h y 1 = l + g + h y 2 = 2 L G sin(g + h) y 3 = 2 G H sin h No need of recomputing the theory: use Brouwer s results (g + δg) + (h + δh) = g + h + (δg + δh) cancels terms in 1/s 1 c cos(h + δh) = 1 c cos h ( 1 c δh) sin h + O(J 2 2 ) avoid trouble with low I (l + δl) + (g + δg) + (h + δh) = l + g + h + δl + δg 1 + δg 2 + δh cancels terms in 1/s and in 1/e No singularities, but more involved series 9
25 Simplification: Polar-nodal variables Further simplifications to Brouwer s formulas reformulate Brouwer s corrections in polar variables l = l(r, θ, R, Θ, N), g =..., h =... trig. functions reduced to one half! cos 2θ, sin 2θ, r = p/(1 + e cos f), R = (G/p)e sin f still remain singularities for low e and I Izsak (1963) approach: removes singularity for zero e orbits theory correctly computed in action-angle variables, but corrections to polar-nodal variables instead of Delaunay r = r(r, θ, R, Θ, N), θ =..., ν =..., R =... dramatically reduces size of the short-period corrections Extended by Aksnes (1978) to the long-period corrections 10-1
26 [ (2 r = γp ) ( 3s κ 1 + η + 2η ) ] s 2 cos 2θ 1 + κ { [ ] θ = γ 3 ( 4 5s 2) φ s2 + ( 4 6s 2) κ sin 2θ [ ( +2σ 5 6s ) 2 + κ 2 s2 1 + η + ( ]} 1 2s 2) cos 2θ [ ] ν = γc 6(σ + φ) (3 + 4κ) sin 2θ + 2σ cos 2θ R = γ Θ { p (1 + κ)2 2s 2 sin 2θ ( [ ] } 2 3s 2) η + η σ (1 + κ) ] 2 Θ = γθs [( κ) cos 2θ + 2σ sin 2θ N = 0 φ eq. of the center α 2 γ = 1 4 C 2,0 p, σ = pr 2 Θ, κ = p r 1, p = Θ2 µ, 10-2
27 Third-Body perturbations Mass-point approx.: V 3 = µ 3 r 3 Expand in Legendre poly: V 3 = β n2 3 a3 3 r 3 ( r 3 x r 3 x r 3 j 2 r3 2 ( r r 3 ) ) j P j (cos ψ) ψ (x, r 3 ), β = m 3 /(m 3 + m) 3rd-body reduced mass Moon perturbations: slow convergence of the expansion Integral may require up to P 6, SIMOLX up to P 8... very long correction series! Perturbation theory by Lie transforms (Hori 67, Deprit 69) theory properly computed in Delaunay variables transformation to mean variables from generating func. closed form: dl = (1 e sin u)du 11-1
28 Example: P 2 only Delaunay var.: generating function (trouble for low e, i) W 1 = L n2 3 n 2 β a3 3 r j=0 l= 2m= 3 m 0 C j,m (e) Q j,l (I) [ T l (ˆr 3 ) cos(2jω + lω + mu) S l (ˆr 3 ) sin(2jω + lω + mu) u eccentric anomaly, ˆr 3 moon direction (ephemeris) Fourier series with 81 different trig. terms Polar-nodal var.: (no trouble for low e, trouble for low i) V 1 = Θ n2 3 n 2 β a3 3 r 3 3 r 3 p j=0 l= 2 Q j,l [ ] (κ j T l σ j S l ) cos(2jθ + lν) (κ j S l + σ j T l ) sin(2jθ + lν) κ j κ j (r, R, Θ), σ j σ j (r, R, Θ), e-type coefficients polar-nodal: only 14 trigonometric terms! Time in evaluating the corrections reduces to < 1 3!! 11-2 ]
29 The case of low inclinations Polar-nodal: singularities remain for low inclinations Alternative: treat separately that case Aksnes (1972): ephemeris in spherical coord. (r, ϕ, ϑ) sin ϕ = s sin θ, sin(ϑ ν) = c sin θ, cos(ϑ ν) = cos θ s = O(J 2 ) compute ϑ = θ + ν, ϕ = (s sin θ) Alternative: reformulate corrections in nonsingular variables non-canonical elements based on polar-nodal variables ψ = θ + ν, ξ = s sin θ, χ = s cos θ, r, R, Θ evaluation: similar performances as polar-nodal variables may be used only for the case of low I dramatic simplifications of the corrections 11-3
30 Conclusions S/A theories still useful for current engineering problems properly constructed in Delaunay variables virtual singularities for circular and equatorial orbits Reformulate the theory in Poincaré nonsingular variables make the theory as much general as possible, but... evaluation of long Fourier series degrade performance Polar-nodal var.: compact form, straightforward evaluation almost equatorial orbits must be treated separately lack of universality fixed with a conditional statement compensated by the simplif. of the theory for low I New nonsingular variables based on the polar-nodal set similar evaluation performances as polar-nodal 12
Previous Lecture. The Von Zeipel Method. Application 1: The Brouwer model. Application 2: The Cid-Lahulla model. Simplified Brouwer transformation.
2 / 36 Previous Lecture The Von Zeipel Method. Application 1: The Brouwer model. Application 2: The Cid-Lahulla model. Simplified Brouwer transformation. Review of Analytic Models 3 / 36 4 / 36 Review:
More informationA SEMI-ANALYTICAL ORBIT PROPAGATOR PROGRAM FOR HIGHLY ELLIPTICAL ORBITS
A SEMI-ANALYTICAL ORBIT PROPAGATOR PROGRAM FOR HIGHLY ELLIPTICAL ORBITS M. Lara, J. F. San Juan and D. Hautesserres Scientific Computing Group and Centre National d Études Spatiales 6th International Conference
More informationCeres Rotation Solution under the Gravitational Torque of the Sun
Ceres Rotation Solution under the Gravitational Torque of the Sun Martin Lara, Toshio Fukushima, Sebastián Ferrer (*) Real Observatorio de la Armada, San Fernando, Spain ( ) National Astronomical Observatory,
More informationNATURAL INTERMEDIARIES AS ONBOARD ORBIT PROPAGATORS
(Preprint) IAA-AAS-DyCoSS2-05-02 NATURAL INTERMEDIARIES AS ONBOARD ORBIT PROPAGATORS Pini Gurfil and Martin Lara INTRODUCTION Short-term satellite onboard orbit propagation is required when GPS position
More informationEFFICIENT DESIGN OF LOW LUNAR ORBITS BASED ON KAULA RECURSIONS. GRUCACI University of La Rioja C/ Madre de Dios 53, Logroño, La Rioja, Spain
EFFICIENT DESIGN OF LOW LUNAR ORBITS BASED ON KAULA RECURSIONS Martin Lara, Rosario López, Iván Pérez, Juan F. San-Juan GRUCACI University of La Rioja C/ Madre de Dios 53, 6006 Logroño, La Rioja, Spain
More informationarxiv: v1 [math.ds] 10 Jun 2014
On Inclination Resonances in Artificial Satellite Theory arxiv:1406.2634v1 [math.ds] 10 Jun 2014 Abstract Martin Lara 1 Columnas de Hércules 1, ES-11000 San Fernando, Spain The frozen-perigee behavior
More informationOrbital and Celestial Mechanics
Orbital and Celestial Mechanics John P. Vinti Edited by Gim J. Der TRW Los Angeles, California Nino L. Bonavito NASA Goddard Space Flight Center Greenbelt, Maryland Volume 177 PROGRESS IN ASTRONAUTICS
More informationA SEMI-ANALYTICAL ORBIT PROPAGATOR PROGRAM FOR HIGHLY ELLIPTICAL ORBITS
A SEMI-ANALYTICAL ORBIT PROPAGATOR PROGRAM FOR HIGHLY ELLIPTICAL ORBITS Martin Lara & Juan F. San-Juan University of La Rioja GRUCACI Scientific Computation Group 264 Logroño, Spain Denis Hautesserres
More informationThe two body problem involves a pair of particles with masses m 1 and m 2 described by a Lagrangian of the form:
Physics 3550, Fall 2011 Two Body, Central-Force Problem Relevant Sections in Text: 8.1 8.7 Two Body, Central-Force Problem Introduction. I have already mentioned the two body central force problem several
More informationSymbolic Solution of Kepler s Generalized Equation
Symbolic Solution of Kepler s Generalized Equation Juan Félix San-Juan 1 and Alberto Abad 1 Universidad de La Rioja, 6004 Logroño, Spain juanfelix.sanjuan@dmc.unirioja.es, Grupo de Mecánica Espacial, Universidad
More informationMassachusetts Institute of Technology Department of Physics. Final Examination December 17, 2004
Massachusetts Institute of Technology Department of Physics Course: 8.09 Classical Mechanics Term: Fall 004 Final Examination December 17, 004 Instructions Do not start until you are told to do so. Solve
More informationPhysical Dynamics (PHY-304)
Physical Dynamics (PHY-304) Gabriele Travaglini March 31, 2012 1 Review of Newtonian Mechanics 1.1 One particle Lectures 1-2. Frame, velocity, acceleration, number of degrees of freedom, generalised coordinates.
More informationLecture 1: Oscillatory motions in the restricted three body problem
Lecture 1: Oscillatory motions in the restricted three body problem Marcel Guardia Universitat Politècnica de Catalunya February 6, 2017 M. Guardia (UPC) Lecture 1 February 6, 2017 1 / 31 Outline of the
More informationarxiv: v1 [math.ds] 27 Oct 2018
Celestial Mechanics and Dynamical Astronomy manuscript No. (will be inserted by the editor) Element sets for high-order Poincaré mapping of perturbed Keplerian motion David J. Gondelach Roberto Armellin
More informationANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS
ANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS A1.1. Kepler s laws Johannes Kepler (1571-1630) discovered the laws of orbital motion, now called Kepler's laws.
More informationM2A2 Problem Sheet 3 - Hamiltonian Mechanics
MA Problem Sheet 3 - Hamiltonian Mechanics. The particle in a cone. A particle slides under gravity, inside a smooth circular cone with a vertical axis, z = k x + y. Write down its Lagrangian in a) Cartesian,
More informationTwo-Body Problem. Central Potential. 1D Motion
Two-Body Problem. Central Potential. D Motion The simplest non-trivial dynamical problem is the problem of two particles. The equations of motion read. m r = F 2, () We already know that the center of
More informationThird Body Perturbation
Third Body Perturbation p. 1/30 Third Body Perturbation Modeling the Space Environment Manuel Ruiz Delgado European Masters in Aeronautics and Space E.T.S.I. Aeronáuticos Universidad Politécnica de Madrid
More informationPhysical Dynamics (SPA5304) Lecture Plan 2018
Physical Dynamics (SPA5304) Lecture Plan 2018 The numbers on the left margin are approximate lecture numbers. Items in gray are not covered this year 1 Advanced Review of Newtonian Mechanics 1.1 One Particle
More informationSolution Set Two. 1 Problem #1: Projectile Motion Cartesian Coordinates Polar Coordinates... 3
: Solution Set Two Northwestern University, Classical Mechanics Classical Mechanics, Third Ed.- Goldstein October 7, 2015 Contents 1 Problem #1: Projectile Motion. 2 1.1 Cartesian Coordinates....................................
More informationThe two-body Kepler problem
The two-body Kepler problem set center of mass at the origin (X = 0) ignore all multipole moments (spherical bodies or point masses) define r := r 1 r 2,r:= r,m:= m 1 + m 2,µ:= m 1 m 2 /m reduces to effective
More informationP321(b), Assignement 1
P31(b), Assignement 1 1 Exercise 3.1 (Fetter and Walecka) a) The problem is that of a point mass rotating along a circle of radius a, rotating with a constant angular velocity Ω. Generally, 3 coordinates
More informationLong-term passive distance-bounded relative motion in the presence of J 2 perturbations
Celest Mech Dyn Astr (015) 11:385 413 DOI 10.1007/s10569-015-9603-x ORIGINAL ARTICLE Long-term passive distance-bounded relative motion in the presence of J perturbations Jing Chu Jian Guo EberhardK.A.Gill
More informationCentral force motion/kepler problem. 1 Reducing 2-body motion to effective 1-body, that too with 2 d.o.f and 1st order differential equations
Central force motion/kepler problem This short note summarizes our discussion in the lectures of various aspects of the motion under central force, in particular, the Kepler problem of inverse square-law
More informationThe Restricted 3-Body Problem
The Restricted 3-Body Problem John Bremseth and John Grasel 12/10/2010 Abstract Though the 3-body problem is difficult to solve, it can be modeled if one mass is so small that its effect on the other two
More informationThe first order formalism and the transition to the
Problem 1. Hamiltonian The first order formalism and the transition to the This problem uses the notion of Lagrange multipliers and Legendre transforms to understand the action in the Hamiltonian formalism.
More informationThe restricted, circular, planar three-body problem
The restricted, circular, planar three-body problem Luigi Chierchia Dipartimento di Matematica Università Roma Tre Largo S L Murialdo 1, I-00146 Roma (Italy) (luigi@matuniroma3it) March, 2005 1 The restricted
More informationRestricted three body problems in the Solar System: simulations
Author:. Facultat de Física, Universitat de Barcelona, Diagonal 645, 0808 Barcelona, Spain. Advisor: Antoni Benseny i Ardiaca. Facultat de Matemàtiques, Universitat de Barcelona, Gran Via de les Corts
More informationAstromechanics. 10. The Kepler Problem
Astromechanics 10. The Kepler Problem One of the fundamental problems in astromechanics is the Kepler problem The Kepler problem is stated as follows: Given the current position a velocity vectors and
More informationTwo dimensional oscillator and central forces
Two dimensional oscillator and central forces September 4, 04 Hooke s law in two dimensions Consider a radial Hooke s law force in -dimensions, F = kr where the force is along the radial unit vector and
More informationAnalytical Method for Space Debris propagation under perturbations in the geostationary ring
Analytical Method for Space Debris propagation under perturbations in the geostationary ring July 21-23, 2016 Berlin, Germany 2nd International Conference and Exhibition on Satellite & Space Missions Daniel
More informationA SYMPLECTIC KEPLERIAN MAP FOR PERTURBED TWO-BODY DYNAMICS. Oier Peñagaricano Muñoa 1, Daniel J. Scheeres 2
AIAA/AAS Astrodynamics Specialist Conference and Exhibit 8 - August 008, Honolulu, Hawaii AIAA 008-7068 AIAA 008-7068 A SYMPLECTIC KEPLERIAN MAP FOR PERTURBED TWO-BODY DYNAMICS Oier Peñagaricano Muñoa,
More informationTHE exploration of planetary satellites is currently an active area
JOURNAL OF GUIDANCE, CONTROL, AND DYNAMICS Vol. 9, No. 5, September October 6 Design of Science Orbits About Planetary Satellites: Application to Europa Marci E. Paskowitz and Daniel J. Scheeres University
More informationSOLUTIONS, PROBLEM SET 11
SOLUTIONS, PROBLEM SET 11 1 In this problem we investigate the Lagrangian formulation of dynamics in a rotating frame. Consider a frame of reference which we will consider to be inertial. Suppose that
More informationOrbital Mechanics! Space System Design, MAE 342, Princeton University! Robert Stengel
Orbital Mechanics Space System Design, MAE 342, Princeton University Robert Stengel Conic section orbits Equations of motion Momentum and energy Kepler s Equation Position and velocity in orbit Copyright
More informationLecture 21 Gravitational and Central Forces
Lecture 21 Gravitational and Central Forces 21.1 Newton s Law of Universal Gravitation According to Newton s Law of Universal Graviation, the force on a particle i of mass m i exerted by a particle j of
More informationISIMA lectures on celestial mechanics. 1
ISIMA lectures on celestial mechanics. 1 Scott Tremaine, Institute for Advanced Study July 2014 The roots of solar system dynamics can be traced to two fundamental discoveries by Isaac Newton: first, that
More informationAnalytical theory of a lunar artificial satellite with third body perturbations
Celestial Mechanics and Dynamical Astronomy (6) 95:47 43 DOI.7/s569-6-99-6 ORIGINAL ARTICLE Analytical theory of a lunar artificial satellite with third body perturbations Bernard De Saedeleer Received:
More informationCelestial Mechanics II. Orbital energy and angular momentum Elliptic, parabolic and hyperbolic orbits Position in the orbit versus time
Celestial Mechanics II Orbital energy and angular momentum Elliptic, parabolic and hyperbolic orbits Position in the orbit versus time Orbital Energy KINETIC per unit mass POTENTIAL The orbital energy
More informationAnalysis of frozen orbits for solar sails
Trabalho apresentado no XXXV CNMAC, Natal-RN, 2014. Analysis of frozen orbits for solar sails J. P. S. Carvalho, R. Vilhena de Moraes, Instituto de Ciência e Tecnologia, UNIFESP, São José dos Campos -
More informationChapter 2 Introduction to Binary Systems
Chapter 2 Introduction to Binary Systems In order to model stars, we must first have a knowledge of their physical properties. In this chapter, we describe how we know the stellar properties that stellar
More informationLecture 2c: Satellite Orbits
Lecture 2c: Satellite Orbits Outline 1. Newton s Laws of Mo3on 2. Newton s Law of Universal Gravita3on 3. Kepler s Laws 4. Pu>ng Newton and Kepler s Laws together and applying them to the Earth-satellite
More informationSolution of Liouville s Equation for Uncertainty Characterization of the Main Problem in Satellite Theory 1,2
Copyright 016 Tech Science Press CMES, vol.111, no.3, pp.69-304, 016 Solution of Liouville s Equation for Uncertainty Characterization of the Main Problem in Satellite Theory 1, Ryan Weisman 3, Manoranjan
More informationQuestion 1: Spherical Pendulum
Question 1: Spherical Pendulum Consider a two-dimensional pendulum of length l with mass M at its end. It is easiest to use spherical coordinates centered at the pivot since the magnitude of the position
More information06. Lagrangian Mechanics II
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 06. Lagrangian Mechanics II Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License
More informationOrbital dynamics in the tidally-perturbed Kepler and Schwarzschild systems. Sam R. Dolan
Orbital dynamics in the tidally-perturbed Kepler and Schwarzschild systems Sam R. Dolan Gravity @ All Scales, Nottingham, 24th Aug 2015 Sam Dolan (Sheffield) Perturbed dynamics Nottingham 1 / 67 work in
More informationSATELLITE RELATIVE MOTION PROPAGATION AND CONTROL. A Thesis PRASENJIT SENGUPTA
SATELLITE RELATIVE MOTION PROPAGATION AND CONTROL IN THE PRESENCE OF J 2 PERTURBATIONS A Thesis by PRASENJIT SENGUPTA Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment
More informationChapter 5 - Part 1. Orbit Perturbations. D.Mortari - AERO-423
Chapter 5 - Part 1 Orbit Perturbations D.Mortari - AERO-43 Orbital Elements Orbit normal i North Orbit plane Equatorial plane ϕ P O ω Ω i Vernal equinox Ascending node D. Mortari - AERO-43 Introduction
More informationPHY411 Lecture notes Part 2
PHY411 Lecture notes Part 2 Alice Quillen April 6, 2017 Contents 1 Canonical Transformations 2 1.1 Poisson Brackets................................. 2 1.2 Canonical transformations............................
More informationSpacecraft Dynamics and Control
Spacecraft Dynamics and Control Matthew M. Peet Arizona State University Lecture 5: Hyperbolic Orbits Introduction In this Lecture, you will learn: Hyperbolic orbits Hyperbolic Anomaly Kepler s Equation,
More informationA645/A445: Exercise #1. The Kepler Problem
A645/A445: Exercise #1 The Kepler Problem Due: 2017 Sep 26 1 Numerical solution to the one-degree-of-freedom implicit equation The one-dimensional implicit solution is: t = t o + x x o 2(E U(x)). (1) The
More informationContinuum Polarization Induced by Tidal Distortion in Binary Stars
Continuum Polarization Induced by Tidal Distortion in Binary Stars J. Patrick Harrington 1 1. On the Roche Potential of Close Binary Stars Let Ψ be the potential of a particle due to the gravitational
More informationHamiltonian Lecture notes Part 3
Hamiltonian Lecture notes Part 3 Alice Quillen March 1, 017 Contents 1 What is a resonance? 1 1.1 Dangers of low order approximations...................... 1. A resonance is a commensurability.......................
More informationUse conserved quantities to reduce number of variables and the equation of motion (EOM)
Physics 106a, Caltech 5 October, 018 Lecture 8: Central Forces Bound States Today we discuss the Kepler problem of the orbital motion of planets and other objects in the gravitational field of the sun.
More informationAn Introduction to Celestial Mechanics
An Introduction to Celestial Mechanics This accessible text on classical celestial mechanics the principles governing the motions of bodies in the solar system provides a clear and concise treatment of
More informationLunar Mission Analysis for a Wallops Flight Facility Launch
Lunar Mission Analysis for a Wallops Flight Facility Launch John M. Dolan III Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University in partial fulfillment of the requirements
More informationREVIEW. Hamilton s principle. based on FW-18. Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws!
Hamilton s principle Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws! based on FW-18 REVIEW the particle takes the path that minimizes the integrated difference
More informationA qualitative analysis of bifurcations to halo orbits
1/28 spazio A qualitative analysis of bifurcations to halo orbits Dr. Ceccaroni Marta ceccaron@mat.uniroma2.it University of Roma Tor Vergata Work in collaboration with S. Bucciarelli, A. Celletti, G.
More informationCALCULATION OF POSITION AND VELOCITY OF GLONASS SATELLITE BASED ON ANALYTICAL THEORY OF MOTION
ARTIFICIAL SATELLITES, Vol. 50, No. 3 2015 DOI: 10.1515/arsa-2015-0008 CALCULATION OF POSITION AND VELOCITY OF GLONASS SATELLITE BASED ON ANALYTICAL THEORY OF MOTION W. Góral, B. Skorupa AGH University
More informationLecture XIX: Particle motion exterior to a spherical star
Lecture XIX: Particle motion exterior to a spherical star Christopher M. Hirata Caltech M/C 350-7, Pasadena CA 95, USA Dated: January 8, 0 I. OVERVIEW Our next objective is to consider the motion of test
More informationTHIRD-BODY PERTURBATION USING A SINGLE AVERAGED MODEL
INPE-1183-PRE/67 THIRD-BODY PERTURBATION USING A SINGLE AVERAGED MODEL Carlos Renato Huaura Solórzano Antonio Fernando Bertachini de Almeida Prado ADVANCES IN SPACE DYNAMICS : CELESTIAL MECHANICS AND ASTRONAUTICS,
More informationPerturbations to the Lunar Orbit
Perturbations to the Lunar Orbit th January 006 Abstract In this paper, a general approach to performing perturbation analysis on a two-dimensional orbit is presented. The specific examples of the solar
More informationTowards stability results for planetary problems with more than three bodies
Towards stability results for planetary problems with more than three bodies Ugo Locatelli [a] and Marco Sansottera [b] [a] Math. Dep. of Università degli Studi di Roma Tor Vergata [b] Math. Dep. of Università
More informationAPPENDIX B SUMMARY OF ORBITAL MECHANICS RELEVANT TO REMOTE SENSING
APPENDIX B SUMMARY OF ORBITAL MECHANICS RELEVANT TO REMOTE SENSING Orbit selection and sensor characteristics are closely related to the strategy required to achieve the desired results. Different types
More informationFundamentals of Astrodynamics and Applications
Fundamentals of Astrodynamics and Applications Third Edition David A. Vallado with technical contributions by Wayne D. McClain Space Technology Library Published Jointly by Microcosm Press Hawthorne, CA
More informationMarion and Thornton. Tyler Shendruk October 1, Hamilton s Principle - Lagrangian and Hamiltonian dynamics.
Marion and Thornton Tyler Shendruk October 1, 2010 1 Marion and Thornton Chapter 7 Hamilton s Principle - Lagrangian and Hamiltonian dynamics. 1.1 Problem 6.4 s r z θ Figure 1: Geodesic on circular cylinder
More informationHAMILTON S PRINCIPLE
HAMILTON S PRINCIPLE In our previous derivation of Lagrange s equations we started from the Newtonian vector equations of motion and via D Alembert s Principle changed coordinates to generalised coordinates
More informationFormation flying in elliptic orbits with the J 2 perturbation
Research in Astron. Astrophys. 2012 Vol. 12 No. 11, 1563 1575 http://www.raa-journal.org http://www.iop.org/journals/raa Research in Astronomy and Astrophysics Formation flying in elliptic orbits with
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationIdentifying Safe Zones for Planetary Satellite Orbiters
AIAA/AAS Astrodynamics Specialist Conference and Exhibit 16-19 August 2004, Providence, Rhode Island AIAA 2004-4862 Identifying Safe Zones for Planetary Satellite Orbiters M.E. Paskowitz and D.J. Scheeres
More informationThe Two -Body Central Force Problem
The Two -Body Central Force Problem Physics W3003 March 6, 2015 1 The setup 1.1 Formulation of problem The two-body central potential problem is defined by the (conserved) total energy E = 1 2 m 1Ṙ2 1
More informationPHY321 Homework Set 10
PHY321 Homework Set 10 1. [5 pts] A small block of mass m slides without friction down a wedge-shaped block of mass M and of opening angle α. Thetriangular block itself slides along a horizontal floor,
More informationHamilton-Jacobi Modelling of Stellar Dynamics
Hamilton-Jacobi Modelling of Stellar Dynamics Pini Gurfil Faculty of Aerospace Engineering, Technion - Israel Institute of Technology Haifa 32000, Israel N. Jeremy Kasdin and Egemen Kolemen Mechanical
More informationOrbital Motion in Schwarzschild Geometry
Physics 4 Lecture 29 Orbital Motion in Schwarzschild Geometry Lecture 29 Physics 4 Classical Mechanics II November 9th, 2007 We have seen, through the study of the weak field solutions of Einstein s equation
More informationBeyond Janus & Epimetheus: Momentum Trading Among Co-Orbiting Satellite Groups
Beyond Janus & Epimetheus: Momentum Trading Among Co-Orbiting Satellite Groups DOUG BALCOM U NIVERSITY OF WASHINGTON APPLIED MATHEMATICS Special Thanks to Sasha Malinsky Janus and Epimetheus: Momentum
More informationAnalysis of Lunisolar Resonances. in an Artificial Satellite Orbits
Applied Mathematical Sciences, Vol., 008, no., 0 0 Analysis of Lunisolar Resonances in an Artificial Satellite Orbits F. A. Abd El-Salam, Yehia A. Abdel-Aziz,*, M. El-Saftawy, and M. Radwan Cairo university,
More informationResearch Article An Economic Hybrid J 2 Analytical Orbit Propagator Program Based on SARIMA Models
Mathematical Problems in Engineering Volume 212, Article ID 27381, 15 pages doi:1.1155/212/27381 Research Article An Economic Hybrid J 2 Analytical Orbit Propagator Program Based on SARIMA Models Juan
More informationI ve Got a Three-Body Problem
I ve Got a Three-Body Problem Gareth E. Roberts Department of Mathematics and Computer Science College of the Holy Cross Mathematics Colloquium Fitchburg State College November 13, 2008 Roberts (Holy Cross)
More informationOrbit Characteristics
Orbit Characteristics We have shown that the in the two body problem, the orbit of the satellite about the primary (or vice-versa) is a conic section, with the primary located at the focus of the conic
More informationFrom the Earth to the Moon: the weak stability boundary and invariant manifolds -
From the Earth to the Moon: the weak stability boundary and invariant manifolds - Priscilla A. Sousa Silva MAiA-UB - - - Seminari Informal de Matemàtiques de Barcelona 05-06-2012 P.A. Sousa Silva (MAiA-UB)
More information= 0. = q i., q i = E
Summary of the Above Newton s second law: d 2 r dt 2 = Φ( r) Complicated vector arithmetic & coordinate system dependence Lagrangian Formalism: L q i d dt ( L q i ) = 0 n second-order differential equations
More informationLong-Term Evolution of High Earth Orbits: Effects of Direct Solar Radiation Pressure and Comparison of Trajectory Propagators
Long-Term Evolution of High Earth Orbits: Effects of Direct Solar Radiation Pressure and Comparison of Trajectory Propagators by L. Anselmo and C. Pardini (Luciano.Anselmo@isti.cnr.it & Carmen.Pardini@isti.cnr.it)
More informationfor changing independent variables. Most simply for a function f(x) the Legendre transformation f(x) B(s) takes the form B(s) = xs f(x) with s = df
Physics 106a, Caltech 1 November, 2018 Lecture 10: Hamiltonian Mechanics I The Hamiltonian In the Hamiltonian formulation of dynamics each second order ODE given by the Euler- Lagrange equation in terms
More information[#1] R 3 bracket for the spherical pendulum
.. Holm Tuesday 11 January 2011 Solutions to MSc Enhanced Coursework for MA16 1 M3/4A16 MSc Enhanced Coursework arryl Holm Solutions Tuesday 11 January 2011 [#1] R 3 bracket for the spherical pendulum
More informationQuiz # 2. Answer ALL Questions. 1. Find the Fourier series for the periodic extension of. 0, 2 <t<0 f(t) = sin(πt/2), 0 <t<2
Quiz # Answer ALL Questions 1 Find the Fourier series for the periodic extension of ½,
More informationSecular Evolution of Extrasolar Planetary Systems:
Secular Evolution of Extrasolar Planetary Systems: an Extension of the Laplace-Lagrange Secular Theory Marco Sansottera [a] [a] Namur Center for Complex Systems (naxys) Based on a research work in collaboration
More informationProblem Set 5 Solution
Problem Set 5 Solution Friday, 14 October 011 Physics 111 Proble The Setup For a pair of point masses, and m, interacting via a conservative force directed along the line joining the masses, show that
More informationDynamical properties of the Solar System. Second Kepler s Law. Dynamics of planetary orbits. ν: true anomaly
First Kepler s Law The secondary body moves in an elliptical orbit, with the primary body at the focus Valid for bound orbits with E < 0 The conservation of the total energy E yields a constant semi-major
More informationMotion under the Influence of a Central Force
Copyright 004 5 Motion under the Influence of a Central Force The fundamental forces of nature depend only on the distance from the source. All the complex interactions that occur in the real world arise
More informationCopyright. Ashley Darius Biria
Copyright by Ashley Darius Biria 2017 The Dissertation Committee for Ashley Darius Biria certifies that this is the approved version of the following dissertation: Revisiting Vinti Theory: Generalized
More informationCentral Forces. c 2007 Corinne Manogue. March 1999, last revised February 2009
Central Forces Corinne A. Manogue with Tevian Dray, Kenneth S. Krane, Jason Janesky Department of Physics Oregon State University Corvallis, OR 97331, USA corinne@physics.oregonstate.edu c 2007 Corinne
More informationCurves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations,
Physics 6010, Fall 2010 Hamiltonian Formalism: Hamilton s equations. Conservation laws. Reduction. Poisson Brackets. Relevant Sections in Text: 8.1 8.3, 9.5 The Hamiltonian Formalism We now return to formal
More informationTheorem of the keplerian kinematics
1 Theorem of the keplerian kinematics Hervé Le Cornec, France, herve.le.cornec@free.fr Abstract : Any mobile having a velocity which is the addition of a rotation velocity and a translation velocity, both
More informationCentrifugal force in Kerr geometry
Centrifugal force in Kerr geometry Sai Iyer and A R Prasanna Physical Research Laboratory Ahmedabad 380009 INDIA Abstract We have obtained the correct expression for the centrifugal force acting on a particle
More informationResearch Article Variation of the Equator due to a Highly Inclined and Eccentric Disturber
Hindawi Publishing Corporation Mathematical Problems in Engineering Volume 009, Article ID 467865, 10 pages doi:10.1155/009/467865 Research Article Variation of the Equator due to a Highly Inclined and
More informationTheory of mean motion resonances.
Theory of mean motion resonances. Mean motion resonances are ubiquitous in space. They can be found between planets and asteroids, planets and rings in gaseous disks or satellites and planetary rings.
More informationComputational Problem: Keplerian Orbits
Computational Problem: Keplerian Orbits April 10, 2006 1 Part 1 1.1 Problem For the case of an infinite central mass and an orbiting test mass, integrate a circular orbit and an eccentric orbit. Carry
More informationHamiltonian. March 30, 2013
Hamiltonian March 3, 213 Contents 1 Variational problem as a constrained problem 1 1.1 Differential constaint......................... 1 1.2 Canonic form............................. 2 1.3 Hamiltonian..............................
More informationAs a starting point of our derivation of the equations of motion for a rigid body, we employ d Alembert s principle:
MEG6007: Advanced Dynamics -Principles and Computational Methods (Fall, 017) Lecture 10: Equations of Motion for a Rigid Body 10.1 D Alembert s Principle for A Free Rigid Body As a starting point of our
More informationNonlinear Control of Electrodynamic Tether in Equatorial or Somewhat Inclined Orbits
Proceedings of the 15th Mediterranean Conference on Control & Automation, July 7-9, 7, Athens - Greece T-5 Nonlinear Control of Electrodynamic Tether in Equatorial or Somewhat Inclined Martin B. Larsen
More information