# Triple configurations of pursuit shock waves in conditions of ambiguity of the solution

### Abstract

#### For citation:

Chernyshov M.V., Kapralova A.S. Triple configurations of pursuit shock waves in conditions of ambiguity of the solution. *Journal of «Almaz – Antey» Air and Space Defence Corporation*. 2017;(4):46-52.
https://doi.org/10.38013/2542-0542-2017-4-46-52

## Introduction

Triple configurations of shock waves, present in stream and nozzle gas flows implemented in jet aviation and rocketry technologies, affect the performance of supersonic air intakes and other equipment based on jet flow technologies.

At the present time, developers keep on searching for effective solutions for ramjet, rotating, and pulse detonation engines, so the problem of analysing interaction between compression shocks, shock and explosive waves is especially relevant. To solve this problem, it is necessary to analyse all triple configurations that may form in given conditions, depending on device parameters. It is also important to analyse a variety of possible solutions to the problems regarding development of explosive-proof equipment, determination of damage effect caused by condensed substance explosion associated with irregular interaction of air shock waves and their Mach reflection.

This paper briefly reviews the properties of optimal triple configurations that correspond to the maximum variations in the parameters of flows after them and are inherent to the basic and alternative solutions within the framework of the local triple shock theory. In this respect, when searching for the optimal flow conditions in triple configurations, we must consider the ambiguity in solution to the constitutive system of equations.

## General information on triple configurations

A triple configuration of compression shocks is a shock-wave structure consisting of three shock waves with a common triple point (point T in Fig. 1). Triple configurations of shock waves, being stationary within a selected coordinate system (compression shocks), are present in stream and nozzle gas flows implemented in jet aviation and rocketry technologies [1–3]. They affect the performance of supersonic air intakes and other equipment based on jet flow technologies [4, 5]. Triple configurations of moving (travelling) shock waves appear under their Mach reflection and irregular interaction [6–10], affecting the efficiency of the mechanical impact of an explosion, as well as the performance of explosion-proof devices intended to suppress the high-explosive effect [11–13], in particular, in case of multiple interaction of shock waves in confined volumes [14–16]. Gas flows having passed through different compression shock wave systems (sequence of shocks 1 and 2 or single shock 3) are separated by tangential discontinuity τ. The parameters of shocks are connected by the conditions of consistency on tangential discontinuity, written in the form [17–19]:

J_{1} J_{2} = J_{3} ; (1)

β_{1} + β_{2} = β_{3} . (2)

**Fig. 1.** Types of triple configurations: а – TC-1; b – TC-2; c – TC-3; d – SMC; e – TTC-2-3

Here, J_{i} (i=1...3 ) – intensity of the i-th shock (ratio of static pressures after and before the shock);

β_{i} – flow turn angle on the surface of the i-th shock.

Angles β_{i} and Mach numbers M_{i} after the i-th shock are associated with shock intensity and Mach number M_{k} before the shock by the known [1] classic relationships.

Depending on the direction of flow turn on shocks 1–3, three types of triple configurations are distinguished. In configurations of the first type (TC-1, Fig. 1, a), flow turn on shock 1 occurs in a direction different from that on shocks 2 and 3. For example, at β_{1} < 0 , angles β_{2} > 0 , β_{3} > 0. In configurations of the second type (TC-2, Fig. 1, b), the direction of turn on shock 2 is different from the others, and in configurations of the third type (TC-3, Fig. 1, c), flow turn occurs in the same direction on all the shocks. The stationary Mach configuration (SMC, Fig. 1, d) with direct main shock (β_{3} = 0) and configuration TTC-2-3 (Fig. 1, e) with direct shock 2 (β_{2} = 0) are transient.

Setting adiabatic index γ, Mach number M of flow before the triple point, and branching shock intensity J_{1} does not always allow to explicitly define the properties of other shocks in the system of equations (1)–(2). The same parameters γ , M , and J_{1} are matched by up to three physically based solutions with different values of β_{2} and β_{3} . The basic solution to the system of equations (1)–(2) is defined in the widest area of parameter space (γ, M, J_{1} ), and two alternative solutions – only at subsets of the region of the basic solution definition. Triple configurations corresponding to the basic solution may belong to all three types, as well as to the transient configurations SMC and TTC-2-3. Alternative triple configurations (ATC) relate to the third type (see Fig. 1, c), and the flow after shock 2 is supersonic at that. They are formed as a result of interaction between pursuit shocks.

Many parameters of gas flows after triple configurations have substantial differences. Those of interest are the differences in stagnation pressures p_{0}, velocities V, flow rate functions q = ρV , flow strength d = ρV^{2} , flow pulses j = ρ + ρV^{2} after the triple point. The measure of difference here are their ratios on the tangential discontinuity sides. Triple configurations with extremal values of such ratios are called optimal configurations. Investigation of the optimal configurations may have practical importance in analysing the origin of self-oscillation regimes of flows in free and impact supersonic jets [20], when designing the equipment generating pulsating gas flows.

Further, we shall analyse the properties of optimal triple configurations corresponding to both the basic and alternative solutions. The numerical results are given for γ = 1, 4.

## Optimal configurations corresponding to the basic solution

The properties of triple configurations of compression shock waves are analysed on the plane of parameters M and σ_{1} (Fig. 2), where σ_{1} – angle of shock 1 to the direction of flow before the shock. Angle σ_{1} correlates with shock intensity J_{1} as

J_{1} = (1 + ε) M^{2} sin^{2} σ_{1} - ε,

where ε = (γ − 1) / (γ + 1).

**Fig. 2.** Parameters of triple configurations

The range of angles σ_{1} variation is limited from below by curve 1, which corresponds to shock transformation into a weak discontinuity (σ_{1} = α(Μ) = arcsin (1/M), J_{1} = 1). The values σ_{1} are limited from above as well, at least by the requirement that shock 2 must exist in the supersonic flow after the branching shock. This requirement is conformed to by the region under curve 2, which is plotted proceeding from condition M_{1} = 1 after shock 1.

For the existence of triple configuration, presence of a supersonic flow after shock 1 is insufficient. A solution to the system (1)–(2) exists only in the region between curves 1 и f_{1} therefore curve f_{1} is the exact upper boundary of the region under consideration and is defined by the equation common for curves f_{i} (i = 1, 2):

Curve f_{1} starts in point F_{1} on curve 1, where M_{F1} = 1.245, M_{F2} = 2.54.

Solutions to the system (1)–(2) in the region under consideration may correspond to configurations of different types. In the subregion between curves 1 and 3, configurations of TC-1 type are implemented; between curves 3 and 4 – TC-2, and between curves 4 and f_{1} – TC-3. Curve 3 corresponds to the stationary Mach configuration and is plotted from solution to the equation

where a = (1− ε)(1+ ε J_{1}) ;

J_{m} – intensity of direct shock, formed in the flow with a given Mach number,

J_{m} =(1 + ε)Μ^{2} -ε;

b = — [(1 + ε — ε^{2} + ε^{3} )J_{1}^{2} + ε(1 + ε) J_{1} + (1 — ε)];

с = J_{1} ((1 - ε^{2}) J_{1}^{2} - (1 + ε^{2}) J_{1} - 2ε).

Transient configurations TTC-2-3 (curve 4) are determined analytically as well [18, 19].

Intensities and other parameters of compression shock waves across the entire region of basic solution existence are changing continuously. The parameters of individual shocks take extremal and specific values (e.g., shocks 2 and 3 may correspond to the maximum deviation points, Crocco points, constant pressure points, and sonic point [19]).

The properties of flows after the triple configuration are determined from the system (1)–(2) and ratios on the compression shocks. For example, the ratios between total pressures p_{0}, velocities V, flow rate functions q, flow strengths d, flow pulses j at tangential discontinuity are as follows:

where – inverse ratio of densities on the shock.

The lower boundary of the region of solution existence (curve 1) corresponds to transformation of shock 1 and the upper boundary (curve f_{1}), of shock 2 – into a weak discontinuity. In these cases, all the considered parameter ratios after the triple point are equal to unity. At fixed Mach number M , the only extremum point of the considered functions in a range between the definition region boundaries is the maximum point. The configurations corresponding to these maxima are optimal at fixed Mach number.

The parameters of configurations optimal with respect to the target functions (4) are shown in Fig. 2 by curves 5–9, respectively. At low Mach numbers, these configurations belong to the third type. The intersections of curves 5–9 with curve 4 corresponds to the optimal transient configurations. In this case, the optimal ratios of parameters (I_{po} = 1,076; I_{v} = 1,085; I_{q} = 1,107; I_{d} = 1,201; I_{j} = 1,090) are low, and the ratios of Mach numbers (M = 1.596; M = 1.567; M = 1.571; M = 1.569; M = 1.584) are very close.

With an increase of Mach numbers, the optimal curves 5–9 come close and intersect in one point a, which corresponds to the stationary Mach configuration (SMC) with Mach number . The intensities of incident shock 1 and reflected shock 2 of the compression shock waves in such SMC are equal to:. It is proved [21, 22] that equality of the shock wave intensities leads to the total pressure maximum after the shock-wave system if the product of those intensities is a fixed value. It can be shown that in the SMC such a product ( J_{3} intensity), even though it is not a fixed value, obeys the above theorem, therefore it is exactly the Mach configuration with equal shock intensities that is the optimal one. The parameter ratios after the optimal SMC:

At greater Mach numbers, configurations of the first type are optimal. The optimal values of target functions increase monotonically but in a limited way, while the optimal intensities of shocks 1 and 3 at M → ∞ tend to infinity. The configurations optimal in terms of I _{p0} have the following finite limits:

The Mach number after shock 1 tends to the infinite limit (order √M) and after shocks 2 and 3 – to finite limits. Ratio I _{p0} itself tends to the value

The limit values of other functions in configurations optimal in terms of I _{p0} , are as follows:

and, as a rule, they are close to the optimal values reached on curves 6–9 (see Fig. 2): I_{V} → 5.261, I_{d} → 155.8, I_{q} → 30.41, I_{j} → 30.22; therefore, optimisation of configurations with respect to these parameters is sometimes substituted for optimisation as per I _{p0} [18]. In configurations optimal with respect to these four parameters, intensity J_{1} has order M_{2} , and values M_{1} and J_{2} tend to high finite values. The angle of shock 1 tends to a low finite value at that, rather than to zero.

The optimal values (especially, the total pressure ratios) tend to their limits slowly: at M = 8 , optimal I _{p0} = 19.36, and at M = 200 , I _{p0} = 439.2. The optimisation of configurations leads to notable increase of the target functions. Thus, at M → ∞ , optimal I _{p0} → 529.1, while I _{p0} → 69.72 in the SMC and I _{p0} → 1 in TTC-2-3.

## Alternative triple configurations

Starting from certain Mach number (M = 2,542 in optimisation as per the ratio of total pressures), parameters ( M , σ_{1} ) of the optimal basic configurations determine two more solutions, and at M > 2.61 – one solution that describes ATC of 2 6, 1 the third type corresponding to one of the alternative solutions, which exist along with the basic one at the same Mach numbers of the intensity flow of shock 1 (branching) and gas adiabat.

The alternative solutions to the system (1)–(2) appear on curve bc (see Fig. 2) as a result of decomposition of shock isomachs [13]. There are two different ATCs in curvilinear triangle F_{2}cb . One of the solutions in segment F_{2}c of curve f_{2} corresponds to value J_{1} < 1, due to which it stops being implemented. At the same time, a new and the only possible solution for ATC appears on curve f_{2} , after point c . Curve f_{2} and point F_{2} are defined by relationship (3), and points b ( ) (M_{b} = 2.089) and c ( M_{c} = 3.117) – by high-degree (for point b – eighth) algebraic equations.

The maxima of equations (4) can be achieved in the ATCs corresponding to the solution which is continuous across the entire region beyond curves bc and f_{2} (curves 10–14). At M → ∞ , the optimal value of I _{p0} tends to the limit (5) and can be achieved at J_{3}/M^{2} → C_{1}, J_{1}/M → √C_{1}, J_{2}/M → √C_{1} . The flow turn angle on shock 3 in an optimal asymptotic ATC is opposite to its value in the “basic” configuration.

The limits of other parameter ratios in the optimal ATCs are at least comparable to the “basic” configurations: in the ATC, at M = 199.3, maximal I_{V} = 4.858, I_{d} = 133.1, I_{q} = 27.47, I_{j} = 28 and in the “basic” configurations, I_{V} = 5,257, I_{d} = 151, I_{q} = 29.23, I_{j} = 28.56. The relative position (from bottom to top) of the optimal curves 10–14 is opposite to the position of curves 5–9 at high Mach numbers.

With parameter γ increased, Mach numbers at which the ATCs are formed increase substantially and tend to infinity at γ → 5/3. At γ ≥ 5/3, the system of equations (1)–(2) has no more than one physically based solution.

## Conclusion

The conducted calculation and parametric analysis of triple configurations forming under all theoretically feasible flow parameters before them serve for optimisation of systems and devices that employ the effects of interaction and reflection of compression shock waves, blast shock waves, and detonation waves.

The study demonstrates that triple configurations corresponding to different physically feasible solutions can be optimal: after such configurations, the maximum and quite high ratios of total pressures, velocities, flow strengths, and other flow parameters can be achieved on different sides of tangential discontinuity originating from the triple point. This statement holds true both for the basic (traditionally considered) and additional (alternative) solutions defining triple configurations, therefore, when searching for the optimal flow conditions in triple configurations, it is necessary to consider the ambiguity of solution to the defining system of equations.

The results obtained using theoretical and numerical methods can be used in various applications of gas dynamics. For instance, high total differential pressures in a supersonic gas jet initiate self-oscillation regimes when a jet interacts with obstacles, and lead to extreme acoustic and force loads when executing starting tasks. The different translational (transferred) impact of blast shock waves on bodies located above and below the triple point is achieved due to a considerable difference between flow strengths on the opposite sides of tangential discontinuity. This phenomenon can be used in design of explosion-proof devices and in an analysis of blast effect (especially in confined spaces with inevitable multiple reflection of shock waves and their irregular interaction). Moreover, high values of flow parameters after triple configurations hamper initiation of detonation in aircraft and rocket engines of appropriate type and shall be eliminated at the development phase of such devices.

## References

1. Дулов В. Г., Лукьянов Г. А. Газодинамика процессов истечения. Новосибирск: Наука, 1984. 234 с.

2. Омельченко А. В., Усков В. Н., Чернышов М. В. Об одной приближенной аналитической модели течения в первой бочке перерасширенной струи // Письма в журнал технической физики. 2003. Т. 29. Вып. 6. С. 56–62.

3. Silnikov M. V., Chernyshov M. V., Uskov V. N. Two-dimensional over-expanded jet flow parameters in supersonic nozzle lip vicinity // Acta Astronautica. 2014. Vol. 97. Pp. 38–41.

4. Oswatitsch K. Gas Dynamics. New York: Academic Press, 1956. 610 p.

5. Герман Р. Сверхзвуковые входные диффузоры / под ред. Г. Н. Абрамовича. М.: ФИЗМАТГИЗ,1960. 290 с.

6. Баженова Т. В., Гвоздева Л. Г. Нестационарные взаимодействия ударных волн. М.: Наука, 1977. 274 с.

7. Усков В. Н., Мостовых П. С. Тройные конфигурации бегущих ударных волн в потоках невязкого газа // Прикладная механика и техническая физика. 2008. Т. 49. № 3. С. 3–10.

8. Гельфанд Б. Е., Сильников М. В. Химические и физические взрывы. Параметры и контроль. СПб.: Полигон, 2003. 416 с.

9. Гельфанд Б. Е., Сильников М. В. Баротермическое действие взрывов. СПб.: Астерион, 2006. 658 с.

10. Гельфанд Б. Е., Сильников М. В. Фугасное действие взрывов. СПб.: Астерион, 2007. 252 с.

11. Гельфанд Б. Е., Сильников М. В. Выбор оптимальной схемы подавления воздушных ударных волн при взрыве ВВ // Доклады Академии наук. 2002. Т. 383. № 1. С. 37.

12. Гельфанд Б. Е., Сильников М. В., Компан Ф. М., Чернышов М. В. К вопросу об эффективности конструктивных схем зарубежных локализаторов взрыва // Вопросы оборонной техники. Сер. 16. Технические средства противодействия терроризму. 2010. Вып. 9–10. С. 3–10.

13. Gelfand B. E., Silnikov M. V., Chernyshov M. V. On the efficiency of semi-closed blast inhibitors // Shock Waves. 2010. Vol. 20. No. 4. Pp. 317–321.

14. Гельфанд Б. Е., Сильников М. В., Михайлин А. И., Чернышов М. В. Защита широкофюзеляжного самолета от взрывных нагрузок // Проблемы управления рисками в техносфере. 2009. Т. 9–10. № 1–2. С. 21–31.

15. Сильников М. В., Михайлин А. И., Чернышов М. В., Шишкин В. Н. Защита узкофюзеляжного воздушного судна от поражающего действия внутреннего взрыва // Известия Российской академии ракетных и артиллерийских наук. 2011. № 1 (67). С. 18–27.

16. Silnikov M. V., Mikhaylin A. I. Protection of flying vehicles against blast loads // Acta Astronautica. 2014. Vol. 97. Pp. 30–37.

17. Адрианов А. Л., Старых А. Л., Усков В. Н. Интерференция стационарных газодинамических разрывов. Новосибирск: Наука, 1995. 180 с.

18. Усков В. Н., Чернышов М. В. Теоретический анализ особенностей тройных конфигураций скачков уплотнения // Современные проблемы неравновесной газо- и термодинамики. СПб.: Изд-во БГТУ «ВОЕНМЕХ», 2002. С. 75–99.

19. Усков В. Н., Чернышов М. В. Особые и экстремальные тройные конфигурации скачков уплотнения // Прикладная механика и техническая физика. 2006. Т. 47. № 4. С. 39–53.

20. Горшков Г. Ф., Усков В. Н. Автоколебания в сверхзвуковых перерасширенных импактных струях // Прикладная механика и техническая физика. 2002. Т. 43. № 5. С. 49–54.

21. Омельченко А. В., Усков В. Н. Оптимальные ударно-волновые системы // Известия РАН. Сер. Механика жидкости и газа. 1995. № 5. С. 118–126.

22. Омельченко А. В., Усков В. Н. Оптимальные ударно-волновые системы при ограничениях на суммарный угол поворота потока // Известия РАН. Сер. Механика жидкости и газа. 1996. № 4. С. 142–150.

### About the Authors

**M. V. Chernyshov**Russian Federation

Chernyshov Mihail Viktorovich – Doctor of Engineering Sciences, Professor, Department of Plasmosdynamics and Heat Engineering, Faculty of Rocket and Space Equipment, deputy vice-rector. Science research interests: gas dynamics, shock and blast waves, interaction of gas-dynamic discontinuities, explosion protection.

Saint Petersburg

**A. S. Kapralova**Russian Federation

Kapralova Anna Sergeevna – post-graduate student, Department of Aircraft Engines and Power Plants, Faculty of Rocket and Space Equipment. Science research interests: shock and blast waves, explosion protection.

Saint Petersburg

### Review

#### For citation:

Chernyshov M.V., Kapralova A.S. Triple configurations of pursuit shock waves in conditions of ambiguity of the solution. *Journal of «Almaz – Antey» Air and Space Defence Corporation*. 2017;(4):46-52.
https://doi.org/10.38013/2542-0542-2017-4-46-52