Page 1 of 1

Precessionrate calculations

Posted: 14.12.2006, 12:59
by rthorvald
I am unable to get this calculation to work. Either i am misunderstanding something, or there are inaccuracies i am unaware of... Please comment:

(I have simplified the numbers):

- Planet RotationPeriod is 18 hours
- Satellite Period is 0.75 days, with no RotationPeriod (tidally locked)
- Satellite Obliquity of 45 degrees
... To make the satellite stand still (not change orientation relative to the planet), as i understand it, the Precessionrate should be:

360 divided by 0.75 = 480.

... This works short term. But in a thousand years, there is a considerable drift.

What am i doing wrong?

- rthorvald

Posted: 14.12.2006, 18:16
by Cham
The problem with PrecessionRate is we can't change its axis. I asked several times to have the ability to change the PrecesionRate axis with an SSC command like "PrecessionAxis [x x x]". I don't know why Chris is hesitant to add this.

Posted: 14.12.2006, 18:28
by rthorvald
Cham wrote:The problem with PrecessionRate is we can't change its axis.


Thanks, but is that neccecary to solve my problem?

- rthorvald

Posted: 14.12.2006, 18:30
by Cham
I don't know. Your number should be ok, if the values are exact. If you rounded some numbers, then it may explain the drift.

Posted: 14.12.2006, 18:38
by rthorvald
Cham wrote:I don't know. Your number should be ok, if the values are exact. If you rounded some numbers, then it may explain the drift.


I have simplified the numbers to get the concept right, since the real values gave me too many decimals for my calculator to handle... And i thought the problem was that the result was not accurate enough. So, i re-stated the problem with simpler values to test it, but it comes out wrong still.

- rthorvald

Posted: 14.12.2006, 18:46
by Cham
The drift is most probably caused by your rounded values. Try some EXACT values (don't do approximations). Try some simple integers and select some so the division could give another integer.

EDIT : I've experienced similar problems before, and I had to fine-tune my precession rate value by trial and error, since the exact calculations wasn't giving the right value for a long period of time. I suspect the cause may be an approximation in Celestia's code. There is probably a Pi somewhere in the code, rounded to something like Pi = 3.1415927, or something like that, since a true rotation rate is given by

Code: Select all

Omega = 2 Pi / Period

Posted: 14.12.2006, 19:10
by rthorvald
Cham wrote:The drift is most probably caused by your rounded values. Try some EXACT values (don't do approximations). Try some simple integers and select some so the division could give another integer.
This is exactly what i have done.

Cham wrote:EDIT : I've experienced similar problems before, and I had to fine-tunned my precession rate value by trial and error, since the exact calculations wasn't giving the right value for a long period of time.

Yes... I am playing with it right now. I need the thing to stay in place for at least three thousand years, but after 7-800, it is noticeably off track.

- rthorvald