Closed huhell closed 4 years ago
Hi,
Unfortunately it's not well defined to just set the eccentricity or semimajor axis (you have to choose what to do with the remaining elements), so REBOUND doesn't let you do that.
If you just want exponential damping of the eccentricity and semimajor axis, you can just use modify_orbits_forces that's already implemented, and just set the appropriate tau_a and tau_e using your tidal parameters (tau_e = 1/de/dt).
If that's not good enough, you're doing the right thing by implementing an operator, but you have to set xyzzy and vxvyvz for an operator. One option would be to follow appendix A of Lee and Peale 2002.
Hi, Thank you very much for the clarification !
Hello, I am trying to add tidal dissipation as a custom operator, which would update the eccentricity, and semimajor axis as e += dedt dt and a += dadt dt. I use Equations (1) and (2) from Jackson et al. (2008), ApJ, 678, 1396. Here is my code
But when I integrate using
sim.integrate()
, I get the following error:Is it impossible to update the eccentricity and semimajor axis ?
Thanks a lot !