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.

Αντικείμενο της παρούσας διδακτορικής διατριβής αποτελεί η θεωρητική και υπολογιστική προσομοίωση της πτήσης ατμοσφαιρικής καθόδου υψηλών ταχυτήτων, με ιδιαίτερη έμφαση στις υπερυπερηχητικές ταχύτητες και την αεροδυναμική και θερμοδυναμική ανάλυσή της. Η μεθοδολογία που ακολουθείται έχει ως κύριο σκοπό την πρόβλεψη του ίχνους επιστροφής υπερυπερηχητικών οχημάτων στη Γη, με δεδομένο οι ταχύτητες πτήσης να φτάνουν τα 15000 m/s (αριθμός Mach ~ 40) και οι θερμοκρασίες ανακοπής τους 15000 Κ. Στην εργασία αυτή ακολουθείται μεθοδολογική προσέγγιση και συστηματική βιβλιογραφική έρευνα στον επιστημονικοτεχνικό τομέα της αεροδυναμικής υπερθέρμανσης, που υφίστανται Υπερυπερηχητικά Οχήματα (Υ/Ο), κατά την πτήση τους στην ατμόσφαιρα. Η ανάλυση της υπερυπερηχητικής ροής, των θερμικών φορτίων και οι μέθοδοι θερμοθωράκισης των Υ/Ο αναπτύσσονται τις τελευταίες δεκαετίες και εξελίσσονται συνεχώς, αλλά μέχρι και σήμερα (2008), προκύπτουν αστοχίες, κυρίως κατά την επανείσοδό τους στην ατμόσφαιρα από το διάστημα. Για την προσομοίωση της υπερυπερηχητικής πτήσης καθόδου στην ατμόσφαιρα γίνεται μελέτη και ανάλυση των θερμοχημικών ιδιοτήτων και παραμέτρων, οι οποίες χαρακτηρίζουν τα ροϊκά πεδία στις υψηλές ταχύτητες και υψηλές θερμοκρασίες. Ο προσδιορισμός των μεταβαλλόμενων ροϊκών συνθηκών, των θερμοδυναμικών ιδιοτήτων και μεγεθών μεταφοράς θα συμβάλλει στις θεωρίες μοντελοποίησης και υπολογισμού της υψηλής θερμικής καταπόνησης, της σύστασης του υπέρθερμου αέρα και της υψηλής θερμοκρασίας του στην εξωτερική επιφάνεια του θερμομονωτικού τοιχώματος οχημάτων υψηλών ταχυτήτων, κατά την πτήση τους στην ατμόσφαιρα της Γης. Στα πλαίσια αυτά αναπτύσσονται υπολογιστικοί αλγόριθμοι προσομοίωσης και επιτυχούς πρόβλεψης του ατμοσφαιρικού ίχνους καθόδου οχημάτων υπερυπερηχητικής μεταφοράς. Τα βασικά αεροδυναμικά χαρακτηριστικά των ανωστικών οχημάτων στις υψηλές ταχύτητες υπολογίζονται με κατάλληλη προσαρμογή της κλασσικής Νευτώνειας Ροής στις υψηλές ταχύτητες και επιβεβαιώνονται με πειραματικά δεδομένα αεροδυναμικών μετρήσεων σε αεροσήραγγες υψηλών ταχυτήτων της NASA. Επιπλέον στη διατριβή ενσωματώνονται μοντέλα προσομοίωσης των ιδιοτήτων συμπεριφοράς του υπέρθερμου αέρα ως πραγματικού αερίου σε συνθήκες θερμοδυναμικής ισορροπίας σε ακραίες ροϊκές συνθήκες (υψηλές ταχύτητες, υψηλές πιέσεις και χαμηλές πυκνότητες). Πιο συγκεκριμένα γίνεται ανάλυση της μεταβαλλόμενης σύστασης και των ιδιοτήτων του αέρα (θερμοδυναμικά μεγέθη και μεγέθη μεταφοράς) στις πολύ υψηλές θερμοκρασίες του υπέρθερμου ροϊκού πεδίου. Η μελέτη αυτή θα οδηγήσει στην εξαγωγή ενός νέου μοντέλου υπολογισμού της μεταφερόμενης θερμικής ροής στην περιοχή ανακοπής υπερυπερηχητικών οχημάτων, το οποίο επιβεβαιώνεται τόσο με διεθνώς αναγνωρισμένα θερμικά μοντέλα (όπως αυτό των Fay-Riddell) όσο και με υπολογιστικές προσομοιώσεις άλλων ερευνητών. Η μεθοδολογία που ακολουθείται θα αποτελέσει τη βάση για την εξαγωγή ενός νέου θεωρητικού μοντέλου υπολογισμού της θερμοκρασίας του αέρα στην επιφάνεια του θερμοθώρακα του υπερυπερηχητικού οχήματος. Η αποτελεσματικότητα της προτεινόμενης μεθοδολογίας, με το συνδυασμό κινηματικής και αεροθερμοδυναμικής ανάλυσης σε ενιαίο υπολογιστικό αλγόριθμο, αποδεικνύεται και ενισχύεται μέσω της εφαρμογής των ανωτέρω μοντέλων προσομοίωσης σε βαλλιστικό (διαστημικός θαλαμίσκος Απόλλων) και ανωστικό (Σύστημα Διαστημικής Μεταφοράς, γνωστό ως Διαστημικό Λεωφορείο Space Shuttle) υπερυπερηχητικό όχημα. Τα αποτελέσματα της προτεινόμενης μεθόδου συγκρίνονται, αξιολογούνται και πιστοποιούνται με παρόμοια από διεθνώς αναγνωρισμένα υπολογιστικά συστήματα και πειραματικά δεδομένα μετρήσεων σε υπερυπερηχητικές αεροσήραγγες οχημάτων υψηλών ταχυτήτων. Μέσα στις επόμενες δεκαετίες, η εξέλιξη των υπερυπερηχητικών συστημάτων μεταφοράς θα είναι ραγδαία και ιδιαίτερα ελκυστική από τεχνολογικής πλευράς με ελπιδοφόρα μηνύματα για "γήϊνες" μετακινήσεις με υψηλές ταχύτητες σε σύντομους χρόνους, αλλά και πιθανές μετοικήσεις ανθρώπων-αστροναυτών σε άλλους πλανήτες του ηλιακού συστήματος. Έτσι η παρούσα διδακτορική διατριβή, αξιοποιώντας τις προόδους που έχουν συντελεστεί τα τελευταία χρόνια, εκτιμάται ότι συμβάλλει ουσιαστικά στο ακόμα υπό εξέλιξη – για τα ελληνικά δεδομένα πρωτόγνωρο – επιστημονικό πεδίο της Υπερυπερηχητικής Αεροθερμοδυναμικής, και στην ανάπτυξη μεθοδολογίας για τον αεροθερμοδυναμικό σχεδιασμό υπερυπερηχητικών οχημάτων. Το προτεινόμενο σύνθετο μοντέλο κινηματικής και θερμοδυναμικής υπολογιστικής ανάλυσης και η μεθοδολογία που ακολουθείται έχει ισχύ και εφαρμογή στην πλειονότητα των περιπτώσεων των υπερυπερηχητικών πτήσεων στη γήινη και, με κατάλληλες προσαρμογές, σε οποιαδήποτε πλανητική ατμόσφαιρα. Βεβαίως, η προτεινόμενη μεθοδολογία έχει περιθώρια εξέλιξης και ανάπτυξης, προκειμένου να βελτιωθεί η απόδοσή της με την εισαγωγή και νέων στοιχείων, αλλά και να διευρυνθεί το φάσμα των δυνατών εφαρμογών της. / The object of the present doctoral thesis constitutes the theoretical and calculating simulation of high-speed atmospheric flightmotion, with particular emphasis in hypersonic speeds and aerodynamic and thermodynamic analysis. The methodology that is followed has as main aim the prediction of return flight trace of hypersonic vehicles in Earth’s ground, assuming the speeds of flight to reach the 15000 m/s (Mach number ~ 40) and stagnation temperatures the 15000 K. In this work are followed methodological approach and systematic bibliographic research in the scientific field of aerodynamic overheating that suffers Hypersonic Vehicles(H/V) at their flight in Earth’s atmosphere. The analysis of hypersonic flow, thermal loads and new methods for thermal protection systems of H/V are developed the last decades and are evolved continuously, but until today (2008) resulting failures mainly at their reentry in atmosphere from the interval. For the simulation of the atmospheric hypersonic flight motion become study and analysis of thermodynamic properties and parameters, that characterize the flow fields in high speeds and high temperatures. The determination of altered flow conditions, thermodynamic properties and transport magnitudes will contribute in the theories of modelling and high thermal strain calculations, in the constitution of overheating air and it’s high temperature determination in the external surface of heat insulation wall of high-speed vehicles at their flight in Earth’s atmosphere. In this work are developed calculating algorithms of simulation for the successful prediction of atmospheric flight motion of hypersonic transport vehicles. The basic aerodynamic characteristics of lifting vehicles in high speeds are calculated with suitable adaptation of classic Newtonian Flow in high speeds and are confirmed with experimental data of aerodynamic measurements in high speed wind-tunnels of NASA. Moreover in this thesis are incorporated models for the simulation of properties of overheating air behaviour as real gas in thermodynamic balance in extreme flow conditions (high speeds, high pressures and low densities). Particularly becomes analysis of the altered constitution and properties of air (thermodynamic magnitudes and transport magnitudes) in very high temperatures on hypersonic flow field. This study will lead to the export of a new model for the transported thermal stagnation region calculations of hypersonic vehicles, which is confirmed so much with internationally recognized thermal models (as that of Fay-Riddell) and with calculating simulations of other researchers. The methodology that is followed will constitute the base for the export of a new theoretical model for high temperature air calculations in thermal protection systems surface of hypersonic vehicles. The effectiveness of proposed methodology, with the combination kinematic and aerothermodynamics analysis in united calculating algorithm, is proved and confirmed via the application of above simulation models on ballistic (Apollo Command Module) and lifting (Space Transportation System, known as Space Shuttle) hypersonic vehicles. The results of proposed method are compared, evaluated and certified with similarly by internationally recognized calculating systems and experimental data of measurements in hypersonic wind-tunnels for high speed vehicles. In the next decades, the development of hypersonic transportation systems will be rapid and particularly attractive from technological side with hopeful messages for “earthy” transfers with high speeds in short times, but also likely persons-astronaut missions in other planets of solar system. Thus the present doctoral thesis, developing the progress that has taken place in the past few years, is appreciated that it contributes substantially in still under development - for the Greek data unusual - scientific field of Hypersonic Aerothermodynamics, and in the growth of methodology for the aerothermodynamic construction of hypersonic vehicles. The proposed complex model of kinematic and thermodynamic calculating analysis and it’s methodology that is followed have influence and application in the majority of cases of hypersonic flights in earthy and, with properly adaptations, in anyone planetary atmosphere. Of course, the proposed methodology has margins of development and growth in order to improve her efficiency with the import of also new elements and also extended the spectrum of her possible applications.

One of the characteristics of the high speed ows over blunt bodies is the detached shock formed in front of the body. The distance of the shock from the stagnation point measured along the stagnation streamline is termed as the shock stand o distance or the shock detachment distance. It is one of the most basic parameters in such ows. The need to know the shock stand o distance arises due to the high temperatures faced in these cases. The biggest challenge faced in high enthalpy ows is the high amounts of heat transfer to the body. The position of the shock is relevant in knowing the temperatures that the body being subjected to such ows will have to face and thus building an efficient system to reduce the heat transfer. Despite being a basic parameter, there is no theoretical means to determine the shock stand o distance which is accepted universally. Deduction of this quantity depends more or less on experimental or computational means until a successful theoretical model for its predictions is developed. The experimental data available in open literature for spherical bodies in high speed ows mostly lies beyond the 2 km/s regime. Experiments were conducted to determine the shock stand o distance in the velocity range of 1-2 km/s. Three different hemispherical bodies of radii 25, 40 and 50 mm were taken as test models. Since the shock stand o distance is known to depend on the density ratio across the shock and hence gamma (ratio of specific heats), two different test gases, air and carbon dioxide were used for the experiments here. Five different test cases were studied with air as the test gas; Mach 5.56 with Reynolds number of 5.71 million/m and enthalpy of 1.08 MJ/kg, Mach 5.39 with Reynolds number of 3.04 million/m and enthalpy of 1.42 MJ/kg Mach 8.42 with Reynolds number of 1.72 million/m and enthalpy of 1.21 MJ/kg, Mach 11.8 with Reynolds number of 1.09 million/m and enthalpy of 2.03 MJ/kg and Mach 11.25 with Reynolds number of 0.90 million/m and enthalpy of 2.88 MJ/kg. For the experiments conducted with carbon dioxide as test gas, typical freestream conditions were: Mach 6.66 with Reynolds number of 1.46 million/m and enthalpy of 1.23 MJ/kg. The shock stand o distance was determined from the images that were obtained through schlieren photography, the ow visualization technique employed here. The results obtained were found to follow the same trend as the existing experimental data in the higher velocity range. The experimental data obtained was compared with two different theoretical models given by Lobb and Olivier and was found to match. Simulations were carried out in HiFUN, an in-house CFD package for Euler and laminar own conditions for Mach 8 own over 50 mm body with air as the test gas. The computational data was found to match well with the experimental and theoretical data

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.

Hypersonic flow is both fascinating and intriguing mainly because of presence of strong entropy and viscous interactions in the flow field. Notwithstanding the tremendous advancements in numerical modeling in the last decade separated hypersonic flow still remains an area where considerable differences are observed between experiments and numerical results. Lack of reliable data base of surface static pressures with good spatial resolution in hypersonic separated flow field is one of the main motivations for the present study. The experiments in hypersonic shock tunnels has an advantage compared to wind tunnels for simulating the total energy content of the flow in addition to the Mach and Reynolds numbers. However the useful test time in shock tunnels is of the order of few milliseconds. Hence in shock tunnel experiments it is essential to have pressure measurement devices which has special features such as small in size, faster response time and the sensors in array form with improved spatial resolutions. Micro Electro Mechanical Systems (MEMS) is an emerging technology, which holds lot of promise in these types of applications. In view of the above requirement, MEMS based pressure sensor array was developed to measure the static pressure distribution. The study is comprised of two parts: one is on the development of MEMS based pressure sensor array, which can be used for hypersonic application and other is on experimental static pressure measurement using MEMS based sensors in separated hypersonic flow over a backward facing step model. Initially a static pressure sensor array with 25 sensors was developed. The static calibration of sensor array was carried out to characterize the sensor array for various characteristic parameters. The preliminary experimental study with cluster of 25 MEMS sensor array mounted on the flat plate did not provide reliable and repeatable results, but gave valuable inputs on the typical problems of using MEMS sensors in short duration hypersonic ground test facilities like shock tunnels. Incidentally, to the best of our knowledge this is first report on use of MEMS based pressure sensors in hypersonic shock tunnel. Later cluster of 5 sensor array was developed with improved electronic packaging and surface finish. The experiments were conducted with flat plate by mounting 5 sensor array shows good agreement in static pressure measurement compared with standard sensors. In the second part of the study a backward facing step model, which simulates the typical gasdynamic flow features associated with hypersonic flow separation is designed. Backward facing step model with step height of 3 mm was mounted with sensor array along the length of model. Just after the step, static pressure measurements were carried out with MEMS sensors. It is important to note that, in the space available in backward facing step model we could mount only one conventional Kulite pressure transducer. The experiments were conducted at Mach number of 6.3 and at stagnation enthalpy of 1.5 MJ/kg in hypersonic shock tunnel (HST-5) at IISc. Based on the static pressure measurement on backward facing step, the location of separation and reattachment points were clearly identified. The static pressure values show that reattachment of flow takes place at about 7 step heights. Numerical simulations were carried out using commercial CFD code, FLUENT for flat plate and backward facing step models to compliment the experiments. The experimental tests results match well with the illustrative numerical simulations results.

Lattice Boltzmann Method has been quite successful for incompressible flows. Its extension for compressible (especially supersonic and hypersonic) flows has attracted lot of attention in recent time. There have been some successful attempts but nearly all of them have either resulted in complex or expensive equilibrium function distributions or in extra energy levels. Thus, an efficient Lattice Boltzmann Method for compressible fluid flows is still a research idea worth pursuing for. In this thesis, a new Lattice Boltzmann Method has been developed for compressible flows, by using the concept of a relaxation system, which is traditionally used as semilinear alternative for non-linear hypebolic systems in CFD. In the relaxation system originally introduced by Jin and Xin (1995), the non-linear flux in a hyperbolic conservation law is replaced by a new variable, together with a relaxation equation for this new variable augmented by a relaxation term in which it relaxes to the original nonlinear flux, in the limit of a vanishing relaxation parameter. The advantage is that instead of one non-linear hyperbolic equation, two linear hyperbolic equations need to be solved, together with a non-linear relaxation term. Based on the interpretation of Natalini (1998) of a relaxation system as a discrete velocity Boltzmann equation, with a new isotropic relaxation system as the basic building block, a Lattice Boltzmann Method is introduced for solving the equations of inviscid compressible flows. Since the associated equilibrium distribution functions of the relaxation system are not based on a low Mach number expansion, this method is not restricted to the incompressible limit. Free slip boundary condition is introduced with this new relaxation system based Lattice Boltzmann method framework. The same scheme is then extended for curved boundaries using the ghost cell method. This new Lattice Boltzmann Relaxation Scheme is successfully tested on various bench-mark test cases for solving the equations of compressible flows such as shock tube problem in 1-D and in 2-D the test cases involving supersonic flow over a forward-facing step, supersonic oblique shock reflection from a flat plate, supersonic and hypersonic flows past half-cylinder.

The numerical modelling of the VKI-Longshot facility remains a challeng-ing task as it requires multi-physical numerical methods in order to simulate all the components. In the current dissertation, numerical tools were developed in order to study each component of the facility separately and a deep investigations of each stage of the shot were performed. This helped to better understand the different processes involved in the flow development inside this hypersonic wind tunnel. However the numerical computation of different regions of the facility treated as independent from each others remains an approximation at best.The accuracy of the rebuilding code for determining the free stream conditions and the total enthalpy in the VKI-Longshot facility was investigated by using a series of unsteady numerical computations of axisymmetric hypersonic flow over a heat flux probe. Good agreement was obtained between the numerical results and the measured data for both the stagnation pressure and the heat flux dur- ing the useful test time.The driver-driven part of the Longshot facility was modelled using the quasi one-dimensional Lagrangian solver L1d2. The three main conditions used for the experiments —low, medium and high Reynolds number —were considered.The chambrage effect due to the junction between the driver and the driven tubes in the VKI-Longshot facility was investigated. The computation showed great ben- efit of the chambrage in increasing the speed of the piston and thus the final compression ratio of the test gas.Two dimensional simulations of the flow in the driver and the driven tube were performed using Arbitrary Lagrangian Eulerian (ALE) solver in COOLFLuiD. A parallel multi-domain strategy was developed in order to integrate the moving piston within the computational domain.The computed pressure in the reservoir is compared to the one provided by the experiment and good agreement was obtained for both con- editions.Finally, an attempt was made to compute the starting process of the flow in the contoured nozzle. The transient computation of the flow showed how the primary shock initiates the flow in the nozzle before reaching the exit plan at time of 1.5 [ms] after the diaphragm rupture. The complex interactions of the reflected shocks in the throat raise the temperature above 9500 [K] which was not expected. Chemical dissociation of Nitrogen was not taken into account during this transient investigation which may play a key role considering the range of temperature reached near the throat. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

The interaction of a shock wave and boundary layer often occurs in high speed flows. For sufficiently strong shock strengths the boundary layer separates, generating shock patterns in the contiguous inviscid flow (termed strong interactions); which may also affect the performances of the systems where they occur, demanding control of the interaction to enhance the performances. The case of impinging shock wave boundary layer interaction is of fundamental importance and can throw light on the physics of the interaction in general. Although various aspects of the interaction are studied at supersonic speeds, the complexities involved in the interaction at hypersonic speeds are not well understood. Of importance is the high total enthalpy associated with hypersonic flows the simulation of which requires shock tunnels. The present experimental study focuses on the interaction between strong impinging shock and boundary layer in hypersonic flows of moderate to high total enthalpies. Experiments are performed in hypersonic shock tunnels HST-2 and FPST (free piston driven shock tunnel), at nominal Mach numbers 6 and 8, with total enthalpy ranging from 1.3 MJ/kg to 6 MJ/kg, and freestream Reynolds number ranging from 0.3 million/m to 4 million/m. The strong impinging shock is generated by a wedge of angle 30.960 to the freestream. The shock is made to impinge on a flat plate (made of Hylem which is adiabatic, except for one case with plate made of aluminium which allows heat transfer). The position of (inviscid) shock impingement may be varied (from 55 mm from the leading edge to 100 mm from the leading edge) by moving the plate back and forth on the fixture which holds the wedge and the plate. Expectedly the strong shock generates a large separation bubble of length comparable to the distance of the location of shock impingement from the leading edge of the plate. Such large separation bubbles are typical of supersonic/hypersonic intakes at off-design operation. The evolution of the flow field- including the evolution of impinging shock and subsequent evolution of the large separation bubble- within the short test duration of the shock tunnels is one of the main concerns addressed in the study. Time resolved schlieren flow visualizations using high speed camera, surface pressure measurements using PCB, kulite and MEMS sensors, surface convective heat transfer measurements using platinum thin film sensors are the flow diagnostics used. From the time resolved visualizations and surface pressure measurements with the fast response sensors, the flow field, even with a separation bubble as large as 75 mm (at Mach 5.96, with shock impingement at 95 mm from the leading edge) was found to be established within the short shock tunnel test time. The effects of various parameters- freestream Mach number, distance of the location of shock impingement, freestream total enthalpy and wall heat transfer- on the interaction are investigated. With increase in Mach number from 5.96 to 8.67, for nearly the same shock impingement locations (95 mm and 100 mm from the leading edge respectively), the separation length decreased from 75 mm to 60 mm despite the fact that the shocks are doubly stronger at the higher Mach number. Inflectional trend in separation length was observed with enthalpy at nominal Mach number 8- separation length increased from 60 mm at 1.6 MJ/kg to 70 mm at 2.4 MJ/kg, and decreased drastically to ~40 mm at 6 MJ/kg (when dissociations are expected). The separation length Lsep for all the experiments, except the experiments at 6 MJ/kg, were found to be large, i.e. comparable with the distance xi of location of shock impingement from the leading edge of the flat plate. The scaled separation length (with Hylem wall) was found to obey the inviscid similarity law proposed from the present study for large separation bubbles with strong impinging shocks, where M∞ is the freestream Mach number, p∞ is the freestream pressure and pr is the measured reattachment pressure; this holds for freestream total enthalpy ranging from 1.3 MJ/kg to 2.4 MJ/kg and Reynolds number (based on location of shock impingement) ranging from 1x105 to 4x105. While the increase in separation length from 1.6 MJ/kg to 2.4 MJ/kg could thus be attributed to the small difference in Mach number between the cases (due to inverse variation with cube of Mach number), the decrease in separation length and the non-confirmation to the proposed similarity law for the 6 MJ/kg case is attributed to the real gas effects. At Mach 6 the flow was observed to separate close to the leading edge, even when the (inviscid) shock impingement was at 95 mm from the leading edge. This prompted the proposal of an approximate inviscid model of the interaction for the Mach 6 case with separation at leading edge, and reattachment at the location of (inviscid) shock impingement; Accordingly, the closer the location of impingement, the more the angle that the separated shear layer makes with the plate and hence more the pressure inside the separation bubble. A small reduction in separation length was also observed with aluminium wall when compared with Hylem wall, emphasizing the importance of wall heat conductivity (especially when concerning separated flows) even within the short test durations of shock tunnels. The free interaction theory over adiabatic wall was found to predict the pressure at the location of separation, but under-predict the plateau pressure (at nominal Mach number 8). Numerical simulations (steady, planar) were also carried out using commercial CFD solver FLUENT to complement the experiments. Simulations using one equation turbulence model (Spalart-Allmaras model) were closer to the experimental results than the laminar simulations, suggesting that the flow field may be transitional or turbulent after separation. Significant reduction of the separation bubble length was demonstrated with the control of the interaction using boundary layer bleed within the short test time of the shock tunnel; with tangential blowing at the separation location20% reduction in separation length was observed, while with suction at separation location the reduction was 13.33 %.

With the age of hypersonic flight imminent just beyond the horizon, researchers are working hard at designing work-arounds for all the major problems as well as the minor quirks associated with it. One such issue, seemingly innocuous but one that could be potentially deadly, is the problem of interference heating due to surface protuberances. Although an ideal design of the external surfaces of a high-speed aircraft dictates complete smoothness to reduce drag, this is not always possible in reality. Control surfaces, sheet joints, cable protection pads etc. generate surface discontinuities of varying geometries, in the form of both protrusions as well as cavities. These discontinuities are most often small in dimension, comparable to the local boundary layer thickness at that location. Such protuberances always experience high rates of heat transfer, and therefore should be appropriately shielded. However, thermal shielding of the protrusions alone is not a full solution to the problem at hand. The interference caused to the boundary layer by the flow causes the generation of local hot spots in the vicinity of the protuberances, which should be properly mapped and adequately addressed. The work presented in this thesis aims at locating and measuring the heat flux values at these hot spots near the protrusions, and possibly formulating empirical correlations to predict the hot spot heat flux for a given set of flow conditions and protrusion geometry. Experimental investigations were conducted on a flat plate model and a cone model, with interchangeable sharp and blunt nose tips, with attached 3D protuberances. Platinum thin-film sensors were placed around the protrusion so that the heat fluxes could be measured in its vicinity and the hottest spot located. These experiments were carried out at five different hypersonic free stream flow conditions generated using two shock tunnels, one of the conventional type, and the other of the free-piston driven type. The geometry of the protrusions, i.e., the height and the deflection angle, was also parametrically varied to study its effect on the hot spot heat flux. The results thus obtained for the flat plate case were compared to existing correlations in the open literature from a similar previous study at a much higher Reynolds number range. Since a mismatch was observed between the results of the current experiments and the existing correlations, a new empirical correlation has been developed to predict the hot spot heat flux, that is valid within the range of flow conditions studied here. A similar attempt was made for the case of the cone model, for which no previous correlations exist in the open literature. However, a global correlation covering the entire range of flow conditions used here could not be formed. A correlation that is valid for just one out of the five flow conditions used here is presented for the cones with sharp and blunt nose tips separately. Schlieren flow visualization was carried out to obtain a better understanding of the shock structures near the protuberances on both models. For most cases, where the protrusion height and deflection angle were large enough to cause flow separation immediately upstream of the protuberance, a separation shock was manifested which deflected some part of the boundary layer above the protuberance, while the rest of the fluid in the boundary layer entered a recirculating region in the separated zone before escaping to the side. Some preliminary computational analysis was conducted which confirmed this qualitatively. However, the quantitative match of surface heat flux between the simulations and experiments were not encouraging. Schlieren visualization revealed that for the flat plate case, the foot of the separation shock was located at a distance of 10.5 to 12 times the protrusion height ahead of it, whereas in the case of the sharp cone, it was at a distance of 9 to 10.5 times the protrusion height. The unsteady nature of the separation shock was also captured and addressed. Some preliminary experiments on boundary layer tripping were also conducted, the results of which have been presented here. From this analysis, it has become evident that a single global correlation cannot be formed which could be used for a wide range of flow conditions to predict the hot spot heat flux in interference interactions. The entire range of conditions that may be encountered during hypersonic flight has to be broken down into sections, and the interference heating pattern should be studied in each of these sections individually. By doing so, a series of different correlations can be formed at the varying flow conditions which will then be available for high-speed aircraft designers.

Shock-shock interactions are among the most basic gas-dynamic problem, and are almost unavoidable in any high speed light, where shock waves generating from different sources crosses each other paths. These interactions when present very close to the solid surface lead to very high pressure and thermal loads on the surface. The related practical problem is that experienced at the cowl lip of a scramjet engine, where the interfering shock waves leads to high heat transfer rates which may also lead to the damage of the material. The classification by Edney (1968) on the shock-shock interaction patterns based on the visualization has since then served the basis for such studies. Though the problem of high heating on the surface in the vicinity of the shock-shock interactions has been studied at length at supersonic Mach numbers, the study on the topic at the hypersonic Mach numbers is little sparse. Even in the studies at hypersonic Mach numbers, the high speeds are not simulated, which is the measure of the kinetic energy of the ow. Very few experimental studies have addressed this problem by simulating the energy content of the ow. Also, some of the numerical studies on the shock-shock interactions suggest the presence of unsteadiness in the shock-shock interaction patterns as observed by Edney (1968), though this observation is not made very clearly in the experimental studies undertaken so far. In the present study, experiments are carried out in a conventional shock tunnel at Mach number of 5.62 (total enthalpy of 1.07 MJ/kg; freestream velocity of 1361 m/s), with the objective of mapping the surface pressure distribution and surface convective heat transfer rate distribution on the hemispherical body in the presence of the shock-shock interactions. A shock generator which is basically a wedge of angle = 25 , is placed at some dis-dance in front of the hemispherical body such that the planar oblique shock wave from the shock generator hits the bow shock wave in front of the hemi-spherical body. The relative distance between the wedge tip and the nose of the hemispherical body is allowed to change in di erent experiments to capture the whole realm of shock-shock interaction by making the planar oblique shock wave interact with the bow shock wave at different locations along its trajectory. The study results in a bulk of data for the surface pressure and heat transfer rates which were obtained by placing 5 kulites pressure transducers, 1 PCB pressure transducer and 21 platinum thin lm gauges along the surface of the hemispherical body in a plane normal to the freestream velocity direction. Along with the measurement of the surface pressure and the surface heat transfer rates, the schlieren visualization is carried out to capture the shock waves, expansion fans, slip lines, present in a certain shock-shock interaction pattern and the measured values were correlated with the captured schlieren images to evaluate the ow build up and steady and useful test time thereby helping in understanding the ow physics in the presence of the shock-shock interactions. From the present study it has been observed that in the presence of Edney Type-I and Edney Type-II interaction, the heat transfer rates on the hemi-spherical body are symmetrical about the centerline of the body, with the peak heating at the centerline which drops towards the shoulder. For Edney Type-III, Edney Type-IV, Edney Type-V and Edney Type-VI interaction pattern, the distribution in not symmetrical and shifts in peak heat transfer rates being on the side of the hemispherical from which planar oblique shock wave is incident. Also, it is observed that for the interactions which appear within the sonic circle, Edney Type-III and Edney Type-IV, the heat transfer rates observe an unsteadiness, such that the gauges located close to the interaction region experiencing varying heat transfer rates during the useful test time of the shock tunnel. Few experiments were conducted at Mach 8.36 (total enthalpy of 1.29 MJ/kg; freestream velocity of 1555.25 m/s) and Mach 10.14 (total enthalpy of 2.67 MJ/kg; freestream velocity of 2258.51 m/s) for the con gurations representing Edney Type-III interaction pattern to further evaluate the unsteady nature observed at Mach 5.62 ows. The unsteadiness was evident in both the cases. It is realized that the short test times in the shock tunnels pose a constraint in the study of unsteady flow fields, and the use of tailored mode operation of shock tunnel can alleviate this constraint. Also, limited number of experiments in the present study, which are carried out in a Free Piston Shock Tunnel, helps to understand the need to conduct such study in high enthalpy test conditions.