Open access peer-reviewed chapter

Investigation of Hypersonic Conic Flows Generated by Magnetoplasma Light-Gas Gun Equipped with Laval Nozzle

Written By

Pavel P. Khramtsov

Reviewed: July 14th, 2021 Published: December 10th, 2021

DOI: 10.5772/intechopen.99457

Chapter metrics overview

77 Chapter Downloads

View Full Metrics


This chapter introduces new approach of hypersonic flow generation and experimental study of hypersonic flows over cones with half- angles τ1 = 3◦ and τ2 = 12◦. Mach number of the of the incident flow was M1 = 18. Visualization of the flow structure was made by the schlieren method. Straight Foucault knife was located in the focal plane of the receiving part of a shadow device. Registration of shadow patterns was carried out using high- speed camera Photron Fastcam (300 000 fps) with an exposure time of 1 μs. The Mach number on the cone was calculated from inclination angle of shock wave in the shadowgraph.


  • hypersonic conic flows
  • light-gas magnetoplasma launcher
  • photometric shadow method

1. Introduction

Currently, light-gas cannons are one of the most effective devices for accelerating shells of relatively large mass to speeds comparable to the orbital. In this area of research, intensive and extensive experimental work is being carried out to study high-speed impact, resistance of materials to high-speed action of solid particles, hypersonic flows around axisymmetric bodies and are closely related to the development of aviation and space technology. In this case, some types of guns with the replacement of the accelerating channel by the Laval nozzle can be effectively used to study hypersonic flows around bodies, since the parameters of the working gas required for high-speed throwing and generation of a hypersonic air flow are identical. High-speed throwing system parameters and technical requirements are described in [1, 2, 3, 4, 5]. Since the 40s of the 20th century the generation of hypersonic flows, development of measurement methods and study of such flows structure over various test bodies were carried out by researchers from all over the world [6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16]. Numerical modeling methods mentioned, for example in [10, 11, 17, 18, 19], were added to the experimental tools of the research with the development of computer technology. However, it should be noted that experimental studies of hypersonic flows around the objects are still preferred. Due to the importance of problems to be solved flight tests are often used for higher authenticity of the experimental data. Usually total and static pressure is measured by means of Pitot-Prandtl probe to determine the Mach number, but in large number of research works measurements are not carried out due to complicated experimental techniques. Researchers [6, 7, 11, 12, 13, 14, 15, 16] used the initial Mach number calculated from Laval nozzle geometrical parameters. Pressure measurement using a Pitot-Prandtl probe and temperature measurement using a thermocouple significantly distort the flow structure, complicate the design due to the need to place measuring elements inside and on the surface of the models under study, which, in the absence of flow visualization, can introduce uncontrollable errors in the experimental results. In hypersonic aerodynamic facilities preferable methods of research are optical methods with high-speed cameras as the reg-istering devices. Methods of optical visualization of supersonic flows include a method of a laser knife, smoke visualization, visualization by means of the spark discharge, etc. [13, 15, 16, 18]. Optical methods based on the relationship between the optical and thermodynamic properties of the medium (shadow methods, interference methods, spectral diagnostics) allow measurements in hypersonic flows with high accuracy. In some cases, to visualize and measure the temperature field on the surface of the model, a thermosensitive coating is used as an alternative to measurements using thermocouples [8, 9]. Research in the field of high-speed impact and diagnostics of hypersonic flow are often interrelated tasks, which makes it necessary to create experimental installations equipped with systems for optical visualization and diagnostics of the flow and measurement of the Mach number, which makes it possible to study hypersonic flow around objects without suppressing the investigated gas flow.


2. Experimental equipment for hypersonic flow generation and optical diagnostics

In most cases, to implement a hypersonic gas flow with a high Mach number, a shock tube equipped at the end with a diaphragm and a Laval nozzle is used [18, 19, 20]. In this case, the hot gas behind the front of the reflected shock wave flows out through the nozzle into the vacuum reservoir at a hypersonic speed. However, due to the large expansion of the gas, its density turns out to be rather low, which greatly complicates the visualization of the flow. Thus, it is possible to estimate the angles of deflection εxand εyof the probe light beam in geometrical optics approach while schlieren method of hypersonic flow visualization is used [21].


At the same time gradients of gas refraction index nxyyand nxyxare very small when the light gas outflows into vacuum. If helium is used as a light gas and the residual pressure in the vacuum chamber is 1 Torr and the characteristic scale of the streamlined body is ∼10−2 m, then the characteristic gradients of the refractive index do not exceed 5×106m1. Under these conditions, the use of shadow methods and interference diagnosis becomes almost impossible. Therefore, in such cases, it is necessary to use the methods of multipass interferometry, multibeam interferometry, interferometry in polarized light, holographic methods or shadow methods with a large focal length and special imaging diaphragms [22]. On the other hand, the use of light-gas launcher where accelerating chan-nel is replaced with Laval nozzle allows to receive a hypersonic flow with high optical density and to apply Tepler method to its visualization. Hypersonic flow facility studied in this research is a modified light-gas launcher we used for ballistic tests and described in details in article [23]. The scheme of the experimental setup is shown in Figure 1. This light-gas gun allows to accelerate balls with a diameter of 2.5 mm to 4 mm made of high-alloy steel to speeds of 2.5–4 km/s.

Figure 1.

Scheme of experimental facility.

During the experiments, the light gas section was filled with helium to a pressure of 40 bar. Experiments on hypersonic flow around cones were carried out in a vacuum chamber with a residual gas pressure of 1 Torr. As a result of a high-current discharge in a magnetoplasma accelerator filled with gunpowder, the piston is set in motion and causes an adiabatic compression of the light gas. The capacitor bank of 1 200 μF was charged to the voltage 4.5 kV. In front of the confuser of Laval nozzle 5 brass diaphragms 100 μm thick each were placed. The rupture of diaphragms occurred when pressure in light-gas section reaches ∼1 600 bars. Compression rate and temperature of working gas can be calculated from the Poisson equation:


where Tis gas temperature, Vis gas volume, Pis gas pressure and γ=53for the monoatomic gas. Subscript 1 corresponds to initial state of gas and subscript 2 is for state of gas at diaphragms rupture. When diaphragms rupture compression rate of the working gas is about 10 and temperature is about 1500 K. With the simultaneous outflow of gas through the critical section of the Laval nozzle, further compression of the gas continues until the piston reaches the end point. The maximum compression ratio is about 50 and can be calculated as the ratio of the volume of the light gas section to the volume of the confuser of the Laval nozzle. The ratio of the throat to the outlet of the diffuser in the Laval nozzle for a Mach number of 18 was calculated using the formula [22, 24]:


where Ais area of Laval nozzle exit section and Ais area of Laval nozzle throat. The test cones were fixed coaxially with the nozzle at a distance of 20 mm from the nozzle outlet. To visualize the hypersonic gas flow, the shadow knife-and-slit method was used. Width of a slit was 0.16 mm. Foucault’s knife was placed in the focal plane of a receiving part of the shadow device and shaded a half of the image of a slit. The 150 W halogen lamp was used as a light source. The focal length of the shadow device was 1 000 mm. Registration of shadow patterns was carried out using high-speed camera Photron Fastcam (300 000 fps) with an exposure time of 1 μs (Figure 2).

Figure 2.

Shadowgraph images of a hypersonic flow over cone with half-angleτ2=12°.


3. Experimental results and discussion

The analysis of shadow pictures sequence shows that at initial time the flow contains some particles formed as a result of diaphragm destruction and about each of which Mach cone (Figure 2) is observed. Further there is an increase of speed of a flow during 700 μs and its subsequent recession during 3 400 μs. At the maximum speed of the flow in shadow pictures the most acute angle of a shock wave inclination is observed. In Figures 3 and 4 shadow pictures of a hypersonic flow over cones with half-angles τ1=3°and τ2=12°are presented at maximum speed of the hypersonic flow.

Figure 3.

Shadowgraph of hypersonic flow over the cone with a half-angleτ1=3°.

Figure 4.

Shadowgraph of hypersonic flow over the cone with a half-angleτ2=12°.

Correlation of Mach number and shock wave inclination angle observed on shadow images of hypersonic flow over sharp cone was described in details by [25] for the first time. Calculation of Mach number was carried out on the basis of measurement of shock wave inclination angle σobserved on shadow pictures (Figures 3 and 4) using formula [26]:




The Mach number calculation results were M1 = 18 for hypersonic flow over cone with half-angle τ1=3°and M2=14.4for the cone with τ2=12°.


4. The analysis of shadowgraphs and calculation of helium density field at a hypersonic flow over a cone

As it is known from the theory of the Schlieren technique [27], relative change of light intensity ΔI/I0in geometrical optics approach in case of small deflection angles of the light εis directly proportional to the magnitude of the deflection angle of the light rays in the optical inhomogeneity and can be defined by a ratio


where d = 0:16 mm is a width of a slit, f = 1 000 mm is a focal length of optical system, I0is light intensity in the absence of an optical inhomogeneity, Iis light intensity in the presence of an optical inhomogeneity. Since the streamlined body has an axial symmetry, when processing shadow pictures the assumption of axisymmetric structure of the flow was used. The coordinate system is chosen in such a way that the x-axis is perpendicular to the edge of the imaging knife, the z-axis is directed along the probing light beam of the shadow device, and the y-axis is along the symmetry axis of the cone and the flowing stream. The origin is set at the apex of the cone. In this coordinate system, the Abel equation can be written in the form (where R is the radius of the inhomogeneity) [21, 27, 28]:


For calculations we use Abel equation transformation where dlnn/dris expressed explicitly:


After repeated integration [21, 27, 28] the distribution of refraction index is received:


where n0is known in advance value of refraction index on the line of integration, for example in the region of undisturbed flow n0=nR1.

The flow area was divided into annular zones, the number of which was selected in accordance with the resolution of the video camera matrix. The deflection angle within the annular zone riri+1is assumed to be constant. Then the expression (9) can be written as a sum of elementary integrals:


The distribution of the refractive index in the flow is calculated as a result of independent integration at each section separately. With such processing of shadow, characteristic distortions arise in the distribution of the refractive index, which are associated with a high level of noise in the experimentally measured light intensity along the y axis. To eliminate random noise, a preliminary smoothing of the luminous intensity function along the y-axis was carried out using the least squares method [29, 30, 31, 32]:


At the same time the smoothed function values should differ from the experimental function values not more than the standard deviation δof the probe light intensity I0noise. Calculation of helium density ρwas carried out with use of Gladstone-Dale Equation [21, 28].


where KHe=0.19607cm/g3is Gladstone-Dale constant for helium. The value for the given gas and given wavelength is considered to be a constant in the wide pressure range. Its value was calculated for gas under normal conditions with use of tabular values of density and refraction index [21]. As far as the initial values n0and ρ0are arbitrarily chosen values, we can calculate Δρfrom the following equation:


Calculation results for helium density fields at a hypersonic flow over the cones with a half-angle τ1=3°and τ2=12°are shown in Figures 5 and 6.

Figure 5.

Results of density field calculations (hypersonic helium flow past cone with the half-angleτ1=3°).

Figure 6.

Results of density field calculations (hypersonic helium flow past cone with the half-angleτ2=12°).


5. Conclusions

The use of light gas ballistic systems, in which the projectile acceleration channel is replaced by a Laval nozzle, makes it possible to generate hypersonic flows with a high Mach number. The gas density of such flows is sufficient for photometric optical measurements by shadow or interference methods. This makes it possible to study hypersonic flow models of any complexity. In addition, due to the high density of the hypersonic flow, it becomes possible to experimentally study the destructive effect of the hypersonic gas flow on the bows of hypersonic missiles and aircraft.


  1. 1. Bulah, B. M. 1970 Nonlinear Conic Gas Flows.Moscow Nauka[In Russian] 343p.
  2. 2. Golubyatnikov, A. N. , Pilyugin, N. N. & Leontiev, N. E. 2003 Methods of Up-grading the Efficiency of Light-gas Guns.Usp. Meh. vol. 2 (1), pp. 97–124.
  3. 3. Smith, F. 1963 Theory of a Two-Stage Hypervelocity Launcher to Give Constant Driving Pressure at the Model.J. Fluid Mech. vol. 17 (1), pp. 113–125.
  4. 4. Putzar, R. & Schaefer, F. 2015 EMI’s TwinGun Concept for a New Light-Gas Gun Type Hypervelocity Accelerator. Procedia Eng. , vol. 103, pp. 421 – 426.
  5. 5. Putzar, R. & Schaefer, F. 2016 Concept for a New Light-Gas Gun Type Hypervelocity Accelerator. Int. J. Impact Eng.88, 118 – 124.
  6. 6. Abramovich, G.N. 1973 Applied Gas Dynamics. Defense Technical Information Center. ed. 3rd.
  7. 7. Rogers, E. W. E., Berryand, C. J. & Davis, B. M. 1971 Experiments with Cones in Low-Density Flows at Mach Numbers near 2. Her Majestys Stationery Office3505.
  8. 8. Hubner, J. P. & Carroll, B. F. & Schanze, K. S. & Ji, H. F. & Holden, M. S. 2001 Temperature- and Pressure-Sensitive Paint Measurements in Short-Duration Hypersonic Flow.AIAA J./vol. 39(4), 564–659.
  9. 9. Kurita, M. & Nakakita, K. & Mitsuo, K. & Watanabe, S. 2006 Temperature Correction of Pressure-Sensitive Paint for Industrial Wind Tunnel Testing.J. Aircr./vol. 43(5), 1499–1505.
  10. 10. Kussoy, M. I. & Brown, J. D. & Brown, J. L. & Lockman, W. K. & Horstman, C. C. 1988 Fluctuations and Massive Separation in Three-Dimentional Shock-Wave/Boundary-Layer Interactions.Tech. Rep. 89224. NASA Tech. Mem.
  11. 11. Brown, J. D.& Brown, J. L. & Kussoy, M. I. 1988 Fluctuations and Massive Separation in Three-dimentional Shock-Wave/Boundary-Layer Interactions.Tech. Rep. 101008. NASA Tech. Mem.
  12. 12. Horstmanand, C. C. & Kussoy, M. I. 1992 Intersecting Shock-Wave/ Turbulent Boundary-Layer Interactions at Mach 8.3.Tech. Rep. 103909. NASA Tech. Mem.
  13. 13. Radzig, A. N. 2004 Experimental Hidroaeromechanics.Moscow MAI[In Russian].
  14. 14. Leutin, P. G. 1976 Preasure Distribution on Sharp Cones at Angle of Sttackα10°in Supersonic Flow.TsAGI Sci. J. vol. 7 (2), 163–166.
  15. 15. Ridyard, H. W. 1954 The Aerodynamic Characteristics of Two Series of Lifting Bodies at Mach Number 6.86.Tech. Rep. L54C15. NASA Tech. Mem.
  16. 16. Warren, W. R. & Kaegi, E. M. & Harris, C. J. & Geiger, R. E. 1962 Shock Tunnel Studies of the Aerodynamics of Atmospheric EntryTech. Rep. R62SD56. General Electric CO.
  17. 17. Zhivotnov, S. D. & Nikolaev, V. S. & Provotvorov, V. P. 1999 Aerodynamic Characteristics of Cones at Attack Angle Under Viscid-Nonviscid Interaction at Supersonic Velocities.TsAGI Sci. J. vol. 30 (3–4), 97–105.
  18. 18. Erdem, E. & Yang, L. & Kontis, K. 2009 Drag Reduction by Energy Deposition in Hypersonic Flows. In16th AIAA/DLR/DGLR Int. Space Planes Hypersonic Syst. Technol. Conf. 19-22 Oct., Bremen, Germany.AIAA Paper 2009-5347.
  19. 19. Saravanan, S. & Jagadeesh, G. & Reddy, K. P. J. 2009 Convective Heat-Transfer Rate Distributions Over a Missile Shaped Body Flying at Hypersonic Speeds.Exp. Therm. Fluid Sci. vol. 33 (4), 782–790.
  20. 20. Gaydon, A. G. & Hurle, I. R. 1963 The Shock Tube in High Temperature Chemical Physics. Verlag Chapman and Hall Ltd.
  21. 21. Vasilev, L. A. 1966 Shadow Methods. [In Russian] Nauka.
  22. 22. Belozerov, A. F. 2007 Optical Methods of Gas Flow Visualization. [In Russian] Kazan State Techn. Univ.
  23. 23. Khramtsov, P. P. & Vasetskiy, v. A. & Grishchenko, V. M. & Makhnach, A. I. and Chernik, M. Yu. and Shikh, I. A. 2015 Experimental and Analytical Study of a Two Stage Light-Gas Magnetoplasma Launcher for Ballistic Tests under Vacuum Conditions.High Temp. Mater. Processesvol. 19 (3–4), 209–219.
  24. 24. Anderson, J. D. 2001 Fundamentals of Aerodynamics. McGraw-Hill.
  25. 25. Liepmann, H. W. & Roshko, A. 1957 Elements of Gas Dynamics. John Wiley & Sons Inc.
  26. 26. Pottsepp, L. 1960 Inviscid Hypersonic Flow Over Unyawed Circular Cones.J. Aerosp. Sci. vol. 27 (7), 558–559.
  27. 27. Taylor, H. S. and Ladenburg, R. W. 1955 Physical Measurements in Gas Dynamics and Combustion. Oxford University Press.
  28. 28. Skotnikov, M. M. 1976 Qualitative Shadow Methods in Gas Dynamics. [In Russian] Nauka.
  29. 29. Savitzky, A. and Golay, M. J. E. 1964 Smoothing and Differentiation of Data by Simplified Least Squares Procedures.Anal. Chem. vol. 36 (8), 1627–1639.
  30. 30. Steinier, J. & Termonia, Y. and Deltour, J. 1972 Comments on Smoothing and Differentiation of Data by Simplified Least Square Procedure.Anal. Chem. vol. 44 (11), 1906-1909.
  31. 31. Bockasten, K. 1961 Transformation of Observed Radiances into Radial Distribution of the Emission of a Plasma.J. Opt. Soc. Am. vol. 51 (9), 943–947.
  32. 32. Gander, W. and Hrebicek, J. 2004 Solving Problems in Scientific Computing Using Maple and MATLAB. Springer-Verlag Berlin Heidelberg.

Written By

Pavel P. Khramtsov

Reviewed: July 14th, 2021 Published: December 10th, 2021