One of the current goals of research in hypersonic, airbreathing propulsion is access to higher Mach numbers. A strong driver of this goal is the desire to integrate a scramjet engine into a transatmospheric vehicle airframe in order to improve performance to low Earth orbit (LEO) or the performance of a semi-global transport. An engine concept designed to access hypervelocity speeds in excess of Mach 10 is the shock-induced combustion ramjet (i.e. shcramjet). This dissertation presents numerical studies simulating the physics of a shcramjet vehicle traveling at hypervelocity speeds with the goal of understanding the physics of fuel injection, wall autoignition mitigation, and combustion instability in this flow regime. This research presents several unique contributions to the literature. First, different classes of injection are compared at the same flow conditions to evaluate their suitability for forebody injection. A novel comparison methodology is presented that allows for a technically defensible means of identifying outperforming concepts. Second, potential wall cooling schemes are identified and simulated in a parametric manner in order to identify promising autoignition mitigation methods. Finally, the presence of instabilities in the shock-induced combustion zone of the flowpath are assessed and the analysis of fundamental physics of blunt-body premixed, shock-induced combustion is accelerated through the reformulation of the Navier Stokes equations into a rapid analysis framework. The usefulness of such a framework for conducting parametric studies is demonstrated.

High Mach number, low-Reynolds number (high-altitude) reentry flowfield predictions are an important problem area in computational aerothermodynamics. Available numerical tools for handling such flows are very few and significantly limited in their applicability. A new implicit fully-iterative Parabolized Navier-Stokes (PNS) scheme is developed to accurately predict such low-Reynolds number flows. In this new approach the differential equations governing the conservation of mass, momentum and energy, and the algebraic equation of state for a perfect gas are solved simultaneously in a coupled manner. The idea is presented that by treating the governing equations in this manner (rather than eliminating the pressure terms in the governing equations by using appropriate differentiated forms of the equation of state) it may be possible to have an unconditionally time-like numerical scheme. The stability of a simplified version of this new PNS scheme is also studied, and it is demonstrated that these simplified equations are unconditionally time-like in the subsonic as well as the supersonic flow regions. A pseudo-time integration approach is used in addition to a new second-order accurate fully-implicit smoothing, to improve the efficiency of the solution algorithm. The new PNS scheme is used to predict the flowfield around a seven-deg sphere-cone vehicle under high- and low-Reynolds number conditions. Two test case, Case A and Case B, are chosen such that Case A has a large freestream Reynolds number (2.92x10⁵), whereas Case B has a freestream Reynolds number of 1.72x10³, which is smaller than the usual limit of applicability of the non-iterative PNS schemes (Re~10⁴ or larger). Comparisons are made with other available numerical schemes, and the results substantiate the stability, accuracy and efficiency claims of the new Parabolized Navier-Stokes scheme. / Ph. D.

To advance the design of hypersonic vehicles, high-fidelity multi-physics CFD is used to characterize 3-D scramjet flow-fields in two novel streamline traced configurations. The two inlets, Jaws and Scoop, are analyzed and compared to a traditional rectangular inlet used as a baseline for on/off-design conditions. The flight trajectory conditions selected are Mach 6 and a dynamic pressure of 1,500 psf (71.82 kPa). Analysis of these hypersonic inlets is performed to investigate distortion effects downstream with multiple single cavity combustors acting as flame holders, and several fuel injection strategies. The best integrated scramjet inlet/combustor design is identified. The flow physics is investigated and the integrated performance impact of the two innovative scramjet inlet designs is quantified. Frozen and finite rate chemistry is simulated with 13 gaseous species and 20 reactions for an Ethylene/air finite-rate chemical model. In addition, URANS and LES modeling are compared to explore overall flow structure and to contrast individual numerical methods. The flow distortion in Jaws and Scoop is similar to some of the distortion in the traditional rectangular inlet, despite design differences. The baseline and Jaws performance attributes are stronger than Scoop, but Jaws accomplishes this while eradicating the cowl lip interaction, and lessening the total drag and spillage penalties. The innovative inlets work best on-design, whereas for off-design, the traditional inlet is best. Early pressure losses and flow distortions in the isolator aid the mixing of air and fuel, and improve the overall efficiency of the system. Although the trends observed with and without chemical reactions are similar, the former yields roughly 10% higher mixing efficiency and upstream reactions are present. These show a significant impact on downstream development. Unsteadiness in the combustor increases the mixing efficiency, varying the flame anchoring and combustion pressure effects upstream of the step.

Traditionally, the design of complex aerospace systems requires iteration among segregated disciplines such as aerodynamic modeling and trajectory optimization. Multidisciplinary design optimization algorithms have been developed to efficiently orchestrate the interaction among these disciplines during the design process. For example, vehicle capability is generally obtained through sequential iteration among vehicle shape, aerodynamic performance, and trajectory optimization routines in which aerodynamic performance is obtained from large pre-computed tables that are a function of angle of attack, sideslip, and flight conditions. This numerical approach segregates advancements in vehicle shape design from advancements in trajectory optimization. This investigation advances the state-of-the-art in conceptual hypersonic aerodynamic analysis and trajectory optimization by removing the source of iteration between aerodynamic and trajectory analyses and capitalizing on fundamental linkages across hypersonic solutions. Analytic aerodynamic relations, like those derived in this investigation, are possible in any flow regime in which the flowfield can be accurately described analytically. These relations eliminate the large aerodynamic tables that contribute to the segregation of disciplinary advancements. Within the limits of Newtonian flow theory, many of the analytic expressions derived in this investigation provide exact solutions that eliminate the computational error of approximate methods widely used today while simultaneously improving computational performance. To address the mathematical limit of analytic solutions, additional relations are developed that fundamentally alter the manner in which Newtonian aerodynamics are calculated. The resulting aerodynamic expressions provide an analytic mapping of vehicle shape to trajectory performance. This analytic mapping collapses the traditional, segregated design environment into a single, unified, mathematical framework which enables fast, specialized trajectory optimization methods to be extended to also include vehicle shape. A rapid trajectory optimization methodology suitable for this new, mathematically integrated design environment is also developed by relying on the continuation of solutions found via indirect methods. Examples demonstrate that families of optimal hypersonic trajectories can be quickly constructed for varying trajectory parameters, vehicle shapes, atmospheric properties, and gravity models to support design space exploration, trade studies, and vehicle requirements definition. These results validate the hypothesis that many hypersonic trajectory solutions are connected through fast indirect optimization methods. The extension of this trajectory optimization methodology to include vehicle shape through the development of analytic hypersonic aerodynamic relations enables the construction of a unified mathematical framework to perform rapid, simultaneous hypersonic aerodynamic and trajectory optimization. Performance comparisons relative to state-of-the-art methodologies illustrate the computational advantages of this new, unified design environment.

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<p>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.<p><p>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.<p><p>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.<p><p>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. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

Laminar turbulent transition is known for a long time as a critical phenomenon influencing the thermal load encountered by hypersonic vehicle during their planetary re-entry trajectory. Despite the efforts made by several research laboratories all over the world, the prediction of transition remains inaccurate, leading to oversized thermal protection system and dramatic limitations of hypersonic vehicles performances. One of the reasons explaining the difficulties encountered in predicting transition is the wide variety of parameters playing a role in the phenomenon. Among these parameters, surface roughness is known to play a major role and has been investigated in the present thesis.<p><p>A wide bibliographic review describing the main parameters affecting transition and their coupling is proposed. The most popular roughness-induced transition predictions correlations are presented, insisting on the lack of physics included in these methods and the difficulties encountered in performing ground hypersonic transition experiments representative of real flight characteristics. This bibliographic review shows the importance of a better understanding of the physical phenomenon and of a wider experimental database, including real flight data, for the development of accurate prediction methods.<p><p>Based on the above conclusions, a hypersonic experimental test campaign is realized for the characterization of the flow field structure in the vicinity and in the wake of 3D roughness elements. This fundamental flat plate study is associated with numerical simulations for supporting the interpretation of experimental results and thus a better understanding of transition physics. Finally, a model is proposed in agreement with the wind tunnel observations and the bibliographic survey.<p><p>The second principal axis of the present study is the development of a hypersonic in-flight roughness-induced transition experiment in the frame of the European EXPERT program. These flight data, together with various wind tunnel measurements are very important for the development of a wide experimental database supporting the elaboration of future transition prediction methods. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

<p align="justify">Cette thèse traite de l'étude numérique et expérimentale du transport et de la déposition d'aérosols dans les poumons. La partie numérique du travail porte sur des simulations uni-, bi- et tridimensionnelles du comportement des aérosols dans la structure pulmonaire. Les simulations unidimensionnelles (1D) sont effectuées dans des modèles trompettes et multibranche similaires à ceux utilisés dans les études de transport et de mélange gazeux dans les poumons. Le dépôt total, le profil des dépôts le long des différentes générations de l'arbre bronchique ainsi que la dispersion de boli d'aérosols sont calculés en fonction de la taille des particules et du protocole respiratoire. Un bolus consiste en un faible volume d'aérosols inhalé sous la forme d'un pic de concentration au cours d'une inspiration d'air pur. Les résultats montrent les limitations intrinsèques liées aux modèles 1D quant à la description du transport des aérosols dans les poumons et suggèrent l'utilisation d'équations multidimensionnelles pour décrire le transport de particules. Des simulations bidimensionnelles (2D) sont alors développées pour décrire le comportement des aérosols dans un modèle représentatif de la zone alvéolaire du poumon humain. Les simulations montrent que les particules ne se déposent pas uniformément sur les parois alvéolaires des conduits mais qu'elles sont principalement localisées près de l'entrée des alvéoles et ceci principalement dans le cas de petites particules (diamètre inférieure à 0.5 mm). De plus, les résultats montrent que le traditionnel coefficient de dispersion utilisé dans l'approche unidimensionnelle ne peut pas être extrapolé dans la zone alvéolaire du poumon.</p><p><p align="justify">Finalement, des simulations tridimensionnelles (3D) sont réalisées dans un modèle d'un conduit pulmonaire entouré d'alvéoles et confirment la déposition largement hétérogène des aérosols calculée dans l'étude bidimensionnelle suggérant que les concentrations locales et moyennes en aérosols peuvent être substantiellement différentes.</p><p><p align="justify">Parallèlement, des données expérimentales de déposition totale et de dispersion de boli d'aérosols sont obtenues et comparées aux résultats numériques. Des indices tels que la dispersion du bolus expiré, la déposition totale ou le déplacement du mode entre les courbes de concentration des boli inspiré et expiré mesurés au niveau de la bouche ont été évalués. Des simulations numériques similaire aux tests expérimentaux sont également effectuées. Bien qu'une approche relativement simplifiée soit utilisée, il apparaît que les simulations décrivent raisonnablement bien les résultats expérimentaux.</p><p> / Doctorat en sciences appliquées / info:eu-repo/semantics/nonPublished

<p>  </p> <p>As space agencies look to conduct more scientific missions beyond Earth orbit, low-cost access to space becomes indispensable. Small satellites (smallsats) fulfill this need as they can be developed at a fraction of the cost of traditional large satellites. Consequently, smallsats are being envisioned for planetary science missions at several destinations including Mars. However, a significant challenge for interplanetary smallsats is performing fully-propulsive orbit insertion because modern smallsat propulsion technologies have limited total velocity change capabilities. At destinations with significant atmospheres, this challenge can be circumvented via <em>aerocapture</em>, a technique that uses a single atmospheric pass to convert a hyperbolic approach trajectory into a captured elliptical orbit. Aerocapture has been shown to enable significant propellant mass savings as compared to fully-propulsive orbit insertion, making it an attractive choice for smallsats. Performing aerocapture with smallsats requires a suitable vehicle design that satisfies the associated control requirements and volumetric constraints. To address this requirement, this dissertation proposes the <em>morphable entry system </em>(MES), a conceptual deployable entry vehicle that utilizes shape morphing to follow a desired atmospheric flight profile during aerocapture. The aerocapture performance of the MES at Mars is investigated using a six degree-of-freedom aerocapture simulation environment. The shape morphing strategy employed by the MES is shown to be feasible for targeting desired angle of attack and sideslip angle profiles that lead to successful orbit captures. Furthermore, the robustness of the MES to simulated day-of-flight uncertainties while employing angle of attack control is demonstrated through a Monte Carlo dispersion analysis. The major contributions of this research as well as areas of future work are described.</p>