Combustion research is at the forefront of development of clean and efficient IC engines, gas turbines, rocket propulsion systems etc. With the advent of faster computers and parallel programming, computational studies of turbulent combustion is increasing rapidly. Many turbulent combustion models have been previously developed based on certain underlying assumptions. One of the major assumptions of the models is the regime it can be used for: either premixed or non-premixed combustion. However in reality, combustion systems are multi-regime in nature, i.e.,\ co-existence of premixed and non-premixed modes. Thus, there is a need for development of multi-regime combustion models which closely follows the physics of combustion phenomena. Much of previous modeling efforts for multi-regime combustion was done using flamelet-type models. As a first, the current study uses the highly robust transported Probability Density Function (PDF) method coupled with Large Eddy Simulation (LES) to develop a multi-regime model. The model performance is tested for Sydney Flame L, a piloted methane-air turbulent flame. The concept of flame index is used to detect the extent of premixed and non-premixed combustion modes. The drawbacks of using the traditional flame index definition in the context of PDF method are identified. Necessary refinements to this definition, which are based on the species gradient magnitudes, are proposed for the multi-regime model development. This results in identifying a new model parameter beta which defines a gradient threshold for the calculation of flame index. A parametric study is done to determine a suitable value for beta, using which the multi-regime model performance is assessed for Flame L by comparing it against the widely used non-premixed PDF model for three mixing models: Modified Curl (MCurl), Interaction by Exchange with Mean (IEM) and Euclidean Minimum Spanning Trees (EMST). The multi-regime model shows a significant improvement in prediction of mean scalar quantities compared to the non-premixed PDF model when MCurl mixing model is used. Similar improvements are observed in the multi-regime model when IEM and EMST mixing models are used. The results show potential foundation for further multi-regime model development using PDF model.

<p>Within the gas turbine industry, turbofan engines are widely implemented to enhance engine efficiency, specific thrust, and specific fuel consumption. However, these turbofans have yet to be widely implemented into microgas turbine engines. As turbofans become implemented into smaller engines, the need to design engine intakes for high-speed mission becomes more vital. In this work, a design procedure for compact, highly diffusive engine intakes for high subsonic speed applications is set about. The aerodynamic tradeoffs between cruise and takeoff flights are discussed and methods to enhance takeoff performance without negatively impacting high-speed cruise performance is discussed. Intake performance is integrated into overall engine analysis to help guide future mission analyses. Finally, an experimental model for engine intakes is developed for application to linear wind tunnels; allowing future designers to effectively validate numerical results.<br></p> <p>A multi-objective optimization routine is performed for compact engine intakes at a Mach number of 0.9. This optimization routine yielded a family of related curves that maximize intake diffusive capability and minimize intake pressure losses. Design recommendations to create such optimal intakes are discussed in this work so that future designers do not need to perform an optimization. Due to high diffusion rate of the intake, the intake performance at takeoff suffers greatly (as measured by massflow ingestion). Methods to enhance takeoff performance, from designing a variable geometry intake, to creating slots, to sliding intake components are evaluated and ranked for future designers to get an order of magnitude understanding of the types of massflow enhancements possible. Then, off-design performance of the intake is considered: with different Mach number flights, non-axial flow conditions, various altitudes, and unsteady engine operation considered. These off-design effects are evaluated to generate an intake map across a wide engine operational envelope. This map is then inputted into an engine model to generate a performance map of an engine; which allows for mission planning analysis. Finally, various methods to replicate intake flow physics in a linear wind tunnel are considered. It is shown that replicating diffuser curvature in a linear wind tunnel allows for best replication of flow physics. Additionally, a method to non-dimesnsionalize intake performance for application to a wind tunnel is developed. </p> <p>This work can be utilized by future engine intake designers in a variety of ways. The results shown here can help guide future designers create highly compact diffuser technology, capable of operating across a wide breadth of conditions. Methods to assess intake performance effects on overall engine performance are demonstrated; and an experimental approach to intake analysis is developed.</p>

In the aerospace industry, gas discharges have gained importance with the exploration of their performance and capabilities for flow control and combustion. Tunable properties of plasma make gas discharges efficient tools for various purposes. Since the scales of plasma and the available technology limit the knowledge gained from experimental studies, computational studies are essential to understand the results of experimental studies. The temporal and spatial scales of plasma also restrict the numerical studies. It is a necessity to use an idealized model, in which enough physics is captured, while the computational costs are acceptable.<br><br>In this work, numerical simulations of different low-pressure gas discharges are presented with a detailed analysis of the numerical approach. A one moment model is employed for DC glow discharges and nanosecond-pulse discharges. The cheap-est method regarding the modeling and simulation costs is chosen by checking the requirements of the fundamental processes of gas discharges. The verification of one-moment 1-D glow discharges with constant electron temperature variation is achieved by comparing other computational results.<br><br>The one moment model for pulse discharge simulation aims to capture the information from the experimental data for low-pressure argon discharges. Since the constant temperature assumption is crude, the local field approximation is investigated to obtain the data for electron temperature. It was observed that experimental data and computational data do not match because of the stagnant decay of electron number densities and temperatures. At the suggestion of the experimental group, water vapor was added as an impurity to the plasma chemistry. Although there was an improvement with the addition of water vapor, the results were still not in good agreement with experiment.<br><br>The applicability of the local field approximation was investigated, and non-local effects were included in the context of an averaged energy equation. A 0-D electron temperature equation was employed with the collision frequencies obtained from the local field approximation. It was observed that the shape of the decay profiles matched with the experimental data. The number densities; however, are less almost an order of magnitude.<br><br>As a final step, the two-moment model, one-moment model plus thermal electron energy equation, was solved to involve non-local effects. The two-moment model allows capturing of non-local effects and improves agreement with the experimental data. Overall, it was observed that non-local regions dominate low-pressure pulsed discharges. The local field approximation is not adequate to solve these types of discharges.

A body exposed to hypersonic flow is subjected to extremely high wall heating rates, owing to the conversion of the kinetic energy of the oncoming flow into heat through the formation of shock waves and viscous dissipation in the boundary layer and this is one of the main concerns in the design of any hypersonic vehicle. The conventional way of tackling this problem is to use a blunt fore-body, but it also results in an increase in wave drag and puts the penalty of excessive load on the propulsion system. An alternative approach is to alter the flow field using external means without changing the shape of the body; and several such methods are reported in the literature. The superiority of such methods lie in the fact that the effective shape of the body can be altered to meet the requirements of low wave drag, without having to pay the penalty of an increased wall heat transfer rate. Among these techniques, the use of local energy addition in the freestream to alter the flow field is particularly promising due to the flexibility it offers. By the suitable placement of the energy source relative to the body, this method can be effectively used to reduce the wave drag, to generate control forces and to optimise the performance of inlets. Although substantial number of numerical investigations on this topic is reported in the literature, there is no experimental evidence available, especially under hypersonic flow conditions, to support the feasibility of this concept. The purpose of this thesis is to experimentally investigate the effect of energy deposition on the flow-field of a 120� apex angle blunt cone in a hypersonic shock tunnel. Energy deposition is done using an electric arc discharge generated between two electrodes placed in the free stream and various parameters influencing the effectiveness of this technique are studied. The effect of energy deposition on aerodynamic parameters such as the drag force acting on the model and the wall heat flux has been investigated. In addition, the unsteady flow field is visualised using a standard Z-type schlieren flow visualisation setup. The experimental studies have shown a maximum reduction in drag of 50% and a reduction in stagnation point heating rate of 84% with the deposition of 0.3 kW of energy. The investigations also show that the location of energy deposition has a vital role in determining the flow structure; with no noticeable effects being produced in the flow field when the discharge source is located close to the body (0.416 times body diameter). In addition, the type of the test gas used is also found to have a major influence on the effectiveness of energy deposition, suggesting that thermal effects of energy deposition govern the flow field alteration mechanism. The freestream mass flux is also identified as an important parameter. These findings were also confirmed by surface pressure measurements. The experimental evidence also indicates that relaxation of the internal degrees of freedom play a major role in the determination of the flow structure. For the present experimental conditions, it has been observed that the flow field alteration is a result of the interaction of the heated region behind the energy spot with the blunt body shock wave. In addition to the experimental studies, numerical simulations of the flow field with energy deposition are also carried out and the experimentally measured aerodynamic drag with energy deposition is found to match reasonably well with the computed values.

The emerging and competitive environment in the space technology requires the improvements in the capability of aerodynamic vehicles. This leads to the analysis in drag reduction of the vehicle along with the minimized heat transfer rate. Using forward facing solid aerospike is the simplest way among the existing drag reduction methodologies for hypersonic blunt cone bodies. But the flow oscillations associated with this aerospike makes it difficult to implement. When analyzing this flow, it can be understood that this oscillating flow can be compared to conical cavity flow. Therefore in the spiked flows, it is decided to implement the technique used in reducing the flow oscillation of the cavities. Based on this method the shallow conical cavity flow generated by the aerospike fixed ahead of a 120o blunt cone body is fissured as multiple cavities by so many disks formed from 10o cone. Now the deep conical cavities had the length to mean depth ratio of unity; this suppresses the unnecessary oscillations of the shallow cavity. The total length of the telescopic aerospike is fixed as 100mm. And one another conical tip plain aerospike of same length is designed for comparing the telescopic spike’s performance at hypersonic flow Mach numbers of 5.75 and 7.9. A three component force balance system capable of measuring drag, lift and pitching moment is designed and mounted internally into the skirt of the model. Drag measurement is done for without spike, conical tip plain spiked and telescopic spiked blunt cone body. The three configurations are tested at different angles of attack from 0 to 10 degree with a step of 2. A discrete iterative deconvolution methodology is implemented in this research work for obtaining the clean drag history from the noisy drag accelerometer signal. The drag results showed the drag reduction when compared to the without spike blunt cone body. When comparing to the plain spiked, the telescopic spiked blunt cone body has lesser drag at higher angles of attack. Heat transfer measurements are done over the blunt cone surface using the Platinum thin film gauges formed over the Macor substrate. These results and the flow visualization give better understanding of the flow and the heat flux rate caused by the flow. The enhancement in the heat flux rate over the blunt cone surface is due to the shock interaction. And in recirculation region the heat flux rate is very much lesser when compared to without spike blunt cone body. It is observed that the shock interaction in the windward side is coming closer towards the nose of the blunt cone as the angle of attack increases and the oscillation of the oblique shock also decreases. Schlieren visualization showed that there is dispersion in the oblique shock, particularly in the leeward side. In the telescopic spike there are multiple shocks generated from each and every disk which coalesces together to form a single oblique shock. And the effect of the shock generated by the telescopic spike is stronger than the effect of the shock generated by the conical tip plain spike.

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.

A body in a hypersonic flow field will experience very high heating especially during re-entry. Conventionally this problem is tackled to some extent by the use of large angle blunt cones. At the cost of increased drag, the heat transfer rate is lower over most parts of the blunt body, except in a region around the stagnation point. Thus even with blunt cones, management of heat transfer rates and drag on bodies at hypersonic speeds continues to be an interesting research area. Various thermal protection systems have been proposed in the past, like heat sink cooling, ablation cooling and aerospikes. The ablative cooling system becomes extremely costly when reusability is the major concern. Also the shape change due to ablation can lead to issues with the vehicle control. The aerospikes themselves may become hot and ablate at hypersonic speeds. Hence an alternate form of cooling system is necessary for hypersonic flows, which is more feasible, cost effective and efficient than the conventional cooling systems. Injection of a mass of cold fluid into the boundary layer through the surface is one of the potential cooling techniques in the hypersonic flight corridors. These kinds of thermal protection systems are called mass transfer cooling systems. The injection of the mass may be through discrete slots or through a porous media. When the coolant is injected through a porous media over the entire surface, the coolant comes out as a continuous mass. Such a cooling system is also referred as “transpiration cooling system”. When the fluid is injected through discrete slots, the system is called as “film cooling system”. In either case, the coolant absorbs the incoming heat through its rise in enthalpy and thus modifies the boundary layer characteristics in such a way that the heat flow rate to the surface is less. Injection of a forward facing jet (opposite to the freestream direction) from the stagnation point of a blunt body can be used for mitigating both the aerodynamic drag and heat transfer rates at hypersonic Mach numbers. If the jet has enough momentum it can push the bow shock forward, resulting in reduced drag. This will also reduce heat transfer rate over most part of the body except around the jet re-attachment region. A reattachment shock impinging on the blunt body invariably increases the local heat flux. At lower momentum fluxes the forward facing jet cannot push the bow shock ahead of the blunt body and spreads easily over the boundary layer, resulting in reduced heat transfer rates. While the film cooling performance improves with mass flow rate of the jet, higher momentum flow rates can lead to a stronger reattachment leading to higher heat transfer rate at the reattachment zone. If we are able to reduce the momentum flux of the coolant for the same mass flow rate, the gas coming out can easily spread over the boundary layer and it is possible to improve the film cooling performance. In all the reported literature, the mass flow rate and the momentum flux are not varied independently. This means, if the mass flow rate is increased, there is a corresponding increase in the momentum flux. This is because the injection (from a particular orifice and for a particular coolant gas) is controlled only by the total pressure of injection and free stream conditions. The present investigation is mainly aimed at demonstrating the effect of reduction in momentum of the coolant (injected opposing a hypersonic freestream from the stagnation point of a blunt cone), keeping the mass flow rate the same, on the film cooling performance. This is achieved by splitting a single jet into a number of smaller jets of same injection area (for same injection total pressure and same free stream conditions). To the best of our knowledge there is no report on the use of forward facing micro-jet array for film cooling at hypersonic Mach numbers. In this backdrop the main objectives of the present study are: • To experimentally demonstrate the effect of splitting a single jet into an array of closely spaced smaller micro-jets of same effective area of injection (injected opposite to a hypersonic freestream from the stagnation zone), on the reduction in surface heat transfer rates on a large angle blunt cone. · Identifying various parameters that affect the flow phenomenon and doing a systematic investigation of the effect of the different parameters on the surface heat transfer rates and drag. Experimental investigations are carried out in the IISc hypersonic shock tunnel on the film cooling effectiveness. Coolant gas (nitrogen and helium) is injected opposing hypersonic freestream as a single jet (diameter 2 mm and 0.9 mm), and as an array of iv micro jets (diameter 300 micron each) of same effective area (corresponding to the respective single jet). The coolant gas is injected from the stagnation zone of a blunt cone model (58o apex angle and nose radius of 35 mm). Experiments are performed at a flow freestream Mach number of 5.9 at 0o angle of attack, with a stagnation enthalpy of 1.84 MJ/Kg, with and without injections. The ratios of the jet stagnation pressure to the pitot pressure (stagnation pressure ratio) used in the present study are 1.2 and 1.45. Surface convective heat transfer measurements using platinum thin film sensors, time resolved schlieren flow visualization and aerodynamic drag measurements using accelerometer force balance are used as flow diagnostics in the present study. The theoretical stagnation point heat transfer rate without injection for the given freestream conditions for the test model is 79 W/cm2 and the corresponding aerodynamic drag from Newtonian theory is 143 N. The measured drag value without injection (125 N) shows a reasonable match with theory. As the injection is from stagnation zone it is not possible to measure the surface heat transfer rates at the stagnation point. The sensors thus are placed from the nearest possible location from the stagnation point (from 16 mm from stagnation point on the surface). The sensors near the stagnation point measures a heat transfer rate of 65 W/cm2 on an average without any injection. Some of the important conclusions from the study are: • Up to 40% reduction in surface heat transfer rate has been measured near the stagnation point with the array of micro jets, nitrogen being the coolant, while the corresponding reduction was up to 30% for helium injection. Considering the single jet injection, near the stagnation point there is either no reduction in heat transfer rate or a slight increase up to 10%. · Far away from stagnation point the reduction in heat transfer with array of micro-jets is only slightly higher than corresponding single jet for the same pressure ratio. Thus the cooling performance of the array of closely spaced micro jets is better than the corresponding single jet almost over the entire surface. • The time resolved flow visualization studies show no major change in the shock standoff distance with the low momentum gas injection, indicating no major changes in other aerodynamic aspects such as drag. · The drag measurements also indicate that there is virtually no change in the overall aerodynamic drag with gas injection from the micro-orifice array. · The spreading of the jets injected from the closely spaced micro-orifice array over the surface is also seen in the visualization, indicating the absence of a region of strong reattachment. · The reduction in momentum flux of the injected mass due to the interaction between individual jets in the case of closely spaced micro-jet array appears to be the main reason for better performance when compared to a single jet. The thesis is organized in six chapters. The importance of film cooling at hypersonic speeds and the objectives of the investigation are concisely presented in Chapter 1. From the knowledge of the flow field with counter-flow injection obtained from the literature, the important variables governing the flow phenomena are organized as non-dimensional parameters using dimensional analysis in Chapter 2. The description of the shock tunnel facility, diagnostics and the test model used in the present study is given in Chapter 3. Chapter 4 describes the results of drag measurements and flow visualization studies. The heat transfer measurements and the observed trends in heat transfer rates with and without coolant injection are then discussed in detail in Chapter 5. Based on the obtained results the possible physical picture of the flow field is discussed in Chapter 6, followed by the important conclusions of the investigation.

In the flying field of space transportation domain, the increased efforts involving design and development of hypersonic flight for space missions is on toe to provide the optimum aerothermodynamic design data to satisfy mission requirements. Aerothermodynamics is the basis for designing and development of hypersonic space transportation flight vehicles such as X 51 a, and other programmes like planetary probes for Moon and Mars, and Earth re-entry vehicles such as SRE and space shuttle. It enables safe flying of aerospace vehicles, keeping other parameters optimum for structural and materials with thermal protection systems. In this context, the experimental investigations on hypersonic waverider are carried out at design Mach 6. The hypersonic waverider has high lift to drag ratio at design Mach number even at zero degree angle of incidence, and this seems to be one of the special characteristics for its shape at hypersonic flight regime. The heat transfer rates are measured using 30 thin film platinum gauges sputtered on a Macor material that are embedded on the test model. The waverider has 16 sensors on top surface and 14 on bottom surface of a model. The surface temperature history is directly converted to heat transfer rates. The heat transfer data are measured for design (Mach 6) and off-design Mach numbers (8) in the hypersonic shock tunnel, HST2. The results are obtained at stagnation enthalpy of ~ 2 MJ/kg, and Reynolds number range from 0.578 x 106 m-1 to 1.461 x 106 m-1. In addition, flow visualization is carried out by using Schlieren technique to obtain the shock structures and flow evolution around the Waverider. Some preliminary computational analyses are conducted using FLUENT 6.3 and HiFUN, which gave quantitative results. Experimentally measured surface heat flux data are compared with the computed one and both the data agree well. These detailed results are presented in the thesis.

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.

An approach for designing a high-enthalpy wind tunnel driven by exothermic chemical reactions was developed. Nozzle contours were designed using CONTUR, a program implementing the method of characteristics, to design nozzle contours at various flow conditions. A reacting mixture including nitrous oxide has been identified as the best candidate for providing clean air at high temperatures. The nitrous oxide has a few performance factors that were considered, specifically the combustion of the gas. Initial CFD simulations were performed on the nozzle and test region to validate flow characteristics and possible issues. Initial results show a fairly uniform exit velocity and ability to perform testing. In a second phase of the work, two generic, high-speed missile configurations were explored using numerical simulation. The mean flow was computed on both geometries at 0 and 45 roll and 0, 1, and 10 angle of attack. The computations identified complex flow structures, including three-dimensional shock/boundary-layer interactions, that varied considerably with angle of attack.