The rotation of Io with a liquid core
|
|
- Monica Williams
- 5 years ago
- Views:
Transcription
1 Celest Mech Dyn Astr (008) 101:1 1 DOI /s ORIGINAL ARTICLE The rotation of Io with a liquid core Jacques Henrard Received: 1 October 007 / Accepted: 0 March 008 / Published online: 14 May 008 Springer Science+Business Media B.V. 008 Abstract In a previous paper (Henrard, Celest. Mech. Dyn. Astron. 178, , 005c) we have developed an analytical theory of the rotation of the Galilean satellite Io, considered as a rigid body and based on a synthetic theory of its orbital motion due to Lainey (Théorie Dynamique des Satellites Galiléens. PhD dissertation, Observatoire de Paris, 00) (see also Lainey et al., A&A, 40, a; A&A, 47, , 004b). One of the most important causes of departure of the actual rotation from the rigid theory is thought to be the existence of a liquid core, the size of which is unknown but would be an important piece of information concerning the structure of the interior of the satellite. In this contribution we develop the analytical theory of a liquid core contained in a cavity filled by an inviscid fluid of constant uniform density and vorticity. The theory is based on Poincaré (Bull. Astron. 7, , 1910) model and is developed by a Lie transform perturbation method, very much like in our previous contribution. Our main conclusion is that the addition of a degree of freedom (the spin of the core) with a frequency close to the orbital frequency multiplies the possibility of resonances and that for some particular size of the core one may expect a (possibly small) region of chaotic behaviour in the vicinity of the Cassini state. Keywords Rotational dynamics Liquid core Jovian planets Io 1 Introduction Due to its proximity of Jupiter and the complexity of its orbits captured in the Laplacian resonance with Europe and Ganymede, the rotational motion of Io is far from an uniform rotation around a Cassini s equilibrium. We have developed in a previous contribution (Henrard 005c) a theory based on the assumption of a rigid satellite and a model of the orbital motion given by the synthetic theory of Lainey (Lainey 00). J. Henrard (B) University of Namur, 8, Rempart de la Vierge, 5000 Namur, Belgium jhenrard@math.fundp.ac.be
2 J. Henrard The volcanic activity of Io is a clear evidence that it cannot be considered as a rigid body as the moon might be, but that the existence of a core, more or less fluid must be considered. In order to investigate the influence of such a core on the rotational motion of the satellite, we have developed a theory of the rotation of Io modeled as a rigid ellipsoidal mantle containing a cavity filled with an inviscid fluid of constant uniform density and vorticity. Poincaré (1910) provided the foundations of a conservative model of the problem. Hence our previous work can be naturally extended to include its effects. This is the aim of the present contribution. Hamiltonian formulation The differential equations describing the motion of a rigid mantle with a traxial ellipsoidal cavity filled with inviscid fluid of constant uniform density and vorticity, have been derived by Hough (1895), and Poincaré (1910). More recently, Touma and Wisdom (001) have derived a Hamiltonian formulation on the base of Poincaré s formalism. We will adopt this formalism in what follows. The components (v 1,v,v 3 ) of the velocity field at the location x i inside the liquid core, in the frame of the principal axes of inertia of the mantle, are assumed to be: v 1 =ω + (a/c)ν x 3 ω 3 + (a/b)ν 3 x v =ω 3 + (b/a)ν 3 x 1 ω 1 + (b/c)ν 1 x 3 v 3 =ω 1 + (c/b)ν 1 x ω + (c/a)ν x 1, (1) where (a b c) are the principal elliptical radii of the cavity and where (ω 1,ω,ω 3 )are the components of the angular velocity of the mantle with respect to an inertial frame. The vector of coordinates (ν 1,ν,ν 3 ) specify the velocity field of the core with respect to the moving mantle. For a spherical cavity they reduce to the components of the angular velocity of the core, moving as a rigid body. Both vectors are expressed in the frame of the principal axes of inertia of the mantle. This velocity field (1) satisfy Helmholtz s equation. The angular momentum N c of the core is obtained by integration of ( x v) ρdw: ( ( ) c b c N c = Mc 5 + Mc 5 + Mc 5 ) b ν 1 + ω 1 b + ( c ) a ν + ω ( b a ν 3 + ω 3 a + ) f 1 c ν 1 + ω 1 ( a ) ac ν + ω c f ( a ) a + b ν 3 + ω 3 b f 3, () where f i are the unit vectors along the principal moments of inertia and M c is the mass of the core. Considering that the moments of inertia of the core are: A c = Mc 5 (b + c ); B c = Mc 5 (a + c ); C c = Mc 5 (a + b ), (3) and defining the following quantities: we obtain: D 1 = Mc 5 bc; D = Mc 5 ac; D 3 = Mc ab, (4) 5 N c =A c ω 1 + D 1 ν 1 f 1 +B c ω + D ν f +C c ω 3 + D 3 ν 3 f 3. (5)
3 The rotation of Io with a liquid core 3 The total angular momentum of the system N = N m + N c is then: N =Aω 1 + D 1 ν 1 f 1 +Bω + D ν f +Cω 3 + D 3 ν 3 f 3, (6) where (A B C) are the moments of inertia of the total system. Consider the kinetic energy of the system: T = Aω1 + Bω + Cω 3 + Ac ν1 + Bc ν + Cc ν3 + D 1 ω 1 ν 1 + D ω ν + D 3 ω 3 ν 3, (7) and its derivatives: T T T = Aω 1 + D 1 ν 1 ; = Bω + D ν ; = Cω 3 + D 3 ν 3. (8) ω 1 ω ω 3 We check that ( T/ ω i ) = N i, the components of the total angular momentum. On the other hand the derivatives: T = D 1 ω 1 + A c T ν 1 ; = D ω + B c T ν ; = D 3 ω 3 + C c ν 3, (9) ν 1 ν ν 3 are not the components of the angular momentum of the core but are close to it for a cavity close to spherical. Indeed, if we note ( T/ ν i ) by Ni c,wehave: N1 c N 1 c = (Ac D 1 )(ω 1 ν 1 ) = Mc 5 (c b) (ω 1 ν 1 ). (10) Similar equations hold for the second and third components. The difference is thus of the second order in the departure from sphericity. With these notations, The Poincaré-Hough s equations of motion, for the system mantlecore in the absence of external torque, are (see Touma and Wisdom 001): d N = N dt N H; d N c = N c dt N c H (11) The function H is the kinetic energy expressed in terms of the components of the vectors N and N c : H = 1 α (Ac N1 + A(N 1 c ) D 1 N 1 N1 c ) + 1 β (Bc N + B(N c ) D N N c ) + 1 γ (Cc N3 + C(N 3 c ) D 3 N 3 N3 c ), (1) where α = AA c D1,β = BBc D,andγ = CCc D3. 3 Additional approximations In order to define a simplified approximation, we will assume that the ellipsoid of inertia of the core and of the mantle are aligned and proportional, i.e A c = δ A, B c = δb, and C c = δc. Actually it may be more reallistic to assume that the core is axi-symmetric, as in (Touma and Wisdom 001), but they assume also that the mantle is axi-symmetric which is not compatible with the spin-orbit resonance assumption. In any case we have no indication concerning the flattening of such an axi-symmetric shape and this parameter, the flattening, is probably more significant than the (small) departure from axi-symmetry. Because of the lack of knowledge of the physical parameters, our aim cannot be to provide a simulation of
4 4 J. Henrard the true motion but to explore the possible effects of a liquid core. A simplified model will do that. We will also linearize the Hamiltonian (1) with respect to the small parameters: ɛ 1 = C A B ; ɛ = B A C C. (13) With these notations, we have A = C(1 ɛ 1 ɛ ), B = C(1 ɛ 1 + ɛ ), D 1 = A c, D = B c, D 3 = C c, and the Hamiltonian reads: H = 3 ( ) / Ni N i Ni c + (Ni c C(1 δ) ) δ i=1 i=1 i=1 ɛ 1 + (N C(1 δ) 1 + N ) (N 1N1 c + N N c ) + ((N 1 c ) + (N c ) )/δ ɛ + (N C(1 δ) 1 N ) (N 1N1 c N N c ) + ((N 1 c ) (N c ) )/δ. (14) 4 Andoyer variables Let us introduce the two sets of Andoyer s variables (Andoyer 196), (l, g, h, L, G, H) and (l c, g c, h c, L c, G c, H c ). The angles (h, K, g) are the Euler angles of the vector n, node of the equatorial plane over the plane perpendicular to the angular momentun N; the angles (J,l) position the axis of least inertia ( f 1 ) with respect to n (see Fig. 1 and Deprit 1967). Correspondingly the angles (h c, K c, g c ) are the Euler angles of the vector n c, and (J c,l c ) position the axis of least inertia with respect to n c. The variables are (h, g,l) and (h c, g c,l c ) and the corresponding momenta (H = N cos K, G = N, L = N cos J) and (H c = N c cos K c, G c = N c, L c = N c cos J c ). Expressed in Andoyer variables the components of N and N c are: N 1 = G L sin l, N1 c = G c L c sin l c N = G L cos l, N c = G c L c cos l c (15) N 3 = L, N3 c = L c. In order to take into account the minus sign in the differential equation for N c (11), we change the sign of the G c L c and of L c in (15) and obtain: H = 1 G + G c C(1 δ) /δ + LL c + (G L )(G c L c ) cos(l l c) + ɛ 1 C(1 δ) ɛ C(1 δ) (G L ) + (G c Lc )/δ + (G L )(G c L c ) cos(l l c) (G L ) cos l + (G c L c )/δ cos l c + (G L )(G c L c ) cos(l + l c). (16) Let us introduce the canonical change of variables:
5 The rotation of Io with a liquid core 5 Fig. 1 The Andoyer s variables (h, K ) locate the plane perpendicular to the angular momentum with respect to the inertial frame ( e 1, e, e 3 ) (the vector n 1 being along the ascending node of this plane) and the angles (g, J,l)locate the first vector of the body frame ( f 1, f, f 3 ) (the vector n being along the ascending node of the equatorial plane on the plane perpendicular to the angular momentum. The variables (J c,l c ) locate the pseudo angular momentum of the core with respect to the mantle body frame p = l + g + h, P = G/N 1, G = (1 + P)N q = l, Q = (G L)/N, L = (1 + P Q)N r = h, R = (G H)/N, H = (1 + P R)N p c = l c + g c + h c, P c = G c /N δ, G c = (δ + P c )N q c = l c, Q c = (G c + L c )/N, L c = (δ + P c Q c )N r c = h c, R c = (G c H c )/N, H c = (δ + P c R c )N, (17) of multiplier 1/N. The constant N (resp. δn) is the norm G (resp. G c ) of the angular momentum (resp. of the pseudo-angular momentum of the core). The Hamiltonian becomes: ω H = P + P (1 + P Q)(δ + P c Q c ) + (Pc (1 δ) + P c)/δ + ωɛ 1 (1 δ) ωɛ (1 δ) + QQ c ( + P Q)(δ + P c Q c ) cos (q + q c ) ( + P Q)Q + (δ + P c Q c )Q c /δ + QQ c ( + P Q)(δ + P c Q c ) cos (q + q c ) ( + P Q)Q cos q + (δ + P c Q c )(Q c /δ) cos q c + QQ c ( + P Q)(δ + P c Q c ) cos (q q c ). (18) We have used the fact that N/C = ω, the norm of the angular velocity of the full body.
6 6 J. Henrard 5 Librations around the equilibrium Let us introduce the cartesian coordinates x = Q sin q, y = Q cos q and x c = Qc sin q c, y c = Q c cos q c, and let us retain only the quadratic terms. We set P = P c = 0, and obtain: H = ω (1 δ) + ωɛ 1 (1 δ) + ωɛ (1 δ) δ(x + y ) + (xc + y c ) + δ(yy c xx c ) (x + y ) + (xc + y c ) + δ(yy c xx c ) (x y ) + (xc y c ) δ(yy c + xx c ). The two degrees of freedom (x, y) and (x c, y c ) are not independent of each other. In order to separate them, we use the canonical untangling transformation (see Henrard and Lemaître 005): (19) x = ξ 1 cos β ξ sin β, y = ρη 1 cos β + η sin β x c = ρ ξ 1 sin β + ξ cos β, y c = η 1 sin β + η cos β (0) The multiplier of this canonical transformation is ρ cos β. Let us write the Hamiltonian (19), as: ω a11 x + a 1 xx c + a xc (1 δ) + b 11y + b 1 yy c + b yc, (1) with: a 11 = δ + ɛ 1 + ɛ, a 1 = δ(1 + ɛ 1 + ɛ ), a = 1 + ɛ 1 + ɛ b 11 = δ + ɛ 1 ɛ, b 1 = δ(1 + ɛ 1 ɛ ), b = 1 + ɛ 1 ɛ. () After transformation, the cross product terms are: ξ 1 ξ a 1 ρ (a 11 + ρ a ) sin β+η 1 η b 1 ρ + (b 11 ρ + b ) sin β. (3) They disappears if sin β is defined by: a 1 ρ sin β = (a 11 + a ρ ) = b 1 ρ (b + b 11 ρ ) = δ 1 + δ The above equalities are compatible if: (1 δ) 1 (1 + δ) ɛ 1. (4) ρ = a 11b 1 + a 1 b = 1 ɛ. (5) a 1 b 11 + a b 1 After transformation the Hamiltonian becomes: H = ω { ɛ1 1 δ + ɛ ξ1 1 δ + ɛ1 1 δ ɛ } η1 1 δ + ω { 1 + ɛ 1 1 δ δɛ ξ 1 δ ɛ 1 1 δ + δɛ } η. (6) 1 δ A change of scale (ξ i ξ i α i )and(η i η i /α i ), with α1 4 = ɛ 1 ɛ ; α 4 = 1, (7) ɛ 1 + ɛ
7 The rotation of Io with a liquid core 7 leads to: H = ω 1 (ξ 1 + η 1 ) + ω (ξ + η ), (8) where the frequencies ω i are: ω 1 ω = ɛ1 ɛ (1 δ) ; ω ω = 1 + ɛ 1 1 δ δɛ 1 δ. (9) When ɛ = 0 (i.e. when the body and the core are axisymmetric), we recover the result of Touma and Wisdom (001). Notice that when δ 0, i.e. when the core vanishes, the angle β 0, hence the degree of freedom (ξ 1,η 1 ) corresponds to the wobble of the full body and the degree of freedom (ξ,η ) corresponds to the motion of the core, which is irrelevant as the core has disappeared. In that case and when ɛ = 0wehaveω 1 ɛ 1 ω and ω (1 ɛ 1 )ω. 6 Perturbation of Io by Jupiter The gravitational pull of Jupiter on the rotation of the mantle of a satellite is modelized by (see Henrard 005c): 3 V = C m n d0 d0 δ1 1 + δ s (x + y ) + δ (x y ), (30) d d where d is the distance between the planet and the center of mass of the satellite, d 0 the mean value of this distance, C m = (1 δ)c the largest moment of inertia of the mantle, δ 1 = (3/)(n /n) (C A B)/C and δ = (3/)(n /n) (B A)/C. The constant n = GM J /d0 3 is very close to the orbital frequency n = ω. Dividing the potential by ωc as we did for the Hamiltonian of the free rotation, we obtain: 3 Cm d0 δ1 (x + y ) + δ (x y ). (31) V ωc = ω C d 1 + δ s d0 d The components (x, y) of the unit vector pointing to Jupiter and the distance d are computed on the basis of the synthetic theory of the motion of Io developed by Lainey (Lainey et al. 004a,b)asinHenrard (005c) (Tables 1, ). 7 Cassini states and librations As for most regular satellites in the Solar System, the spin period of Io is equal to its orbital period. In such a case, as Cassini pointed out, the node of the orbit and the node of the equator have, in the mean, the same period. There is thus a resonance between the angles λ 1 and p Table 1 Constants used in the theory (Sohl et al. 00; Zhang 003) J = C = C /MR J = R J /d 0 = J (Jup) = C /M J R J = 0.6 δ 1 = δ =
8 8 J. Henrard Table The angular variables in the description of Io orbit (Lainey 00) Angle Frequency (in rd/day) Phase (in rd) λ 1 n 1 = λ n = λ 3 n 3 = λ 4 n 4 = ϖ 1 ϖ 1 = ϖ ϖ = ϖ 3 ϖ 3 = ϖ 4 ϖ 4 = = = = = l = Laplace lib. l = λ Sun n S = Table 3 Obliquity and frequencies (in unit of the orbital frequency) of the free librations R c K (radian) Latitude (ω v ) Wobble (ω w ) Core(ω c ) δ = δ = δ = δ = δ = on one hand, and the angles 1 and r. To take this into account,we define the two resonance angular variables σ and ρ and introduce a new set of canonical variables: σ = p L 1 + π, P ρ = r + 1, R. (3) The resonant angle σ is defined as the difference between the angle p and the mean longitude of Jupiter as seen from Io. Similarly ρ is defined as the difference between the angle h = r and the ascending node of the orbit of Jupiter. The variable L 1 in Lainey s theory is the mean longitude of Io as seen from Jupiter. The mean longitude of Jupiter as seen from Io is thus L 1 + π. The transformation being time dependent we ought to add its remainder function n 1 P + 1 R to the Hamiltonian. As in (Henrard 005a, b) we isolate K 0 the part of the Hamiltonian which do not depend on the angular variables describing the perturbed motion of Io. It corresponds to a problem for which Io is on a circular Keplerian orbit around Jupiter. Equilibria (the Cassini equilibria) are found for (σ = ρ = x = y = x c = y c = 0) and P = P and R = R such that the derivatives of K 0 with respect to P and R vanish. The values of P is extremely small (10 10 ) and the value of K = cos 1 (1 R ) for the main Cassini equilibrium is given in Table 3 for various values of δ. The value of the mean radius of the core (considering that core and mantle have the same density) is also given.
9 The rotation of Io with a liquid core 9 Expanding K 0 in powers of the small canonical departure from the main Cassini state: ξ σ = σ, ξ ρ = (R R ) sin r, η σ = P P η ρ = ηρ + (R R ) cos r ξ w = x, η w = y ξ c = x c, η c = y c, (33) where η ρ = R, and keeping only the quadratic terms, we define a linear (and thus integrable) system for which we introduce angle-action variables as in (Henrard 005c) for the first two degrees of freedom and as in (0) forthelasttwo: ξ σ = UU sin u γ VV sin v η σ = (1 αγ) U/U cos u α U/U cos u ξ ρ = α UU sin u + (1 αγ) VV sin v η ρ = γ U/U cos u + U/U cos u ξ w = cos β WW sin w sin β ZZ sin z η w =cos β W/W cos w sin β Z/Z cos z/ cos β ξ c = sin β WW sin w + cos β ZZ sin z/ cos β η c = sin β W/W cos w + cos β Z/Z cos z. (34) The parameters α, β, γ and the scaling factors U, V, W, Z are function of δ and adjusted in such a way that, after transformation, the quadratic Hamiltonian reduces to N = ω u U + ω v V + ω w W + ω z Z. (35) The frequencies of the free libration in longitude (ω u ), of the free libration in latitude (ω v ), of the wobble (ω w ) and of the free libration of the core (ω c ) are functions of δ and tabulated in Table 3. The frequency of the free libration in longitude (ω u ) does not depend on the value of δ, and is equal to Perturbation theory We impose on the perturbation function P = H N the transformations described in the previous sections and expand it, up to the fourth degree, in powers of the amplitudes U, V, W, Z of the oscillations around the Cassini s equilibrium. A fourth order perturbation theory by Lie transform is now straightforward except that, because of the introduction of an extra frequency ω c, close to the orbital frequency, several harmonics are close to be resonant. An attempt at averaging all the quasi periodic contributions is bound to fail. Instead we decided to eliminate from the Hamiltonian only the terms the frequency of which were much larger than their amplitude. This partial averaging left us with a Hamiltonian still containing quasi periodic terms such as the ones listed in Table 3. Thefirsttwoare not really dangerous. The amplitudes are small and the frequencies do not go to zero. They can be eliminated in a further transformation with a rather small effect on the solution. The other terms (and some smaller ones not listed) are of the second order in the amplitudes U, V, W, Z, which means that the origin ( U = V = W = Z = 0) is still an equilibrium of the partially averaged Hamiltonian. Nevertheless, the existence of those terms, the frequencies of which goes to zero for some values of the parameter δ
10 10 J. Henrard Table 4 Frequencies (in unit of 10 ) of some resonant terms Angular variable Amp. δ = 10 3 δ = 0.1 δ = 0. δ = 0.3 δ = 0.4 λ 1 1 z λ 1 z λ v z λ ϖ 1 + v z λ 1 λ + 1 v w λ 1 λ + w z Fig. Two of the resonant frequencies: above the frequency of the angular variable λ v z, below the frequency of λ ϖ v z (as shown in Table 4) is an indication that the averaging process may be flawed and that the vicinity of the Cassini equilibrium may be chaotic for values of the parameter δ around δ = 0.1, where the frequency of the perturbing term in cos(λ ϖ 1 + v z) is close to zero, or around δ = 0., where it is the frequency of the rather large perturbing term cos(λ v z) which vanishes, or around δ = 0.4, where it is the frequency of cos(λ 1 λ + w z) (Fig. ). Even if, for particular values of the parameter δ, the possibility of the existence of resonances and chaotic motions in the vicinity of the Cassini equilibrium, may cast a doubt on the validity of the perturbative approach we have taken, it may be useful to report the effects of the partial averaging. We list in Table 5 the non-linear corrections to the frequencies of the free libration (in unit of the orbital frequency). We list also in Tables 6 and 7 the expressions of P 1 and P, the first two component of the unit vector perpendicular to Jupiter equator (in the frame of the principal moment of inertia of Io), and, in Tables 8 and 9, the principal terms of the free librations of these two components. Tables 6 and 7 show that in general the presence of a liquid core does not modify considerably the behaviour of the spin of Io except for some particular values of the parameter δ which are such that some combination of frequencies goes to zero implying a second order resonance. Close to these particular values (for instance for δ = 0. and δ = 0.4) new combination of angular variables appear in the Fourier expansions.
11 The rotation of Io with a liquid core 11 Table 5 Non-linear corrections to the frequencies of the free librations Longitude (ω u ) Latitude (ω v ) Wobble (ω w ) Core(ω c ) δ = δ = δ = δ = δ = Table 6 Expressions (truncated at 10 )of10 4 P 1 for different values of δ δ = 0 δ = 0.1 δ = 0. δ = 0.3 δ = 0.4 sin(λ 1 1 ) sin(λ 1 ) sin(λ 1 3 ) sin(λ 1 4 ) sin(λ 1 λ + ) sin(3λ 1 4λ + ) cos(λ 1 λ + ϖ 3 ) Table 7 Expressions (truncated at 10 )of10 4 P for different values of δ δ = 0 δ = 0.1 δ = 0. δ = 0.3 δ = 0.4 cos(λ 1 1 ) cos(λ 1 ) cos(λ 1 3 ) cos(λ 1 4 ) cos(λ 1 λ + ) sin(3λ 1 4λ + ) Table 8 Free librations of P 1 (truncated at ) δ = δ = 0.1 δ = 0. δ = 0.3 δ = 0.4 V sin(λ1 1 + v) W sin w Z sin z Table 9 Free librations of P (truncated at 10 1 ) δ = 0 δ = 0.1 δ = 0. δ = 0.3 δ = 0.4 V cos(λ1 1 + v) W cos w Z cos z
12 1 J. Henrard 9Conclusion This preliminary investigation shows that the consideration of a liquid core introduces in the problem of the rotation of Io a new frequency (the frequency of the spin of the core) close to the orbital frequency. The orbital motion and the rotation of Io is already characterized by long period contributions. The introduction of a new frequency close to the main frequency multiplies the possibility of commensurability and resonance. We have detected three possible such resonances for values of δ (ratio of the core inertia to the full inertia of the body) close to 0.1, 0. and 0.4. These are second order resonances which do not displace the Cassini state but might generate a chaotic vicinity around it. Of course, this is a very preliminary exploration of the effect of a liquid core on the rotation of Io. Our assumptions on the shape of the cavity containing the core may not be optimal but we do not think this should alter significantly the results. We plan to test an axi-symmetric model of the cavity in a further publication. What could be more important is our neglect of the effects of the possible dissipation occuring at core-mantle boundary. Such effects will also be investigated in further publications. References Andoyer, H.: Mécanique Céleste. Gauthier-Villars, Paris (196) Deprit, A.: Free rotation of an rigid body studied in the phase space. Am. J. Phys. 35, (1967) Goldreich, P., Peale, S.J.: Spin-orbit coupling in the solar system. Astron. J. 71, (1966) Henrard, J.: The rotation of Europa. Celest. Mech. Dyn. Astron. 91, (005a) Henrard, J.: Complements to the rotation of Europa. Celest. Mech. Dyn. Astron. 93, (005b) Henrard, J.: The rotation of Io. Icarus. 178, (005c) Henrard, J., Lemaître, A.: The untangling transformation. Astron. J. 130(5), (005) Hough, SS.: The oscillations of a rotating ellipsoidal shell containing fluid. Philos. Trans. R. Soc. London A 186, (1895) Lainey, V. Théorie Dynamique des Satellites Galiléens. PhD dissertation, Observatoire de Paris (00) Lainey, V., Arlot, J.E., Vienne, A.: New accurate ephemerides for the Galilean satellites of Jupiter: I Numerical integration of elaborated equation of motion. A&A 40, (004a) Lainey, V., Arlot, J.E., Vienne, A.: New accurate ephemerides for the Galilean satellites of Jupiter: II Fitting the observations. A&A 47, (004b) Poincaré, H.: Sur la précession des corps déformables. Bull. Astron. 7, (1910) Sohl, F., Spohn, T., Breuer, D., Nagel, K.: Implication from Galileo observations on the interior structure and chemistry of the Galilean satellites. Icarus 157, (00) Touma, J., Wisdom, J.: Nonlinear core-mantle coupling. AJ 1, (001) Zhang, H.: Internal structure models and dynamical parameters of the Galilean satellites. Celest. Mech. Dyn. Astron. 87, (003)
THE LUNAR LIBRATION: COMPARISONS BETWEEN VARIOUS MODELS- A MODEL FITTED TO LLR OBSERVATIONS
THE LUNAR LIBRATION: COMPARISONS BETWEEN VARIOUS MODELS- A MODEL FITTED TO LLR OBSERVATIONS J. CHAPRONT, G. FRANCOU SYRTE - Observatoire de Paris - UMR 8630/CNRS 61, avenue de l Observatoire 75014 Paris
More informationANALYTICAL MODELING OF THE RIGID INTERNAL MOTIONS OF A THREE-LAYER CELESTIAL BODY THROUGH HAMILTON S PRINCIPLE
ANALYTICAL MODELING OF THE RIGID INTERNAL MOTIONS OF A THREE-LAYER CELESTIAL BODY THROUGH HAMILTON S PRINCIPLE A. ESCAPA Department of Applied Mathematics, University of Alicante PO Box 99, E-03080 Alicante,
More informationarxiv: v1 [physics.geo-ph] 16 May 2012
Obliquity of the Galilean satellites: The influence of a global internal liquid layer arxiv:125.3628v1 [physics.geo-ph] 16 May 212 Abstract Rose-Marie Baland a,b, Marie Yseboodt b, Tim Van Hoolst b a Université
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 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 informationDYNAMICS OF THE LUNAR SPIN AXIS
The Astronomical Journal, 131:1864 1871, 2006 March # 2006. The American Astronomical Society. All rights reserved. Printed in U.S.A. DYNAMICS OF THE LUNAR SPIN AXIS Jack Wisdom Massachusetts Institute
More informationarxiv: v4 [astro-ph] 12 Mar 2008
Celestial Mechanics and Dynamical Astronomy manuscript No. (will be inserted by the editor) Benoît Noyelles arxiv:0709.2293v4 [astro-ph] 12 Mar 2008 Titan s rotational state The effects of a forced free
More informationProbing planetary interiors by spacecraft orbital observations
Probing planetary interiors by spacecraft orbital observations Alexander Stark, Jürgen Oberst, Frank Preusker, Klaus Gwinner, Gregor Steinbrügge, Hauke Hussmann Funded by Deutsche Forschungsgemeinschaft
More informationAstronomy 6570 Physics of the Planets. Precession: Free and Forced
Astronomy 6570 Physics of the Planets Precession: Free and Forced Planetary Precession We have seen above how information concerning the distribution of density within a planet (in particular, the polar
More informationPhysics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top
Physics 106a, Caltech 4 December, 2018 Lecture 18: Examples on Rigid Body Dynamics I go through a number of examples illustrating the methods of solving rigid body dynamics. In most cases, the problem
More informationEvolution of Mercury s obliquity
Icarus 181 2006) 327 337 www.elsevier.com/locate/icarus Evolution of Mercury s obliquity Marie Yseboodt, Jean-Luc Margot Department of Astronomy, Cornell University, Ithaca, NY 14853, USA Received 24 June
More informationSection B. Electromagnetism
Prelims EM Spring 2014 1 Section B. Electromagnetism Problem 0, Page 1. An infinite cylinder of radius R oriented parallel to the z-axis has uniform magnetization parallel to the x-axis, M = m 0ˆx. Calculate
More informationHAMILTONIAN STABILITY OF SPIN ORBIT RESONANCES IN CELESTIAL MECHANICS. 1. Introduction
HAMILTONIAN STABILITY OF SPIN ORBIT RESONANCES IN CELESTIAL MECHANICS ALESSANDRA CELLETTI 1 and LUIGI CHIERCHIA 2 1 Dipartimento di Matematica, Università di Roma Tor Vergata, Via della Ricerca Scientifica,
More informationarxiv: v2 [astro-ph.ep] 28 Sep 2017
Noname manuscript No. (will be inserted by the editor) Tilting Styx and Nix but not Uranus with a Spin-Precession-Mean-motion resonance Alice C. Quillen 1, Yuan-Yuan Chen 1,2, Benoît Noyelles 3, & Santiago
More information2 so that the (time-dependent) torque, as measured in space-xed coordinates, is ~N = ma2! 2 (^x cos!t +^y sin!t) 10 = a p 5 2 ^z ~ F 2 ~F 1 = a p 5^z
1 PHYS 321 Homework Assignment #9 Due: Friday, 22 November 2002 (4 probs) 1. Problem 7-14, B&O p. 280. The situation is as shown above. We take the body-xed axes ^x 0 and ^y 0 to lie in the plane of the
More informationTides and Lagrange Points
Ast111, Lecture 3a Tides and Lagrange Points Arial view of Tidal surge in Alaska (S. Sharnoff) Tides Tidal disruption Lagrange Points Tadpole Orbits and Trojans Tidal Bulges Tides Tidal Force The force
More informationResonance Capture. Alice Quillen. University of Rochester. Isaac Newton Institute Dec 2009
Resonance Capture Alice Quillen University of Rochester Isaac Newton Institute Dec 2009 This Talk Resonance Capture extension to nonadiabatic regime Astrophysical settings: Dust spiraling inward via radiation
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 informationResonance and chaos ASTRONOMY AND ASTROPHYSICS. II. Exterior resonances and asymmetric libration
Astron. Astrophys. 328, 399 48 (1997) ASTRONOMY AND ASTROPHYSICS Resonance and chaos II. Exterior resonances and asymmetric libration O.C. Winter 1,2 and C.D. Murray 2 1 Grupo de Dinâmica Orbital e Planetologia,
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 informationPh1a: Solution to the Final Exam Alejandro Jenkins, Fall 2004
Ph1a: Solution to the Final Exam Alejandro Jenkins, Fall 2004 Problem 1 (10 points) - The Delivery A crate of mass M, which contains an expensive piece of scientific equipment, is being delivered to Caltech.
More informationInformation on internal structure from shape, gravity field and rotation
Information on internal structure from shape, gravity field and rotation Seismological information is available only for the Earth and in limited amounts for the Moon. Various geodetic data put constraints
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 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 informationAstronomy 6570 Physics of the Planets
Astronomy 6570 Physics of the Planets Planetary Rotation, Figures, and Gravity Fields Topics to be covered: 1. Rotational distortion & oblateness 2. Gravity field of an oblate planet 3. Free & forced planetary
More informationSPACECRAFT DYNAMICS NEAR A BINARY ASTEROID. F. Gabern, W.S. Koon and J.E. Marsden
PROCEEDINGS OF THE FIFTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS June 16 19, 2004, Pomona, CA, USA pp. 1 10 SPACECRAFT DYNAMICS NEAR A BINARY ASTEROID F. Gabern, W.S. Koon
More informationEulerian equilibria of a rigid body in the three body problem
Proceedings of the nd WSEAS Int. Conference on Applied and Theoretical Mechanics, Venice, Italy, November 0-, 006 68 Eulerian equilibria of a rigid body in the three body problem J.A. VERA Universidad
More informationOn the Secular Evolution of Extrasolar Planetary Systems
On the Secular Evolution of Extrasolar Planetary Systems Marco Sansottera [a] [a] Namur Center for Complex Systems (naxys) Based on a research work in collaboration with Anne-Sophie Libert [a] JISD 2012
More informationTorque and Rotation Lecture 7
Torque and Rotation Lecture 7 ˆ In this lecture we finally move beyond a simple particle in our mechanical analysis of motion. ˆ Now we consider the so-called rigid body. Essentially, a particle with extension
More informationDynamics of the conservative and dissipative. spin orbit problem
Dynamics of the conservative and dissipative spin orbit problem Alessandra Celletti Dipartimento di Matematica, Università di Roma Tor Vergata, Via della Ricerca Scientifica, I-0033 Roma (Italy), celletti@mat.uniroma2.it
More informationEART162: PLANETARY INTERIORS
EART162: PLANETARY INTERIORS Francis Nimmo Last Week Applications of fluid dynamics to geophysical problems Navier-Stokes equation describes fluid flow: Convection requires solving the coupled equations
More informationRotational Mechanics and Relativity --- Summary sheet 1
Rotational Mechanics and Relativity --- Summary sheet 1 Centre of Mass 1 1 For discrete masses: R m r For continuous bodies: R dm i i M M r body i Static equilibrium: the two conditions for a body in static
More informationarxiv: v1 [astro-ph.ep] 12 Dec 2016
Astronomy & Astrophysics manuscript no. manuscript_arxiv c ESO 2018 December 20, 2018 Frequency-dependent tidal dissipation in a viscoelastic Saturnian core and expansion of Mimas semi-major axis arxiv:1612.03664v1
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 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 informationGeodesy and interior structure of Mercury
Geodesy and interior structure of Mercury Tim Van Hoolst and the ROB planets team Royal Observatory of Belgium February 26, 2007 Outline 1 Tidal potential Potential Response 2 Interior structure models
More informationIf the symmetry axes of a uniform symmetric body coincide with the coordinate axes, the products of inertia (Ixy etc.
Prof. O. B. Wright, Autumn 007 Mechanics Lecture 9 More on rigid bodies, coupled vibrations Principal axes of the inertia tensor If the symmetry axes of a uniform symmetric body coincide with the coordinate
More informationarxiv:astro-ph/ v1 17 Jul 2003
Astronomy & Astrophysics manuscript no. RB-letter September 25, 217 (DOI: will be inserted by hand later) Mercury s spin-orbit model and signature of its dynamical parameters arxiv:astro-ph/37328v1 17
More informationTranscendental cases in stability problem. Hamiltonian systems
of Hamiltonian systems Boris S. Bardin Moscow Aviation Institute (Technical University) Faculty of Applied Mathematics and Physics Department of Theoretical Mechanics Hamiltonian Dynamics and Celestial
More informationModels for locations and motions in the solar system
Journal of Physics and Astronomy Research Vol. 2(1), pp. 054-058, March, 2015. www.premierpublishers.org, ISSN: 2123-503x JPAR Review Models for locations and motions in the solar system Lena Strömberg
More informationPeriodic Orbits in Rotating Second Degree and Order Gravity Fields
Chin. J. Astron. Astrophys. Vol. 8 (28), No. 1, 18 118 (http://www.chjaa.org) Chinese Journal of Astronomy and Astrophysics Periodic Orbits in Rotating Second Degree and Order Gravity Fields Wei-Duo Hu
More informationAST111 PROBLEM SET 2 SOLUTIONS. RA=02h23m35.65s, DEC=+25d18m42.3s (Epoch J2000).
AST111 PROBLEM SET 2 SOLUTIONS Home work problems 1. Angles on the sky and asteroid motion An asteroid is observed at two different times. The asteroid is located at RA=02h23m35.65s, DEC=+25d18m42.3s (Epoch
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 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 informationPC 1141 : AY 2012 /13
NUS Physics Society Past Year Paper Solutions PC 1141 : AY 2012 /13 Compiled by: NUS Physics Society Past Year Solution Team Yeo Zhen Yuan Ryan Goh Published on: November 17, 2015 1. An egg of mass 0.050
More informationarxiv: v1 [nucl-th] 8 Sep 2011
Tidal Waves a non-adiabatic microscopic description of the yrast states in near-spherical nuclei S. Frauendorf, Y. Gu, and J. Sun Department of Physics, University of Notre Dame, Notre Dame, IN 6556, USA
More informationPDFlib PLOP: PDF Linearization, Optimization, Protection. Page inserted by evaluation version
PDFlib PLOP: PDF Linearization, Optimization, Protection Page inserted by evaluation version www.pdflib.com sales@pdflib.com Relative Equilibria for General Gravity Fields in the Sphere-Restricted Full
More informationMoment of inertia. Contents. 1 Introduction and simple cases. January 15, Introduction. 1.2 Examples
Moment of inertia January 15, 016 A systematic account is given of the concept and the properties of the moment of inertia. Contents 1 Introduction and simple cases 1 1.1 Introduction.............. 1 1.
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 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 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 informationAnalytical Solution of the Perturbed Oribt-Attitude Motion of a Charged Spacecraft in the Geomagnetic Field
(IJACSA) International Journal of Advanced Computer Science Applications, Vol. 4, No., Analytical Solution of the Perturbed Oribt-Attitude Motion of a Charged Spacecraft in the Geomagnetic Field Hani M.
More informationDetermination of the rotation of Mercury from satellite gravimetry
Mon. Not. R. Astron. Soc. 427, 468 482 (2012) doi:10.1111/j.1365-2966.2012.21919.x Determination of the rotation of Mercury from satellite gravimetry S. Cicalò and A. Milani Università degli Studi di Pisa,
More informationRigid bodies - general theory
Rigid bodies - general theory Kinetic Energy: based on FW-26 Consider a system on N particles with all their relative separations fixed: it has 3 translational and 3 rotational degrees of freedom. Motion
More informationPart 8: Rigid Body Dynamics
Document that contains homework problems. Comment out the solutions when printing off for students. Part 8: Rigid Body Dynamics Problem 1. Inertia review Find the moment of inertia for a thin uniform rod
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 informationLecture AC-1. Aircraft Dynamics. Copy right 2003 by Jon at h an H ow
Lecture AC-1 Aircraft Dynamics Copy right 23 by Jon at h an H ow 1 Spring 23 16.61 AC 1 2 Aircraft Dynamics First note that it is possible to develop a very good approximation of a key motion of an aircraft
More informationRotational Motion. Chapter 4. P. J. Grandinetti. Sep. 1, Chem P. J. Grandinetti (Chem. 4300) Rotational Motion Sep.
Rotational Motion Chapter 4 P. J. Grandinetti Chem. 4300 Sep. 1, 2017 P. J. Grandinetti (Chem. 4300) Rotational Motion Sep. 1, 2017 1 / 76 Angular Momentum The angular momentum of a particle with respect
More informationarxiv: v3 [astro-ph.ep] 20 Aug 2015
Astronomy & Astrophysics manuscript no. nix hydra c ESO 8 August 9, 8 Letter to the Editor Spin-orbit coupling and chaotic rotation for circumbinary bodies Application to the small satellites of the Pluto-Charon
More informationarxiv: v1 [astro-ph.ep] 30 Aug 2018
Astronomy & Astrophysics manuscript no. article c ESO 218 August 31, 218 Long-term orbital and rotational motions of Ceres and Vesta T. Vaillant, J. Laskar, N. Rambaux, and M. Gastineau ASD/IMCCE, Observatoire
More informationAS3010: Introduction to Space Technology
AS3010: Introduction to Space Technology L E C T U R E S 8-9 Part B, Lectures 8-9 23 March, 2017 C O N T E N T S In this lecture, we will look at factors that cause an orbit to change over time orbital
More informationORBITS WRITTEN Q.E. (June 2012) Each of the five problems is valued at 20 points. (Total for exam: 100 points)
ORBITS WRITTEN Q.E. (June 2012) Each of the five problems is valued at 20 points. (Total for exam: 100 points) PROBLEM 1 A) Summarize the content of the three Kepler s Laws. B) Derive any two of the Kepler
More informationarxiv: v1 [astro-ph.ep] 11 Aug 2016
Noname manuscript No. (will be inserted by the editor) Secular and tidal evolution of circumbinary systems. Alexandre C. M. Correia Gwenaël Boué Jacques Laskar arxiv:168.3484v1 [astro-ph.ep] 11 Aug 216
More informationexpression that describes these corrections with the accuracy of the order of 4. frame usually connected with extragalactic objects.
RUSSIAN JOURNAL OF EARTH SCIENCES, English Translation, VOL, NO, DECEMBER 998 Russian Edition: JULY 998 On the eects of the inertia ellipsoid triaxiality in the theory of nutation S. M. Molodensky Joint
More informationOn Sun-Synchronous Orbits and Associated Constellations
On Sun-Synchronous Orbits and Associated Constellations Daniele Mortari, Matthew P. Wilkins, and Christian Bruccoleri Department of Aerospace Engineering, Texas A&M University, College Station, TX 77843,
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 informationLecture 23: Jupiter. Solar System. Jupiter s Orbit. The semi-major axis of Jupiter s orbit is a = 5.2 AU
Lecture 23: Jupiter Solar System Jupiter s Orbit The semi-major axis of Jupiter s orbit is a = 5.2 AU Jupiter Sun a Kepler s third law relates the semi-major axis to the orbital period 1 Jupiter s Orbit
More informationA REGION VOID OF IRREGULAR SATELLITES AROUND JUPITER
The Astronomical Journal, 136:909 918, 2008 September c 2008. The American Astronomical Society. All rights reserved. Printed in the U.S.A. doi:10.1088/0004-6256/136/3/909 A REGION VOID OF IRREGULAR SATELLITES
More informationSimple and Physical Pendulums Challenge Problem Solutions
Simple and Physical Pendulums Challenge Problem Solutions Problem 1 Solutions: For this problem, the answers to parts a) through d) will rely on an analysis of the pendulum motion. There are two conventional
More informationM04M.1 Particles on a Line
Part I Mechanics M04M.1 Particles on a Line M04M.1 Particles on a Line Two elastic spherical particles with masses m and M (m M) are constrained to move along a straight line with an elastically reflecting
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 informationarxiv: v1 [astro-ph.ep] 23 Nov 2011
Draft version September 3, 18 Preprint typeset using L A TEX style emulateapj v. 11/1/9 PUMPING THE ECCENTRICITY OF EXOPLANETS BY TIDAL EFFECT Alexandre C.M. Correia Department of Physics, I3N, University
More informationOrbital Stability Regions for Hypothetical Natural Satellites
Orbital Stability Regions for Hypothetical Natural Satellites By Samantha RIEGER, 1) Daniel SCHEERES, 1) 1) Ann and H.J. Smead Department of Aerospace Engineering Sciences, University of Colorado Boulder,
More informationLong-term evolution of orbits about a precessing oblate planet. 2. The case of variable precession
Celestial Mech Dyn Astr 006 96:59 88 DOI 101007/s10569-006-9046-5 ORIGINAL ARTICLE Long-term evolution of orbits about a precessing oblate planet The case of variable precession Michael Efroimsky Received:
More informationClassical Mechanics. Luis Anchordoqui
1 Rigid Body Motion Inertia Tensor Rotational Kinetic Energy Principal Axes of Rotation Steiner s Theorem Euler s Equations for a Rigid Body Eulerian Angles Review of Fundamental Equations 2 Rigid body
More informationPrecise modeling of the precession-nutation of Mars
Astron. Astrophys. 345, 8 97 (999) ASTRONOMY AND ASTROPHYSICS Precise modeling of the precession-nutation of Mars S. Bouquillon and J. Souchay Observatoire de Paris/DANOF, URA 5 du CNRS, 6 Avenue de l
More informationThe in-plane Motion of a Geosynchronous Satellite under the Gravitational Attraction of the Sun, the Moon and the Oblate Earth
J Astrophys Astr (1990) 11, 1 10 The in-plane Motion of a Geosynchronous Satellite under the Gravitational Attraction of the Sun, the Moon and the Oblate Earth Κ Β Bhatnagar Department of Mathematics,
More informationSecular dynamics of extrasolar-systems
Secular dynamics of extrasolar-systems Marco Sansottera [1] [1] Università degli Studi di Milano [2] Namur Center for Complex Systems (naxys) Based on a research project in collaboration with A.-S. Libert
More informationSymplectic Correctors for Canonical Heliocentric N-Body Maps
Symplectic Correctors for Canonical Heliocentric N-Body Maps J. Wisdom Massachusetts Institute of Technology, Cambridge, MA 02139 wisdom@poincare.mit.edu Received ; accepted 2 ABSTRACT Symplectic correctors
More informationPerturbed Earth rotation. L. D. Akulenko, S. A. Kumakshev, and A. M. Shmatkov
Perturbed Earth rotation L. D. Akulenko, S. A. Kumakshev, and A. M. Shmatkov 1 Introduction Numerous astrometric studies are based on the dynamic theory of the Earth s rotation with respect to the center
More informationTidal Effects on Earth s Surface
Tidal Effects on Earth s Surface Tom Murphy February, 1 This treatment follows the conventions of F. D. Stacey s Physics of the Earth, and is largely an elaboration on this work. 1 Tidal Potential The
More informationFor a rigid body that is constrained to rotate about a fixed axis, the gravitational torque about the axis is
Experiment 14 The Physical Pendulum The period of oscillation of a physical pendulum is found to a high degree of accuracy by two methods: theory and experiment. The values are then compared. Theory For
More informationLecture 38: Equations of Rigid-Body Motion
Lecture 38: Equations of Rigid-Body Motion It s going to be easiest to find the equations of motion for the object in the body frame i.e., the frame where the axes are principal axes In general, we can
More informationPhys 214. Planets and Life
Phys 214. Planets and Life Dr. Cristina Buzea Department of Physics Room 259 E-mail: cristi@physics.queensu.ca (Please use PHYS214 in e-mail subject) Lecture 28. Search for life on jovian moons. March
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 informationRotational & Rigid-Body Mechanics. Lectures 3+4
Rotational & Rigid-Body Mechanics Lectures 3+4 Rotational Motion So far: point objects moving through a trajectory. Next: moving actual dimensional objects and rotating them. 2 Circular Motion - Definitions
More informationISIMA lectures on celestial mechanics. 3
ISIMA lectures on celestial mechanics. 3 Scott Tremaine, Institute for Advanced Study July 2014 1. The stability of planetary systems To understand the formation and evolution of exoplanet systems, we
More informationEffective resonant stability of Mercury
Mon. Not. R. Astron. Soc. 452, 4145 4152 (15) Printed October 22, 15 (MN LATEX style file v2.2) Effective resonant stability of Mercury M. Sansottera,C. Lhotka,A. Lemaître, Dipartimento di Matematica,
More informationThe Three Body Problem
The Three Body Problem Joakim Hirvonen Grützelius Karlstad University December 26, 2004 Department of Engineeringsciences, Physics and Mathematics 5p Examinator: Prof Jürgen Füchs Abstract The main topic
More informationThe rotation of Mimas
Astronomy & Astrophysics manuscript no. rotation mimas c ESO 2018 October 2, 2018 The rotation of Mimas B. Noyelles 1,2, Ö. Karatekin 3, and N. Rambaux 4,2 1 FUNDP-University of Namur - Department of Mathematics
More information" = Y(#,$) % R(r) = 1 4& % " = Y(#,$) % R(r) = Recitation Problems: Week 4. a. 5 B, b. 6. , Ne Mg + 15 P 2+ c. 23 V,
Recitation Problems: Week 4 1. Which of the following combinations of quantum numbers are allowed for an electron in a one-electron atom: n l m l m s 2 2 1! 3 1 0 -! 5 1 2! 4-1 0! 3 2 1 0 2 0 0 -! 7 2-2!
More informationOn the extension of the Laplace-Lagrange secular theory to order two in the masses for extrasolar systems
Celestial Mechanics and Dynamical Astronomy DOI 1007/s10569-013-9501-z On the extension of the Laplace-Lagrange secular theory to order two in the masses for extrasolar systems Anne-Sophie Libert Marco
More information1/30. Rigid Body Rotations. Dave Frank
. 1/3 Rigid Body Rotations Dave Frank A Point Particle and Fundamental Quantities z 2/3 m v ω r y x Angular Velocity v = dr dt = ω r Kinetic Energy K = 1 2 mv2 Momentum p = mv Rigid Bodies We treat a rigid
More informationFormation of the Solar System
Formation of the Solar System V.I. Ferronsky S.V. Ferronsky Formation of the Solar System A New Theory of the Creation and Decay of the Celestial Bodies 123 V.I. Ferronsky Water Problems Institute of
More informationAstronomy 241: Review Questions #2 Distributed: November 7, 2013
Astronomy 241: Review Questions #2 Distributed: November 7, 2013 Review the questions below, and be prepared to discuss them in class. For each question, list (a) the general topic, and (b) the key laws
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 informationMean-Motion Resonance and Formation of Kirkwood Gaps
Yan Wang Project 1 PHYS 527 October 13, 2008 Mean-Motion Resonance and Formation of Kirkwood Gaps Introduction A histogram of the number of asteroids versus their distance from the Sun shows some distinct
More informationResonancias orbitales
www.fisica.edu.uy/ gallardo Facultad de Ciencias Universidad de la República Uruguay Curso Dinámica Orbital 2018 Orbital motion Oscillating planetary orbits Orbital Resonances Commensurability between
More information5.1. Accelerated Coordinate Systems:
5.1. Accelerated Coordinate Systems: Recall: Uniformly moving reference frames (e.g. those considered at 'rest' or moving with constant velocity in a straight line) are called inertial reference frames.
More information