Effect of Spacecraft Aerodynamics and Heat Shield Characteristics on Optimal Aeroassisted Transfer ()
1. Introduction
One of the future objectives in space activities is to use aeroassisted orbital maneuvers, i.e., maneuvers that are carried out with the aid of the atmosphere, to satisfy the increasingly stringent constraints of the cost of space missions. Depending on the goal of the mission, aeroassisted maneuvers are performed in different ways: aerogravity assist, aerobraking, and aerocapture. To these maneuvers, one must add the re-entry maneuvers and some possible variants such as the skip re-entry technique. These maneuvers are very attractive, although they are only theoretical at this point (except for rare uses in aerobraking), because they can substantially reduce the requirements of a space mission in terms of propulsion and flight time in favor of, among other things, the possibility of housing a larger payload.
Underlying this approach is the possibility of exploiting the presence of the atmosphere of the celestial body around which one wants to operate to lower the overall energy required. In practice, one tries to execute the complete operation with the help of aerodynamics because orbital maneuvers are quite expensive in terms of propulsion, especially those outside the orbital plan. Thus, the design of a spacecraft with specific reference to atmospheric portion of flight must present an efficient aerodynamic configuration; however, the priority constraint of the overall cost of space missions must be satisfied.
From this brief introduction, it is already clear that once shown the technical feasibility of the aeroassisted maneuver, one must evaluate its convenience compared with alternative hypotheses, such as a classical purely propulsive maneuver. In fact, optimizing a spacecraft and its mission, or more specifically its trajectory for a given mission, is always a compromise between the interests of performance, security, and economics, which are almost always in mutual conflict.
Usually, minimizing the mass of the heat shield while respecting the limits of safety is a primary requirement. Indeed, any savings in terms of the mass of the thermal protection system (TPS) can be translated into an increase of the “useful” mass, namely, into an opportunity to accommodate a greater payload or extend other spacecraft subsystems. The problem becomes more complex when considering an aeroassisted maneuver and taking into account the mass of propellant still required during some of the various phases of the maneuver itself.
In this regard, a basic scheme is usually adopted in the literature to perform an aeroassisted maneuver—the use of propulsion only in space together with one or more atmospheric segments of pure aerodynamic flight—which is also considered in this work. Maximizing the benefits at the propulsive level for an aeroassisted maneuver normally requires a more intensive use of the atmospheric phase of flight. Thus, the resulting trajectory touches the denser layers of the atmosphere with longer crossing times. Therefore, a larger TPS with its relative higher mass fraction must be adopted.
A legitimate question that arises at this point is whether the resulting increase in the mass of the TPS, together with the mass of propellant required to enter the atmosphere and then to achieve the final orbit, may override the convenience of the aeroassisted maneuver compared with a classical operation, which is a purely propulsive extra-atmospheric maneuver. It seems evident that this is an optimization problem.
All of the above considerations concerning both the evaluation of the convenience of the aeroassisted approach and the search for the optimal solution and the methodological sphere itself are the basis of the research questions that motivated this work.
An original procedure was carried out jointly with the implementation of software developed by the authors to optimize the aeroassisted orbital maneuvers using a genetic algorithm (GA) with simultaneous evaluation of the optimal configuration of the associated heat shields and the coupling of the dynamic and thermal analyses. The tool was verified by comparing its results with those found in the literature. This tool, because of its level of implementation details, is suitable for the conceptual development stages of a spacecraft and its mission.
The initial intention was to evaluate the influence of various parameters—orbital, aerodynamic, and dimensional—on the feasibility and convenience of the mission. For an introductory analysis of the problem and to evaluate the importance of various factors, coplanar transfer (aerobraking) from a high Earth orbit (HEO) to a low Earth orbit (LEO), both circular, is proposed as a case study. The vehicle is a delta wing spacecraft that is protected by an ablative TPS with uniform thickness. The problem is studied in the absence of constraints on the maximum allowable entering heat flux.
Beyond the presentation of the problem in Section 1, Section 2 describes the model and the optimization procedure. The case study is presented in Section 3, and the relevant results and analyses are discussed in Section 4. Finally, Section 5 offers a summary, conclusions, and recommendations for future improvements.
2. Models and Optimization
The description of the models, the governing equations, and the relevant assumptions for the problem—i.e., thermal models, aerodynamic heating, atmospheric flight mechanics, and heat shield configuration—are thoroughly presented in [1] and [2]. Reference [2] can be consulted for full details of the original optimization procedure.
All of the analysis presented here was performed using the cited software developed by the authors called ATHSHO (Aeroassisted Trajectory and Heat SHield Optimization). It is important to recall that thermal analysis is performed with a one-dimensional plane model and that the adopted GA refers to a mixed one-point/two-point crossover operator together with a reproduction plan that provides a full generational replacement with elitism [3-5].
3. Case Study
A parametric analysis was performed by varying the following characteristics of the vehicle: aerodynamic efficiency, mass-to-surface ratio, deorbit impulse, and the altitude of the initial orbit. The influence of the characteristics of the TPS was examined by assessing the type of ablative material adopted, the thermal safety factor, and the allowable temperature for the bond-line, i.e., the adhesive layer on the substructure.
3.1. Hypotheses
The main hypotheses considered are the following:
• The initial total mass of the vehicle is given.
• The attitude control is accomplished only through the angle of attack.
• The entering heat flux is unconstrained.
• The TPS is fully ablative with uniform thickness.
3.2. Vehicle
The vehicle model is a delta wing shuttle with a high ratio, which is comparable in the first instance to the configuration and dimensions of the Boeing X-37A vehicle (Figure 1).
The dimensions, sizes, and aerodynamic characteristics of the vehicle used were taken in part from [6-8]. Other data were reasonable assumptions made by the authors. The main dimensions and characteristics of the vehicle are listed in Table 1, and the principal aerodynamic and propulsive parameters are listed in Table 2.
The values shown are those considered for the nominal reference case. Some of them may vary depending on the purpose of the study (Section 4).
Figure 1. Boeing X-37A (source NASA.gov).
Table 1. Vehicle dimensions and characteristics.
Table 2. Vehicle’s aerodynamic and propulsive characteristics.
3.3. Mission and Maneuver
Table 3 lists the values of the altitudes for the initial Geostationary Earth Orbit (GEO) and final LEO for the aerobraking maneuver as well as the conventional height assumed for the atmosphere.
Even in this case, the indicated values are those used for the reference case. In particular, in the parametric analysis performed, the value of the initial HEO altitude was varied. The physical properties of the atmosphere were derived from the model 1976 US Standard Atmosphere [9].
Figure 2 shows a classic schematic for a HEO-LEO aerobraking maneuver. The strategy involves the combined use of aerodynamic maneuvering in the atmosphere and some extra-atmospheric propulsion phases. More precisely, one assumes that the propulsive phases are concentrated in three impulses in space and that the portion of atmospheric flight is performed without the
Figure 2. Schematic of a HEO-LEO aerobraking maneuver.
use of the propulsion system. The first propulsive impulse (deorbit) at the altitude of the initial HEO varies the vehicle’s speed by to enter the atmosphere along an elliptic orbit segment. The second impulse (boost) is applied upon the exit from the atmosphere to achieve the final LEO by ascending once more along an elliptic orbit segment. The boost is expressed by a speed variation of the vehicle,.
The third and final impulse (circularizing) varies the speed by to circularize the vehicle’s path within the altitude of the final LEO.
It is useful at this point to describe the various phases of the maneuver in more detail as follows. Initially, the vehicle is moving on a circular orbit of radius with a speed around the Earth, which has radius. The expression of the circular speed is the following:
(1)
where
(2)
The deorbit is accomplished by applying the first impulse in the opposite direction of the spacecraft’s speed. This impulse puts the vehicle along an elliptic orbit with the perigee inside the dense layers of the atmosphere.
The atmospheric region below the altitude in which the aerodynamic effects are considered to be conventionally present is denoted the sensible atmosphere. One can determine the speed and the flight path angle of the vehicle’s trajectory at the atmospheric entry point () through the following relationships, which are obtained according to energetic considerations:
(3)
(4)
where
(5)
Clearly, it is necessary that the applied be greater than the minimum for which there would be only a tangential trajectory to the edge of the sensible atmosphere. This is presented in Figure 3 and given by the following expression:
(6)
During the atmospheric portion of flight, the vehicle performs the required maneuver, which is optimally controlled by modulations of the angle of attack (leaving the bank angle fixed) subject to the heating constraints, if any. During this phase, the vehicle’s speed decreases because of aerodynamic drag, and because of this loss of energy, a new impulse is necessary to achieve the final altitude.
At the end of atmospheric flight, the vehicle is situated at an altitude again, is driven at a speed, and has a flight path angle equal to. At this moment, a boost impulse is applied to enter an ascending elliptic orbit with the apogee equal to the radius of the final circular orbit.
The required, as a function of and, can be found from the following expression:
(7)
where
(8)
Figure 3. ΔVl,min vs initial HEO radius.
Once the final altitude is reached, the third impulse is applied to circularize the final orbit. The expression for the third is the following:
(9)
3.4. Fitness Function, Objective Function, and Constraints
The total initial mass of the vehicle is the sum of the propellant mass, the TPS mass, and the structural and payload masses. The goal of the current optimization problem is to perform the assigned orbital transfer while minimizing the sum of the mass of the propellant and the mass of the TPS needed. The objective function, which must be maximized, is then given by the final mass of the vehicle, which can be defined as the performance index of the problem. There are some events that cause variation in the mass of the vehicle. The first is the consumption of fuel due to the deorbit impulse; thus, one can calculate the mass of the vehicle entering the atmosphere. The latter can be derived directly using the Tsiolkovsky equation for the impulse in question:
(10)
This mass is further reduced by during the passage through the atmosphere because of TPS ablation. The mass loss due to ablation comes from both the surface recession and the material density change due to pyrolysis. Thus, the mass of the vehicle at the atmosphere exit is the following:
(11)
At this point, the boost and circularization impulses are applied in sequence, and the final vehicle mass is obtained by two successive applications of the Tsiolkovsky equation:
(12)
Compliance with the eventual constraint on the entering heat flux is ensured through a reward factor that is added to the objective function () by means of an appropriate multiplicative weight, which is chosen specifically to achieve rapid convergence. Likewise, the objective function is multiplied by its weight. Thus, the final expression for the fitness function for the genetic optimizer is as follows:
(13)
In this case study, the heat flux is considered unconstrained, and assuming wm = 1, the previous expression reduces to the following:
(14)
That is, the fitness function coincides with the objective function.
3.5. Assessment of the Convenience
The convenience of the aeroassisted maneuver can be assessed with respect to an equivalent Hohmann transfer. Such a transfer is completely “all propulsive” outside the atmosphere and based on two impulsive velocity changes at two points with radii RA and RB, respectively:
(15)
(16)
Thus, the total “all propulsive” impulse is:
(17)
Consequently, the expression for the final mass of the vehicle in the case of the Hohmann transfer is as follows:
(18)
4. Results and Analysis
The case study for the HEO-LEO coplanar transfer through aerobraking, which was described in the previous Section 3, has been addressed by parameterising the mass-tosurface ratio of the vehicle, the corresponding to the deorbit impulse applied, the aerodynamic efficiency (with fixed), the ablative material used, the maximum allowable temperature for the bondline, and the thermal safety factor. In addition, the above parametric analysis was performed for different values of the initial HEO altitude, which are expressed as fractions of the GEO.
More specifically, the following values were adopted:
• HEO Altitude: 35,786 km (=GEO); 53,679 km (HEOH = 1.5 GEO); 17,893 km (HEOL = 0.5 GEO). The altitude for the final LEO is the same for all cases and equal to 480 km.
• : 200 kg/m2, 300 kg/m2 and 450 kg/m2, corresponding, respectively, to low, medium and high load values. Incidentally, 300 kg/m2 is the value found most frequently in the literature. These load values correspond to the initial mass of the vehicle: 2338.0 kg, 3507.0 kg, and 5260.5 kg.
• range: the range of the parametric first variation of (in steps of 10 m/s) between the minimum value (rounded to the nearest ten higher) given by Equation (6) and the maximum value that allows a solution to complete the aeroassisted maneuver. The minimum values of the deorbit impulse for the three HEO altitudes analyzed are (Figure 3) 1490 m/s for GEO, 1440 m/s for HEOH, and 1420 m/s for HEOL.
• : 1.5 and 3 correspond, respectively, to a medium and a high aerodynamic efficiency. These values were obtained by varying and.
• : 1 and 2 correspond, respectively, to a minimum nominal thickness of the TPS and a “safer” double thickness.
• : 176.85˚C (450 K) and 400˚C (673.15 K) were both acquired from data in the literature; the first is used for standard adhesives, and the second is used for high-performance adhesives.
• Existing adhesives with enhanced characteristics allow the threshold temperature to be increased to 400˚C, which results in greater heat shield thickness savings [10].
• Ablative material: PICA-15 (hereinafter referred to as PICA), AVCOAT 5026-HCG (hereinafter referred to as AVCOAT), and FM 5055 CP in its Reduced Density version (hereinafter referred to as RDCP).
Among all the possible combinations of these parameters with their associated values, 12 cases were analyzed for each value of the ratio for a total of 36 different scenarios. The combinations analyzed are presented in Tables 4-6, where the index (a, b, and c) marks the adopted value for the load. Moreover, because the range of variation of is analyzed by steps of 10 m/s, 390 different cases form the database of results for the analysis of the case study in question.
As a reference, the nominal case is chosen for the GEOHEO transfer that provides a TPS made of PICA, with = 450 K, = 1, and = 1. The twelve combinations chosen and the corresponding purpose of
Table 4. HEO-LEO transfer, parameters and values (cases 1, 2, 3, and 4).
Table 5. HEO-LEO transfer, parameters and values (cases 5, 6, 7, and 8).
Table 6. HEO-LEO transfer, parameters and values (cases 9, 10, 11, and 12).
the evaluation are summarized in Table 7. Figures 4 to 15 report the trends of the mass gain compared with the “all propulsive” case with varied applied and parameterized as a function of the three chosen values of, which were analyzed for each scenario.
The graph on the left of each figure shows the cited mass gain in absolute terms, whereas the graph on the right presents its performance in terms of the percentage of the initial vehicle’s mass.
The mass gain is therefore a measure of the convenience of the aeroassisted maneuver compared with the corresponding Hohmann transfer. When reading these graphics, the definitions of the “range of feasibility” and the “range of convenience” are introduced.
The first one is the interval of in which the aeroassisted maneuver is feasible, whereas the second interval is the range in which the aeroassisted operation is more convenient than the “all propulsive” one. As mentioned above, having made the run for all values of leading to a solution, the range of feasibility is directly readable as the x-axis interval of the definition of the curves themselves.
The range of convenience is by definition the interval in which the curve takes positive values. The range of convenience is always within the range of feasibility. Table 8 shows, for the value = 200 kg/m2 (case “a” of Tables 4-6), the values of the three propulsive impulses applied in each scenario, which in each case give the highest convenience for the maneuver (maxima of the “blue” curves in the graphs of Figures 4 to 15).
It is evident that the utmost convenience for the aeroassisted maneuver corresponds to a value of that is very close to zero, i.e., in the absence of the boost impulse (case 5a has no feasible solutions for the considered value for the load).
This outcome is valid in general for the other values of. Thus, the trajectory optimization process tends to reduce the vehicle’s speed while exiting the atmosphere to the “right” value to reach the final LEO with minimal energy expenditure. Saving this fuel for the second impulse is the basis for the achievement of the highest convenience.
Comparing Figures 4 to 15, one can draw some interesting conclusions:
• The increase of in all cases involves a decrease in the range of feasibility; in particular, the maximum possible diminishes (the right bound of the range).
• The increase of (Figure 5) from the standard value to the high-performance value due to the reduced thickness of the TPS improves the maximum achievable convenience (by approximately 2%).
Table 7. HEO-LEO transfer, scenarios analyzed.
Table 8. Propulsive impulses in the more convenient maneuvers for “a” cases.
Figure 4. HEO-LEO scenario 1: reference nominal case.
• The increase in the aerodynamic efficiency of the vehicle (Figure 6) involves a significant extension of the range of feasibility, but a more moderate increase in the range of convenience. Cases with higher benefit most. Even here, the maximum achievable convenience increases (approximately 2% more).
• The doubling of the thermal safety factor (Figure 7), which increases the mass of the TPS, represents a significant performance penalty for the aeroassisted maneuver, which remains affordable for a short interval of only for the case with highest.
• The use of other ablative materials (Figures 8 and 9),
Figure 5. HEO-LEO scenario 2: TBL,lim influence.
Figure 6. HEO-LEO scenario 3: e influence.
Figure 7. HEO-LEO scenario 4: TSF influence.
Figure 8. HEO-LEO scenario 5: material (RDCP) influence.
Figure 9. HEO-LEO scenario 6: material (AVCOAT) influence.
compared with the reference material PICA, appears to be inappropriate for a GEO initial altitude. Using the RDCP, the sole case with the highest is feasible, but still not convenient. Conversely, using the AVCOAT increases the feasibility compared with RDCP but limits convenience compared with PICA. Thus, the relatively low density of the PICA is a discriminating element in the selection of the ablative material. Actually, these results are conservative because of the assumption that the vehicle is covered with a TPS with a uniform initial thickness.
In Figure 8, all runs with load values of 200 kg/m2 and 300 kg/m2 give results, but the curves are not represented because the final “useful” masses (i.e., all the masses that are not propellant and TPS) are negative. Thus, the vehicle is able to get out of the atmosphere, but without enough fuel to complete the mission.
• Figures 10 (HEO high altitude) and 14 (HEO low altitude) show that the aeroassisted maneuver is much more convenient if the HEO altitude is higher. Lowering the initial altitude reduces the convenience but considerably extends the range of feasibility. The reduced convenience at low altitudes is due to the fact that the Hohmann transfer requires a lower fuel quantity with the lowered the initial altitude for the same final LEO. In contrast, for the corresponding aeroassisted maneuver, even if the TPS is less consistent, the propulsive contribution remains significant.
• Figures 10-13, considered in pairs (the first two and the second two), provide a performance comparison between PICA and AVCOAT, respectively, for high initial orbits in the case of low and high aerodynamic efficiency. The comparison is in favor of the use of PICA, albeit less marked than for GEO orbits. For high orbits for a vehicle with high, the performance of PICA remains considerable; however, the performance of AVCOAT, although lower than PICA, is notable. In this context, it seems appropriate to comment on NASA’s recent evaluation concerning the choice of material for ablative TPS of the Orion spacecraft. The final decision was between PICA and AVCOAT, and NASA ultimately chose the latter. PICA has been used to date only on a re-entry capsule of the Stardust probe, while AVCOAT has a long history of reliability (Apollo missions). Though the results of the present work provide evidence of the tangible prevalence of the PICA performance, the reliability of AVCOAT most likely outweighs these performance considerations.
• The increase in aerodynamic efficiency has less impact in the case of lower orbits (Figures 14 and 15).
Figure 10. HEO-LEO scenario 7: HEO altitude influence; reference case for HEOH.
Figure 11. HEO-LEO scenario 8: material (AVCOAT) influence w.r.t. HEOH.
Figure 12. HEO-LEO scenario 9: e influence w.r.t. HEOH.
Figure 13. HEO-LEO scenario 10: Material (AVCOAT) + e influence w.r.t. HEOH.
Figure 14. HEO-LEO scenario 11: HEO altitude influence; reference case for HEOL.
Figure 15. HEO-LEO scenario 12: e influence w.r.t. HEOL.
5. Summary and Conclusions
Some aeroassisted transfers between coplanar orbits were analyzed to examine the influence of the aerodynamic characteristics of a spacecraft and the characteristics of its heat shield on the maneuver’s final performance.
The study was conducted with an original tool developed by the authors comprising highly representative models for thermal analysis, atmospheric flight dynamics, and an optimizer based on a genetic algorithm. The convenience of the maneuver was evaluated by comparison with the corresponding extra-atmospheric Hohmann transfer with two impulses.
The results of the case study indicated that there is a range of feasibility for the maneuver, depending on the first propulsive impulse of deorbit, within which there is, in turn, an interval in which the maneuver is more convenient than the homologous Hohmann transfer in terms of the total mass savings achieved (TPS and propellant). In particular, the results show that in terms of the convenience of the aeroassisted maneuver, the influence of the ablative material type used (low-density advantage) and of the aerodynamic efficiency is the most relevant, whereas other characteristics play a minor role. Significant benefits can be obtained from a careful thermal design, with a suitable thermal safety factor and high-performance adhesives for the bond-line.
For all cases, the highest convenience, and consequently the corresponding optimal trajectories, is obtained when it is not necessary to apply the second propulsion impulse. The most important future development for this type of analysis is increasing the number of cases analyzed in terms of materials, types of maneuvers and spacecraft characteristics.
Nomenclature
zero-lift drag coefficient
derivative of w.r.t. the angle of attack
maximum lift coefficient
fitness function
gravitational acceleration at sea level (m/s2)
altitude of the initial LEO (m)
altitude of the sensible atmosphere (m)
altitude of the final LEO (m)
propellant specific impulse (s)
induced drag factor
vehicle length (m)
TPS mass lost during atmospheric pass (kg)
final vehicle mass (kg)
“all propulsive” case final vehicle mass (kg)
vehicle mass at atmospheric entry (kg)
initial vehicle mass (kg)
vehicle mass at atmospheric exit (kg)
initial HEO radius (m)
radius of the sensible atmosphere (m)
final LEO radius (m)
reward factor for component “j”
Earth’s radius (m)
vehicle body radius (m)
vehicle reference surface (m2)
vehicle TPS total surface (m2)
bond-line temperature (K)
bond-line limit temperature (K)
thermal safety factor
velocity modulus (m/s)
circular orbit speed in initial HEO (m/s)
circular orbit speed in final LEO (m/s)
speed at atmospheric entry (m/s)
speed at atmospheric exit (m/s)
multiplicative weight of component “j”
vehicle wing cord (m)
vehicle wing span (m)
speed variation, impulse (m/s)
first impulse for Hohmann transfer (m/s)
second impulse for Hohmann transfer (m/s)
“all propulsive” (impulse) speed variation (m/s)
deorbit (impulse) speed variation (m/s)
minimum deorbit impulse (m/s)
boost (impulse) speed variation (m/s)
circularizing (impulse) speed variation (m/s)
angle of attack (rad)
aerodynamic efficiency
flight path angle at atmospheric entry (rad)
flight path angle at atmospheric exit (rad)
gravitational parameter (m3/s2)
bank angle (rad)
Subscripts
value at initial HEO
“all propulsive”
value at final LEO
bond-line
heat flux
value at atmospheric entry
lost
mass
Thermal Protection System
value at atmospheric exit
vehicle