EmelineBolmont / mercury-90

Automatically exported from code.google.com/p/mercury-90
3 stars 2 forks source link

angular velocity #3

Closed GoogleCodeExporter closed 9 years ago

GoogleCodeExporter commented 9 years ago
The angular velocity (included in the calculation of the scaleheight) is 
defined how?

Omega = sqrt(GM/a^3) is in what condition? with e=0? Does we need an averaged 
angular velocity or an instant one for torque calculation?

Original issue reported on code.google.com by ccos...@gmail.com on 29 Jul 2011 at 2:36

GoogleCodeExporter commented 9 years ago
We need to use the averaged angular velocity on the orbit. Instant angular 
velocity makes no sense knowing the fact that torque of the disk are calculated 
over at least one orbit. 

Thus, the angular velocity does not depend on the eccentricity because the 
period of the planet does not depend on the eccentricity either. 

Original comment by ccos...@gmail.com on 8 Aug 2011 at 8:57

GoogleCodeExporter commented 9 years ago
[deleted comment]
GoogleCodeExporter commented 9 years ago
In revision 5f1e8da7fb64a8dcc69f4b8e8063bb2262c00b1a

Changed again because during close encounters, semi major axis can loose any 
physical meaning du to graviational perturbations of the other body. As a 
consequence I can have sometimes a>1e4 AU while r \sim 4AU. 

To solve this issue I try to use r instead of a. With low eccentricities, this 
should not change results while it should avoid problems during close 
encounters. 

Original comment by ccos...@gmail.com on 26 Apr 2013 at 1:58