The purpose of this thesis is to perform numerical solutions of hypersonic, high temperature, perfect gas flows over various geometries. Three dimensional, thin layer, compressible, Navier-Stokes equations are solved. An upwind finite difference approach with Lower Upper-Alternating Direction Implicit (LU-ADI) decomposition is used. Solutions of laminar, hypersonic, high temperature, perfect gas flows over flat plate and compression corners (qw=5&deg / , 10&deg / , 14&deg / , 15&deg / , 16&deg / , 18&deg / and 24&deg / ) with eight different free-stream and wall conditions are presented and discussed. During the analysis, air viscosity is calculated from the Sutherland formula up to 1000&deg / K, for the temperature range between 1000 &ordm / K and 5000 &ordm / K a curve fit to the estimations of Svehla is applied. The effects of Tw/T0 on heat transfer rates, surface pressure distributions and boundary layer characteristics are studied. The effects of corner angle (&amp / #952 / w) on strong shock wave/boundary layer interactions with extended separated regions are investigated. The obtained results are compared with the available experimental data, computational results, and theory.

This study is comprised of two separate parts: (1) modeling thermochemical nonequilibrium processes, and (2) flow field simulations of spark-induced plasma. In the first part, the methodology and literature for modeling thermochemical nonequilibrium processes in partially ionized air is presented and implemented in a zero-dimensional solver, termed as NEQZD. The solver was verified for a purely reacting flow case as well as two thermochemical nonequilibrium flow cases. A three-temperature electron-electronic model for thermochemical nonequilibrium partially ionizing air mixture was implemented and demonstrated the ability to capture additional physics compared to the legacy two-temperature model through the inclusion of electronic energy nonequilibrium. In the second part of this work, full scale axisymmetric simulations of the flow field produced by the abrupt heat release of spark-induced plasma were presented and analyzed for two electrode configurations. The heat release was modeled based on data from experiments and assumed that all electrical power supplied to the electrodes is converted to thermal energy. It was found that steeper electrode walls lead to a greater region of hot gas, a stronger shock front, and slightly larger vortices.

Thermochemical nonequilibrium widely exists in supersonic combustion, cold plasma and hypersonic flight. The effect can influence heat transfer, surface ablation and aerodynamic loads. One distinct feature of it is the coupling between internal energy excitation and chemical reactions, particularly the vibration-dissociation coupling. The widely used models are empirical and calibrated based on limited experimental data. Advances in theories and computational power have made the first-principle calculation of thermal nonequilibrium reaction rates by methods like quasi-classical trajectory (QCT) almost a routine today. However, the approach is limited by the uncertainties and availability of potential energy surfaces. To the best of our knowledge, there is no study of thermal nonequilibrium transport properties with this approach. Most importantly, non-trivial effort is required to process the QCT data and implement it in flow simulation methods. In this context, the first part of this work establishes the approach to compute transport properties by the QCT method and studies the influence of thermal nonequilibrium on transport properties for N₂-O molecules. The preponderance of the work is the second part, a comprehensive study of the development of a new thermal nonequilibrium reaction model based on reasonable assumptions and approximations. The new model is as convenient as empirical models. By validating against recent QCT data and experimental results, we found the new model can predict nonequilibrium characteristics of dissociation reactions with nearly the same accuracy as QCT calculations do. In general, the results show the potential of the new model to be used as the standard dissociation model for the simulation of thermochemical nonequilibrium flows.

Master’s thesis deals with aerodynamic heating of launch vehicle during ascent phase by using CFD simulation. Ascent trajectory and payload fairing geometry is design using data of existing small launch vehicles. Critical flight regimes are identified using 2D calculations, and in these regimes analysis is performed by axially symmetric simulations. Simulation results are compared to values obtained from theoretical and semi-empirical calculations.

No description available.

No description available.

Boundary-layer transition from a laminar to a turbulent regime is a critical driver in the design of high-speed vehicles. The aerothermodynamic loads associated with transitional or fully turbulent hypersonic boundary layers are several times higher than those associated with laminar flow. The presence of isolated roughness elements on the surface of a body can accelerate the growth of incoming disturbances and introduce additional instability mechanisms in the flow field, eventually leading to a premature occurrence of transition. This dissertation studies the instabilities induced by three-dimensional discrete roughness elements located inside a high-speed boundary layer developing on a flat plate. Two-dimensional local linear stability theory (2D-LST) is employed to identify the instabilities evolving in the three-dimensional flow field that characterizes the wake induced by the roughness elements and to investigate their evolution downstream. A formulation of the disturbance energy evolution equation available for base flows depending on a single spatial direction is generalized for the first time to base flows featuring two inhomogeneous directions and perturbations depending on three spatial directions. This generalization allows to obtain a decomposition of the temporal growth rate of 2D-LST instabilities into the different contributions that lead to the production and dissipation of the total disturbance energy. This novel extension of the formulation provides an additional layer of information for understanding the energy exchange mechanisms between a three-dimensional base flow and the perturbations resulting from 2D-LST. Stability computations for a calorically perfect gas illustrate that the wake induced by the roughness elements supports the growth of different sinuous and varicose instabilities which coexist together with the Mack-mode perturbations that evolve in the flat-plate boundary layer, and which become modulated by the roughness-element wake. A single pair of sinuous and varicose disturbances is found to dominate the wake instability in the vicinity of the obstacles. The application of the newly developed decomposition of the temporal growth rate reveals that the roughness-induced wake modes extract most of their potential energy from the transport of entropy fluctuations across the base-flow temperature gradients and most of their kinetic energy from the work of the disturbance Reynolds stresses against the base-flow velocity gradients. Further downstream, the growth rate of the wake instabilities is found to be influenced by the presence of Mack-mode disturbances developing on the flat plate. Strong evidence is observed of a continuous synchronization mechanism between the wake instabilities and the Mack-mode perturbations. This phenomenon leads to an enhancement of the amplification rate of the wake modes far downstream of the roughness element, ultimately increasing the associated integrated amplification factors for some of the investigated conditions. The effects of vibrational molecular excitation and chemical non-equilibrium on the instabilities induced by a roughness element are studied for the case of a high-temperature boundary layer developing on a sharp wedge configuration. For this purpose, a 2D-LST solver for chemical non-equilibrium flows is developed for the first time, featuring a fully consistent implementation of the thermal and transport models employed for the base flow and the perturbation fields. This is achieved thanks to the automatic derivation and implementation tool (ADIT) available within the von Karman Institute extensible stability and transition analysis (VESTA) tool-kit, which enables an automatic derivation and implementation of the 2D-LST governing equations for different thermodynamic flow assumptions and models. The stability computations for this configuration show that sinuous and varicose disturbances also dominate the wake instability in the presence of vibrational molecular energy mode excitation and chemical reactions. The resulting base-flow cooling associated with the modeling of such high-temperature phenomena is found to have opposite stabilizing and destabilizing effects on the streamwise evolution of the sinuous and varicose instabilities. The modeling of vibrational excitation and chemical non-equilibrium acting exclusively on the perturbations is found to have a stabilizing influence in all cases. / Doctorat en Sciences de l'ingénieur et technologie / info:eu-repo/semantics/nonPublished

Aerodynamic cavities are common features on hypersonic vehicles which are caused in both large and small scale features like surface defects, pitting, gap in joints etc. In the hypersonic regime, the presence of such cavities alters the flow phenomenon considerably and heating rates adjacent to the discontinuities can be greatly enhanced due to the diversion of flow. Since the 1960s, a great deal of theoretical and experimental research has been carried out on cavity flow physics and heating. However, most of the studies have been done to characterize the effect downstream and within the cavity. In the present study, a series of were carried out in the shock tunnel to investigate the heating characteristics, upstream and on the lateral side of the cavity. Heat flux measurement has been done using indigenously developed high resistance platinum thin film gauges. High resistance gauges, as contrary to the conventionally used low resistance gauges were showing good response to the extremely low heat flux values on a flat plate with sharp leading edge. The experimental measurements of heat done on a flat plate with sharp leading edge using these gauges show good match with theoretical relation by Crabtree et al. Flow visualization using high speed camera with the cavity model and shock structures visualized were similar to reported in supersonic cavity flow. This also goes to state that in spite of the fluctuating shear layer-the main feature of hypersonic flow over a cavity ,reasonable studies can be done within the short test time of shock tunnel. Numerical Simulations by solving the Navier-Stokes equation, using the commercially available CFD package FLUENT 13.0.0 has been done to complement the experimental studies.

Despite years of research in high speed boundary layer flow, there is still a need for insightful experiments to realize key features of the flow like boundary layer response to different conditions and related transition mechanisms. Volumes of data on the these problems point to the fact that there is still much to be understood about the nature of boundary layer instability causing transition and growth of boundary layer in different conditions. Boundary layer stability experiments have been found to be more useful, in which the boundary layer is perturbed and its behavior observed to infer useful conclusions. Also, apart from the stability part, the effect of various changes in boundary layer due to the perturbation makes interesting observation to gain more insight into the understood and the not so understood facets of the same. In view of the above, the effect of a steady axisymmetric thermal bump is investigated on a hypersonic boundary layer over a 60º sharp cone cylinder model. The thermal bump, placed near tip of the cone, perturbs the boundary layer, the behavior of which is observed by recording the wall heat flux on the cone and cylinder surface using platinum thin film sensors. The state of the boundary layer is qualitatively assessed by the wall heat flux comparisons between laminar and turbulent values. The same thermal bump also acts as a heat addition source to boundary layer in which case this recorded data provides a look into the effect of the heat addition to the wall heat flux. To gain a larger view of heat addition causing changes to the flow, effects of change in enthalpy are also considered. Experiments are performed in the IISc HST2 shock tunnel facility at 2MJkg−1 stag-nation enthalpy and Mach number of 8,with and without the thermal bump to form comparisons. Some experiments are also performed in the IISc HST3 free piston driven shock tunnel facility at 6MJkg−1, to investigate the effect of change in stagnation enthalpy on the wall heat flux. To support the experimental results theoretical comparisons and computational studies have also been carried out. The results of experiments show that the laminar boundary layer over the whole model remains laminar even when perturbed by the thermal bump. The wall heat flux measurements show change on the cone part where there seems to be fluctuation in the temperature gradients caused by the thermal bump, which decrease at first and then show an increase towards the base of the cone. The cylinder part remains the same with and without the thermal bump, indicating heavy damping effects by the expansion fan at cone cylinder junction. A local peak in wall heat flux is observed at the junction which is reduced by 64% by the action of the thermal bump. The possible reason for this is attributed to the increased temperature gradients at the wall due to delayed dissipation of heat that is accumulated in the boundary layer as a result of the thermal bump action. The comparison of data for enthalpies of 2MJkg−1 and 6MJkg−1 show that there are negligible real gas effects in the higher enthalpy case and they do not affect the wall heat flux much. Also it is found that the thermal bump fails to dump heat into the flow directly though it creates heat addition virtually by mere discontinuity in the surface temperature and causes temperature gradients fluctuation in the boundary layer. Considering the thermal bump action and the change in stagnation enthalpy of the flow, there seems to be no change in both cases that can be attributed to a common observation resulting from the factor of change in heat inside the boundary layer.

The growing requirements to achieve hypersonic flights, as in the case of reentry vehicles, pose a serious challenge to the designers. This demands an understanding of the features of hypersonic flow and its effect on hypersonic vehicles. Hypersonic shock tunnels are one of the most widely used facilities for the purpose of obtaining valuable design data by conducting experiments on scaled down models. They are operated by conventional shock tubes by rupturing metal diaphragms placed between the driver and driven sections of the shock tube. Shock tunnels are being extensively used in spite of some of the drawbacks they possess. Due to the varying nature of metal diaphragm rupture, reproducibility of the experiment results is difficult to obtain. Damage to model and inner surface of the shock tube can happen when the diaphragm petal breaks away from the diaphragm. Lastly the time consuming diaphragm replacement process is not desired in applications which require quick loading of shock waves on the specimen. All these disadvantages call for the replacement of the diaphragm mode of operation with a diaphragmless mode of operation for the generation of shock waves. The main objective of the present study is to design and demonstrate the working of a diaphragmless hypersonic shock tunnel. The motivation for the present study comes from the fact that the diaphragmless operation of a shock tunnel has not been reported so far in the open literature. All the research works carried out deal with diaphragmless drivers operating only a shock tube. In the present work, the conventional metal diaphragm is substituted by fast acting pneumatic valves which serve the purpose of quickly opening the driven section of the shock tube to allow the driver gas to rush in, resulting in the formation of a shock wave. To design a diaphragmless driver, a detailed study of the shock formation process is accomplished which helps in understanding the effect of valve opening time on the shock formation distance. Also the theoretical basis for the design of a pneumatic cylinder is understood. Following the theoretical studies, three types of diaphragmless drivers are designed and tested. The first setup incorporates a rubber membrane, which acts as a valve. The rubber membrane when bulged closes the mouth of the driven section and on retraction the driven section is opened to the driver gas. The second and the third setups utilise two different types of double acting pneumatic cylinders. Experimental results of the three diaphragmless drivers operating a shock tube are analysed and compared with the ideal shock tube theory. Better repeatability in terms of shock Mach number is shown with all three diaphragmless shock tubes when compared with a conventionally operated shock tube. Finally, the best among the three systems is identified to operate the hypersonic shock tunnel 2 (HST2) facility of the Shock Waves laboratory, IISc. Demonstration of the working of the diaphragmless shock tunnel is shown by performing heat transfer measurements on a 3 mm backward facing step flat plate model. The experimental results are compared with those obtained in a conventional shock tunnel. CFD studies on diaphragmless shock tube model are done to have an idea on the flow in the shock tube there by identifying the shock formation distance. ANSYS-CFX package is used for this purpose. Further, results from the numerical simulation of hypersonic flow over the backward facing step model are compared with the experimental results thus validating the code.