tides v4
DOCX · 76.7 KB
Open DOCX file
Short physics note by Phil dated 8.13.12, a version 4 outline of the facing-planet tides problem. It sets up the Earth-Moon triangle geometry, balances gravity and centrifugal force at the center of object 1 to fix ω, then sums real and fictitious forces at a surface point. It evaluates the far-distance limit at the near point A, far point B and the top, finding the centrifugal term dominates and the top may be lighter than A or B. Phil flags the result as possibly in error.
AI-written summary; may contain errors.
Extracted text (machine-read; may contain errors)
Facing Planet Problem PhL 8.13.12
This is a basic outline of how you might do this problem. This is tides on the Earth if the earth faced the moon during its rotation.
Consider this picture
Before doing anything, we have a lot of geometry to figure out. Assume that r12, M1, M2, R2 and θ are given quantities. Then here is how we figure out everything else shown:
θ' = π-θ // determines θ'
R = r12 [M1/(M1+ M2)] // determines R
r1c2 = r122 + R22 - 2r12R2cosθ' // determines r1c
sinβ = sinθ(R2/r1c) // determines sinβ
R22 = r1c2 + r122 - 2r1cr12cosβ // determines cosβ
R32 = R2+R22- 2RR2cosθ' // determines R3
sinα = (R2/R3)sinθ' // determines sinα
R22 = R2+R32- 2RR3cosα // determines cosα
So now we know all quantities in the picture.
Next, we consider a particle at the center of object 1. We will have there this obvious force balance
Fgrav1(center) = -(mM1G/r122) real
Fcent(center) = mω2R fict
Since a particle at the center has a' = 0, these must cancel and we get
M1G/r122 = ω2R // determines ω
Next, we compute the total force on the Particle at C. The magnitude of the upper vector is ω2R3 and this is a fictitious force, the angle is α as shown. The grav1 vector has magnitude mM1G/r1c2 and it is at angle β as shown. Now stare at the picture to conclude that, at point C,
Fgrav1(C) = (mM1G/r1c2) [ cos(π/2-θ'- β) + sin(π/2-θ'- β) ] real
Fcent(C) = (mω2R3) [ -cos(α - π/2 +θ') + sin(α - π/2 +θ') ] real
Fgrav2(C) = -(mM2G/R22) fict
Fsurf(C) = Fsurf real
The total force F'eff acting on Particle C is the sum of these four forces
F'eff = [ (mM1G/r1c2) cos(π/2-θ'- β) - (mω2R3) cos(α - π/2 +θ')]
+ [(mM1G/r1c2) sin(π/2-θ'- β) + (mω2R3) sin(α - π/2 +θ') – (mM2G/R22) + Fsurf ]
I think this is the tidal force at point C. Now simplify by using
cos(π/2-x) = sin(x)
sin(π/2-x) = cos(x)
Then
cos(π/2-θ'- β) = cos(π/2 - [θ'+β]) = sin(θ'+β)
cos(α - π/2 +θ') = cos(π/2 -α - θ') = cos(π/2-[α+θ']) = sin(α+θ')
sin(π/2-θ'- β) = sin(π/2 - [θ'+β]) = cos(θ'+β)
sin(α - π/2 +θ') = - sin(π/2-α - θ') = -sin(π/2- [α+θ']) = -cos(α+θ')
F'eff = [ (mM1G/r1c2) sin(θ'+ β) - (mω2R3) sin(α+θ')]
+ [(mM1G/r1c2) cos(θ'+β) – (mω2R3) cos(α+θ') – (mM2G/R22) + Fsurf]
Now we would like to know the direction of this tidal force.
We finally simplify by assuming that object 1 is far away so r12 >> R2. Then our list above becomes
θ' = π-θ // determines θ'
R = r12 [M1/(M1+ M2)] // determines R
r1c2 = r122 // determines r1c
sinβ = sinθ(R2/r12) // determines sinβ
0 = r122 + r122 - 2r12r12cosβ // determines cosβ ≈ 1 so β ≈ 0 but see above
R32 = R2+R22- 2RR2cosθ' // determines R3 , no simplification
sinα = (R2/R3)sinθ' // determines sinα, no simplification
R22 = R2+R32- 2RR3cosα // determines cosα
M1G/r122 = ω2R // fact also known
Then we have
F'eff ≈ [ (mM1G/r122) sin(θ'+ β) - (mω2R3) sin(α+θ')]
+ [(mM1G/r122) cos(θ'+β) – (mω2R3) cos(α+θ') – (mM2G/R22) + Fsurf]
F'eff ≈ [ (m ω2R) sin(θ'+ β) - (mω2R3) sin(α+θ')]
+ [(m ω2R) cos(θ'+β) – (mω2R3) cos(α+θ') – (mM2G/R22) + Fsurf]
What happens at θ' = 0 ? At that point α = π and R3 = R2-R and β = 0
F'eff ≈ [ (m ω2R) sin(0) - (mω2R3) sin(π)]
+ [(m ω2R) cos(0) – (mω2R3) cos(π) – (mM2G/R22) + Fsurf]
F'eff ≈ + [(m ω2R)+ (mω2R3) – (mM2G/R22) + Fsurf]
F'eff ≈ + [(m ω2R)+ (mω2[R2-R]) – (mM2G/R22) + Fsurf]
F'eff(A) ≈ + [ mω2R2 – (mM2G/R22) + Fsurf]
and since A is at rest, this adds up to 0.
Normally we would have Fsurf = (mM2G/R22) for something at rest on the surface, But now we have'
Fsurf = (mM2G/R22) - mω2R2
so the earth presses less than normal, so something is pulling to the left, and
Ftid(A) = mω2R2 to the left
This is the main term I got in my version 2 document. My approx here killed off the other term which recall is 41x smaller for moon-earth. The centrifugal effect dominates.
What about point B? There θ' = π, α = 0, R3 = R+R2, β= 0. So
F'eff ≈ [ (m ω2R) sin(θ'+ β) - (mω2R3) sin(α+θ')]
+ [(m ω2R) cos(θ'+β) – (mω2R3) cos(α+θ') – (mM2G/R22) + Fsurf]
F'eff ≈ [ (m ω2R) sin(π+ 0) - (mω2R3) sin(0+π)]
+ [(m ω2R) cos(π+0) – (mω2R3) cos(0+π) – (mM2G/R22) + Fsurf]
F'eff ≈ + [- (m ω2R) + (mω2R3) – (mM2G/R22) + Fsurf]
F'eff ≈ + [- (m ω2R) + (mω2[R+R2]) – (mM2G/R22) + Fsurf]
F'eff(B) ≈ + [mω2R2 – (mM2G/R22) + Fsurf]
Now we have
Fsurf(B) = (mM2G/R22) - mω2R2
and again our particle is "lighter" by the same amount mω2R2
Now let's try a point at the top, so θ' = π/2. Then we have
sinα = R2/R sinβ ≈ R2/r12 = very small R3 = [ R2 + R22]1/2
cosα = R/R3
so:
F'eff ≈ [ (m ω2R) sin(θ'+ β) - (mω2R3) sin(α+θ')]
+ [(m ω2R) cos(θ'+β) – (mω2R3) cos(α+θ') – (mM2G/R22) + Fsurf]
F'eff ≈ [ (m ω2R) sin(π/2) - (mω2R3) sin(α+π/2)]
+ [(m ω2R) cos(π/2) – (mω2R3) cos(α+π/2) – (mM2G/R22) + Fsurf]
F'eff ≈ [ (m ω2R) 1 - (mω2R3) cos(α)]
+ [(m ω2R) 0 + (mω2R3)sin(α) – (mM2G/R22) + Fsurf]
F'eff ≈ [ (m ω2R) - (mω2R3) cos(α)]
+ [(mω2R3)sin(α) – (mM2G/R22) + Fsurf]
F'eff ≈ [ (m ω2R) - (mω2R3) (R/R3)]
+ [(mω2R3)(R2/R) – (mM2G/R22) + Fsurf]
F'eff ≈ + [(mω2R3)(R2/R) – (mM2G/R22) + Fsurf]
Now we find that
Fsurf = (mM2G/R22) – (mω2R3)(R2/R)
and again we are lighter. The amount lighter is mω2R2(R3/R) so we have an extra factor. This is a fact I have never computed before. the tidal force is outward at the top and it is larger than at A and B ! This suggests to me that the bulge will be around the belt, and we have dips at the A and B. Maybe I have made an error.
OK, that is enough for now.