A unique implementation of commonly used engineering methods for estimating supersonic and hypersonic aerodynamics of arbitrary shapes is discussed. A validation case comparing the aerodynamic coefficients generated by this technique with both wind tunnel and published computational results for the NASA HL-20 configuration from Mach 1.2 to Mach 10 is presented. These results indicate that the method can generate aerodynamic databases where the lift and drag coefficients are within 15% of experimental data in minutes on a typical workstation. Therefore, the benefits of this rapid aerodynamic analysis approach are further investigated by integrating it into a multidisciplinary design optimization framework for waverider-derived crew reentry vehicles. A shape design code, mass estimating model, and trajectory/aeroheating analysis are linked using a genetic algorithm optimization process to generate a pareto front comparing vehicle downrange versus landed mass. The results of this analysis indicate that a higher lift-to-drag ratio (L/D) typically results in a longer range, but also correlates with a larger vehicle mass. Finally, a comparison with the HL-20 is performed to highlight how the waverider-derived design approach can produce configurations exhibiting a 20% increase in maximum L/D with only a 4% increase in vehicle length.
Nomenclatureܽ = speed of sound ܥ = pressure coefficient ℎ = altitude ܦ/ܮ = Lift-to-drag ratio ܯ = Mach number ݊ ො = surface normal vector = pressure ݎ = radial coordinate in Taylor-Maccoll equation ܴ ௌ = waverider shock profile curve radius for cone region in base plane ܴ݁ = Length-based Reynolds number ݑ = radial velocity in Taylor-Maccoll equation ݒ = angular velocity in Taylor-Maccoll equation ̅ݔ = axial center of gravity location from vehicle nose (non-dimensionalized to vehicle length) ݕ = waverider lower surface base curve control point (vertical) in base plane ݖ = waverider lower surface base curve control point (horizontal) in base plane ܸ ሬԦ = velocity vector ߚ = shock wave angle ߛ = ratio of specific heats 1 Senior Project Engineer, Launch Strike and Range Department, 2310 E. El Segundo Blvd., Mail Stop M1-557; AIAA Member. Downloaded by UNIVERSITY OF NEW SOUTH WALES on August 23, 2015 | http://arc.aiaa.org | 2 ߥ = Prandtl-Meyer function value ߩ = density ߠ = surface inclination angle; angular coordinate in Taylor-Maccoll equation Subscripts 02 = station number for stagnation point behind a normal shock wave ݅ = surface panel index = streamline segment index ∞ = freestream