Next Article in Journal
Comparative Study of a Powerplant Life Consumption Rate When Installed in Two Different Aircraft Variants
Next Article in Special Issue
Spillage-Adaptive Fixed-Geometry Bump Inlet of Wide Speed Range
Previous Article in Journal
Multimodal Estimation of Sine Dwell Vibrational Responses from Aeroelastic Flutter Flight Tests
Previous Article in Special Issue
Aerodynamic Characteristics of Re-Entry Capsules with Hyperbolic Contours
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Redistribution of Energy during Interaction of a Shock Wave with a Temperature Layered Plasma Region at Hypersonic Speeds

1
Federal Research Center “Computer Science and Control” of the Russian Academy of Sciences, Vavilova St. 40, 119333 Moscow, Russia
2
Ioffe Institute, Politekhnicheskaya Str. 26, 194021 St. Petersburg, Russia
3
Faculty of Mechanics and Mathematics, Leninskiye Gory, Main Building, Lomonosov Moscow State University, 119991 Moscow, Russia
*
Author to whom correspondence should be addressed.
Aerospace 2021, 8(11), 326; https://doi.org/10.3390/aerospace8110326
Submission received: 29 September 2021 / Revised: 27 October 2021 / Accepted: 28 October 2021 / Published: 1 November 2021
(This article belongs to the Special Issue Hypersonics: Emerging Research)

Abstract

:
The paper is devoted to the problem of the interaction between a shock wave and a thermally stratified energy source for the purpose of supersonic/hypersonic flow control realization. The effect of the thermally stratified energy source on a shock wave with the Mach number in the range of 6–12 is researched numerically based on the Navier-Stokes system of equations. Redistribution of specific internal energy and volume density of kinetic energy behind the wave front is investigated. Multiple manifestations of the Richtmyer-Meshkov instability has been obtained which has caused the blurring and disappearance of shock wave and contact discontinuity fronts in density fields. A study of the efficiency of using a stratified energy source instead of a homogeneous one with the same value of the full energy is carried out. The agreement with the available experimental data for the shock wave Mach number 6 has been obtained.

1. Introduction

The problems in the interaction of a shock wave with plasma structures can be considered as part of the direction of supersonic hypersonic flow control, since it simulates the interaction of energy sources with a bow shock wave created by an aerodynamic (AD) body in the flow. Currently, the research of flow/flight control via energy deposition is at the forefront of investigations in aerospace engineering. To solve the problems of controlling supersonic/hypersonic flows and shock wave configurations, improving the flow process, reducing aerodynamic drag, controlling the bow shock wave and boundary layers, it was proposed to supply energy both to the incoming flow, and on the surfaces of AD bodies.
In [1], D. Knight presented a review of CFD methods for the issue of studying the shock interactions in nonequilibrium laminar flows. In an overview [2], M. Ahmed and N. Qin presented a comparison of efficiency of various types of devices, mechanical devices (spikes), fluidic devices (opposing jets), and energy deposition devices to control the flowfield structure, drag and lift forces, surface pressure, heating and flow instabilities. A. Russel et al. presented the works on energy deposition into supersonic flows to control shock-wave configurations due to electro-gasdynamic action and Joule heating [3]. In [4], S. Leonov et al. reviewed the plasma base methods for control of flow and shock wave configurations by the creation of different types of gas discharge near surfaces of AD bodies.
For the purpose of flow control, bow shock wave control in particular, various types of gas discharges, as well as the imposition of external electric and magnetic fields can be used. Experimental results of the effect of longitudinal energy supply by a microwave discharge on the standoff and propagation of a bow shock wave are presented in the research of Yu. Kolesnichenko et al. [5]. The laser pulse effect was investigated by P. Tretyakov et al. [6] and by L. Myrabo et al. [7]. Pulsed-periodic high-voltage discharges for this purpose was used by V. Bityurin et al. [8] and D. Roupassov et al. [9]. In [10], K. Kourtzanidis et al. investigated the plasma influence on bow shock wave configuration by means of gas discharge, created near the surface of AD body, and showed the possibility of influence to the shock position by plasma effects.
The possibility of controlling shock waves by magnetohydrodynamic methods was experimentally shown by T. Lapushkina et al. in [11]. In [12], Q. Sun et al. investigated the effect of an arc gas discharge on the position of an oblique shock wave and on the distribution of the flow pressure. The results are explained by thermal effects and modeled in [13] by H. Yan, et al. The thermal effect on a shock wave propagating in the wake of a nanosecond high-energy discharge, called a plasma sheet, is considered by I. Znamenskaya, et al. [14]. O. Azarova obtained the vortex structure of the shock layer region formed under the action of energy source, and proposed a vortex mechanism of action on the body, leading to a decrease in its wave drag. Additionally, it was shown that the initialization of the vortex occurs as a result of the manifestation of the Richtmyer-Meshkov instability [15,16].
The creation of plasma temperature inhomogeneities along the path of propagation of shock waves can lead to a change in its shape and velocity. This has been shown by experiments on the interaction of shock waves with plasma structures. In [17], N. Apazidis et al. used a set of heated thin wires to create thermal and density inhomogeneities of various geometric configurations in air, which led to the generation of the Richtmyer-Meshkov instability and a train of vortices generated by the Kelvin-Helmholtz instability. In [18], V. Bityurin et al. observed the curvature of the shock wave as it passed through the region of a longitudinal pulsed discharge. In [19,20,21] S. Leonov et al. investigated the influence of an inhomogeneously layered plasma on the position of the reflected shock wave. In the experiments, the blurring of the front of the reflected shock wave was obtained due to the organization of a system of plasma filaments created by a high-frequency discharge [20,21]. In [22] T. Gan et al. used an array of sixteen surface arc plasma actuators to impact on the shock wave in the boundary layer which provided the disappearance of the foot part of the separation shock wave.
The curvature of the shock wave and generation of a contact discontinuity caused by the interaction of the shock wave with the glow discharge region were investigated by T. Lapushkina, et al. [23,24]. The effect of an ionization-unstable discharge on a plane shock wave was obtained and numerically studied [24]. Due to this interaction, blurring, curvature and the disappearance of the shock front were obtained. In our works the possibility of the formation of flow regions behind the wave front with significantly increased values of the gas temperature due to the interaction of a thermally stratified energy source with the shock wave front has been shown [25]. In addition, the influence of a thermally stratified energy source on the supersonic flow around an AD body was investigated and a new multi-vortex mechanism of the action of the energy source on the body surface has been established [26]. In our previous article, the redistribution of energy types in a curved shock wave under the conditions of the experiment [24] was estimated for M = 2 and 5 [27].
Our current paper focuses on the interaction of a preliminary plane shock wave with an area of thermally stratified energy source, which is accompanied by generation of instabilities and vortices. The study is based on the same ideas used in [27], but here we consider the redistribution of different types of energy for hypersonic shock wave velocities. The purpose of this work is to establish how this redistribution depends on the shock wave Mach number, which varies from 6 to 12, and on the degree of gas rarefaction in the energy source layers. In addition, we aimed to obtain and study the complicated vortex shock-wave structures, characterized by multiple manifestation of the Richtmyer-Meshkov instabilities. The research can be important for achievement of high-speed flow control and improvement of the conditions for the processes of mixing and ignition in jets and combustion chambers.

2. Numerical Method and Statement of the Problem

The effect of an energy source in the form of a thermally stratified plasma region on an initially flat shock wave is simulated based on the Navier-Stokes system of equations for viscous and heat conductive gas (air) [28] (p. 329).
U ¯ t + ( F ¯ + F v ¯ ) x + ( G ¯ + G v ¯ ) y = 0 ¯
U ¯ = ( ρ ρ u ρ v E ) , F ¯ = ( ρ u p + ρ u 2 ρ u v u ( E + p ) ) , G ¯ = ( ρ v ρ u v p + ρ v 2 v ( E + p ) )
F v ¯ = ( 0 μ / Re ( 4 / 3 u x 2 / 3 v y ) μ / Re ( v x + u y ) μ π 1 / Re + ( 1 / N ) k T x ) , G v ¯ = ( 0 μ / Re ( v x + u y ) μ / Re ( 4 / 3 v y 2 / 3 u x ) μ π 2 / Re + ( 1 / N ) k T y )
π 1 = u ( 4 / 3 u x 2 / 3 v y ) + v ( v x + u y ) , π 2 = v ( 4 / 3 v y 2 / 3 u x ) + u ( v x + u y )
E = ρ ( ε + 0.5 ( u 2 + v 2 ) ) , N = RePr ( γ 1 ) / γ .
Here Re is the Reynolds number (Re = 9500) and Pr is the Prandtl number (Pr = 0.7), ρ, p, u, v are the gas density, pressure, x- and y- are the velocity components; ε is the specific internal energy,
ε = p ρ ( γ 1 ) ,
which corresponds to the state equation of the perfect gas. For the dependence of dynamic viscosity µ on temperature the Sutherland law,
μ = T 1.5 ( 1 + s 1 ) / ( T + s 1 ) ,
with the constant s1 = 0.41 (120K) was used. It was assumed that the coefficient of thermal conductivity k depends on temperature as
k = T0.5.
To describe the thermodynamic properties of the plasma region the “effective value” of the ratio of specific heats γ = 1.2 was accepted. This value of γ corresponds to the degree of non-equilibrium of 0.015 and to the ionization degree of the gas medium of 0.00015 [10].
The problem is solved in the dimensionless variables. Dimensionless parameters for time t, spatial variables x, y, velocity components u, v, sound velocity c, gas density ρ, pressure p, and temperature T are expressed throw the dimensional ones (marked with the index “dim”) as follows:
t = t d i m t n ,     ,   x = x d i m l n ,     , y = y d i m l n ,     , u = u d i m u n ,     , v = v d i m u n ,     , c = t d i m t n ,     ,
ρ = ρ d i m ρ n     , p = p d i m p n     ,   T = T d i m T n     .
Here the following scales for the parameters are accepted:
ρ n = ρ , p n = p , l n = D , T n = T , u n = ( p / ρ ) 0.5 , t n = l n / u n ,
where the index ∞ denotes the parameters of undisturbed flow. Normalizing values for density and pressure, ρn = 0.01205 kg/m3 and pn = 1013.25 Pa, and the Reynolds number Re = 9500 corresponds to the conditions of the experiments in [24].
The effect of the energy source was modeled by a stationary region of heated gas layers of the same width with zero velocities and the reduced density, ρs = αρ. Here, the parameter α characterizes the medium rarefaction, α < 1. The pressure in the source layers was equal to its undisturbed value, so, the temperature was increased in these layers. The gaps between the layers had a width equal to a half of a layer width. Below, for brevity, such a heated area is referred to as an energy source. This heated area is set at the initial time instant at some distance from the shock wave (SW), xsw is the initial shock wave coordinate, and xs is the stratified energy source boundary coordinate (Figure 1a).
Horizontal walls were supposed to be cold with the value of undisturbed flow temperature, which is in accordance with the experimental conditions. As was shown in [23], these boundary conditions provide the curvature of the shock wave front, which is adequate to that obtained in experiments. Additionally, the adhesion conditions and the condition of absence of the flows normal to the walls for the according parameters were employed.
The domestic numerical code was employed based on the conservative difference schemes of the second order of approximation in space and in time [29]. The TVD variant of the scheme was used, with the limiter introducing to the scheme guided by the min-mod approach. The schemes are constructed on staggered grids and use a stencil of the Lax’s scheme. For the second approximation order obtaining the systems for the space derivatives were employed. There were 2000 nodes on the transverse size of the computation domain (1.6 × 106 nodes in the calculation area, counting the central node). Location of the stratified energy source on difference grid in enlarged form is presented in Figure 1b,c.

3. Analysis of the Grid Convergence

Grid convergence analysis was conducted for four difference grids. Characteristics of the grids are presented in Table 1.
Comparison of the calculations on these four different grids of the field of internal energy is presented in Figure 2. Despite the fact that the grids differ considerably, Figure 2a shows that the main elements of the flowfield (source layers, shock wave front, details of the unstable shock-wave structure) appear in the calculations in a similar way. Behind the shock-wave structure, the fields are also almost identical. Dynamics of εmax, calculated on these four difference grids, are also quite close (Figure 2b). These facts indicate the presence of grid convergence for the considered grids. Note that we used the finest grid (Grid1) for the simulations.

4. An Example of the Thermally Layered Gas Media Formation

The organization of a temperature-layered structure of the medium is possible by creating an ionization-unstable gas discharge [24]. In this work, the interaction of a shock wave (with a Mach number of 5–6) with the afterglow zone of a glow gas discharge characterized by the appearance of ionization instability is considered. The ionization instability is accompanied by arising the ionization waves—ionization strata [30] (Figure 3). The electron temperature in the strata is an order of magnitude higher than the gas temperature of the medium and, due to heat transfer, contributes to an increase in the gas temperature in the region of the striation. This is because the discharge with strata exists for quite a long time (several minutes) before the arrival of the shock wave.
As was shown in the schlieren images obtained in the experiments for M = 5–6 [24], when a shock wave passes through an afterglow medium layered in gas temperature, a complete blurring of the shock wave front occurs.

5. Results of the Simulation

In the simulation specific internal energy, ε = p / ( ρ ( γ 1 ) ) , and volume density of kinetic energy, E = ρ ( ε + 0.5 ( u 2 + v 2 ) ) was investigated for the interaction of a thermally stratified energy source with a preliminary plane shock wave of the Mach numbers M = 6, 7, 8, 10, 12. The fields of the parameters are forming as the solution of the 2D Riemann problem describing the interaction of the shock wave with the boundaries of heated layers forming the thermally stratified source. This solution includes a shock wave and contact discontinuity, moving to the left, and a rarefaction wave to the right of them. Here the initial shock wave coordinate was xsw = 1.5 and the stratified energy source boundary coordinate xs = 1.4.
In Figure 4 typical flow dynamics resulting from stratified energy source—shock wave interaction is presented. It can be seen that the fronts of the shock wave and contact discontinuity become wavy and complicated flow structure is generated including the manifestation of instabilities. Dynamics of density, pressure and two types of energy fields for the interaction are presented in Figure 5. Time moments are indicated in the right lower angle in the images. In Figure 6, Figure 7, Figure 8 and Figure 9 the fields of considered parameters are presented for M = 6, 8, 10, 12, accordingly. The fields for M = 6 correspond qualitatively to the situation of the complete disappearance of the shock front at the schlieren pictures in the experiments described in Part 3 [24].
The multiple manifestations of the Richtmyer-Meshkov instabilities which are exhibited as local peaks in the density, pressure and kinetic energy fields can be seen (see, also, [15]). One can see that in the density field the shock wave front is practically disappeared under the action of the multiple generation of the Richtmyer-Meshkov instabilities. This explains the effect of the disappearance of the shock fronts at the obtained schlieren photographs in the case of small-scaled ionization strata—shock wave interaction [24]. In this experiment, the shock wave completely disappeared in the very beginning of the interaction with the afterglow region, which had a thermally stratified nature. It should be noted that, as calculations have shown, for the Mach numbers under consideration, instabilities manifest themselves in the first microseconds after the start of the interaction process, and at this time the shock wave already becomes wavy and blurred.
It is also important that, at sufficiently high Mach numbers, the density changes behind the shock front are limited. Consequently, the density perturbations arising during the development of instability are also limited.
Nevertheless, it should be underlined that the front of the shock wave doesn’t disappear at all in the fields of other parameters despite the presence of the instabilities. At the same time, this instability does not manifest itself in the internal energy fields. This can be explained by the fact that entropy perturbations, whose propagation velocity is less than the speed of sound, make a significant contribution to the variation of ε near the shock front. Furthermore, one can see a series of acoustic waves emitted by the resulting shock-wave structures with instabilities.
Note that the curvature of the shock wave front is due to the boundary conditions on the horizontal walls and for M = 6 it is in accordance with the front curvature at the schlieren pictures for the case of large-scaled ionization strata—shock wave interaction [24].
In Figure 10 the dynamics of maximum value εmax in the obtained shock-wave structures is presented for different rarefaction parameters α for all considered Mach numbers of the shock wave. Here, the maximum value of the internal energy is calculated over the coordinates in the region of the obtained shock-wave structure with instabilities. This value has a sense of the greatest local increase in the internal energy due to the formation of such a structure as a result of the action of a stratified energy source. It can be seen that these dependences have maxima on t at the beginning stage of the interaction. These maximum values are larger for smaller α and for larger M. Figure 11 shows these dependencies for different M for all considered α. It can be seen that at the first stage, the dynamics of εmax is changing stronger for smaller α and larger M. Note that fluctuations in internal energy are directly connected with fluctuations in the temperature of the gas in the flow, therefore, from the results shown in Figure 10 and Figure 11, similar conclusions follow for the gas temperature [25].
The same analysis is presented in Figure 12 and Figure 13 for the dynamics of maximum value Emax in the obtained unstable shock-wave structures. This value has a sense of the greatest local increase in the kinetic energy due to the formation of a shock-wave structure as a result of the use of a stratified energy source. In Figure 12, the dependences Emax(t) are presented for different rarefaction parameter α for the considered shock wave’s Mach numbers. One can see that the dependences have maxima on t which are larger for larger α and larger M. Figure 13 shows these dependencies for different M for all considered α It can be seen from these dependences that the growth of the kinetic energy occurs quicker for larger α and larger M.
Figure 10 and Figure 11 show that the density, pressure and internal energy fluctuations are stronger at the initial stage of the interaction while the kinetic energy maxima are generated at the middle stage in time (Figure 12 and Figure 13). This can be explained by the fact that the kinetic energy maximum is defined not only by the peaks, connected with the unstable nature of the gas density, but by the behavior of the flow velocities. Thus, the kinetic energy has the maximal attitudes of fluctuations when vortices are generating in the flow and the formation of the vortices takes some additional time.
It should be noted that this paper is devoted to the redistribution of energy and both considered types of energy are characterized by possibly achievable maximum values. The histories of their behavior over time are given for these quantities. At the same time, for the internal energy, the considered time interval up to t = 0.12 turned out to be sufficient (since its maximum value is reached at the initial stage in time). Here, the initial shock wave coordinate xsw = 1.5 and the source boundary coordinate xs = 1.4. For the study of kinetic energy, this time interval turned out to be not enough (because its maximum value is reached at the middle stage in time). Here, we used a different geometry of the calculation domain with the position of the shock wave xsw = 2.25 and the boundary of the energy source xs = 2.15 with the same vertical dimensions. This made it possible to study the time history of kinetic energy in the time interval up to t = 0.25 for the considered shock wave Mach numbers.
To evaluate the “energy” efficiency of using a thermally stratified energy source, the fields of parameters for stratified and homogeneous sources were compared. The homogeneous source was defined at the initial time in the form of a homogeneous region of heated gas with the sizes that coincide with the sizes of the stratified energy source. The parameters of a homogeneous energy source were chosen in such a way that the average values of internal energy in the stratified and homogeneous sources were equal:
ε h   a v e r a g e d = ε a v e r a g e d = ( N i N j ) 1 i , j ε i , j .
Here Ni and Nj—the amounts of grid nodes in i- and j-directions, εaveraged and εh averaged are the averaged values of internal energy in the stratified and homogeneous sources, accordingly.
It is easy to conclude that in this case, the values of full energies for these energy sources are also equal (since the velocity components in the energy sources are equal to zero). Thus, as the result, we can evaluate the transformation of different types of energy only due to the redistribution of a source energy into layers.
Figure 14 demonstrates typical fields of internal (Figure 14a) and kinetic (Figure 14b) energy for the homogeneous energy source (compare with the images for ε and E in Figure 7 for the stratified energy source). The fields in Figure 14 are forming as the solution of the 2D Riemann problem describing the interaction of a preliminary plane shock wave with a heated area produced by the homogeneous source. It is seen that in the region of the shock wave structure under consideration the fields of εh and Eh are close to constant.
Efficiency in the internal energy ηε and in the kinetic energy ηE are evaluated as follows:
η ε = max t ε m a x ε h ε h ,   η E = max t E m a x E h E h .
Here εh and Eh are the values inside the region between the shock wave and contact discontinuity. The values ηε and ηE have the meaning of the maximum relative difference of the considered types of energy over time between their values for stratified and homogeneous energy sources.
In Figure 15, the dependences of ηε and ηE on the rarefaction parameter α are presented for the considered Mach numbers. It is seen that the relative deviation in the local peaks of ε can achieve up to 29% (for M = 6, α = 0.1) and increases with decreasing α and M (Figure 15a). The relative deviation in the local peaks of E can achieve up to 8.3 times (for M = 6, α = 0.1) and also increases with decreasing α and M (Figure 15b). It should be pointed out that this value depends on the temperature of the layers in the energy source, more so than on the Mach number of the shock wave. Note that the difference in the behavior of ηε for M = 6 and for M = 8–10 is connected with a different location of the maximum value εmax in the considered shock-wave structures.

6. Manifestation of the Richtmyer-Meshkov Instability

It should be noted that the effect of the manifestation of the Richtmyer-Meshkov instability is clearly seen for wider layers in the energy source (Figure 16). Here, the mushroom-like structures consisting of two oppositely rotating vortices were obtained which indicate the formation of this type of instability. In smaller scaled structures these mushroom structures are not clearly seen but it can be stated that the considered shock wave structures contain vortices. Manifestation of the Richtmyer-Meshkov instabilities in problems of supersonic flow control via external energy deposition was obtained and studied in [15].

7. Conclusions

Interaction of a thermally stratified energy source with an initially plane shock wave is simulated for the Mach numbers of the shock wave lying in the hypersonic range, M = 6–12. Shock-wave structures characterized by manifestation of the Richtmyer-Meshkov instabilities were obtained. For the Mach number 6, an almost complete blurring of the shock wave front in the density field is obtained, which is in agreement with the available experimental data. Redistribution of the specific internal energy and the volume density of kinetic energy behind the shock wave front was investigated. The studies were carried out for the rarefaction parameter α lying in the range 0.08–0.4. Analysis of the energy efficiency has also been conducted. It was shown that:
  • The use of stratified energy sources gives the possibility to redistribute the energy behind the shock wave front in such a way that results in the distortion or complete disappearance of the shock wave fronts (in density fields).
  • The mechanism of shock wave blurring is associated with multiple manifestations of the Richtmyer-Meshkov instabilities.
  • The use of stratified energy sources gives the possibility to obtain the local zones of specific internal energy and volume kinetic energy behind the shock wave with the relative difference exceeding up to 29% and 8.3 times, accordingly, these values for the homogeneous energy source with the same total energy; these relative differences increase with decreasing α and M.
  • The dependences of specific internal energy have maxima on t at the beginning stage of the interaction which are larger for smaller α and for larger M, and the internal energy is changing more significantly for these α and M. The dependences of kinetic energy have maxima on t which are larger for larger α and larger M, and the growth of the kinetic energy occurs more quickly for these α and M.
  • The values of the considered types of energy in the obtained shock-wave structures can be controlled at hypersonic speeds via the rarefaction parameter (or temperature) in the layers of the stratified energy source.

Author Contributions

Conceptualization, O.A.A., T.A.L. and K.V.K.; Data curation, O.A.A., T.A.L., K.V.K. and O.V.K.; Formal analysis, K.V.K. and O.V.K.; Investigation, O.A.A. and T.A.L.; Methodology, O.A.A., T.A.L. and K.V.K.; Software, O.A.A. and O.V.K.; Validation, O.A.A. and K.V.K.; Visualization, O.A.A. and T.A.L.; Writing—original draft, O.A.A. and T.A.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data supporting reported results can be provided by O.A. Azarova by the according request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Knight, D.D. Energy Deposition for High-Speed Flow Control; Cambridge University Press: Cambridge, UK, 2019; p. 450. [Google Scholar]
  2. Ahmed, M.Y.M.; Qin, N. Forebody shock control devices for drag and aero-heating reduction: A comprehensive survey with a practical perspective. Prog. Aerosp. Sci. 2020, 112, 100585. [Google Scholar] [CrossRef]
  3. Russel, A.; Zare-Bentash, H.; Kontis, K. Joule heating flow control methods for high-speed flows. J. Electrost. 2016, 80, 34–68. [Google Scholar] [CrossRef] [Green Version]
  4. Leonov, S.B. Review of plasma-based methods for high-speed flow control. AIP Conf. Proc. 2011, 1376, 498–502. [Google Scholar] [CrossRef]
  5. Kolesnichenko, Y.; Brovkin, V.; Azarova, O.; Grudnitsky, V.; Lashkov, V.; Mashek, I. Microwave energy release regimes for drag reduction in supersonic flows. In Proceedings of the 40th AIAA Aerospace Sciences Meeting & Exhibit, Reno, NV, USA, 14–17 January 2002; Paper AIAA-2002-353. pp. 1–13. [Google Scholar] [CrossRef]
  6. Tretyakov, P.K.; Fomin, V.M.; Yakovlev, V.I. New principles of control of aerophysical processes—research development. In Proceedings of the International Conference on the Methods of Aerophysical Research, Novosibirsk, Russia, 29 June–3 July 1996; p. 210. [Google Scholar]
  7. Myrabo, L.N.; Raizer, Y.P. Laser-induced air spike for advanced transatmospheric vehicles. In Proceedings of the 25th Plasmadynamics and Lasers Conference, Colorado Springs, CO, USA, 20–23 June 1994. Paper AIAA-1994-2451. [Google Scholar] [CrossRef]
  8. Bityurin, V.A.; Klimov, A.I.; Leonov, S.B.; Potebnya, V.G. On interaction of longitudinal pulse discharge with bow shock. In Proceedings of the Workshop Perspectives of MHD and Plasma Technologies in Aerospace Applications, IVTAN, Moscow, Russia, 24–25 March 1999; pp. 114–119. [Google Scholar]
  9. Roupassov, D.; Nikipelov, A.; Nudnova, M.; Starikovsky, A. Flow separation control by plasma actuator with nanosecond pulsed-periodic discharge. AIAA J. 2009, 47, 168–185. [Google Scholar] [CrossRef]
  10. Kourtzanidis, K.; Raja, L.; Coumar, S.; Lago, V. Numerical simulation of DC glow discharges for shock wave modification. In Proceedings of the 54th AIAA Aerospace Sciences Meeting, San Diego, CA, USA, 4–8 January 2016; Paper AIAA-2016-2157. pp. 1–11. [Google Scholar] [CrossRef]
  11. Lapushkina, T.A.; Erofeev, A.V. Supersonic flow control via plasma, electric and magnetic impacts. Aerosp. Sci. Technol. 2017, 69, 313–320. [Google Scholar] [CrossRef]
  12. Quan, S.; Bangqin, C.; Yonggui, Y.; Yiwen, L.; Di, J. A study of variation patterns of shock wave control by different plasma aerodynamic actuations. Plasma Sci. Technol. 2010, 12, 708–714. [Google Scholar]
  13. Yan, H.; Liu, F.; Xu, J.; Xue, Y. Study of oblique shock wave control by surface arc discharge plasma. AIAA J. 2018, 56, 532–541. [Google Scholar] [CrossRef]
  14. Jin, J.; Znamenskaya, I.A.; Sysoev, N.N. Two regimes of pulsed volume discharge action upon a shock wave. Tech. Phys. Lett. 2013, 39, 418–442. [Google Scholar] [CrossRef]
  15. Azarova, O.A. Generation of Richtmyer-Meshkov and secondary instabilities during the interaction of an energy release with a cylinder shock layer. Aerosp. Sci. Technol. 2015, 42, 376–383. [Google Scholar] [CrossRef]
  16. Azarova, O.A. Supersonic flow control using combined energy deposition. Aerospace 2015, 2, 118–134. [Google Scholar] [CrossRef] [Green Version]
  17. Apazidis, N.; Sembian, S.; Liverts, M. Blast wave interaction with thermal and density inhomogeneities in air. In Proceedings of the 32nd International Symposium on Shock Waves (ISSW32), Singapore, 14–19 July 2019; pp. 1657–1669. [Google Scholar] [CrossRef]
  18. Bityurin, V.; Klimov, A.; Leonov, S.; Brovkin, V.; Kolesnichenko, Y.; Popov, N.; Van Wie, D. Shock waves structure and velocity at propagation through non-homogeneous plasma. In Proceedings of the 31st Plasmadynamics and Lasers Conference, Denver, CO, USA, 19–22 June 2000. Paper AIAA-2000-2571. [Google Scholar] [CrossRef]
  19. Leonov, S.B.; Carter, C.D.; Hedlund, B.E.; Houpt, A.W.; Ombrello, T.; Firsov, A.A. Control of amplitude and position of reflected shock wave by stripwise plasma. In Proceedings of the 2018 AIAA Aerospace Sciences Meeting, Kissimmee, FL, USA, 8–12 January 2018. Paper AIAA-2018-0683. [Google Scholar] [CrossRef]
  20. Leonov, S.; Carter, C.; Houpt, A.; Ombrello, T. Mitigation of reflected shock wave by streamwise plasma array. In Proceedings of the 7th European Conference for Aeronautics and Space Sciences (EUCASS), Milan, Italy, 3–6 July 2017; pp. 1–11. [Google Scholar] [CrossRef]
  21. Leonov, S.B.; Firsov, A.A.; Houpt, A.W. Suppression of reflected oblique shock wave by multifilamentary plasma. J. Phys. Conf. Ser. 2018, 1112, 012005. [Google Scholar] [CrossRef]
  22. Gan, T.; Wu, Y.; Sun, Z.; Jin, D.; Song, H.; Jia, M. Shock wave boundary layer interaction controlled by surface arc plasma actuators. Phys. Fluids 2018, 30, 055107. [Google Scholar] [CrossRef] [Green Version]
  23. Lapushkina, T.A.; Erofeev, A.V.; Azarova, O.A.; Kravchenko, O.V. Passage of a plane shock wave through the region of a glow gas discharge. Tech. Phys. 2019, 64, 34–41. [Google Scholar] [CrossRef]
  24. Lapushkina, T.A.; Erofeev, A.V.; Azarova, O.A.; Kravchenko, O.V. Interaction of a plane shock wave with an area of ionization instability of discharge plasma in air. Aerosp. Sci. Technol. 2019, 85, 347–358. [Google Scholar] [CrossRef]
  25. Azarova, O.A.; Kravchenko, O.V.; Lapushkina, T.A.; Erofeev, A.V. Density and temperature fluctuations behind a shock wave under the influence of a stratified energy source. Tech. Phys. Lett. 2020, 46, 649–652. [Google Scholar] [CrossRef]
  26. Azarova, O.A.; Kravchenko, O.V. Impact of a thermally stratified energy source on the bow shock wave and aerodynamic characteristics of a body. J. Phys. Conf. Ser. 2021, 1891, 012025. [Google Scholar] [CrossRef]
  27. Azarova, O.A.; Krasnobaev, K.V.; Kravchenko, O.V.; Lapushkina, T.A.; Erofeev, A.V. Redistribution of energy in a viscous heat-conductive medium during the interaction of a shock wave with a temperature layered plasma region. J. Phys. Conf. Ser. 2020, 1698, 012004. [Google Scholar] [CrossRef]
  28. Roache, P.J. Computational Fluid Dynamics; Mir: Moscow, Russia, 1980; p. 612. [Google Scholar]
  29. Azarova, O.A. Complex conservative difference schemes for computing supersonic flows past simple aerodynamic forms. J. Comp. Math. Math. Phys. 2015, 55, 2025–2049. [Google Scholar] [CrossRef]
  30. Raizer, Y.P. Gas Discharge Physics; Springer: Berlin, Germany, 2001; p. 449. [Google Scholar]
Figure 1. Schematic of the problem statement (a); location of the stratified energy source on difference grid (enlarged) (b); location of one layer of on difference grid (enlarged) (c).
Figure 1. Schematic of the problem statement (a); location of the stratified energy source on difference grid (enlarged) (b); location of one layer of on difference grid (enlarged) (c).
Aerospace 08 00326 g001
Figure 2. Comparison of the calculations on four different grids, M = 8, α = 0.4: (a) isolines of ε, t = 0.12 (b) dynamics of εmax.
Figure 2. Comparison of the calculations on four different grids, M = 8, α = 0.4: (a) isolines of ε, t = 0.12 (b) dynamics of εmax.
Aerospace 08 00326 g002
Figure 3. Examples of large-scale (a) and small-scale (b) strata of an ionization-unstable discharge.
Figure 3. Examples of large-scale (a) and small-scale (b) strata of an ionization-unstable discharge.
Aerospace 08 00326 g003
Figure 4. Dynamics of the density field in isochores: M = 8, α = 0.1, t = 0.017, 0.07 and 0.12 (superposed).
Figure 4. Dynamics of the density field in isochores: M = 8, α = 0.1, t = 0.017, 0.07 and 0.12 (superposed).
Aerospace 08 00326 g004
Figure 5. Dynamics of the considered parameters (ρ, p, ε and E): M = 10, α = 0.1.
Figure 5. Dynamics of the considered parameters (ρ, p, ε and E): M = 10, α = 0.1.
Aerospace 08 00326 g005
Figure 6. Fields of the considered parameters (ρ, p, ε and E): M = 6, α = 0.1, t = 0.15.
Figure 6. Fields of the considered parameters (ρ, p, ε and E): M = 6, α = 0.1, t = 0.15.
Aerospace 08 00326 g006
Figure 7. Fields of the considered parameters (ρ, p, ε and E): M = 8, α = 0.1, t = 0.12.
Figure 7. Fields of the considered parameters (ρ, p, ε and E): M = 8, α = 0.1, t = 0.12.
Aerospace 08 00326 g007
Figure 8. Fields of the considered parameters (ρ, p, ε and E): M = 10, α = 0.1, t = 0.1.
Figure 8. Fields of the considered parameters (ρ, p, ε and E): M = 10, α = 0.1, t = 0.1.
Aerospace 08 00326 g008
Figure 9. Fields of the considered parameters (ρ, p, ε and E): M = 12, α = 0.1, t = 0.08.
Figure 9. Fields of the considered parameters (ρ, p, ε and E): M = 12, α = 0.1, t = 0.08.
Aerospace 08 00326 g009
Figure 10. Dynamics of εmax for different α: (a) M = 6; (b) M = 8; (c) M = 10; (d) M = 12.
Figure 10. Dynamics of εmax for different α: (a) M = 6; (b) M = 8; (c) M = 10; (d) M = 12.
Aerospace 08 00326 g010
Figure 11. Dynamics of εmax for different M: (a) α = 0.1; (b) α = 0.2; (c) α = 0.3; (d) α = 0.4.
Figure 11. Dynamics of εmax for different M: (a) α = 0.1; (b) α = 0.2; (c) α = 0.3; (d) α = 0.4.
Aerospace 08 00326 g011
Figure 12. Dynamics of Emax for different α: (a) M = 6; (b) M = 7; (c) M = 8.
Figure 12. Dynamics of Emax for different α: (a) M = 6; (b) M = 7; (c) M = 8.
Aerospace 08 00326 g012
Figure 13. Dynamics of Emax for different M: (a) α = 0.1; (b) α = 0.2; (c) α = 0.3; (d) α = 0.4.
Figure 13. Dynamics of Emax for different M: (a) α = 0.1; (b) α = 0.2; (c) α = 0.3; (d) α = 0.4.
Aerospace 08 00326 g013
Figure 14. Fields of εh (a) and Eh (b), M = 8, α = 0.1, t = 0.12.
Figure 14. Fields of εh (a) and Eh (b), M = 8, α = 0.1, t = 0.12.
Aerospace 08 00326 g014
Figure 15. Efficiency on ε (a) and E (b) for different M and α.
Figure 15. Efficiency on ε (a) and E (b) for different M and α.
Aerospace 08 00326 g015
Figure 16. Manifestation of the Richtmyer-Meshkov instabilities in a case of wider layers, density, M = 3, α = 0.4.
Figure 16. Manifestation of the Richtmyer-Meshkov instabilities in a case of wider layers, density, M = 3, α = 0.4.
Aerospace 08 00326 g016
Table 1. Characteristics of the using difference grids.
Table 1. Characteristics of the using difference grids.
Difference
Grid
A Number of Working Nodes in the Calculation AreaThe Space Steps Values
Grid 11.6 × 106hx = hy = 0.0005
Grid 21.2 × 106hx = hy = 0.000571
Grid 30.9 × 106hx = hy = 0.000667
Grid 40.63 × 106hx = hy = 0.0008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Azarova, O.A.; Lapushkina, T.A.; Krasnobaev, K.V.; Kravchenko, O.V. Redistribution of Energy during Interaction of a Shock Wave with a Temperature Layered Plasma Region at Hypersonic Speeds. Aerospace 2021, 8, 326. https://doi.org/10.3390/aerospace8110326

AMA Style

Azarova OA, Lapushkina TA, Krasnobaev KV, Kravchenko OV. Redistribution of Energy during Interaction of a Shock Wave with a Temperature Layered Plasma Region at Hypersonic Speeds. Aerospace. 2021; 8(11):326. https://doi.org/10.3390/aerospace8110326

Chicago/Turabian Style

Azarova, O. A., T. A. Lapushkina, K. V. Krasnobaev, and O. V. Kravchenko. 2021. "Redistribution of Energy during Interaction of a Shock Wave with a Temperature Layered Plasma Region at Hypersonic Speeds" Aerospace 8, no. 11: 326. https://doi.org/10.3390/aerospace8110326

APA Style

Azarova, O. A., Lapushkina, T. A., Krasnobaev, K. V., & Kravchenko, O. V. (2021). Redistribution of Energy during Interaction of a Shock Wave with a Temperature Layered Plasma Region at Hypersonic Speeds. Aerospace, 8(11), 326. https://doi.org/10.3390/aerospace8110326

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop