Scramjet engines are key for sustained hypersonic flights. Analytic models play a critical role in the preliminary design of a scramjet engine configuration. The objective of this research is to develop and validate a quasi-1D model for the scramjet engine encompassing inlet, isolator and combustor, to evaluate the impact of flight conditions and design parameters on the engine functionality. The model is developed assuming isentropic flow in the inlet with a single turn; modified Fanno-flow equations in the isolator that account for the area change of the core flow; and the combustor is modeled using Heiser-Pratt equations accounting for the fuel mixing efficiency. The isolator and combustor models are validated against experimental results. The model accounts for twelve parameters allowing for a decent range of possible configurations. Finally, the model was applied to five sets of parametric studies to evaluate the effect of multiple parameters on the engine functionality.

A computational study of sonic light-gas jet injection into a supersonic cross flow was conducted. The scope of the numerical analysis encompassed many studies that affect how the flow-field is numerically modeled and the behavior, specifically mixing, of the flow-field itself. A single, round injector was used for the Baseline design. Simulated conditions involved sonic injection of helium heated to 313 K into a Mach 4 air cross-stream with average Reynolds number 5.77 e+7 per meter and a freestream momentum flux ratio of 2.1. Experiments at these conditions were available for comparison. The primary numerical flow solver employed was GASP v. 4.2. The Menter Shear Stress Transport (SST) turbulence model was used, since the algorithm has good capability of solving both wall-bounded and free-shear flows. The SST model was able to capture the mixing behavior of the complex flow-field. Important numerical parameters that affect the capabilities of the numerical solver were studied for the Baseline injector. These sensitivity studies varied the choice of turbulent Prandtl number, Schmidt number, freestream turbulence intensity, boundary layer size, steady and unsteady approaches and computational software packages. A decrease in the turbulent Prandtl number resulted in better mixing behavior of the prediction and better agreement with the experiment. An increase in the turbulent Schmidt number had a small adverse effect on the predictions. The mixing characteristics remained constant with an increase in freestream turbulence intensity. The best Baseline prediction was then compared to three different injector configurations: an aerodynamic ramp consisting of four injectors in an array, a diamond injector both aligned and yawed 15° to the oncoming flow. The Computational Fluid Dynamics (CFD) tools were more accurate compared to experiment in the prediction of the aeroramp injector than the diamond-shaped injectors. The aeroramp injector slightly improved mixing efficiency over the Baseline injector at these conditions. Both of the diamond-shaped injectors had similar mixing as the Baseline injector but did not predict significant improvement in penetration for the analyzed conditions. Additional studies involving the interaction of transverse injection with impinging oblique shock waves were performed. The impingement of a shock upon light gas jet injection increased mixing. The closer the shock is to the injection point, the larger the effect on mixing and vorticity. The last analyses involved a numerical comparison of a non-reacting model to a reacting hydrogen-air model. The reacting analysis prediction had an improved spreading rate and larger counter-rotating vortex pair with downstream distance over the non-reacting analysis. The mixing was not significantly altered by the addition of hydrogen-air reactions to the numerical equations. The numerical tools used are capable of reasonable accuracy in predicting the complex flow-field of jet injection into a supersonic freestream with proper choice of models and parameters. Numerical modeling offers a way to study the entire flow-field thoroughly in a cost and time efficient manner. / Ph. D.

One of the major technological innovations that are necessary for faster and cheaper access-to-space will be the commercial realization of supersonic combustion jet engines (SCRAMJET). The establishment of the flow through the inlet is one the prime requirement for the success of a SCRAMJET engine. The flow through a SCRAMJET inlet is dominated by inviscid /viscous coupling, transition, shock-shock interaction, shock boundary layer interaction, blunt leading edge effects and flow profile effects. Although the literature is exhaustive on various aspects of flow features associated with SCRAMJET engines, very little is known on the fundamental gasdynamic features dictating the flow establishment in the SCRAMJET inlet. On one hand we need the reduction of flight Mach number to manageable supersonic values inside the SCRAMJET combustor, but on the other hand we have to achieve this with minimum total pressure loss. Hence the dynamics of ramp/cowl shock interaction process ahead of the inlet has a direct bearing on the quality and type of flow inside the SCRAMJET engine. There is virtually no data base in the open literature focusing specifically on the cowl/ramp shock interactions at hypersonic Mach numbers. Hence in this backdrop, the main aim of the present investigation is to systematically understand the ramp/cowl shock interaction processes in front of a generic inlet model. Since we are primarily concerned with the shock interaction process ahead of the cowl all the investigations are carried out without any fuel injection. Variable geometry is necessary if we want to operate the inlet for a wide range of Mach numbers in actual flight. The investigation mainly comprises of three variable geometry configurations; namely, variation of contraction ratios at 00 cowl (CR 8.4, 5.0 and 4.3), variation of cowl length for a given chamber height (four lengths of cowls at 10 mm chamber height) and variation of cowl angle (three angles cowl each for two chamber heights). The change in cowl configuration results in different ramp/cowl shock interaction processes affecting the performance of the inlet. Experiments are performed in IISc hypersonic shock tunnel HST 2 (test time ~ 1 ms) at two nominal Mach numbers 8.0 and 5.74 for design and off-design testing conditions. Exhaustive numerical simulations are also performed to compliment the experiments. Further the effect of concentrated energy deposition on forebody /cowl shock interactions has also been investigated. A 2D, planar, single ramp scramjet inlet model has been designed and fabricated along with various cowl geometries and tested in a hypersonic shock tunnel to characterize the forebody/cowl shock interaction process for different inlet configurations. Further a DC plasma power unit and a plasma torch have been designed, developed and fabricated to serve as energy source for conducting flow-alteration experiments in the inlet model. The V-I characteristics of the plasma torch is studied and an estimation of plasma temperature is also performed as a part of characterizing the plasma flame. Initial standardization experiments of blunt body flow field alteration using the plasma torch and hence its drag reduction, are performed to check the torch’s suitability to be used as a flow-altering device in a shock tunnel. The plasma torch is integrated successfully with the inlet model in a shock tunnel to perform experiments with plasma jet as the energy source. The above experiments are first of its kind to be conducted in a shock tunnel. They are performed at various pressure ratios and supply currents. Time resolved schlieren flow visualization using Phantom 7.1 (Ms Vision Research USA) high speed camera, surface static pressure measurements inside a generic inlet using miniature kulite transducer and surface convective heat transfer rate measurements inside a generic inlet using platinum thin film sensors deposited on Macor substrate are some of the shock tunnel flow diagnostics that have been used in this study. Some of the important conclusions from the study are: • Experiments performed at different contraction ratios show different shock patterns. At CR 8.4, the SOL condition is satisfied, but the flow gets choked due to over contraction and flow through inlet is not established. For CR 5.0, formation of a small Mach stem before the chamber is observed with the reflection point on the cowl and the weak reflected shock entering inside the chamber. The Mach stem grows with time. For CR 4.3, the forebody/cowl shock interference created an Edney’s Type II shock interaction pattern. However, at off-design conditions, for CR 5 the shock reflection is regular and at CR 4.3, the Edney’s Type II pattern lasts for a short time. • For all lengths of cowl tested, 131mm and 141mm showed Edney’s Type II shock interference where as 151mm showed Edney’s Type I pattern at design condition. In all cases the flow is choked for high contraction ratio. At off-design condition these shock patterns do not last for the entire test time but rather it becomes a lambda pattern with the normal shock before the inlet. • For inlet configurations with cowl angle other than 00, the flow is found to be established for all cases at designed condition and except for 100 cowl at off-design condition. • For CR 8.4 the peak value of pressure (~1.7x104 Pa) occurs at a location of 151mm, where as for CR 5.0 and 4.3 they occur at 188mm and 206mm having values ~1.6x104 Pa and ~1.4x104 Pa respectively. These locations indicate the likely locations of shock impingements inside the chamber. • For cowl angle of 00 for a 10 mm chamber the maximum pressure value recorded is ~1.7x104 Pa whereas for 100 and 200 cowl it is ~1.1x104 Pa and 1.2 x104 Pa respectively. This is because in the first case the inlet is choked because of over contraction whereas in the other two cases the CR is less and flow is established inside the inlet. • The average heat transfer rates of last four heat transfer gauges (180 mm, 190 mm, 200 mm and 210 mm from the forebody tip) for all lengths of cowls tested are found to be almost same (~ 20 W/cm2). This is because the flow is choked in all these cases. The numerical simulation also shows uniform distribution here, consistent with the experimental findings. • The locations of heat transfer peaks for 100 cowl at design condition can be observed to be occurring at 170 mm and 200 mm from the forebody tip having values ~44 W/cm2 and ~39 W/cm2 respectively. For a 200 cowl they seem to be occurring at 170 mm and 180 mm from the forebody tip having values ~50 W/cm2 and ~30 W/cm2. These locations indicate the likely locations of shock impingements inside the chamber. With the evolution of concept of upstream fuel injection in recent times these may the most appropriate locations for fuel injection. • At higher jet pressure ratios the plasma jet/ramp shock interaction results in a lambda shock pattern with the triple point forming vertically above the cowl level. This means the normal shock stands in front of the inlet making a part of the flow entering the inlet subsonic. The reflected shock from the triple point also separates the ramp boundary layer. • At lower jet pressure ratios the triple point is formed below the cowl level and the flow entering inside the inlet is supersonic. The reflected shock interacts with the cowl shock and a weak separation shock is seen. • Experiments are performed with concentrated DC electric discharge as energy source. Even though the amount of energy dumped here is less than 0.25% of the total energy it creates a perceptible disturbance in the flow. • Experiments are also performed to see the effect of electric discharge as energy source on height of Mach stem for a given inlet configuration. Deposition of energy in the present location does not seem to alter the Mach stem height. However more experiments need to be performed by varying the energy location to see its effect. Non-intrusive energy sources like microwave and lasers can be thought of as options for depositing energy to study its effect on Mach stem height. Since they provide more flexibility on varying the location of energy the optimum location of energy can be found out for highest reduction of Mach stem height.

Rocket-Based Combined Cycle (RBCC) engines combine the best performance characteristics of air-breathing systems such as ramjets and scramjets with rockets with the goal of increasing payload/structure and propellant performance and thus making LEO more readily accessible. The idea of using RBCC engines for Single-Stage-To-Orbit (SSTO) trans-atmospheric acceleration is not new, but has been known for decades. Unfortunately, the availability of detailed models of RBCC engines is scarce. This thesis addresses the issue through the construction of an analytical performance model of an ejector rocket in a dual combustion propulsion system (ERIDANUS) RBCC engine. This performance model along with an atmospheric model, created using MATLAB was designed to be a preliminary `proof-of-concept' which provides details on the performance and behavior of an RBCC engine in the context of use during trans-atmospheric acceleration, and also to investigate the possibility of improving propellant performance above that of conventional rocket powered systems. ERIDANUS behaves as a thrust augmented rocket in low speed flight, as a ramjet in supersonic flight, a scramjet in hypersonic flight, and as a pure rocket near orbital speeds and altitudes. A simulation of the ERIDANUS RBCC engine's flight through the atmosphere in the presence of changing atmospheric conditions was performed. The performance code solves one-dimensional compressible flow equations while using the stream thrust control volume method at each station component (e.g. diffuser, burner, and nozzle) in all modes of operation to analyze the performance of the ERIDANUS RBCC engine. Plots of the performance metrics of interest including specific impulse, specific thrust, thrust specific fuel consumption, and overall efficiency were produced. These plots are used as a gage to measure the behavior of the ERIDANUS propulsion system as it accelerates towards LEO. A mission averaged specific impulse of 1080 seconds was calculated from the ERIDANUS code, reducing the required propellant mass to 65% of the gross lift off weight (GLOW), thus increasing the mass available for the payload and structure to 35% of the GLOW. Validation of the ERIDANUS RBCC concept was performed by comparing it with other known RBCC propulsion models. Good correlation exists between the ERIDANUS model and the other models. This indicates that the ERIDANUS RBCC is a viable candidate propulsion system for a one-stage trans-atmospheric accelerator.

The main objective of this dissertation is to develop large-eddy simulation (LES) based computational tools for supersonic inlet and combustor design. In the recent past, LES methodology has emerged as a viable tool for modeling turbulent combustion. LES computes the large scale mixing process accurately, thereby providing a better starting point for small-scale models that describe the combustion process. In fact, combustion models developed in the context of Reynolds-averaged Navier Stokes (RANS) equations exhibit better predictive capability when used in the LES framework. The development of a predictive computational tool based on LES will provide a significant boost to the design of scramjet engines. Although LES has been used widely in the simulation of subsonic turbulent flows, its application to high-speed flows has been hampered by a variety of modeling and numerical issues. In this work, we develop a comprehensive LES methodology for supersonic flows, focusing on the simulation of scramjet engine components. This work is divided into three sections. First, a robust compressible flow solver for a generalized high-speed flow configuration is developed. By using carefully designed numerical schemes, dissipative errors associated with discretization methods for high-speed flows are minimized. Multiblock and immersed boundary method are used to handle scramjet-specific geometries. Second, a new combustion model for compressible reactive flows is developed. Subsonic combustion models are not directly applicable in high-speed flows due to the coupling between the energy and velocity fields. Here, a probability density function (PDF) approach is developed for high-speed combustion. This method requires solution to a high dimensional PDF transport equation, which is achieved through a novel direct quadrature method of moments (DQMOM). The combustion model is validated using experiments on supersonic reacting flows. Finally, the LES methodology is used to study the inlet-isolator component of a dual-mode scramjet. The isolator is a critical component that maintains the compression shock structures required for stable combustor operation in ramjet mode. We simulate unsteady dynamics inside an experimental isolator, including the propagation of an unstart event that leads to loss of compression. Using a suite of simulations, the sensitivity of the results to LES models and numerical implementation is studied. / text

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.

No description available.

No description available.

The recent resurgence in hypersonics is centered around the development of SCRAMJET engine technology to power future hypersonic vehicles. Successful flight trials by Australian and American scientists have created interest in the scramjet engine research across the globe. To develop scramjet engine, it is important to study heat transfer effects on the engine performance and aerodynamic forces acting on the body. Hence, the main aim of present investigation is the design of scramjet engine configuration and measurement of aerodynamic forces acting on the model and heat transfer rates along the length of the combustor. The model is a two-dimensional single ramp model and is designed based on shock-on-lip (SOL) condition. Experiments are performed in IISc hypersonic shock tunnel HST2 at two different Mach numbers of 8 and 7 for different angles of attack. Aerodynamic forces measurements using three-component accelerometer force balance and heat transfer rates measurements using platinum thin film sensors deposited on Macor substrate are some of the shock tunnel flow diagnostics that have been used in this study.

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.