11/19/18 at 17:04:24 News: Registration for new users has been disabled to discourage spam. If you would like to join the forum please send me an email with your desired screen name to tony at gravitysimulator dot com.

 Pages: 1 2
 Eccentricity in binaries - simulation 2 (Read 15492 times)
 Tony YaBB Administrator Posts: 1058 Gender: Re: Eccentricity in binaries - simulation 2 Reply #15 - 10/30/06 at 23:44:41   There's nothing in the code that tells it to do inclination cycles.  The code only tells the objects to follow Newton's laws.  Inclination cycles, just like ecc cycles should be a natural consequence of obeying Newton's laws.   Gravity Simulator is giving you all your orbital elements based on the ecliptic, or more specifically, the x-y plane.   The inclination you're reading for the planet is from the ecliptic, not from the plane of the secondary's orbit.  And the cycle you're expecting to see is probably based on the plane of the secondary.  But you can't just subtract one from the other.  Converting coordinates from one plane to the other requires some tricky geometry.  (It's probably not that tricky, but off the top of my head I don't know how to do it).   What if you tried putting the planet in a 75 degree inclination orbit and the secondary in a 0 degree orbit.  Then you wouldn't have to do any conversions since the secondary would be in the ecliptic, and the inclination numbers you plot would now be the planet's plane with respect to the secondary's plane.  I'd bet you'd get the cycles you were looking for. Back to top IP Logged
 EDG Ultimate Member oh, crumbs!!! Posts: 611 Gender: Re: Eccentricity in binaries - simulation 2 Reply #16 - 10/30/06 at 23:52:37   Quote from Tony   on 10/30/06 at 23:44:41:There's nothing in the code that tells it to do inclination cycles. The code only tells the objects to follow Newton's laws. Inclination cycles, just like ecc cycles should be a natural consequence of obeying Newton's laws.   Ok then.     Quote:Gravity Simulator is giving you all your orbital elements based on the ecliptic, or more specifically, the x-y plane. The inclination you're reading for the planet is from the ecliptic, not from the plane of the secondary's orbit. And the cycle you're expecting to see is probably based on the plane of the secondary. But you can't just subtract one from the other. Converting coordinates from one plane to the other requires some tricky geometry. (It's probably not that tricky, but off the top of my head I don't know how to do it).   I didn't get this last time you mentioned it, but now you say that it's relative to the x-y plane it makes more sense. But still... in this sim the planet starts in the x-y plane and the companion is inclined at 75° to it. As the sim progresses the planet inclination starts increasing, but that's still measured from the x-y plane isn't it? So if at one point the planet inclination is 35° and the companion inclination is 75° doesn't that mean the relative inclination is 40°, since they're both measured relative to the same x-y plane?     Quote:What if you tried putting the planet in a 75 degree inclination orbit and the secondary in a 0 degree orbit. Then you wouldn't have to do any conversions since the secondary would be in the ecliptic, and the inclination numbers you plot would now be the planet's plane with respect to the secondary's plane. I'd bet you'd get the cycles you were looking for.   Well I can try that next, I'll leave the current sim running for tonight to see where it goes. It's going to stop recording at 500,000 years anyway. Back to top (formerly known as Mal)   IP Logged
 Tony YaBB Administrator Posts: 1058 Gender: Re: Eccentricity in binaries - simulation 2 Reply #17 - 10/31/06 at 00:06:39   Only if your longitude of the ascending nodes are equal in both orbits can you simply add or subtract inclinations and get a useful number.  Consider the following:   All three orbits have inclinations of 30 degrees.  The purple and green orbits have the same longitude of ascending node.  Therefore you can do 30-30=0 and conclude that the purple orbit is inclined 0 degrees to the green orbit.     But the 30-30 trick doesn't work for the orange orbit.   It is clearly 60 degrees inclined from the green and purple orbits. Instead, you need 30+30 = 60.   The green and purple orbits have a longitude of ascending node of 90 degrees and the orange orbit has a longitude of 270, 180 degrees away from the purple and green.   You may be able to express the inclinations relative to each other in excel with this formula:   i1+cos(deltaLAN)*i2   Excel likes radians instead of degrees, so... =E2+cos((G2-O2)/180 * 3.14159)*M2   Assuming my conversion logic is correct. Back to top « Last Edit: 10/31/06 at 01:11:57 by Tony »     IP Logged
 EDG Ultimate Member oh, crumbs!!! Posts: 611 Gender: Re: Eccentricity in binaries - simulation 2 Reply #18 - 10/31/06 at 07:30:01   420,000 years (planet):     So I guess Arg of Peri, LAN, and Inclination are doing weird looking things because they're not in the same plane as the companion? Maybe it'd look more sensible if they were calculated relative to the companion instead?     Here's what the companion's been doing:   LAN and Arg of Peri aren't shown, because they're oscilating slightly around 0 degrees - so they're going from 1 to 359 degrees and messing up the graph. Back to top (formerly known as Mal)   IP Logged
 Tony YaBB Administrator Posts: 1058 Gender: Re: Eccentricity in binaries - simulation 2 Reply #19 - 10/31/06 at 15:33:09   Quote from Tony   on 10/31/06 at 00:06:39: You may be able to express the inclinations relative to each other in excel with this formula: i1+cos(deltaLAN)*i2 Excel likes radians instead of degrees, so... =E2+cos((G2-O2)/180 * 3.14159)*M2 This is wrong.  I think the correct way is this:   i3 = acos(sin(i1)sin(i2)cos(deltaLAN)+cos(i1)cos(i2))   and for Excel:   =ACOS(SIN(E2/180*PI())*SIN(M2/180*PI())*COS((G2-O2)/180*PI())+COS(E2/180*PI())*COS(M2/180*PI()))/PI()*180     Assuming: i1 is in the E column i2 is in the M column LAN1 is in the G column LAN2 is in the O column.   I'll run a sim to verify. Back to top « Last Edit: 11/01/06 at 16:34:28 by Tony »     IP Logged
 EDG Ultimate Member oh, crumbs!!! Posts: 611 Gender: Re: Eccentricity in binaries - simulation 2 Reply #20 - 10/31/06 at 18:03:48   536,700 years:     Two peaks now, but the odd thing is that they're about 218,000 years apart - I thought the calculation said that the cycle should be 335,795 years?! Back to top (formerly known as Mal)   IP Logged
 Tony YaBB Administrator Posts: 1058 Gender: Re: Eccentricity in binaries - simulation 2 Reply #21 - 10/31/06 at 18:22:13   If instead of measuring peak-to-peak, your measure start to first peak and multiply by 2, the accuracy improves a little.  I think in ideal conditions (ie, much slower timestep), the eccentricity should have gone all the way back down to 0.   Also keep in mind that the formulas in the paper don't use the equal sign.  They use the =~ sign. Back to top IP Logged
 Tony YaBB Administrator Posts: 1058 Gender: Re: Eccentricity in binaries - simulation 2 Reply #22 - 11/01/06 at 16:40:59   I verified the formula I gave above by running 2 simulations.:   1:  planet 2.5 AU, 1 Mj      Secondary 200 AU 0.9 Ms, 75 degrees   2:  planet 2.5 AU, 1 Mj, 75 degrees      secondary 200 AU 0.9 Ms   Time step 32K for 2 million years.   With a small modification the formula works flawlessly.  The inclination of the planes of the planet and the secondary are identical when compared to each other instead of against the ecliptic by using the Excel Formula.   I edited the formula in the above post by placing a few missing multiplication symbols (*) and converting the answer to degrees.  Again, the Excel formula is:   =ACOS(SIN(E2/180*PI())*SIN(M2/180*PI())*COS((G2-O2)/180*PI())+COS(E2/180*PI())*COS(M2/180*PI()))/PI()*180   Assuming:   i1 is in the E column   i2 is in the M column   LAN1 is in the G column   LAN2 is in the O column. Back to top IP Logged
 EDG Ultimate Member oh, crumbs!!! Posts: 611 Gender: Re: Eccentricity in binaries - simulation 2 Reply #23 - 11/01/06 at 17:19:58   789,800 years   I'm calling it quits here. Top graph is the planet, bottom is the companion. I got no idea what all the properties of the planets' orbit that aren't eccentricity are doing...         I'll try those formulas when I get a chance... Back to top (formerly known as Mal)   IP Logged