Interesting... I was plotting a maneuver to fully reverse the orbital plane (of a circular orbit) from 180º to 0º by raising the apoapsis, reverse the planes at the new apoapsis and finally lower the apoapsis to its original height and I wanted to find the minimum apoapsis where the maneuver would be more efficient than a straight retrograde burn on the original orbit. Turns out, in KSP any higher apoapsis seems to do. I also found the minimum amount of delta-v required for this maneuver to be roughly 41.4% of the delta-v for the simple retrograde burn. In other words, you save at most 58.6% delta-v doing this maneuver, the higher the better. (Or in KSP's case, aim for an apoapsis close to SoI.) Furthermore, with an apoapsis at least 8 times higher than the circular orbit, you'll save at least 50% delta-v.
Example, if you're orbiting Eve (Radius 0.7 Mm) at 1 Mm (1 Mm + 0.7 Mm = 1.7 Mm), you can save 50% delta-v with a "work" apoapsis at 12.9 Mm (1.7 Mm * 8 - 0.7 Mm).
Assume circular orbit.
Initial velocity (circular orbit): v
Initial orbital height (calculated from planet center): R
"Work" apoapsis: h
Apoapsis speed: u
Planet gravitational parameter: μ
Delta-v for initial transfer burn:
v = sqrt(μ/R)
v1 = sqrt(μ/R)sqrt(2h/(R+h))
|DV1| = v1 - v = sqrt(μ/R)(sqrt(2h/(R+h)) - 1)
Delta-v for "reverse" burn:
|DV2| = 2u = 2*sqrt(μ/h)sqrt(2R/(R+h))
Delta-v for circularizing:
|DV3| = |DV1| = sqrt(μ/R)sqrt(2h/(R+h)) |Initial orbit with reverse direction
Total delta-v:
DVtot = |DV1| + |DV2| + |DV3| = 2*|DV1| + |DV2|
= 2*[sqrt(μ/R)(sqrt(2h/(R+h)) - 1)] + 2*sqrt(μ/h)sqrt(2R/(R+h)) = 2*sqrt(μ)[sqrt(1/R)(sqrt(2h/(R+h)) - 1) + sqrt(1/h)sqrt(2R/(R+h))]
Delta-v for "simple" retrograde burn:
|DVrev| = 2v = 2*sqrt(μ/R)
Minimum height h:
DVtot < |DVrev|
2*sqrt(μ)[sqrt(1/R)(sqrt(2h/(R+h)) - 1) + sqrt(1/h)sqrt(2R/(R+h))] < 2*sqrt(μ/R) |/[2*sqrt(μ)]
sqrt(1/R)(sqrt(2h/(R+h)) - 1) + sqrt(1/h)sqrt(2R/(R+h)) < sqrt(1/R) |*sqrt(R)
sqrt(2h/(R+h)) + R*sqrt(1/h)sqrt(2/(R+h)) < 2 |*sqrt(h)
h*sqrt(2/(R+h)) + R*sqrt(2/(R+h)) < 2*sqrt(h) |*sqrt((R+h)/2)
h + R < 2*sqrt(h)sqrt((R+h)/2) |Both sides positive, square both sides
(h + R)^2 < 2h(R+h)
(h + R)^2 - 2h(R+h) < 0
(R+h)(R-h) < 0 |R+h > 0
h > R
Q.E.D.
Assume h>R for a given R>0. (Specifically above atmosphere...)
DVtot(h) = 2*sqrt(μ)[sqrt(1/R)(sqrt(2h/(R+h)) - 1) + sqrt(1/h)sqrt(2R/(R+h))]
Note that for R > 0, h > 0, the function g(h) = sqrt(2h/(R+h)) - 1, g(h) > 0, h>0 is always increasing, which is easily seen from a first order derivative test.
Thus the function DVtot(h) should also always increase as its second part, sqrt(1/h)sqrt(2R/(R+h)), is never negative.
Therefore, an extrema can only be found in the end intervals of h according to Fermat's theorem. For h = R, the DVtot(R) = 2*sqrt(μ/R) which is the same as DVrev. The upper interval is at infinity: limh→∞ DVtot(h) = limh→∞ 2*sqrt(μ)[sqrt(1/R)(sqrt(2h/(R+h)) - 1) + sqrt(1/h)sqrt(2R/(R+h))] = 2*sqrt(μ/R)(sqrt(2) - 1).
Thus the maximum efficiency is 100% - [2*sqrt(μ/R)(sqrt(2) - 1)/2*sqrt(μ/R)]*100% = 100%*(2-sqrt(2)) ~ 58.6 %
Q.E.D.
2*sqrt(μ)[sqrt(1/R)(sqrt(2h/(R+h)) - 1) + sqrt(1/h)sqrt(2R/(R+h))] = 0.5*2*sqrt(μ/R) |/[2*sqrt(μ)
sqrt(1/R)(sqrt(2h/(R+h)) - 1) + sqrt(1/h)sqrt(2R/(R+h)) = 0.5*sqrt(1/R) |*sqrt(R*h)
h*sqrt(2/(R+h)) + R*sqrt(2/(R+h)) = 1.5*sqrt(h) |*sqrt((R+h)/2)
h + R = 1.5*sqrt(h)*sqrt((R+h)/2) |Square both sides
(h+R)^2 = (9/8)*h(R+h) |/(R+h)
h+R = (9/8)h
h/8 = R
h = 8R
Q.E.D.