Shock Reduction through Opposing Jets—Aerodynamic Performance and Flight Stability Perspectives

: In this research paper, investigations of counter ﬂow (opposing) jet on the aerodynamic performance, and ﬂight stability characteristics of an airfoil with blunt leading-edge in supersonic regime are performed. Unsteady Reynolds-Averaged Navier-Stokes ( URANS ) based solver is used to model the ﬂow ﬁeld. The effect of angle of attack ( α ), free-stream Mach number ( M ∞ ), and pressure ratio ( PR ) on aerodynamic performance of airfoil with and without jet are compared. The results indicate that the opposing jet reduces drag from 30% to 70%, improves the maximum lift-to-drag ratio from 2.5 to 4.0, and increases shock stand-off distance from 15% to 35% depending on ﬂow conditions. The effect of opposing jet on longitudinal ﬂight stability characteristics, studied for the ﬁrst time, indicate improvement in dynamic stability coefﬁcients ( C m q + C m ˙ α ) at low angles of attack. It is concluded that the opposing jet can help mitigate ﬂight disturbances in supersonic regime.


Introduction
The typical objectives of commercial aeronautics research during the twentieth century were centered around achieving long range, high altitude and speed. Efforts in the form of Concorde and Tupolev-144 were made for supersonic commercial travels but due to technical problems, these endeavors were abandoned. Therefore, the speed of commercial aviation was curtailed to subsonic and transonic flow regimes because of multi-faceted concerns [1]. Since 2003, there is no supersonic passenger aircraft in service [2]. Supersonic flight faces multiple problems such as undesirable acoustic signature, high drag, increased aerodynamic heating and flight instabilities. The wave drag is a specific problem of a supersonic flight, which contributes more than half of the total drag and increases rapidly with increasing free stream Mach number. This increase in drag value builds up high temperature and pressure that may damage the leading edge of the blunt body and it must be prevented [3].
This motivates researchers around the globe to explore technologies that may be used to minimize the drag for blunt bodies in supersonic flow regimes. In 1950s, the opposing jet gained immense popularity in the research community due to its durability and relatively economic implementation as compared to other methods [4,5]. The jet pushes the shock-wave away from the stagnation point and develops an interface between the shock-wave and blunt body. The interface is composed of a complex flow structure that include formation of expansion and reflection waves within the jet, barrel shocks to cover these waves, and an enlarged vortical structure named re-circulation region [6]. The re-circulation region ahead of the blunt body plays a prominent role in the reduction of drag. The effective shock reduction through opposing jets is also dependent on flow conditions and geometric shapes [7,8].
The parameters of opposing flow jets have been investigated experimentally, numerically, as well as analytically. Experimentally, researchers ascertained an encapsulation of body due to formation of bow and re-compression shocks, called as shock caps. The ejection of secondary flow field (opposing jet) opposes the free stream flow and forms a contact surface that separate these two shocks. The encapsulation due to shocks provide a protecting blanket to reduce the body temperature [9,10]. Opposing jet also reduces the drag. The ejection of jet distributes the pressure over the surface and generate different low pressure zones (re-circulation or dead air regions and stagnation point or leading edge pressure after ejection of jet) as shown in Figure 1a,b [11,12].  Jet to body diameter ratio, shape of jet orifice and pressure ratio also effect the performance of opposing jet. The exposed front side of blunt body helps in the strength and positioning of the re-circulation regions [13]. The change in pressure ratios between free stream flow and opposing jet develop long penetration mode (LPM) and short penetration mode (SPM) [14].
The under-expanded jet (LPM) attained at low pressure ratio penetrates the bow shock and displaces it to the maximum position for the formation of conical shape [4,15]. An unstable diamond shaped pattern can be observed due to the formation of weak expansion and compression waves (reflected waves) which makes the LPM mode unstable in Figure 1b [16]. The size of re-circulation regions in LPM mode are thinner and remain at the leading edge of the blunt body to cover the jet boundary. The displaced shock-wave reattaches to the surface of blunt body after passing over the re-circulation region and shows changes in pressure [17,18]. These changes are named turning shock-waves. A similar kind of trend is observed in SPM mode, which manifests at higher pressure ratios. The size and positioning of re-circulation regions are different than the LPM mode. Re-circulation regions are enlarged in width and placed at some distance from the vicinity of opposing jet. The reattachment of shockwave also triggers similar turning shock-waves due to change in pressure regions (upstream and downstream pressure differences due to shock-wave) as shown in Figure 1a [19].
Analytical models provide the theoretical basis of opposing jet for the mitigation of shock-wave adverse effects. Finley [20] developed an analytical model of jet flow-field, which primarily predicted the opposing force coefficient (C F ) of jet. Later on, researchers called this force coefficient static pressure of jet. The pressure ratio of jet and body was controlled by using the concept of jet exit Mach number (Mj). Numerically it was found that the interface at front of Mach disc is blunt for large PR values (due to decrease in Mach of jet) as shown in Figure 1a. The simulated result suggested that the Mach disk of unstable mode is not clearly visible (due to absence of barrel shocks). The increase in free stream Mach (due to increase in pressure ratio) convert an unstable LPM into stable SPM.
The development of opposing jet technology led researchers to study various other aspects such as effect of different jet ejection materials (water, gases and plasma). By injection of cold vaporized water in airflow, drag reduced by up to 50% at the cost of very high pressure ratio [21]. Therefore, instead of using vaporized water, different gases (He, H 2 , CO 2 and N 2 ) were used to investigate the reduction of drag [10,22]. Helium (He) is found to have better performance for drag reduction as it is lighter and inert gas (as compared to CO 2 and N 2 ) [23]. Although H 2 is lightest, its reactive property can cause increase in heat transfer to the body [24]. A similar methodology was adopted by Fomin [25], in which plasma materials were used to reduce the shock. To get much better results for drag reduction [26], combination of passive (spike) and active (opposing jet) were tested [27][28][29].
The active and passive combination of flow control technologies on blunt bodies attain some additional factors (separated re-circulation regions) which produce unsteadiness in the flow and effects the stability of body. The ejection of jet reduces the unsteadiness in the flow by displacing the separation region and vortices on the surface of blunt body and stabilize the flow. Previously, researchers used wind tunnel tests [30][31][32][33][34][35][36], empirical and computational fluid dynamics (CFD) tools for the prediction of stability behavior. Empirical methods like DATCOM and ESDU were developed for the correct prediction of stability derivatives [37]. These methods were fast enough for the wind tunnel tests. However, these methodologies do not model the active flow control features. Therefore, a numerical approach like CFD can be used for the capturing the contributions of flow control techniques to stability characteristics [38,39].
Studies conducted to date indicate that extensive investigations have been performed on the use of the opposing jet technique for the reduction of drag and heat transfer from the surface of body. However, nearly all work includes blunt bodies (like sphere, cylinder, truncated cone or aeroshells) for the ejection of jet. A limited work exists in which elliptical bodies (airfoil shaped) were used for the reduction of drag using flow control techniques, especially opposing jet [5]. For a supersonic configuration, such as Unmanned Combat Air Vehicle, a long spike or fuselage ahead of the wings will reduce its stealth features. Therefore, usage of opposing jets should be tested on the wing geometries as well. In this research, the effect of opposing jet on overall aerodynamic performance and the flight stability characteristics of a two-dimensional wing is explored. Specifically, the interaction of opposing jet with the supersonic free-stream flow is investigated using a numerical approach. The purpose is to reduce the drag from the surface of airfoil (NACA 0012) and to improve the understanding of some fundamental phenomena like effects of pressure ratios, free stream flow Mach number and their combined effects with increasing angle of attack. These phenomena affect the ejection of opposing jet and stability parameters of airfoil. The stability derivatives (statically and dynamically) are calculated using a reduced frequency approach. It is pertinent to mention that the effects of acoustic and heating are not considered in this study.
The problem formulation and validation cases are briefly presented in Section 2. The detailed results are discussed in Section 3, which include the influence of free stream Mach number on the ejection of opposing jet in accordance with its required PR. The next section entails the effect of angle of attack on the opposing jet and its stability during perturbation. At the end, concluding remarks are presented in Section 5.

Problem Formulation
The computational setup is created by using two dimensional Unsteady Reynolds Averaged Navier Stokes (URANS) and turbulence model equations. A density based implicit solver is used as it works effectively for supersonic flows. The validation process for all geometries (a, b and c) (as shown in Figure 2a) followed in this research, is a three-step procedure. Like Itsariyapinyo [40], a similar grid formation strategy is adopted for all three validation steps as shown in Figure 2b. In the first step, experimental findings of Graham [41] are used to validate the aerodynamic characteristics and flow-field around a bi-convex airfoil with blunt leading edge. The second step include assessment of computational setup for adequate modeling of the opposing jet. The experimental results of Eugene [5] are reproduced numerically for with and without jet scenarios around elliptical body in supersonic regime. Finally, in the third step, experimental results of Landon and Davis [42] are compared with the computational setup to capture the dynamic stability characteristics.

Computational Setup
In this research, N ACA 0012 airfoil is selected as a reference aerodynamic surface. The airfoil has a blunt leading edge and moderate thickness to chord ratio. Moreover, it is one of the most generic airfoil shapes studied in aerodynamic literature. The airfoil has 1 m chord with 0.004 m jet orifice at leading edge ( Figure 2a).
A fully structured two-dimensional mesh was generated. The boundary layer resolution was made such that Y + is less than one on the surface. Pressure inlet (jet), pressure outlet (static pressure at outlet) and the pressure far field (to model free-stream compressible flow up to an optimum distance) boundary conditions were selected for numerical analysis. Isentropic relations, as in Equation (1), were used to calculate the static pressure and temperature (P ∞ , T ∞ ) at sea level stagnation conditions (P o ∞ , T o ∞ ) [12].
where γ is the specific heat ratio. The discretized domain used for all simulations is shown in Figure 2c,d.
The compressible flow around the airfoil is simulated by using ANSYS R FLUENT. Two dimensional Unsteady Reynolds Averaged Navier Stokes (URANS) and turbulence model equations were solved using a density based implicit solver, as it works effectively for supersonic flows. The ideal gas law was used to model the flow properties in the far field. However, the viscosity of air was defined by Sutherland's viscosity law to capture the complex flow features like shock wave and re-circulation regions accurately. The k-ω SST turbulence model was used to calculate the complex flow field which is expressed as however, ρ is density, ω is specific turbulent dissipation rate and k is the turbulence kinetic energy. Where A 1 in Equation (3) represents the blending function F 1 which ensures the smooth transition between the k-ω and k-. This transition can be expressed in Equation (4) as The value of blending function (F 1 ) becomes one near wall but at far position its value become zero. The building of turbulence in stagnation regions is prevented by the implementation of production limiter (P k ). So, the generation of turbulent viscosity can be expressed as where F 2 is a second blending function with strain rate, S. This model was proved to be superior as it was a combined version of k-ω (for the inner boundary layer) and k-(for the outer region and outside of boundary layer). The models switch for prediction of far field flow characteristic. For the boundary layer (SST inner), constants are defined as However, the constant values for outside of boundary layer (standard k-) are represented as For further details about blending functions (F 1 , F 2 ), production limiter (P k ) and SST turbulence model constants (for inside and outside of boundary layer), readers are referred to Reference [43].

Validation of Computational Setup
The test cases used for validation of computational setup ahead use same boundary conditions, parameters for grid generation, flow solver and turbulence model.

Validation of Supersonic Airfoil with Blunt Leading-Edge
The experimental data generated by Graham [41] on bi-convex airfoil (R3015) is validated numerically. The airfoil of 0.254 m chord length was cut at 15% from the leading edge to make it blunt Figure 3. The flow conditions included Reynolds number of 5 × 10 6 , angle of attack of 2 • , and Mach number of 1.4. The blunt-shaped airfoil is modeled and simulated using the computational setup mentioned above. The comparison of pressure coefficient between our numerical and experimental data show decent agreement in Figure 3. A small pressure drop at point C is observed due to over-expansion. Similarly, a negligible rise in the pressure near trailing edge (point B) is still within the bounds.

Modeling of Opposing Jet Flows
The analysis of with and without opposing jet at supersonic flow regime was carried out using elliptical shaped body with the chord length of 0.1524 m (6 in) [5]. The free stream flow is at Mach number 1.62 with Reynolds number of 7.66 × 10 6 . The opposing jet is ejected from the nozzle of diameter 0.000762 m (0.03 in) at a Mach number of one and its thrust coefficient (C T ) was −0.0156.
The thrust coefficient for without jet case is zero (C T = 0). The comparison of the pressure coefficient in Figure 4a shows that our computational setup generated close agreement with the experimental data. In a similar manner, the changes in the coefficient of pressure due to the ejection of opposing jet is validated with the experimental data as illustrated from Figure 4b.

Problem Setup for NACA 0012
After validating the computational setup, NACA 0012 airfoil, considered in this study, is analyzed in transonic flow regime at free stream Mach number of 0.86 and Reynolds number of (9 × 10 6 ). The numerical results of C p , compared with experimental data of Harris [44] shows excellent agreement in Figure 5a. Parameters like grid sizing, operating conditions and geometric accuracy might have contributed in the small variations of results. The obtained numerical result shows a similar trend as in experimental values. To attain the best agreement between the numerical and experimental data, different grid sizes were used for the analysis. The coefficient of drag ceases to vary after 0.6 mil (Figure 5b) grid. Hence, mesh with 0.6 mil was used for further analysis of opposing jet ejection and to calculate stability parameters.

Validation of Stability Cases
In this research, the longitudinal pitch damping derivatives of NACA-0012 airfoil (with and without opposing jet) were calculated by selecting suitable parameters such as amplitude of oscillation (α A ), mean angles of attack (α ∞ ) and reduced frequency (k). The pivot point of rotation for NACA 0012 is at half of the chord (0.5c) of an airfoil for calculation of lift, drag and moment coefficients which calculate the pitch damping derivatives. The force and moment coefficients against time at different angles of attack were computed. The sinusoidal pitching motion of airfoil was achieved by applying forced pitch oscillation. A user defined function (UDF) generated the oscillatory motion. A circular domain was formed around the airfoil as an interface in Figure 2c and oscillates against its angular frequencies. The pitching oscillation and sinusoidal motion at different angles of attack, can be related as: where α o is the mean angle of attack while α A is oscillatory amplitude and ω is the angular frequency.
The reduced frequency (k c ) approach builds a relationship between the angular frequency (ω) and free stream velocity (V). The forced pitching oscillatory motion (sinusoidal motion) of airfoil continuously change the free stream flow angle that brings the dynamic change in the moment coefficient of airfoil.
Using the Taylor expansion [45], the moment coefficient can be expanded as: neglecting higher order derivatives, Equation (8) will reduce to or where C m α and (C mα + C m q ) express the static and dynamic stability parameters for airfoil respectively. The values of normal force coefficient (C n ) over the range of angle of attack show a close agreement with experimental data [42] (having Mach number 0.57 and Reynolds number of 4 × 10 6 for static case) as shown in Figure 6a. However, Figure 6b shows superimposed numerical and experimental results for dynamic analysis (perturbed airfoil) by the comparison of moment coefficient (C m ) at range of angles of attack.
(a) C n along static angles of attack (b) Pitch damping coefficient (C m ) long angles of attack

Flow Physics of Opposing Jet
In supersonic flow regime, the most significant feature of blunt shape airfoil is formation of bow shock-wave. This work emphasizes on the alleviation of drag by implementing the opposing jet on the leading edge of airfoil. Variation in the angle of attack, pressure ratio, and Mach number is analyzed in detail such that effective shock reduction can be quantified.

Influence of Mach Number
Free stream Mach is an important variable which effects the jet flow structure. A significant amount of drag reduction (up to 70%) is demonstrated in Figure 7d. The percentage of reduction in drag coefficient (C d ) improves with the increase in free stream Mach number (M ∞ = 1.4 to 3.0). Pressure ratio, PR of 1.8 and Mach of jet, M j , equal to one were used to displace the shock wave from the nose of airfoil which primarily helps in the reduction of drag coefficient. However, the effectiveness of opposing jet at these PR and M j decreases beyond M ∞ = 2.5 as shown in Figure 7d.
The Mach contours and the variation of pressure at stagnation point (0 in x-axis) and surface of airfoil (with and without jet) with the increase in free stream Mach number is illustrated in Figure 8a-d. The comparison of Mach contours in Figure 8a-d clearly depicts that the smaller PR ratio is required at higher free stream Mach number to get the visible Mach disc. Fujita [46] also suggested that the increase in the free stream Mach number strengthens the bow shock and, as a result, the total pressure behind the shock decreases. The depletion in total pressure behind the shock minimize the pressure in the recirculating region so that the static pressure ratio of the jet enhances (at M ∞ = 2 & M ∞ = 3), where the pressure behind the shock is considerably low.
M ∞ at 1.4 and 1.8, shows significantly higher drag. As a result, unstable flow structure was formed thereby exhibiting inferior performance of jet. At M ∞ = 2, the jet attains the SPM mode with visible Mach disc and appropriate re-circulation region as well as reduced drag coefficient. But at M ∞ = 3, the opposing jet seems to be more under-expanded as the Mach disc did not maintain its shape.
In Figure 8a-d, the decrease in pressure values due the effects of counter flow jet at the stagnation point of airfoil are plotted. The stagnation point is displaced by means of opposing jet to an optimal distance from the leading edge of airfoil and creates a shock stand-off distance to reduce the drag. However, the graphical comparison of with and without jet raise a point which presumes that there was a sudden rise in pressure after the displacement of shock-wave.
The re-circulation region collapses near (x/c = 0.02) and produce turning shocks along the surface of airfoil [20]. The abrupt change in pressure values take place after the turning shocks as shown in Figure 7b and coalesce along the surface of airfoil. Similarly, these turning shocks are of compressive nature before maximum thickness position and expansive afterwards for no jet case as shown in Figure 7a. However the re-circulation region destroyed and evolved a new free stream boundary layer for different values of M ∞ . This reattachment of free stream flow was another reason for the generation of turning shocks and increase in pressure. The hatched region in Figure 7c explains the difference of pressure between with and without jet curves; which were disturbed due to the generation of turning shocks.
In the similar scenario by increasing the M ∞ , pressure over surface of airfoil continuously decreases as depicted from Figure 7d. The main reason is movement of barrel shocks along the free stream and spreading of re-circulation region (dead-air region) over airfoil surface. So, it significantly decreases the force factors such as drag.

Pressure Ratio Effects
The flow field in the vicinity of leading edge is significantly governed by pressure ratio. The impact of opposing jet on flow field and drag was analyzed for different values of pressure ratios (0. .0) at M ∞ = 1.4 and M j = 1. The variation in pressure ratio (PR) generate two distinct flow modes, LPM and SPM. LPM, unsteady and undesirable, exists at low PR values. SPM, steady and desirable, exist at higher PR, which is more beneficial for drag and heat flux reduction [47]. In between the oscillating/unstable LPM and stable SPM, a transition regime is also observed.
The re-circulation region generated as a result of opposing jet is illustrated in Figure 9c. Interestingly, the pressure on the upper surface of airfoil remains same (while, moving towards the trailing edge) for with and without jet at zero degree angle of attack (as displayed in Figure 9a,b. Points A, B, and C in Figure 9d explain the peak pressure values for flow countered with and without jet. Point B is the stagnation point for no jet case that gets displaced through opposing jet to point A. Point A is the new stagnation point lying on the interface curve between disturbed shock wave and Mach disc of jet. Point C represent the pressure value of the jet at the leading edge of an airfoil that is significantly higher than the free stream total pressure. The modes are shown in Figure 10a-l. For PR = 0.8, the total pressure in the jet reservoir is less than the stagnation pressure. Therefore, the jet ejection from the orifice do not take place. However, when PR is greater than one, the shock standoff distance amplifies with PR. Moreover, the flow pattern of re-circulation zone in the vicinity of opposing jet is also significantly changed. The distinct LPM mode can be seen at PR of 1.9 and 2.4. The oscillatory behavior of pressure can be observed at these values. In LPM, the drag value is recorded enough due to the fluctuations that formulate a series of expansion and compression waves [19] at PR = 2.4 (as in Figure 10f). The pressure of re-circulation regions were significantly higher than the stagnation point in front of Mach disc. From PR of 2.8 and 3.0 onwards, a strong SPM mode is captured. The influence of opposing jet diminishes beyond 10% of the chord length (x/c = 0.1) as depicted in Figure 10a-l.
Similarly, some minor fluctuations were observed for SPM at PR = 3 where pressure of re-circulation region was less than the stagnation point at different PR values. As the free-stream static pressure at the interface curve was greater than the counter flow jet, the stagnation point of interface stopped Mach disc of jet from further protrusion. With the increase in PR value, the Mach disc also got enlarged thereby bending the barrel shocks. The bending of shockwaves along the free-stream flow also moved the re-circulation regions over the airfoil surface. That re-circulation region creates an imaginary body that reduces the skin friction coefficient from the leading edge of airfoil. Figure 11a empirically illustrate the flow field difference between no jet, unsteady LPM and steady SPM modes. The flow manifests the same trend as that of Rebecca [48] who analyzed the flow modes on blunt body at free stream Mach 3.48. There is a transition region that coexists between LPM and SPM modes. At PR = 4, the Mach disc does not maintain its shape as depicted in Figure 11b. The increment in PR increases the mass flow rate [8].
The increase in mass flow rate on small jet to body diameter ratio build more under-expanded jet and less jet stream passes through the barrel shock [49]. Same can be observed in Figure 11b where the increase in mass flow rate by increasing pressure ratio (PR = 4), due to the formation of relatively compressed barrel shock offers more room to the jet's out flow as compared to the lower pressure ratio (PR = 3). Location of Mach Disk for different pressure ratios can be controlled by varying jet diameter.
(a) Regional distribution of Flow modes.
(b) Movment of barrel shockwave due to change in PR value. Figure 11. Under-expansion of opposing jet at different pressure ratio and its regional distribution.

Angle of Attack
The comparison of instantaneous pressure distributions for no jet and jet cases at various angles of attack are illustrated in Figure 12a SPM mode remains steady and stable with the increase in angle of attack (α ∞ ). By increasing α ∞ from 0 • to 6 • , the flow structure of SPM sustained as depicted in Figure 12a-d, though a strong flow asymmetry was introduced at higher values. This asymmetry persists and apparently becomes stronger with increasing angles of attack (8 • and 10 • ) until the counter flowing jet cannot persist the short penetration state. The counter flow jet was detected in one side of an airfoil and shock wave reattached with the surface. The comparison of pressure graphs of with and without jet (Figure 12a-d) reinforces the above explanation.

Aerodynamic Performance and Stability Analysis
The aerodynamic efficiency (lift-to-drag ratio) increases with the increase in angle of attack and maximizes till 8 • . The advantage of reduction in coefficient of drag reduces with the increase in angle of attack as shown in Figure 13b. The overall performance of airfoil increases by using counter flow jet as depicted from Figure 13a. The re-circulation regions tilt in the windward direction and fill the wake areas of airfoil at higher angles of attack. It is evident from the Figure 13a that after the stall of an airfoil, the performance of jet is no longer as efficient as it is for lower angles of attack. Figure 14a show the comparison of lift coefficients of with and without cases at range of angles of attack (α = 0 • -8 • ). It is evident from Figure 14a that by implementing jet the lift increases, drag reduces and static stability remains unchanged at positive angles of attack (Figure 14b).
The opposing jet certainly influence stability parameters as well. These stability parameters were calculated at constant free stream Mach 1.4 for both with and without jet cases with a pivot point at mid (0.5c) of chord of airfoil. The equation of linear line is assumed analogously for the calculation of static stability derivatives.
however, m is the slope of linear line. The adopted analogy for the hysteresis loops of pitching moment coefficient (C m ) is changing as whereas, i and j represents the moment (m) and angle of attack (α) of perturbed airfoil. So, the equation 13 becomes and C m α is the pitching moment slope coefficient of hysteresis loops. To evaluate dynamic stability derivatives, fully converged static cases were initialized with forced oscillations. The sensitivity of dynamic stability derivatives is checked against reduced frequency, mean angle of attack and amplitude of forced oscillations. The reduced frequency contributes a lot in the longitudinal pitch damping derivatives. The dynamic derivatives of with and without opposing jet are shown in Figure 15a,b for several values of reduced frequencies with constant amplitude α A = 4 • . Figure 15a shows that by increasing reduced frequency the body moves towards the unstable behavior (The value of (C m q + C mα ) start increasing). Although it shows a linear behavior along increasing angles of attack. But in the case of ejection of jet these linear behaviors drastically change for low frequencies (k = 0.05, 0.1). However, for higher frequencies in opposing jet shows linear behavior for increasing angles of attack as it was for no jet case. Basically, there are three types of bodies (primary, secondary and tertiary) that affect the moment coefficient and its damping forces as shown in Figure 15a,b. The most important body is the primary body that is an airfoil. As the angles of attack starts increasing from 0 • the secondary body (generation of adverse pressure gradient) starts developing and effects on the overall moment coefficient C m . This adverse pressure gradient (secondary body) brings stall for the body at high angles of attack.
But in case of jet ejection, a jet with recirculation regions represents a third body (tertiary body) generated which reduces drag from the leading edge of an airfoil but also effects the stability derivatives. At low angles of attack, no secondary body is present for low reduced frequency (k = 0.05), but the tertiary body (jet with recirculation regions) still produce effects on the stability derivatives as shown in Figure 15b.
The comparison of without and with jet cases at k = 0.05 shows that the without jet case depicts linear trend whereas the nonlinearity produces for increasing angles of attack by the ejection of jet. As clear from above Figure 12a-d, by increasing angles the jet disturbed from the leading edge causes instability in flow which leads the disturbance in the dynamic stability parameters as illustrated in Figure 15b. By increasing the reduced frequency, a similar linear trend starts building for higher angles of attack as that of without ejection of jet. The reason behind this linearity is reduction of opposing jet (tertiary body) effectiveness on the movement of the airfoil (primary body) but the secondary body (adverse pressure gradient) still exists as they were for without opposing jet.
The other parameter of interest is the influence of the angles of attack on the stability derivatives. The analysis was carried out at several angles of attack ( α o = 0 • , 2 • , 4 • , 6 • & 8 • ) with fixed amplitude of 4 • and reduced frequency of 0.2. By the application of jet, the slope of stability derivatives is damped as compared to without jet cases (Figure 16a,b). With increasing angles of attack, a negative slope starts developing for pitching moment of simple perturbed airfoil (as shown in Figure 16a), while the ejection of opposing jet from the leading edge leads the negative slope towards the positivity as depicts in Figure 16b. At higher angle (8 • ) the effect of opposing jet starts reducing (or vanished) which leads the positive slope towards the negative value.
The number of calculations was made to investigate the influence of mean angles of attack on airfoil stability. The analysis were carried out at several mean angles of attack ( α o = 0 • , 2 • , 4 • , 6 • & 8 • ) with fixed reduced frequency k = 0.2 at mid chord movement. In this analysis the amplitude of motion was also varied and selected as, The comparison in Figure 17a,b shows that by implementing jet on the stagnation point of an airfoil the dynamic stability enhances. It was noticed that for jet case the pitch damping derivatives decreases as the angle of attack increases. But beyond the mean angle 6 • the dynamic pitch coefficient value starts increasing again because at higher mean angle the recirculation regions of jet attached to the free stream flow and destroyed the jet mode, so the value of C m q + C mα starts increasing and reached close to the without jet case.
The reduction at the lower angles are due to the abrupt decrease in the pressure on the upper surface at windward side as shown in Figure 18a-c. This reduction brings the restoring moment in the dynamic system of airfoil, which ultimately decreases the overall dynamic derivatives of stability at lower angles of attack. The positive coefficient of pressure distribution of leading edge is reduced with the ejection of jet. The vicinity on the upper surface contributes a destabilizing to damping coefficient that suppress with introduction of jet.

Conclusions
The current numerical research aimed to understand the phenomenon of opposing jet for drag reduction on a blunt leading edge airfoil (however, discussion related to temperature, density and acoustic were not considered for the analysis of opposing jet). After the numerical validation of baseline experimental cases, a jet is introduced on conventional airfoil. The computational investigation shows that by introducing opposing jet the overall drag of an airfoil reduces up to 70%. However, the performance of supersonic jet is extremely affected by varying operating parameters as follows: • The pressure ratio (PR) is the key parameter that affects the counter flow jet. The LPM exists in a smaller range of pressure ratios (2.4); it produces a long shock standoff distance. As the pressure ratio is beyond the threshold value, the counter flow jet transmits its flow variation from LPM to SPM (PR = 3).

•
Increased in angles of attack influenced the opposing jet efficiency. Beyond α ∞ = 6, the jet destroyed and the drag at the leading edge starts increasing again.
• Free stream Mach is the major influencing parameter. The SPM mode is achieved at M = 2. With further increase in Mach (M = 3) the jet becomes highly under expanded and the shape of the Mach disc changes into a convex form.

•
The static and dynamic stability of airfoil in supersonic flow enhances by the implementation of jet. At higher free stream flow angles (α ∞ = 6) the jet does not work properly as it tilts in the free stream flow direction. So, the stability of an airfoil reduced. Funding: There were no funding sources for this research.