Welcome, Guest

Author Topic: Demonstrate the Kozai Mechanism  (Read 7327 times)

infringement153

  • ****
  • Posts: 48
Demonstrate the Kozai Mechanism
« on: June 22, 2010, 07:21:10 AM »
Does anybody have a ubox/xml file that uses the Kozai mechanism to cause a moon to crash into a planet?  Or, if not, does anybody know how to create the sim myself?  Whenever I try to use the Kozai mechanism to cause crashing by putting a moon into polar orbit, all I get is a gentle oscillation of eccentricity/inclination, not an impact.

Something similar to the earth/moon simulation here.
http://www.orbitsimulator.com/gravity/articles/kozai.html

Thank you.

Laura

  • *****
  • Posts: 123
Re: Demonstrate the Kozai Mechanism
« Reply #1 on: June 22, 2010, 08:59:24 AM »
Do you have the planet orbiting a star in the experiment? I believe it's important that the planet is moving for this effect to show up.

infringement153

  • ****
  • Posts: 48
Re: Demonstrate the Kozai Mechanism
« Reply #2 on: June 22, 2010, 09:38:42 AM »
Yes.  I took the normal "solar system" sim, and added a polar-orbiting moon to earth.

APODman

  • Universe Sandbox 1 Beta Team
  • *****
  • Posts: 225
    • Astronomic Picture of the Day - Brazil
Re: Demonstrate the Kozai Mechanism
« Reply #3 on: June 22, 2010, 09:57:28 AM »
Does anybody have a ubox/xml file that uses the Kozai mechanism to cause a moon to crash into a planet?  Or, if not, does anybody know how to create the sim myself?  Whenever I try to use the Kozai mechanism to cause crashing by putting a moon into polar orbit, all I get is a gentle oscillation of eccentricity/inclination, not an impact.

Something similar to the earth/moon simulation here.
http://www.orbitsimulator.com/gravity/articles/kozai.html

Thank you.

You try convert the kozai simulation file of Gravity Simulator to Universe Sandbox simulation file ?

[ ]´s

atomic7732

  • Global Moderator
  • *****
  • Posts: 3849
  • caught in the river turning blue
    • Paladin of Storms
Re: Demonstrate the Kozai Mechanism
« Reply #4 on: June 22, 2010, 10:31:55 AM »
Or simulate the HD 80606 (Struve 1374?) system.

Dan Dixon

  • Creator of Universe Sandbox
  • Developer
  • *****
  • Posts: 3244
    • Personal Site
Re: Demonstrate the Kozai Mechanism
« Reply #5 on: June 22, 2010, 10:52:47 AM »
Here's the converted XML file of the kozai.gsim simulation.

Gravity Simulator's latest beta version has a Universe Sandbox XML format exporter. I used that to create this file. Although I think there's a bug in the exporter and you need to set your timestep to 1 before you export.

APODman

  • Universe Sandbox 1 Beta Team
  • *****
  • Posts: 225
    • Astronomic Picture of the Day - Brazil
Re: Demonstrate the Kozai Mechanism
« Reply #6 on: June 22, 2010, 10:59:53 AM »
Here's the converted XML file of the kozai.gsim simulation.

Gravity Simulator's latest beta version has a Universe Sandbox XML format exporter. I used that to create this file. Although I think there's a bug in the exporter and you need to set your timestep to 1 before you export.

It works, the Kozai mechanism is demonstrated but an impact don´t occurs.

[ ]´s

frankdepinte

  • **
  • Posts: 17
Re: Demonstrate the Kozai Mechanism
« Reply #7 on: June 22, 2010, 11:40:41 AM »
Has the timestep been checked during the simulation ? if to big the collision doesn't happen

Laura

  • *****
  • Posts: 123
Re: Demonstrate the Kozai Mechanism
« Reply #8 on: June 22, 2010, 10:04:54 PM »
Experiment: Demonstrate Kozai Mechanism

Initial conditions: The kozai.xml file provided by Dan.

Parameters: RK4 accuracy, Timestep 30 minutes (except at the very end, where it was reduced to more accurately check for collision).

Observations: Orbit of moon oscillates around inclination 90, tending more towards positive growth rather than negative. During periods of greatest positive deviation from 90, eccentricity grows markedly, while it comes to a standstill at periods of greatest negative deviation. With each such cycle, the effect becomes more pronounced, and the deviation in inclination greater. Eccentricity never reduces; it only grows exponentially, and as it does so, deviation from inclination 90 becomes increasingly pronounced.
Collision took place near the end of the 5th year, with inclination exceeding 110 degrees, and eccentricity somewhere between 0.99 and 1. Timestep was reduced to 5 minutes for the last few cycles.
Regrettably, I have no visual data for that final year which exhibited the greatest changes - I should not have stuck so doggedly to the policy of taking one picture at the end of each year.
Below is an orbital plot covering the time from the beginning of the simulation to the end of the 4th. year. Each year is color coded to make it easier to see the age of the data.

Conclusion: The Kozai Mechanism can be demonstrated with US2 :)