Parameters in Sutherland law.
As a promising technology, plasma-assisted combustion (PAC) has attracted many researchers to explore the effect of PAC on improving the combustion in propulsion devices, such as scramjet, detonation engines, internal engines, and so on. In this chapter, we aim to exhibit the influence of quasi-DC discharge plasma on the operating performance of scramjet combustor and find the internal mechanisms, which may contribute to the development of PAC technology in supersonic combustion. For case one, a plasma filament is generated upstream of fuel jet through quasi-DC discharge in a scramjet combustor; for case two, the plasma is formed across the backward facing step of a flame holding cavity to improve the flame stabilization of the cavity in the scramjet combustor. The two cases are investigated in detail through three-dimensional numerical simulation based on the dominant thermal blocking mechanism. Important parameters including temperature distribution, separation zone, water production, stagnation pressure loss, combustion efficiency, cavity drag, mass exchange rate, and cavity oscillating characteristics are obtained and analyzed. It shows that the quasi-DC discharge plasma does benefit for the improvement of the combustion in a scramjet combustor.
- quasi-DC discharge
- plasma-assisted combustion
- supersonic combustion
- transverse fuel jet
- numerical simulation
As one of the qualified candidates for hypersonic propulsion system, scramjet combustor has attracted a large amount of attention all over the world. As is now well known, realizing high efficiency and steady combustion in a supersonic combustor always remains as a critical issue for scramjet engines. Numerous researches have indicated that the additions of cavity in combustor and the transverse fuel injection upstream of the cavity can promote the stabilization of combustion and flame [1, 2, 3]. Whereas, it is difficult yet for the fuel to reach properly mixing within the supersonic flow by the mechanical methods [4, 5]. Moreover, it is a great challenge for matching the transverse fuel jet up with the cavity under off-design condition . Furthermore, certain stagnation pressure loss will appear using the approaches. Taken all account, some new methods are imminently needed for keeping stable and highly efficient combustion in combustor with the least penalties adding to the flowfield.
It is rather promising to adopt a discharge plasma for supersonic flow and/or combustion control in aerospace field [6, 7, 8]. It has been widely considered that plasma-assisted combustion is one of the most promising approaches for enhancing ignition and combustion in the environment of scramjet combustion [4, 6, 9]. As a further promising method, there are three advantages, that is: rapid response, less inertia, and flexibility . Former studies clearly show that the quasi-DC discharge plasma can availably modify supersonic flow in a controllable manner among the discharge plasma mentioned above [8, 11], whereas the DBD plasma is commonly used in low-speed flow environment [12, 13]. If we can combine the quasi-DC discharge plasma with cavity and transverse fuel jet together, some new phenomena will surely appear, which may help in the improvement of combustor performance.
There in three primary mechanisms of the plasma effect on flow and combustion can be summarized: (1) momentum transfer, (2) fast local ohmic heating of the medium, and (3) active particles [6, 7, 14]. Ignoring the magnetic field, the mechanism of quasi-DC discharge plasma affecting a supersonic flow is mainly its fast local heating rather than the electrostatic force (i.e., the momentum transfer mechanism) [10, 15, 16].
This chapter aims at investigating the changes of the fuel jet, cavity, and whole scramjet combustor led by the quasi-DC discharge plasma based on the dominant thermal blocking mechanism. Here, a short cavity downstream of a fuel jet orifice is considered in order to simulate the combustor flowfield more realistic. The plasma is set as a controllable heat source. The
2. Case one—effect of plasma on fuel jet
2.1. Numerical methods
2.1.1. Plasma model and its configuration
By applying a high voltage across the anode and cathode, a bright plasma filament appears above the two electrodes. With the impacting of incoming flow, the plasma occurs not only between the two electrodes but also several centimeters downstream on the wall . In consideration of the oscillation character of its electrical parameters and mean temperature in the discharge path, this type of discharge plasma is called “quasi-DC discharge plasma.” However, it is a type of low temperature arc plasma in fact. The main properties of quasi-DC discharge plasma are described in [10, 11, 12, 15].
As depicted in Figure 1, the domain of quasi-DC discharge plasma filament is simplified as a cuboid region presented in red dashed lines for simulation based on its appearance. In Figure 1, the symmetric plotted line indicates that the cathode and anode are symmetrical for the central line of combustor wall and so does the plasma filament.
Because the quasi-DC discharge plasma releases a large amount of heat concentrated in its discharge path, the discharge path goes into very hot. This high temperature domain (i.e., the discharge path) obstructs the supersonic inflow due to a thermal blocking occurred. Hereby, the quasi-DC discharge plasma behaves as a virtual blockage in the high speed flow of scramjet combustor, which is called as “the dominant thermal blocking mechanism” . On the foundation of dominant thermal blocking mechanism, the individual plasma filament can be simulated as a volumetric heat source [10, 12].
Generally speaking, the plasma input power and the effective heat power that transfers into circumambient gas are quite different. Besides, the percentage of loss commonly varies with the power supply and environment, so it is not a suitable way to use the plasma input power as the heat source value of quasi-DC discharge plasma in simulation. In order to straightforwardly describe the plasma heat strength, the average temperature of the plasma zone is acquired by using the numerical simulation here. This way is feasible when the average temperature of the plasma zone is within a reasonable range [10, 16, 18, 19]. Hence, based on the thermal blocking mechanism, a certain temperature which denotes the actuating strength (i.e., the input power) is specified for the plasma filament domain when the actuator works.
The electrodes are flush mounted and do not have influence on the main flow themselves. The length and the section dimensions of individual plasma filament heat source are 20 mm length and 3 × 3 mm, respectively. The plasma filament locates 40 mm upstream of the fuel orifice center. Besides, it is generated near the wall in the center of the combustor, which means the symmetric plane of plasma filament is within the symmetric
2.1.2. Governing equations, physical models, and numerical schemes
The unsteady Reynolds averaged 3D Navier-Stokes equations (URANS) with the
where the coefficients can be found in . The
|H2||273.11||96.67||8.411 × 10−6|
|O2||273.11||138.9||1.919 × 10−5|
|H2O||416.67||861.11||1.703 × 10−5|
|N2||273.11||106.67||1.663 × 10−5|
Since this study focuses on the qualitative effect of plasma on the fuel jet and the scramjet combustor, the finite rate chemical model with the single step H2/O2 combustion kinetics model is applied. The reaction rate constant is derived from the Arrhenius formula. Hence, the computational time can be saved much and the well combustion flowfield can be obtained, too. The relevant parameters of H2/O2 one-step chemical model are shown in Table 2.
|Equation||Pre-exponential factor||Temperature index||Activity energy [J/kg·mol]|
|9.87E + 8||0||3.1E + 7|
In order to capture the shock waves and other complex fluid structures better, the convective fluxes are evaluated using the advection upstream splitting method (AUSM)  with the second-order upwind approach at same time. The viscous fluxes are evaluated by using the second-order central differential scheme. Because the transport process of multispecies and the reaction both exist in the flowfield, a modified LU-SGS implicit method  is adopted for temporal integration.
2.1.3. Computational zone and boundary conditions
Figure 2 shows the entire computational domain which is the half of scramjet combustor together with a short cavity and a fuel jet orifice. The inlet of combustor is located at
The grids number is largely reduced by setting the symmetry face (
The entire computational domain is divided into 32 subdomains, and all assignments are completed by parallel computation on HP senior workstation, which takes about 320 h to obtain a convergence result. The inflow conditions in computation are as follows:
2.2. Results and discussion
Three representative times are selected during one plasma control cycle for comparison after computation converged. The end of the actuator duty time (
2.2.1. Temperature and wall pressure distribution of the jet flowfield
The temperature distribution on the symmetric
The distribution of wall pressure near the fuel orifice on the symmetric plane of the combustor is shown in Figure 4, which is normalized by the value of the static pressure of inflow at the inlet. Compared with no plasma case, the positions that pressure starts to rise at different times move upstream from 56.6 to 52 mm due to the effect of plasma, which indicates the separation shock wave induced by the fuel jet moves upstream. And the first pressure peak decreases from 2.0 to 1.7, because of the weakening of the fuel jet induced shock on the symmetric plane. The results above are similar to the previous studies on nonreaction flow combustor . More details indicate that the first pressure peak of time A is a bit higher than time B and C, but time B equals time C, which is due to the plasma control cycle too short for flow response and the duty cycle ratio also comparatively large in the flow condition here. Hence, the influence of plasma on the shock will be observed a little latter due to an inertia influence. Nevertheless, the effect of plasma on the shock wave for three typical times is highly similar on the whole.
Because the separation zone upstream of the fuel orifice is primarily controlled by the separation shock wave, the size of this separation zone can be regulated by changing the location of the separation shock wave. On the one side, this zone behaves as an main ignition zone which can provide a high temperature and low flow speed environment in the scramjet combustor. On the other side, it will bring in certain pressure loss to the combustor. Hence, we can make use of the separation zone upstream of the fuel orifice by means of using the plasma with proper control parameters.
2.2.2. Characteristics of fuel mixing and combustion
For the sake of determining the fuel mixing quality, the five cross-sectional planes distributed in equivalent ratio along the flow direction are given, as shown in Figure 5. The fuel jet arises with the actuator working, resulting in fuel decreasing near the wall including the bottom wall of the cavity. Whereas, the distribution shape of fuel in the cross-sectional planes varies from a distribution narrow and long profile into a circular profile in the upper space, which shows that the process of fuel spreading into the main supersonic flow is enhanced resulted from the plasma. Furthermore, the fuel distributions at time A are nearly the same as at time C, which is also correlated to the inertial effect mentioned above.
In Figure 6, the distribution of product water is shown. In order to distinguish the extent of product water between time A and C easily, the iso-surfaces of both are plotted by combining the half parts of them together as shown in Figure 6(b). Similar to the changes in Figure 5, more water is generated in the upper space due to plasma. Contrasted with the case without plasma actuator, the iso-surface of water expands much in its center and shrinks near the combustor wall especially downstream of cavity. These above should be attributed to the change of fuel jet spread, as given in Figure 5. Compared time A with C, the only distinction is a little more water formed at time A, which indicates that the combustion becomes weaker in the free time of a plasma control cycle.
It can be seen that the fuel entrainment into the cavity is decreased and the combustion in cavity becomes weaker due to the plasma, as given in Figures 5 and 6.Whereas, the two figures also show that more fuel penetrates into the supersonic air flow, so the mass exchange between the inside cavity and its outside is enhanced due to the plasma. Hence, the plasma improves the fuel mixing above the cavity prominently, which can also be realized as the calculation results of combustion efficiency in Section 2.2.4.
2.2.3. Combustor stagnation pressure
The stagnation pressure recovery coefficient can indicate the pressure loss in a combustor, so it is an important index. When the stagnation pressure recovery coefficient goes up, the capability of combustor outflow will be enhanced. So, it is defined as the ratio of combustor outlet stagnation pressure to inlet stagnation pressure . In fact, calculating the mass flow weighted mean stagnation pressure can gain the stagnation pressure in a cross section.
The calculated data about the stagnation pressure are given in Table 3. In order to analyze the tendency of stagnation pressure loss coefficient further, the average stagnation pressure at the outlet is replaced by the average stagnation pressure at different positions, as shown in Figure 7. At the actuator working time or even at the free time, the stagnation pressure loss can both be increased due to the plasma, as given in Table 3. In (7), the relative change rate of
|Index of performance||Case|
|Plasma off||Plasma on, time A||Plasma on, time B||Plasma on, time C|
In Figure 7, the difference of stagnation pressure loss between plasma cases is shown where no plasma case enlarges along the flow direction (
2.2.4. Combustion efficiency downstream of the fuel orifice
The combustion efficiency
where is the mass flow rate of water at inlet cross section and is
Figure 8 plots the distribution of combustion efficiency in
In order to define the cost to effectiveness of quasi-DC plasma, it is calculated in value
The main results in this section are as follows: (1) The distribution of relatively high temperature zone moves downstream prominently due to the heat release from the quasi-DC discharge plasma. The separation shock wave induced by the fuel jet is partly weakened and moves upstream due to the plasma, which can regulate the size of recirculation zone upstream of the fuel orifice. (2) The fuel jet moves upward integrally resulted from the plasma heating effect. The fuel spread wider along the spanwise and penetrates into the leading flow deeper, resulting in the cross-section shape of the fuel jet varying from a narrow and long profile to a circular profile. Because of the variation of fuel mixing, in the upper space, more water forms while less appears near the wall compared with the case without plasma. (3) The stagnation pressure loss of combustor increases a little as actuator works, but the combustion efficiency in the combustor rises obviously. These above can be summarized as the comprehensive effects of flow structure changes caused by the plasma, including waves induced and heat transfer. Since it is negligible for the relative change of stagnation pressure recovery coefficient in the actuator working cases, and the ratio of deposited plasma energy to the increased combustion heat release is very little, it can be obtained that the quasi-DC discharge plasma can make more benefits than penalties for the scramjet combustor, when proper adopting control parameters of the plasma actuator.
3. Case two—effect of plasma on cavity
3.1. Plasma model
On account of the direction of the quasi-DC discharge path for the inflow, the discharge modes include two types: longitudinal mode and transversal mode . In case one, the quasi-dc discharge operates under transversal mode, while in case two, a longitudinal mode is adopted for the configuration. As shown in Figure 9, there are five anodes set upstream of the backward wall of cavity, while in the bottom wall of cavity the five cathodes are set. The five pairs of identical electrodes are arranged parallel and symmetrically. Besides, the flow will not be disturbed directly resulted from all the electrodes flush mounted in the wall. A filament plasma forms when applied a high voltage (generally 150–1200 V) between anode and cathode, which acts as strongly oscillating and bright and looks like an inverted “L” crossing the backward wall of cavity. Since the major heat energy of the quasi-DC discharge plasma focuses on the bright filament domain [10, 25], the filament plasma domain can be established in such a simplified shape as given in Figure 9.
Just like the way given in case one, every quasi-DC plasma filament is dealt as a volumetric heat source. To represent the plasma heat strength reasonably, the mean temperature of the plasma domain
3.2. Results and discussion
Three representative times are chosen from one plasma cycle, which consists of the later actuator free duration and the actuator working duration, for comparison after computations converged. In the simulation, the duration of one plasma cycle is
3.2.1. Typical parameters distribution of cavity flowfield
The Mach number distribution of local cavity flowfield on the symmetrical
Owing to the periodic heat release from the plasma filaments and the thermal blocking functions, the plasma filaments behave as five knives cut the cavity shear layer. Hence, the mass transportation process will be disturbed by this “cutting” effect, and the moving direction of original shear layer has to be changed also. And then, the turbulence intensity and vorticity magnitude both are increased, so the fuel and air existing around the original shear layer can exchange through the cavity mouth more easily. Furthermore, the combustion enhancement downstream or over the cavity causes the decrease of local flow velocity, which results in the rise of local static pressure and blocks the incoming flow.
Figure 11 presents the distributions of wall pressure near the rear edge of cavity on
As presented in Figure 11b, the wall pressure of cavity rises distinctly on
Because the distributions of combustion products may be impacted by the change of cavity shear layer and waves, the mass fraction iso-surface of water (
The above phenomena can be attributed to the reasons as follows: (1) Plasma filaments existed upstream of the cavity front edge release a large amount of heat, so the local static temperature increases and then the movement of fuel jet is promoted toward upper space. Because of the larger velocity far from the bottom wall (i.e., flow moves faster in the middle height of the combustor), the original product water will surely move downstream quickly. Therefore, the little water is found around the orifice when plasma actuators work, and the symmetry structures are formed in higher space downstream. It can be verified from the “branch pipes” presented in Figure 12b–d. (2) As stated above, the periodic fluctuation of cavity shear layer is related to the plasma filaments “cutting” effect on it, and then the mass transportation process is blocked which leads to the obvious change of combustion over the cavity. Hence, it affects the distribution of product water. Meanwhile, notice that the “branch pipes” curve structures match with the typical spanwise reverse vortex pairs in shape. Generally speaking, the spanwise reverse vortex pairs can extend the contact area between air and fuel by the methods of entraining air into its fuel core, which makes the mixing between fuel and air better and also the local combustion efficiency arisen. Whereas, the periodic disturbance from the plasma filaments impels the air and fuel transporting in the
3.2.2. Analysis of cavity drag and mass exchange rate
The drag and mass exchange rate of cavity are two important cavity performance parameters . The combustor’s drag is mainly generated by the cavity as it has a nearly constant cross section. The cavity drag includes two types: pressure drag and friction drag. Pressure drag is defined as the difference value between the force on cavity front wall and rear wall. Usually, the friction drag is too small compared with pressure drag so it always can be ignored. So the pressure drag is regarded as the cavity drag here. The drag coefficient of flame holding cavity is defined as:
In Table 4, the calculation results of the cavity drag and its drag coefficient are listed. It can be observed that the drag coefficient is 0.060 at the no plasma case, which is smaller than time B and C but larger than time A. In view of the change of cavity shear layer in Figure 10, it should be due to the moving upward of the shear layer as actuator works. As the shear layer moves upward, it deflects to the main flow and no more impacts on the cavity rear wall so the drag decreases. However, at time B and C, because of the variation of the shear layer fluctuation, stronger impact may appear the rear wall of cavity, and the combustion may be boosted in the cavity rear part, which can both lead to the increase of back pressure in the cavity. Considering the drag coefficient values at the three typical times, the time weighted average cavity drag coefficient is calculated to be 0.0668 which means the cavity drag usually increases by the plasma.
|Front wall||Rear wall||Drag||Drag coefficient|
As we know, the cavity drag in reaction flow is closely related to the combustion heat release zone of a cavity. So only reducing the cavity drag may not bring benefit for the whole combustor. It can be seen that the further detail analysis about the combustion distribution in flowfield should be done to understand the effect of plasma on the cavity performance.
The mass exchange rate of cavity
The mass exchange rate of cavity is given in Table 5. Because the species transportation is achieved through the vortices in the cavity shear layer, it is very significant for enhancing the diffusion ability of cavity shear layer. Owing to the break on the original stable structure of the shear layer resulted from the plasma filaments “cutting” effect, the mass exchange rates at plasma existing cases at times A, B, and C are 9.2, 197.2, and 107.8 times, respectively, than the “no plasma” case. In addition, the time-weighted average mass exchange rate is 8.217 g/s, which is 122.6 times than the “no plasma” case. And the variation of the mass exchange rate magnitude from time A to C is in related to fluctuation degree of cavity shear layer. As a result, the quasi-DC plasma does obviously promote the species transportation between external and internal cavity.
|No plasma||Pulse, A||Pulse, B||Pulse, C|
3.2.3. Analysis of pressure oscillation in cavity
In the research fields related to “sound cavity,” “embedded magazine” and “flame holding cavity,” the dynamic study on pressure oscillation in a cavity is always greatly valued around the world . According to the cavity shear layer distribution, four monitor points are chosen to capture the dynamic pressure feature, as shown in Figure 9. F-1 is located on the front wall of cavity with
The SPLs of F-1, F-2, R-1, and R-2 are 168, 193, 188, and 191 dB, respectively. However, at the no plasma case, their SPL values are 120, 118, 124, and 106 dB, respectively. It indicates that the SPL of all the monitor points go up sharply, and the relationship between point location and magnitude of SPL changes prominently. The points are near the mouth of cavity and the effect of plasma filaments on the cavity edges is strongest, so the SPL magnitude increases here. Moreover, through the analysis on cavity drag in Table 4, the conclusion can be obtained that the pressure disturbance on the monitor points is controlled by combustion, and the shear layer strongly affects the change of local combustion.
In Figure 13, the frequency spectrum characteristics of R-1 and R-2 by FFT are depicted. The several dominant frequencies of pressure oscillation are marked with red circles here, which show that the first dominant frequency is almost the double of the plasma actuation frequency 5 kHz. Whereas, it is barely equal to the plasma actuation frequency under nonreaction condition . Furthermore, nearly every dominant frequency is an integer multiple of the plasma actuation frequency. The results above indicate that the pressure oscillation in the cavity is controlled by both local combustion flowfield and the plasma actuation frequency.
3.2.4. Analysis of stagnation pressure loss and combustion efficiency
The combustor stagnation pressure recovery coefficient of time A, B, C, and the no plasma case is 0.75, 0.81, 0.65, and 0.83, respectively. And the calculated time weighted average stagnation pressure recovery coefficient is 0.73. It suggests that the overall stagnation pressure loss is raised due to the plasma filaments, especially for time A and C. The increase of the combustor stagnation pressure loss is caused by several factors, such as the enlargement of separation zones, the rise of back pressure induced by local combustion and the strengthening of shock waves, which is based on the analysis of cavity shear layer and wall pressure distribution.
Considering the interesting variation of product water distribution as shown in Figure 12, the distribution of combustion efficiency in the flow direction may also vary influenced by the plasma filaments. Hence, the combustion efficiency
In a word, the combustion efficiency increases downstream of the front wall of cavity in most of the time due to the plasma. It becomes lower than the no plasma case at several positions for few moments, but the decrease is comparatively small. Furthermore, during one plasma cycle, the combustion efficiency gradually increases from beginning to end, which indicates that the plasma performs better during actuator free period for the pulsed actuation mode.
The ratio of power consumption of the plasma actuator to the increased combustion heat release of the combustor is calculated in value
The main results are as follows: (1) Because of the “cutting” effect of quasi-DC plasma filaments, the cavity shear layer becomes fluctuating and the local combustion changes which also results in the variation of wall pressure distribution near the cavity rear edge. On the mass fraction iso-surface of product water several “branch pipes” curve structures form periodically, which can be put down to the influence of plasma on the local combustion. (2) The effect of plasma on the cavity drag is relatively perplexing and the drag coefficient value is observed to be unsteady. Nevertheless, the mass exchange rate goes up prominently due to plasma, and the magnitude change from time A to C is in agreement with the fluctuation degree of cavity shear layer. (3) The SPLs of four cavity monitor points rise, and the frequency spectrum of monitor points near the rear edge presents that the first dominant frequency is twice the plasma actuation frequency, under the effect of plasma. Hence, the pressure oscillation in the cavity is controlled by both the plasma and local combustion flowfield. (4) On the one side, the combustion efficiency is usually increased due to plasma, and on the other side, the certain pressure loss increases resulted from the changes of combustion and waves in the flowfield by the plasma. Because the ratio of deposited plasma energy to the increased combustion heat release is extraordinarily small, it is considered that if optimal actuation parameters of the actuator are chosen the quasi-DC plasma can bring more benefits than penalties for the cavity in scramjet combustor.
This work is sponsored by the National Natural Science Fund of China, No. 51876219, No. 91441123, and No. 51777214. The authors would like to acknowledge the help of Dr. W. Feng in numerical simulation.
Conflict of interest
The authors declared that they have no conflicts of interest to this work.
Mathur T, Gruber M, Jackson K, et al. Supersonic combustion experiments with a cavity-based fuel injector. Journal of Propulsion and Power. 2011; 17:1305-1312
Ben YA. Experimental investigation of mixing and ignition of transverse jets in supersonic crossflows [thesis]. Palo Alto: Stanford University; 2000
Ding M. Research on the flame holding technology based on cavity in supersonic combustion [thesis]. Changsha: National University of Defense Technology; 2005
Li G, Li H, Yang L, et al. Review of the plasma ignition and assisted combustion in Russia. Science & Technology Review. 2012; 30:66-72
Kim KM, Back SW, Han CY. Numerical study on supersonic combustion with cavity-based fuel injection. International Journal of Heat and Mass Transfer. 2004; 47:271-286
Starikovskiy A, Aleksandrov N. Plasma-assisted ignition and combustion. Progress in Energy and Combustion Science. 2013; 39:61-110
Starikovskaia SM. Plasma assisted ignition and combustion. Journal of Physics D: Applied Physics. 2006; 39:265-299
Do H, Im S, Cappelli MA, et al. Plasma assisted flame ignition of supersonic flows over a flat wall. Combustion and Flame. 2010; 157:2298-2305
Do H. Cappelli MA, Mungal G: Plasma assisted cavity flame ignition in supersonic flows. Combustion and Flame. 2010; 157:1783-1794
Leonov SB, Firsov AA, Yarantsev DA, et al. Plasma effect on shocks configuration in compression ramp. In: 17th AIAA International Space Planes and Hypersonic Systems and Technologies Conference; 11–14 April 2011. San Francisco: AIAA; 2011. pp. 2011-2362
Leonov SB, Savelkin CV, Yarantsev DA, et al. Experiments on plasma-assisted combustion in M=2 hot test-bed PWT-50h. In: 46th AIAA Aerospace Sciences Meeting and Exhibit; 7–10 January 2008. Reno: AIAA; 2008. pp. 2008-1359
Leonov SB, Yarantsev DA, Gromov VG, et al. Mechanisms of flow control by near-surface electrical discharge generation. In: 43rd AIAA Aerospace Sciences Meeting and Exhibit; 10–13 January 2005. Reno: AIAA; 2005. pp. 2005-2780
Enloe CL, McLaughlin TE, VanDyken RD, et al. Mechanisms and responses of a single dielectric barrier plasma. In: 41st AIAA Aerospace Sciences Meeting and Exhibit; 6–9 January 2003. Reno: AIAA; 2003. pp. 2003-1021
Zhou S, Nie W, Che X. Numerical modeling of quasi-dc plasma-assisted combustion for flame holding cavity. Combustion Science and Technology. 2016; 188:1640-1654
Leonov SB, Yarantsev DA, Gromov VG, et al. Mechanisms of flow near surface electrical discharge generation. Vacuum. 2006; 80:1199-1205
Cheng Y, Nie W, Li G. Numerical study of plasma aerodynamic actuation mechanism. Acta Physica Sinica. 2012; 61:060509
Wang J, Li Y, Cheng B, et al. The mechanism investigation on shock wave controlled by plasma aerodynamic actuation. Acta Physica Sinica. 2009; 58:5513-5519
Zhou S, Nie W, Che X. Numerical investigation of influence of quasi-dc discharge plasma on fuel jet in scramjet combustor. IEEE Transactions on Plasma Science. 2015; 43:896-905
Zhou S. Studies on plasma assisted combustion/flame stabilization in scramjet [thesis]. Beijing: Equipment Academy; 2014
Menter FR. Two-equation eddy-viscosity turbulence models for engineering applications. AIAA Journal. 1995; 33:1418-1425
Esch DD, Siripong A, Pike RW. Thermodynamic properties in polynomial form for carbon, hydrogen, nitrogen and oxygen systems from 300 to15000°K. NASA RFL-TR-70-3(NASA CR-111989), Nov. 1970
Liou MS, Steffen CJ. A new flux splitting scheme. Journal of Computational Physics. 1993; 107:23-39
Yoon S, Jameson A. Lower-upper symmetric-gauss-Seidel method for the Euler and Navier-stokes equations. AIAA Journal. 1988; 26:1025-1026
Rajasekaran A, Babu V. Numerical simulation of three-dimensional reacting flow in a model supersonic combustor. Journal of Propulsion and Power. 2006; 22:820-827
Jiang J, Weng J. Cathode Electronics and Gas Discharge. 1st ed. Beijing: Tsinghua University Press; 1980
Sun M. Studies on flow patterns and flameholding mechanisms of cavity flameholders in supersonic flows [thesis]. Changsha: National University of Defense Technology; 2008
Yang D, Li J, Fan Z, et al. Investigation on aerodynamic noise characteristics of cavity flow at high subsonic speeds. Acta Aerodynamica Sinica. 2010; 28:703-707