Abstract The main objective of this thesis was to develop a method than can be used to approximate the pressure forces on air vehicles traveling at hypersonic speed (Mach number > 5). The aerodynamic forces such as lift and drag were calculated from the pressure values on the surface of the airplane. Pitching moment was also tabulated. This work was initiated based on the idea of developing a flow solver proficient and capable of providing aerodynamic data (lift and drag look-up tables) for hypersonic air vehicles that can be fed to a flight simulator (used by the Aviation Systems Department) at the University of Tennessee Space Institute. Several approximation methods are used to solve hypersonic such as shock expansion method. Based on different studies, Computational Fluid Dynamic (CFD) proved to produce very accurate results; however, it is a difficult technique to use. In this thesis work Newtonian Method was adopted as a technique to approximate the aerodynamic forces and hence the performance of hypersonic airplanes, therefore, a computer program (Hyper-N) has been developed for aerodynamic analysis of three dimensional geometries airplane. The program is designed to read in a previously configured list of plates and compute the aerodynamic forces and moments for hypersonic free stream conditions. Programming was completed using MatLab language. The results obtained from the Hyper-N program were for the experimental airplane X-43A which were found to match the results when the shock expansion method is used for the same airplane, [1]. Because of the difficulties involve in using CFD or the complete Navier Stocks equation to obtain the aerodynamic forces on bodies traveling at hypersonic speeds, the Newtonian method is considered to be the most efficient technique to use for preliminary evaluation of the performance of hypersonic airplanes. Modified Newtonian theory and the computational requirement of the code are described. A number of geometric configurations, including the X-43A (experimental hypersonic) airplane, are provided as examples of applications of the Hyper-N program.

Tunable diode laser absorption spectroscopy using an external cavity diode laser operating in the infra-red has been developed to monitor CO2 in the freestream of the Longshot hypervelocity facility at the Von Karman Institute for Fluid Dynamics. The Longshot facility offers a unique European facility for ground testing and numerical validation applications, however, some of the traditional data rebuilding aspects are in question. A non-intrusive absorption sensor could significantly aid in improving the knowledge of freestream static values thereby improving the models used in data rebuilding and numerical simulation. The design of such a sensor also expands the spectroscopic capabilities of the Von Karman Institute. The absorption sensor is designed around the single P12 (00001)-(30013) rovibrational transition near 1.6µm (6218.09cm-1 specifically) which yields relatively weak direct absorption levels at about 3.5% per meter for typical Longshot freestream conditions. However, when handled carefully, adequate signal-to-noise can be acquired to exploit significant flow information. By being able to operate in this range, total sensor cost can be easily an a factor of two or more cheaper than sensors designed for the deeper infrared. All sensor elements were mounted to a compact portable optics bench utilizing single-mode optical fibers to allow for quick installation at different facilities by eliminating tedious optical realigning. Scans at 600Hz were performed over 20ms of the 40ms test time to extract core static temperature, pressure and velocity. These results are compared with the current state of the Longshot data rebuild method. The non-uniform flow properties of the shear layer and test cabin rested gas accumulation was of an initial concern. The temperature and density gradients along with significant radial velocity components could result in DLAS temperature, pressure and velocity that are significantly different than that of the target freestream inviscid core values. Fortunately, with the proper selection of the P12 rotational number, this effect could be more or less ignored as the higher temperature and lower density gas of this region is relatively transparent. Ultimately, acquired temperature and density were moderately accurate when compared to Longshot rebuilt results owing primarily to the baseline extraction which poses issues for such low absorption signals. However, the extracted velocity data are quite accurate. This is a definite puls for the sensor as the freestream enthalpy of cold hypersonic facilities is dictated primarily by the kinetic energy contribution. Being able to compare velocity gives insight to the level of vibration non-equilibrium in the flow. The velocity of the DLAS and the Longshot rebuild are quite close. This adds more weight to the argument that vibrational excitation is very low (if present at all) in the free stream and that the van de derWaals equation of state usage and constant specific heat assumption might be an adequate model for the data rebuild after all.

A parallel implicit adaptive-mesh-refinement scheme is proposed for the solution of the Navier-Stokes equations as applied to two-dimensional steady-state hypersonic laminar flows in conjunction with an equilibrium high-temperature equation of state. A finite-volume discretization is applied to the governing equations. Limited piecewise-linear solution reconstruction and Riemann solvers (Roe and HLLE, both modified for a general equation of state) are used to evaluate the inviscid fluxes. The gradients in the viscous fluxes are calculated using diamond-path reconstruction. The system of non-linear algebraic equations resulting from the finite-volume discretization are solved using an inexact Newton method with GMRES to solve the update step of the Newton method. GMRES is preconditioned with Schwarz preconditioning with local block-fill incomplete lower-upper factorization. Multigrid and pseudo-transient continuation are used for startup. Numerical results, including flows at Mach numbers of 7.0, are discussed and demonstrate the validity and efficiency of the scheme.

A parallel implicit adaptive-mesh-refinement scheme is proposed for the solution of the Navier-Stokes equations as applied to two-dimensional steady-state hypersonic laminar flows in conjunction with an equilibrium high-temperature equation of state. A finite-volume discretization is applied to the governing equations. Limited piecewise-linear solution reconstruction and Riemann solvers (Roe and HLLE, both modified for a general equation of state) are used to evaluate the inviscid fluxes. The gradients in the viscous fluxes are calculated using diamond-path reconstruction. The system of non-linear algebraic equations resulting from the finite-volume discretization are solved using an inexact Newton method with GMRES to solve the update step of the Newton method. GMRES is preconditioned with Schwarz preconditioning with local block-fill incomplete lower-upper factorization. Multigrid and pseudo-transient continuation are used for startup. Numerical results, including flows at Mach numbers of 7.0, are discussed and demonstrate the validity and efficiency of the scheme.

Ballute hypersonic aerodynamic decelerators have been considered for aerocapture since the early 1980's. Recent technology advances in fabric and polymer materials as well as analysis capabilities lend credibility to the potential of ballute aerocapture. The concept of the thin-film ballute for aerocapture shows the potential for large mass savings over propulsive orbit insertion or rigid aeroshell aerocapture. Several technology hurdles have been identified, including the effects of coupled fluid structure interaction on ballute performance and survivability. To date, no aeroelastic solutions of thin-film ballutes in an environment relevant to aerocapture have been published. In this investigation, an aeroelastic solution methodology is presented along with the analysis codes selected for each discipline. Variable-fidelity aerodynamic tools are used due to the long run times for computational fluid dynamics or direct simulation Monte Carlo analyses. The improved serial staggered method is used to couple the disciplinary analyses in a time-accurate manner, and direct node-matching is used for data transfer. In addition, an engineering approximation has been developed as an addition to modified Newtonian analysis to include the first-order effects of damping due to the fluid, providing a rapid dynamic aeroelastic analysis suitable for conceptual design. Static aeroelastic solutions of a clamped ballute on a Titan aerocapture trajectory are presented using non-linear analysis in a representative environment on a flexible structure. Grid convergence is demonstrated for both structural and aerodynamic models used in this analysis. Static deformed shape, drag and stress level are predicted at multiple points along the representative Titan aerocapture trajectory. Results are presented for verification and validation cases of the structural dynamics and simplified aerodynamics tools. Solutions match experiment and other validated codes well. Contributions of this research include the development of a tool for aeroelastic analysis of thin-film ballutes which is used to compute the first high-fidelity aeroelastic solutions of thin-film ballutes using inviscid perfect-gas aerodynamics. Additionally, an aerodynamics tool that implements an engineering estimate of hypersonic aerodynamics with a moving boundary condition is developed and used to determine the flutter point of a thin-film ballute on a Titan aerocapture trajectory.

A numerical study of the flow field generated by injection through diamondshaped orifices into a high-speed flow is presented in this document. Jet interaction flows have a wide range of applications in the field of engineering. These applications include the use of jets for fuel injection in scramjets, for reaction control of high-speed aerodynamic bodies and as cooling jets for skins of high-speed vehicles. A necessary requirement in the use of transverse jets for these and other applications is a thorough understanding of the physics of the interaction between the jet and freestream. This interaction generates numerous flow structures that include multiple shocks, vortices, recirculation regions and shear layers. This study involves diamond-shaped orifices that have the advantage of generating weaker or attached interaction shocks as compared to circular injectors. These injectors also negate the effects due to the recirculation region that is formed upstream of the injector. This study was undertaken in order to gain further understanding of the flow features generated by diamond-shaped injectors in a high-speed flow. Numerical simulations were performed using two different levels of turbulence models. Reynolds™ Averaged Navier-Stokes (RANS) simulations were performed using the GASP flow solver while Detached-Eddy Simulation (DES) runs were performed using the Cobalt flow solver. A total of fifteen diamond injector simulations were performed using the RANS model for a 15 half-angle diamond injector. The fifteen simulations spanned over five different injection angles and three jet total pressures. In addition to these, two circular injector simulations were also performed. In addition, low pressure normal injection through diamond and circular orifices simulations were performed using DES. Results obtained from CFD were compared to available experimental data. The resulting flow structure and the turbulent properties of the flow were examined in detail. The normal injection case through the diamond-shaped orifice at the lowest jet total pressure was defined as the baseline case and is presented in detail. In order to study the effect of different components of the vorticity transport equation, an in-house code was used post-process the results from the RANS runs.

During planetary entry, a vehicle uses drag generated from flight through the planetary atmosphere to decelerate from hyperbolic or orbital velocity. To date, all guided entry systems have utilized lift-modulation trajectory control. Deployable aerodynamic devices enable drag-modulation trajectory control, where a vehicle controls its energy and range during entry by varying drag area. Implementation of conventional lift-modulation systems is challenging for deployable systems. In contrast, drag-modulation trajectory control may be simpler and lower-cost than current state-of-the-art lift-modulation systems. In this investigation, a survey of analytical methods for computing planetary entry trajectories is presented and the approximate analytical solution to the entry equations of motion originally developed by Allen and Eggers is extended to enable flight performance evaluation of drag-modulation trajectory control systems. Results indicate that significant range control authority is available for vehicles with modestly sized decelerators. The extended Allen-Eggers solution is closed-form and enables rapid evaluation of nonlifting entry trajectories. The solution is utilized to develop analytical relationships for discrete-event drag-modulation systems. These relationships have direct application to onboard guidance and targeting systems. Numerical techniques were used to evaluate drag-modulation trajectory control for precision landing and planetary aerocapture missions, including development of prototype real-time guidance and targeting algorithms. Results show that simple, discrete-event drag-modulation trajectory control systems can provide landed accuracies competitive with the current state of the art and a more benign aerothermal environment during entry for robotic-scale exploration missions. For aerocapture, drag-modulation trajectory control is shown to be feasible for missions to Mars and Titan and the required delta-V for periapsis raise is insensitive to the particular method of drag modulation. Overall, results indicate that drag-modulation trajectory control is feasible for a subset of planetary entry and aerocapture missions. To facilitate intelligent system selection, a method is proposed for comparing lift and drag-modulation trajectory control schemes. This method applies nonlinear variational techniques to closed-form analytical solutions of the equations of motion, generating closed-form expressions for variations of arbitrary order. This comparative method is quantitative, performance-based, addresses robustness, and applicable early in the design process. This method is applied to steep planetary entry trajectories and shows that, in general, lift and drag-modulation systems exhibit similar responses to perturbations in environmental and initial state perturbations. However, significant differences are present for aerodynamic perturbations and results demonstrate that drag systems may be more robust to uncertainty in aerodynamic parameters. Finally, the results of these contributions are combined to build a set of guidelines for selecting lift or drag-modulation for a Mars Science Laboratory-class planetary entry mission.

The purpose of the current study was to answer several questions related to hypersonic, low Reynolds number, turbulent boundary layers, of which available data related to turbulence quantities is scarce. To that end, a unique research facility was created, instrumentation was developed to acquire data in the challenging low Reynolds number (low density) domain, and meaningful data was collected and analyzed. The low Reynolds number nature of the boundary layer (Re_theta = 3700) allows for tangible DNS computations/validations using the current geometry and conditions. The boundary layer examined in this experiment resembled other, higher Reynolds number boundary layers, but also exhibited its own unique characteristics. The Van Driest equivalent velocity scaling method was found to perform well, and the log layer of the law of the wall plot matched expected theory. Noticeably absent from the data was an overlap region between the two layers, which suggests a different profile for the velocity profiles at these low Reynolds number, hypersonic conditions. The low density effects near the wall may be having an effect on the turbulence that modifies this region in a manner not currently anticipated. The Crocco-Busemann relation was found to provide satisfactory results under its general assumptions. When compared to available data, the Morkovin scaled velocity fluctuations fell almost an order of magnitude short. Currently, it is not known if this deficit is due to inadequacies with the Strong Reynolds Analogy, or the Morkovin scaling parameters. The trips seem to promote uniformity across the span of the model, and the data seems to generally be in agreement across the spanwise stations. However, additional information is needed to determine if two-dimensional simulations are sufficient for these boundary layers. When the turbulent boundary layer power spectra is analyzed, the result is found to follow the traditional power law. This result verifies that even at low Reynolds numbers, the length scales still follow the behavior described by Kolmogorov. Moving downstream of the trips, the peak RMS disturbance value grows in amplitude until it reaches a critical value. After this point, the peak begins to decrease in amplitude, but the affected region spreads throughout the boundary layer. Once the influenced region covers a significant portion of the boundary layer, transition occurs.

Thesis (Ph. D.)--University of Washington, 1990. / Vita. Includes bibliographical references (leaves [100]-105).

Thesis (Ph. D.)--Virginia Polytechnic Institute and State University, 1993. / Vita. Abstract. Includes bibliographical references (leaves 162-167). Also available via the Internet.