In the study of high speed flow (supersonic and hypersonic) related to aerodynamics of aerospace vehicles, concepts such as heat transfer, non equilibrium flow, ionization, dissociation and other high temperature effects are very important. These effects are very difficult to model theoretically or computationally. It is also impossible to measure the flow conditions which are encountered by the aerospace vehicles during the actual flight; therefore high speed transient flow test facility is an important component in the design process of aeronautics equipment.
The problem of providing a source of air at a sufficiently high temperature and pressure, similar to those encounters by aerospace vehicle to act as the working fluid of wind tunnel has been approached in a variety of ways. Because of the high capital cost and also high operating cost associated with continuous high speed wind tunnels, most of the high speed flow test facilities, especially those of high temperature, are short duration blow down type. For example Hotshot tunnels, Plasma arc tunnels and light gas gun . Furthermore, in high speed fluid flow problems, achieving a high Mach number is less important than achieving the correct flow enthalpy (temperature). This is because the high temperatures which developed around the stagnation region have a profound effect on the gas dynamics and aerodynamics . Continuous facilities cannot provide high enthalpy flows due to design problems associated with heat transfer especially in nozzle and reservoir regions. Transient high speed flow test facilities are very attractive due to their ability to provide high enthalpy flows .
Short duration high speed transient flow test facility is a device that can be used to generate gas flows or gas conditions that are difficult to achieve in other test devices. By its nature, these facilities produce these conditions for a very short duration. The maximum achievable stagnation temperature or enthalpy in a short duration facility is an inverse function of flow duration. For example, expansion tubes can provide stagnation enthalpy in the order of 100 MJ/kg but only for tens of microseconds . The useful flow duration in shock tunnels is longer (in the order of milliseconds) but the maximum stagnation enthalpy decreases to around 20 MJ/kg. Hotshot facilities that produce useful flows for around 100 ms have a maximum stagnation enthalpy of around 4 MJ/kg. Blow down tunnels are near continuous facilities that can provide long flow durations (measured in tens of minutes) but the maximum stagnation enthalpy is generally somewhat less than 4 MJ/kg .
In this article, a new short duration high speed transient flow test facility which has been developed in Universiti Tenaga Nasional – Malaysia is described. The test facility can be operated in four different modes, shock tube, free-piston compressor, shock tunnel and gun tunnel. All utilize the same basic facility and offer alternative modes of operation for the same apparatus. They differ in such away that, in shock tube and shock tunnel, the driver gas is employed to operate a shock tube containing the test gas, whilst in the free-piston compressor the free piston is used to directly compress the test gas in the driven section. The free piston in the gun tunnel compresses the test gas before ejecting it through a high speed flow nozzle.
The developed facility is the first of its kind in Malaysia. In other part of the world, there are currently existing facilities which can reach these high temperatures, but they are limited and rare. Free piston hypersonic shock tunnels have been constructed at California Institute of Technology (Caltech) in the United States and the German Aerospace Center (DLR Deutsches Zentrum fur Luft-und Raumfahrt) in Germany, but most of these high temperature ground facilities were very costly to build and expensive to operate and maintain .
The main reason for development of the test facility is to enable aerodynamics study related to high speed flow. Apart from that, such high speed transient flow test facility is fundamentally important for the development of advanced instrumentation for high speed flow, for example, fiber optic pressure sensors and fast response thermocouples and investigation of flow and heat transfer in high speed low pressure turbine.
Although similar facilities have been built in a few parts of the world , there have been no open literatures available that describe the design procedures and important parameters that need to be considered when building such facilities. Hence, in the present work, a systematic approach has been developed in designing such facilities.
Historically, a large number of methods have been used to improve the performance of shock tubes and shock tunnels . One method is to fill the driver section with a light gas such as helium. Another is to increase the temperature of the driver gas by use of a heater. In both of these methods, the improved performance is achieved by a higher speed of sound than if cold air is used. In the first case, the speed of sound in helium is higher than air because of its lower molecular weight. In the latter, the speed of sound is increased by raising the gas temperature. The higher sound speed results in a lower driver-to-driven tube pressure ratio
Although both of the above mentioned approaches are well established, there are a number of disadvantages associated with them. Both require high pressure pumps and extensive plumbing. Helium is expensive and large amounts must be used to pressurize the driver tube. Helium’s higher specific heat ratio also reduces the advantages somewhat. If the driver tube is to be heated, the tube must be designed to withstand high pressure under high temperatures, adding to the construction costs. The heater also adds to the cost and it must be designed to heat uniformly because hot spots in the driver can cause unpredictable and dangerous driver tube failure. Thermal fatigue may, in fact, limit the useful life of the tube.
For quite sometime, the free-piston driving technique, first proposed by Stalker , has been used to achieve some of the highest enthalpies, culminating in the X3 at the University of Queensland, Brisbane, Australia , the HEG in G¨ottingen, Germany  and the largest known facility, the HIEST in Kakuda, Japan . The free-piston technique involves compressing the driver gas by a heavy piston accelerated to nearly sonic speed. The piston compresses the gas ahead of it to achieve high values of temperature and pressure. The operation of a free-piston shock tunnel, however, appears to be complicated, as a massive piston has to be accelerated rapidly and then must be stopped in a controlled manner.
Many methods of producing strong shock waves have been and are still being analyzed and tried experimentally.
2. Objectives and Methodology
The main objective of this study is to design and develop a short duration high speed transient flow test facility and investigate the parameters that influence the performance of the facility. The performance of this test facility is described by the maximum Mach number that can be generated and the quality of flows that it can produce. The Mach number depends on the maximum temperature and pressure achievable. The facility is aimed to produce a targeted Mach number of 4.0. Different values of Mach number can be obtained by changing the diaphragm pressure ratio, different driver/driven gas combinations and by rising the driver gas temperature or in another word, increasing driver/driven gas temperature ratio
A theoretical model has been developed to determine the shock strength (
In order to investigate the detail flow process inside the shock tube, which will influence its performance, a new two dimensional time accurate Navier-Stokes solver for shock tube applications was developed. The solver uses second order accurate cell-vertex finite volume spatial discretization and fourth order accurate Runge-Kutta temporal integration with Air-Air as working fluids. The solver is validated against experimental measurements in the developed high speed flow test facility. Further investigations were made on the flow process inside the shock tube by using the solver. The shock wave motion, reflection and interaction with the boundary layers were investigated and their influence on the performance of the shock tube was determined.
An in-house made surface junction fast response thermocouple was used to measure the surface temperature change profile during the facility operation. In order to evaluate the heat flux from the surface temperature change history, a MATLAB numerical transient heat transfer model was developed. Extensive experimental measurements were performed at different pressure ratios in order to investigate the performance of the facility. Pressure transient inside the tube were captured at two different locations using fast response transducers.
The last part of the work concerns with further experimental investigation of the performance of the facility using different gas combinations (dissimilar gases). The gas combinations used are He-Air and He-CO2. The experimental results were compared with the theoretical predictions. The discrepancies were discussed based on the factors determined in the previous part.
3. Overview of UNITEN’s Test Facility
There are three important limitations that govern the selection of the shock tube dimensions; cost, available space, and manufacturing process. Taking into consideration all these parameters, the decision was made to choose the dimensions of the test facility as described below. The detail components of the facility are described briefly and shown in Figure 1. More details about the facility are available in .
Driver section:- A high-pressure section (driver) which will contain the high pressure driver gas, which can be either Air, Helium, Hydrogen or other light gases.
Discharge valve:- To discharge the driver section after each run.
Pressure gauge:- To read the pressure inside the driver section, which is also provided with a static pressure transducer to record the exact value of the driver pressure
Vacuum pump:- When the driver gas is not Air (e.g. Helium or Hydrogen) then the driver section should be evacuated and refilled with the required driver gas.
The primary diaphragm:- This is a thin aluminum membrane to isolate the low-pressure test gas from the high-pressure driver gas until the compression process is initiated.
Piston compression section:- A piston is placed in the (driven tube) adjacent to the primary diaphragm so that when the diaphragm ruptures, the piston is propelled through the driven tube, compressing the gas ahead of it. This piston is used in the free-piston compressor and gun tunnel tests.
Discharge valve:- To discharge the driven section after each run.
Vacuum gauge:- To set the pressure inside the driven section to values less than atmospheric value (vacuum).
Driven section:- A shock tube section (smooth bore), to be filled with the required test gas (Air, nitrogen or carbon dioxide).
Driven section extension:- The last half meter of the driven section on which the pressure transducers and thermocouples are mounted.
The secondary diaphragm:- A light plastic diaphragm to separate the low pressure test gas inside the driven section from the test section and dump tank which are initially at a vacuum prior to the run.
Test section:- This section will expand the high temperature test gas through a nozzle to the correct high enthalpy conditions needed to simulate hypersonic flow. A range of Mach numbers is available by changing the diameter of the throat insert.
Vacuum vessel (dump tank):- To be evacuated to about 0.1 mm Hg pressure before running. Prior to a run, the driven section, test section and dump tank are to be evacuated to a low-pressure value.
4. Numerical Formulation of Flow Process
In this section, a two dimensional time accurate Navier-Stokes solver for shock tube applications is described briefly. The solver has been programmed based on the dimensions and configuration of the test facility. The developed solver uses second order accurate cell-vertex finite volume spatial discretization and fourth order accurate Runge-Kutta temporal integration and it is designed to simulate the flow process for similar driver/driven gases (e.g. Air-Air as working fluids). The solver is validated against analytical solution and experimental measurements in the high speed flow test facility. In the next section, further investigations were made on the flow process inside the shock tube by using the solver. The shock wave motion, reflection and interaction with the boundary layers were investigated and their influence on the performance of the shock tube was determined.
It is intended that this CFD solver can be a useful tool in the design process of the test facility. Experimental data which could be achieved from the facility can be verified using this developed numerical tool.
The mathematical model and the fluid flow governing equations are described in ref . This includes the Reynolds-Averaged Navier-Stokes Equations (RANS) and The Mixing Length Turbulence Model. The latter is added to the solver in order to include the viscous effects to the solution. Based on the facility dimensions, mesh generation procedure was performed to transform the flow domain into a finite number of grid points. Finally, using two verification approaches, the code was validated in terms of the ability to capture shocks, rarefaction waves and contact discontinuity and to produce the correct pressure, temperature, density and speed profiles, using two verification approaches. The first one was the validation against a standard analytical solution of the shock tube problem. The second was to compare the code solution with selected experimental measurements for a certain value of diaphragm pressure ratio. Further details about the solver can be found in .
5. Selected Case Studies
In this section, in order to obtain greater understanding of the processes involved, CFD simulations were used for selected cases and the results were analyzed in details. In order to investigate the effect of various parameters on the performance of the facility, various experimental measurements were made. The parameters that were investigated were pressure ratio, gas combination and effect of piston. This section provides detail discussions of the results obtained.
5.1 Inviscid Transient Flow in Shock Tube
CFD solution for an inviscid simulation of a diaphragm pressure ratio
The pressure history for the above mentioned shot is depicted in Figure 3 from which one can follow the physics of the flow inside the shock tube. The first jump represents the shock wave, for which the pressure inside the driven section increases from 100 kPa to around 220 kPa. As the shock wave proceeds to the end of the tube it will reflect and move in the opposite direction, increasing the pressure to about 450 kPa. The shock wave will then interact with the contact surface which is following the shock wave and due to this interaction between the shock wave and the contact surface the pressure will be increased until it reaches its peak pressure value of 530 kPa.
The shock wave speed can be determined from the CFD data obtained in this simulation. As the distance between the two stations is known (0.342 m) and the time of the shock to travel from station 1 to station 2 can be obtained from the pressure history graph (see Figure 4), the shock wave speed in this case was calculated to be 518 m/s. Comparing this value to the theoretical value for this pressure ratio, which is 558 m/s, the percentage difference was found to be around 7%. The difference is probably due to the two-dimensional effects which were not modeled by the theoretical solution. From experimental measurements the shock speed for the same pressure ratio was 450 m/s, which indicate percentage difference of about 13% from CFD results.
Using the same procedure, the reflected shock wave speed can be determined. As the wave reflects from the tube end and moves in the opposite direction (left direction), due to impact with the end wall the wave will lose some of its kinetic energy due to the collision with the end wall and consequently its speed decreases to about 342 m/s.
The same trend can be noted when the temperature history is investigated, as shown in Figure 5. The first jump in the temperature profile represents the shock wave and the second jump is due to the reflected shock wave. The temperature is increased from the initial value 300 K to about 380 K due to shock wave effect. When the shock reflects from the tube end, the temperature rises to 475 K and after interaction between reflected shock wave and the contact surface, the flow temperature becomes about 490 K.
In order to have an overall view of what is happening inside the tube after diaphragm rupture, the
It is interesting to note that after interaction with the reflected shock wave, the contact surface remains at about the same position, indicating achievement of the tailored condition. The presence of the bush is also seen to have prevented the rarefaction wave and the shock wave from passing to the other section. The rarefaction wave and the shock wave are reflected when they reach the bush.
After diaphragm rupture, shock wave travels to the right through the driven section while the expansion wave travels to the left through the driver section. These two waves are captured after 0.005 sec and shown in Figure 9.
The two waves continue their journey towards the tube ends, at time 0.0078 s the shock wave hits the driven section end on the right hand side while the expansion wave reaches the left hand side of the facility as shown in Figure 10.
The reflected shock wave now moves to the left towards the contact surface while the reflected expansion wave moves to the right towards the bush section as shown in Figure 12.
The shock wave interacts with the contact surface and reflects again. This process continues until pressure balance along the whole facility as shown in Figure 13.
5.1.1 Two-Dimensional Effects
The presence of the bush has caused the flow in the facility to be two-dimensional and this requires two-dimensional CFD simulations. The contour plots of the velocity in
The velocity contours at
To investigate further, the velocity profiles at different times at
The formation of the recirculating region in this inviscid simulation is surprising especially considering that the tube is symmetrical . However, it has been previously reported by Xu Fu et al.  that high speed flow tends to become unstable when shock wave interacts with contact discontinuity.
In order to get deeper understanding to what is happening inside the shock tube after diaphragm rupture; velocity, pressure, temperature and density profiles below tube axis (25% of the tube diameter), along the tube axis (50% of the tube diameter) and above tube axis (75% of the tube diameter) were depicted at different times. Figure 20 shows the results at
Figure 21 shows the flow profiles after shock reflection. From the velocity profiles at this time, it can be seen that the flow is still symmetry about the tube longitudinal axis, however, it can be noted that asymmetry has already started. Pressure, temperature and density along the driven section (except the bush section) are all increased due to the effect of the compression waves.
As the reflected shock proceeds to the left, asymmetric flow becomes very clear as shown in Figure 23. The gas flows to the right in the upper side while it flows to the left in the bottom side of the tube, as can be seen from velocity profile shown in Figures 24.
After a series of frequent reflections the gas becomes more quiescent and moves towards a balance state as shown in Figure 25.
5.2 Viscous Transient Flow in Shock Tube
In order to investigate the effect of viscosity on the transient flow in shock tube and how it affects the performance of the facility, a viscous simulation has been conducted using the same boundary conditions as for the inviscid simulation presented in the previous section. The pressure history for the above mentioned shot is depicted in Figure 26. The Figure shows similar trend as for the inviscid flow. The first jump represents the pressure rise due to shock wave, for which the pressure inside the driven section increases from 100 kPa to around 220 kPa. The shock wave then reflects as it hits the end of the tube and moves in the opposite direction, subsequently rising the pressure to about 450 kPa.
The shock wave will then interact with the contact surface (which is following the shock wave) and due to their interaction, the pressure will be increased until it reaches its peak pressure value, which is equal to around 530 kPa in this case.
Using similar procedure as outlined previously, the shock wave speed in this case was calculated to be 456. Comparing to the inviscid value for the same pressure ratio obtained previously (518 m/s) the effect of viscosity becomes obvious. It can be deduced that viscosity decreases the shock wave speed to about 11% due to the boundary layer effects. From experimental measurements the shock speed for the same pressure ratio was 450 m/s, which indicated percentage difference of about 1.3% from the CFD result.
After it hits the tube end, shock wave will be reflected and moves to the left with a slower wave speed of about 311 m/s. Comparing with respect to the reflected shock wave speed for inviscid flow which is 342 m/s, it is apparent that viscosity resists the fluid motion causing slower speed of the shock wave by 9.1%.
Analyzing the temperature history for this simulation, it can be seen that the trend is quite similar to pressure history. The temperature results for this run have been displayed in Figure 28. The first jump represents the shock wave and the second jump is due to the reflected shock wave. The temperature is increased from the initial value of 300 K to about 500 K.
Figure 29 and 30 show
Figure 30 shows the so-called tailored interface, where no disturbance is reflected from the contact surface back towards the rear wall of the shock tube. The “tailored” contact surface configuration offers a number of advantages when applied to the operations of shock tubes, it increases the testing-time and it improves the homogeneity of the working gas parameters (i.e. it decreases possible contamination effects in the test section caused by the driver gas).
As shown in Figure 30, the maximum useful duration time that can be obtained when the prescribed pressure ratio
5.2.1 Shock Wave - Boundary Layer Interaction
Ideally, the reflection of a shock wave from the closed end of a shock tube provides, for laboratory study, a quantity of stationary gas at extremely high temperature. Because of the action of viscosity, however, the flow in the real case is not one-dimensional and a boundary layer grows in the fluid following the initial shock wave. In the flow following the initial shock wave, there is a boundary layer generated near the walls of the shock tube, across which the velocity of the flow decreases from that in the main stream to zero at the walls.
Figure 31 shows the velocity profiles at
In section 5.1.1, it has been shown that the inviscid flow tend to be very unstable in the region close to the bush after shock wave interacts with contact surface and the recirculating flow has developed. In order to investigate the effect of viscosity to this performance, Figure 32 to 34 shows the evolution of the velocity profiles for viscous flow case. It can be seen in Figure 34 that before the shock wave was reflected, after 0.0005 sec from diaphragm rupture, the boundary layer separation has occurred close to the tube wall and the separated region grows.
After the shock wave reflection and subsequent interaction with the contact surface, the separated flow region evolves into a full re-circulating region rotating in the anticlockwise direction. Then the flow returns back to the small separated flow region close to the tube walls.
5.3 The Effect of Diaphragm Pressure Ratio on the Facility Performance
It has been reported by many researchers that increase in diaphragm pressure ratio,
Figures 35 and 36 show that as the pressure ratio increases, the pressure rise and the shock speed increases. The recorded initial rises are 280 kPa, 320 kPa and 420 kPa for pressure ratio 8.8, 15 and 20 respectively. After the shock reflection, the pressure rise further to 580 kPa, 800 kPa and 1000 kPa for pressure ratio 8.8, 15 and 20 respectively.
The corresponding temperature histories are shown in Figure 37. It can be seen that the temperature rise to 500 K, 575 K and 625 after the shock reflection for pressure ratios 8, 15 and 20 respectively. This clearly indicates that the pressure ratio has the effect of rising the temperature and hence the enthalpy of the gas.
The measured shock strength (
Figures 39 and 40 show the measured, theoretical and CFD results for shock speed and shock Mach number respectively. The Mach number is calculated based on the speed of sound of the undisturbed gas. The measured and CFD data are generally lower than the theoretical data and the agreement between CFD and measured data is very good. For the highest allowable pressure ratio of 55, the maximum shock Mach number is 2.116 when Air was allocated in both the driver and driven sections.
Figure 41 shows the measured peak pressure variation with pressure ratio. The peak pressure is defined as the maximum pressure achieved after shock reflection and subsequent interaction with contact surface, which further increases the pressure. It can be seen that the peak pressure rise increase when the pressure ratio is increased but at a slower rate. The maximum peak pressure produced was 2200 kPa.
5.4 The Effect of Gas Combination on the Facility Performance
Theoretical analysis showed that in order to achieve higher shock strength and shock speed, a gas combination with higher ratio of γ4/γ1 should be used. This can be achieved by using lighter gas in the driver section and heavier gas in the driven section. In the present work, combinations of He, CO2 and Air were chosen. Other lighter gas such as Hydrogen was not used due to safety reason.
Figure 42 shows the pressure history at station (1) for pressure ratios of approximately 8, 15 and 20 respectively. The corresponding graphs at station (2) are plotted in Figures 43.
In general, it can be seen from Figures 42 and 43 that the He-Air gas combination provides the best results in terms of shock speed, shock strength and peak pressure. Comparing to Air-Air gas combination, He-CO2 gives better shock strength and peak pressure. The recorded initial rises at station (1) for He-Air are 400 kPa, 450 kPa and 650 kPa for 8, 15 and 20 pressure ratios respectively. After the shock reflection, the pressure rise further and peak pressure of 1500 kPa, 1800 kPa and 2400 kPa for pressure ratio 8, 15 and 20 respectively.
The experimental measurements were performed for pressure ratios from 8 to 55. Figures 44 and 45 show the measured shock speed and shock Mach number, respectively, for different gas combinations at various pressure ratios. Also plotted on the same graphs are the theoretical solutions.
Similar to the Air-Air case, in general the measured data shows lower values compared to theoretical solution. However the trends are very similar. It can be seen that when γ4/γ1 is reduced the shock speed and Mach number increases. The maximum shock Mach number achieved was 3.69 when He-CO2 was used at pressure ratio of 55. The peak pressure values achieved at various pressure ratios and gas combinations are plotted in Figure 46.
5.5 The Effect of Piston on the Facility Performance
Previous researchers have shown that the addition of very light piston immediately after the diaphragm will increase the peak pressure achieved. After diaphragm rupture, the piston velocity will rapidly approach the contact surface velocity of a conventional shock tube. Multiple reflections of the shock wave ahead of the piston between the end of the driven section and the piston will compress the gas non-isentropically and result in a higher peak pressure.
Two experiments have been performed at pressure ratio
It is apparent from Figures 48 and 49 that the shock speed is faster when the piston is not used, this is because the piston slows down the test gas ahead of it and consequently the shock will be slower. However, the peak pressure of that using piston is higher; the piston increases the peak pressure due to the compression process of the test gas. In spite of this, the He-Air as combination shows a tremendous improvement in terms of both shock strength and peak pressure. It is worth to mention here that the He-Air combination can produce a peak pressure of about 70 bar (with piston) as compared to that of Air-Air (33 bar), which has marked percentage of increase of more than 100%.
The shock speed and the shock Mach number can be obtained by comparing the pressure transient at station (1) and station (2). A sample result for Air-Air at pressure ratio 13 is shown in Figure 50. In this particular case the shock speed and shock Mach number is 441.3 m/s and 1.3 respectively.
Similarly, Figure 51 shows the pressure history at stations 1 and 2 for shock tube shot.
As explained previously, due to the piston effects in the free piston compressor test, the shock wave speed and shock Mach number for the shock tube shot are found to be higher, which are 570 m/s and 1.67 respectively.
5.6 Surface Temperature and Heat Flux
The gas temperature increase significantly during the test duration. For example, from the CFD results of Air-Air at a pressure ratio of 20 will produce gas temperature rise of 350 K. This high temperature rise, although for a very short duration will likely cause the tube wall temperature to rise and further cause heat loss that can affect the performance of the facility. In order to ascertain this effect, a surface junction thermocouple was used to measure the wall temperature during the test.
Figure 52 shows the surface temperature change measured experimentally at three different pressure ratio for shock tube shots with Air-Air gas combination. It clearly shows that the difference in surface temperature increases as the pressure ratio increases, which consequently enhance the undesirable heat loss to the surrounding. As explained previously, the first temperature jump is due to shock wave as it compresses the test gas and as the shock wave reflects from the end wall, further temperature rise is achieved.
The corresponding figure for the heat flux is shown in Figure 53 and as depicted, the heat flux increases after shock wave passes through the test gas, which may consequently, influence the facility performance.
Figure 54 presents surface temperature change for three gas combinations at a fixed pressure ratio of 8.4. The results show that maximum heat transfer occurs when He-Air was used.
The same trend can be observed for the free piston compressor shot as shown in Figure 55. It is worth mentioning here that the free piston compressor was able to produce a surface temperature change of about 4 oC (with piston) as compared to about 1.7 oC in the conventional shock tube indicating marked percentage increase of approximately 100%.
6. Discussion and Conclusions
The goal of designing a high speed transient flow test facility with a relatively inexpensive price for research use has been achieved. This design for a multi arrangement high speed fluid flow tunnel will provide high speed fluid flows around test objects for actual flow measurement. Although the control of flow quality in such a tunnel would not be so easy, the actual ability for researchers to get their knowledge in real high speed fluid flow would be an immeasurable benefit.
The experimental and theoretical results agree very well with slight difference in shock speed values which are expected as the theoretical solution disregards the viscous, two dimensional and heat transfer effects which influence the shock strength and shock speed. It can be noted from results that time of shock wave creation for both experimental and CFD data are very much matched. The incident shock wave travels all the way along the driven section until it reflects off of the end wall. Similarly, the expansion wave reflects at the end of the driver section. The reflected shock then interacts with either the contact surface or the reflected expansion wave. It can be said that the general pattern of the experimental and CFD results are quite similar. Experiments showed that the shock wave speed for He/Air driver/driven gases is higher than He-CO2 and Air/Air gas combinations for the same diaphragm pressure ratio.
Shock speed can be increased by raising the diaphragm pressure ratio, or more powerfully, by raising the speed of sound in the driver gas. To achieve high Mach numbers it is essential to raise the speed of sound ratio (
Results when employing different working fluids combination are also presented. The results showed a direct proportional relationship between Mach number and diaphragm pressure ratio and inverse proportion with speed of sound ratio.
The results show that for Air-Air driver/driven gases, the shock speed is decreased when piston is used. However, the peak pressure produced using piston is higher. The same trend can be observed for the He-Air combination. In spite of this, the latter shows a tremendous improvement in terms of both shock strength and peak pressure . It is worth to mention, that He-Air combination can produce a peak pressure of about 70 bar for diaphragm pressure ratio
As the shock wave propagates in the driven section the test gas temperature rises up and consequently heat transfer to the tube wall causes an increase in the wall surface temperature. The same trend can be observed for the free piston compressor shot. It is worth to mention here that the free piston compressor can produce a surface temperature change of about 4 oC (with piston) as compared to that of conventional shock tube (1.7 oC), which has marked a percentage of increase of approximately 100%. The stagnation temperature achieved in free piston compressor is limited by the strength of the light piston used to compress the test gas, and the useful running time is restricted by the large heat losses to the cold walls of the tube.
The maximum testing time is varying between the arrival of the primary shock wave at the tube end, and the arrival of the first reflected wave, which has passed up the low pressure tube, reflected from the high pressure tube and returned to the nozzle station. The useful test time achieved from results was 10 ms which is quite comparable to other existing facilities.
Although this design showed good results, it should be noted that in the future, upgrades could be made to the facility to provide the better flow quality and longer flow times that would be needed for high speed fluid flow research. Detail design calculations, numerical simulations and experimental works have been conducted. The results presented in this chapter show that two-dimensional modeling of the high speed fluid flow test facility is an effective way to obtain facility performance data. Although this thesis focused on the UNITEN’s facility, the CFD code is generic and may be applied to other facilities.