A combined strategy of aerocapture and aerobraking is presented to achieve a near-circular orbit, starting from a hyperbolic trajectory, without requiring an orbital insertion burn. Aerothermodynamic force, moment, and heat flux calculations employ a Maxwellian free-molecular flow model, with Knudsen-number-based interpolations for the transition regime. The six-degree-of-freedom motion model, including quaternion-based attitude dynamics, allows stability and sensitivity analyses for the atmospheric passes. A spacecraft model with two large panels suitable for Earth aerocapture is considered. Minor orbit-correction burns at the apoapsis are provided after each pass for manipulating the periapsis for the next pass to meet the desired aerobraking corridor. It is observed that an initial orbit of eccentricity 1.6 and relative entry velocity of 12 km/s at 300-km altitude can be reduced to an orbit with an eccentricity of 0.02, using a total of six atmospheric passes, without exceeding the peak convective heat flux constraint for the spacecraft or requiring an orbit insertion burn. This result has considerable importance for low-Earth-orbit space-tug captures and Mars missions, wherein the strategy proposed will lead to significant savings in spacecraft propellant mass during the orbit insertion and, subsequently, orbit circularization.
Nomenclature= orbital eccentricity F = external force vector, ∈ 3 × 1 J = inertia tensor of the spacecraft, ∈ 3 × 3 J 11 , J 22 , J 33 = principal moments of inertia of the spacecraft K n = Knudsen number M = external moment vector, ∈ 3 × 1 m = mass of the spacecraft n i = unit normal vector of the ith surface panel p, q, r = angular velocity components about the X , Y , and Z body axes Q = convective heat load on the spacecrafṫ Q = convective heat fluẋ Q max = maximum convective heat flux q 0 , q 1 , q 2 , q 3 = quaternion components used to describe body orientation R = distance of spacecraft center of mass from the planet's center r i = position vector of the ith panel centroid relative to the center of mass of the spacecraft, r i ∈ 3 × 1 S(ω) = skew-symmetric matrix, ∈ 3 × 3 s = molecular speed ratio T W /T ∞ = ratio of wall temperature to freestream temperature V = spacecraft velocity measured in the planet fixed coordinate system X fo , Y fo , Z fo = forces along the wind axes α = angle of attack