To what extent could a single Triton flyby slow down a direct Hohmann transfer to Neptune for NOI?

Question

I can't seem to find any good sources online for this, all I get are documents on how Nice model was used to compute Neptune's gravity capture of Triton via a binary dissociation, possibly because of similar terminology. So to the question;

To what extent, in terms of % delta-v, could Triton's retrograde orbit around Neptune gravitationally assist in slowing down a direct Hohmann transfer from Earth to Neptune for Neptune Orbit Insertion (NOI) during a single close-shave flyby of Triton?

I realize that the transfer would take about good three decades (I calculated 30.61652 years at Hohmann semi-major axis of 15.53545 AU using Kepler's third law and semi-major axes of Earth and Neptune), transfer time is not a concern, but I'm stuck with calculating velocity relative to Neptune at Triton's altitude for this transfer (I'm assuming Triton to be at the right place at the right time) and I didn't yet get to calculate delta-v that a low altitude (100 km above surface) gravity assist flyby of Triton could shed, so I'm not sure I'd know how to do that either. Neptune NOI would be for a Neptune orbit with semi-major axis of roughly 5,000 km above its surface (at 1 atm pressure). There wouldn't be any plane change for the target orbit.

This is not a homework question. It's been about 20 years since I last had to do such calculations back in the college years and I'd appreciate a bit of help in brushing up on this. I could plug this in some software but, call me masochistic, I wanted to do it the old-fashioned way. Back-of-an-envelope calculations would do, preferably discussing any shortcuts there might be in getting first-order approximations faster. No need to be too academic either, but I would like to see some calculations here, if that's not terribly inconvenient, in which case I guess a good reference with a short-ish writeup would do, too.

If you'll opt for the latter option using external references, please also include optimal launch dates, delta-v and Hohmann transfer times.

Answer 1

I'll try to get you started anyway. From the frame of reference of the assisting body, the trajectory of the probe is hyperbolic, with the same $v_\infty$ going out as coming in, but in a different direction. The trajectory is simply bent. The bend angle is:

$$\delta=2\sin^{-1}\left(1\over e\right)$$

where $e$ is the eccentricity of the hyperbola. You can get $e$ from $v_\infty$, the closest approach radius $r$ (measured from the center of the body), and the $GM$ of the body, $\mu$:

$$e={r\,v^2_\infty\over\mu}+1$$

$e$ goes from $1$ to $\infty$, where $1$ is effectively an infinite gravity assist that reverses the direction of the velocity by 180°, and $\infty$ is no assist at all with a bend of 0°. The closer in you can get and the slower you can go, the greater the bend angle is.

With the vector addition of the velocity of the assisting body (Triton) in the frame of the body it is orbiting (Neptune), you can see how the velocity of the probe is changed by the flyby in the Neptune frame, which is what you care about. You will need to do the reverse to get the $v_\infty$ approaching Triton, i.e. subtract the velocity vector of Triton from the approach velocity vector in the frame of Neptune. To get that, you subtracted the velocity vector of Neptune from the velocity in the Sun frame of the Hohmann transfer. Lots of frame changes.

There is a range of geometries for the flyby you can try, depending on where in its orbit Triton is at the time, the altitude of the flyby, and whether the flyby closest approach is on the leading or trailing side of Triton in its orbit.