Spatial-Temporal Instability of an Inviscid Shear Layer
In this work, we explore the transition of absolute instability and convective instability in a compressible inviscid shear layer, through a linear spatial-temporal instability analysis. From linearized governing equations of the shear layer and the ideal-gas equation of state, the dispersion relation for the pressure perturbation was obtained. The eigenvalue problem for the evolution of two-dimensional perturbation was solved by means of shooting method. The zero group velocity is obtained by a saddle point method. The absolute/convective instability characteristics of the flow are determined by the temporal growth rate at the saddle point. The absolute/convective nature of the flow instability has strong dependence on the values of the temperature ratio, the velocity ratio, the oblique angle, and number. A parametric study indicates that, for a great value of velocity ratio, the inviscid shear layer can transit to absolute instability. The increase of temperature ratio decreases the absolute growth rate when the temperature ratio is large; the effect of temperature ratio is opposite when the temperature ratio is relatively small. The obliquity of the perturbations would cause the increase of the absolute growth rate. The effect of number is different when the oblique angle is great and small. Besides, the absolute instability boundary is found in the velocity ratio, temperature ratio, and number space.
A high mixing rate of the fuel and air is desired in scramjet engine for the propulsion of hypersonic aircraft, because the residence time of the fuel and air in the combustion chamber is very short. In the interest of the projected use of the scramjet engine, it is fundamental and also extremely important to understand the stability characteristics of compressible shear/mixing layers. Many experimental studies [1–3] suggested that the mixing rates of shear layers decrease as the Mach number increases from zero. Hence, a major concern in the development of scramjet is the mixing enhancement techniques. Imparting disturbances on the shear layer, which pulsate at some prescribed frequency, is a choice of mixing enhancement techniques. The prescribed frequency can be obtained through linear stability analysis of the compressible shear layer.
There have been numerous literatures on the topic of linear stability analysis of the compressible shear layer, including the earlier studies conducted by Lessen et al. [4, 5] Drazin and Davey  performed a temporal stability analysis of a compressible mixing layer, which has a hyperbolic tangent velocity profile and uniform temperature throughout the layer. Jackson and Grosch  reported the results of the inviscid spatial stability of a parallel compressible mixing layer. All these studies found multiple stability modes. Zhuang et al.  and Ragab  both found a strong stabilization effect on the flow when increasing Mach number. Ho and Huerre  summarized the studies on linear stability analysis of incompressible shear layer. For the studies on linear stability analysis of compressible shear layer, the reader can refer to the relevant literature in the monograph by Criminale et al. .
The studies above are all confined in the scope of temporal or spatial mode. When we study the stability characteristics of shear layers in spatial-temporal mode, which treat both spatial and temporal eigenvalue complex [12, 13], there are two distinct instabilities for spatial-temporal evolving disturbances: convective and absolute instabilities. The concept of absolute and convective instabilities was first put forward by Briggs  in the study of plasma instability and then introduced to classical hydrodynamic stability . A flow is convectively unstable if the unsteady response to an impulsive perturbation grows along some rays that pass away from the forcing location but decays at the forcing location itself. A flow is absolutely unstable if the impulse response grows at the forcing location .
An absolutely unstable flow is not sensitive to external disturbances and initial conditions; consequently, the “flow management” techniques such as forcing the shear layer at some prescribed frequency could be useless. . Thus, if we want to control the downstream evolution of the flow, it is essential to determine whether the shear/mixing layer is convectively or absolutely unstable. Kulikovskii and Shikina [18–20] studied the asymptotic behavior of localized perturbations on the surface of a shear discontinuity separating two homogeneous steady flows of ideal incompressible fluid in the linear approximation. The effect of surface tension, gravity forces, and viscosity is taken into account. Pavithran and Redekopp  investigated the absolute and convective instabilities for a subsonic mixing layer with the hyperbolic-tangent-like profiles in velocity and temperature fields. The study on transition of absolute and convective instability for compressible shear layer is rare. Caillol  analyzed the transition of absolute and convective instabilities for an inviscid mixing layer. The Mach number in his work attained high supersonic values. However, he did not consider the density stratification within the mixing layer. Large density stratification is possible for binary mixing layers at high pressures, which is usually encountered in power and propulsion system.
The aim of this paper is to track the transition between the absolute instability and convective instability in a compressible shear layer. The density stratification is taken into account. The effects of flow parameters on the spatial-temporal stability of a compressible shear layer are examined by observing whether an increase of the value of parameter tends to increase or decrease the value of absolute growth rate.
2. Theoretical Framework
The base flow is a shear layer between two streams with its far-stream condition denoted by the subscripts and , as shown in Figure 1. The thickness of shear layer increases with increasing distance downstream, denoting that the base flow is nonparallel. Here, we made a locally parallel flow assumption, implying that the present results yield the instability characteristics of individual profiles to leading order when the shear layer thickness grows slowly with downstream distance . The streamwise direction is defined as and the cross stream is defined as . In the present study, denotes the velocity, the density, and the temperature. Taking the momentum thickness and average mean velocity as the characteristic length and velocity scale, respectively, the flow properties can be normalized. Superscript denotes a nondimensional quantity.
The linear stability analysis, the features of which are standard and, hence, are not presented in detail here, is performed for a flow governed by the conservation equations of continuity, momentum, and energy. The thermodynamic condition of the layer is assumed far above the critical point, so the ideal-gas equation of state is used to close the formulation. Linearization of these equations starts by splitting the dependent variables into base and perturbation components. A base flow is assumed that is described by the normalized velocity profile and temperature profile specified as (1)-(2) and by a constant thermodynamic pressure :where the velocity ratio is defined as and the temperature ratio is defined as . The density distributions are obtained by employing the ideal-gas equation of state on the shear layer with the temperature distribution defined in (2).
We perturbed the steady-state solution of the flow with a small-amplitude wave in normal modes: where denotes the perturbation components, denotes the perturbation wave amplitude, , , and are the wave number, phase speed, and frequency, respectively. denotes the angle between the wave propagation direction and the streamwise direction of the mean flow. When the disturbance wave propagates in the axial direction, equals zero, and the resultant system is a two-dimensional problem.
To derive the governing equation for , the governing equations and equation of state are linearized by considering small perturbations to the mean flow whose velocity and temperature profiles are and , respectively. By employing a process of elimination, a single equation can be derived for the radial component of the pressure perturbation [25–28]:where , denotes the speed of sound. For simplicity of notation, the subscript stands for the spatial derivative with respect to the coordinate. The governing equation of the two-dimensional incompressible shear layer can be found by setting and .
Solution of (4) poses an eigenvalue problem due to the presence of the wave parameters, , and frequency, . The boundary conditions are obtained by considering that (i) at far stream the base flow derivatives are zero and that (ii) the solution must be bounded. The solution is obtained by treating the frequency as the eigenvalue for an externally specified wavenumber and by using a shooting iterative procedure. ODE45 in MATLAB is employed to perform the numerical integration of the ordinary differential equation, and the Newton-Raphson shooting method is utilized to perform the iterative shooting to find the eigenvalue. Convergence is found to be strongly dependent on a good starting guess for . The detailed process is as follows: to obtain the relations, we begin from an intermediate (e.g., ) and make an initial guess for . Then, the Newton-Raphson shooting method is used, assuming the initial guess for is good enough to obtain the exact . Then we decrease or increase by and use calculated with the old as the initial guess for the new shooting iteration. If the shooting iteration converges again, we will move to the smaller or larger . With this computing process, the curve extends from the middle to both low and high regimes. If, at some point, the shooting method diverges or can not obtain a physical solution, or the resultant curve is not smooth enough, we conclude that the initial guess for at is incorrect and we provide a new guess and start over again. It is by this process that the contours presented in this article are obtained.
This study focuses on the transition between absolute instability (AI) and convective instability (CI) for a compressible inviscid shear layer, which are dominated by modes with zero group velocity. These correspond to complex pairs, , which satisfy (4) and also the further property that .
The saddle point method is used to determine the solution with zero group velocity. In general, there are two distinct spatial branches of solutions of (4) on the complex plane. The two branches will approach each other as the growth rate is decreased from large positive values, and a saddle point on the complex plane will occur at point . It is necessary to identify the physical and nonphysical saddle points. The physical saddle point must satisfy the Briggs-Bers collision criterion [14, 29]. According to the Briggs-Bers collision criterion, the only relevant complex pair is the physical pair, if the two branches come, respectively, from the downstream-propagating branch and upstream-propagating branch on the complex plane. Then if the flow is said to be convectively unstable, while if , the flow is said to be absolutely unstable. Normally, is referred to as the absolute wave number and the absolute growth rate.
3. Results and Discussion
In this section, the effects of dimensionless parameters on the spatial-temporal instability of an inviscid shear layer will be checked by taking a parametric study. The complex frequency is solved as function of the wave number and the parameter set (, , M, and ).
The typical contour plots of the imaginary part of the complex frequency in the complex wave number plane for an inviscid compressible shear layer are shown in Figure 2. Apparently, there is a saddle point in the plot, which is marked with “.” The location of the saddle point is and the value of real part of complex frequency at the saddle point is . The saddle point satisfies the Briggs pinching criterion. In Figure 2, the absolute growth rates are negative, indicating that the inviscid shear layer is convectively unstable with the given dimensionless parameters.
The effect of velocity ratio Λ on the absolute instability of an inviscid shear layer is illustrated in Figure 3. It appears in Figure 3 that the absolute growth rate increases with the increases of velocity ratio . In particular, when the velocity ratio is larger than 1.47, the absolute growth rate becomes positive. The transition velocity ratio 1.47 was close to the value obtained by Huerre and Monkewitz , which was 1.315. Huerre and Monkewitz  stated that the mixing layer is absolutely unstable when velocity ratio is larger than 1.315. That is because Huerre and Monkewitz considered shear layer in incompressible fluid, whereas the present paper deals with perfect gas. The velocity ratio represents the shear force imparted on the shear layer; therefore, it can be said that, for a two-dimensional inviscid shear layer, it is absolutely unstable when the shear force imparted on it is relatively large. Manoharan and Hemchandra  proposed a mechanism through which absolute instability can cause combustion instability in combustor of power generation and propulsion system. Hence, avoiding the occurrence of absolutely unstable flows in combustors, at least at desired operating conditions, represents a possible passive control strategy to mitigate the occurrence of combustion instability. According to Figure 3, to depress absolute instability, a possible smallest velocity ratio is recommended. However, small velocity ratio, denoting weak shear force imparted on the shear layer, shows no advantage on mixing. Thus the velocity ratio should be designed specially in a combustor.
Figure 4 gives the effects of the temperature ratio, , on the absolute growth rate of a two-dimensional inviscid shear layer, by increasing and keeping other parameters constant for every plot. Figure 4 shows that, for different cases, as the temperature ratio increases, the absolute growth rate decreases. Hence the strategy for suppressing absolute instability is to increase temperature ratio.
However, the temperature effects do not always tend to depress the absolute instability. When the temperature ratio is relatively small (), the absolute growth rates shows minor increase with increase of for every case. That is to say when is relatively large, the effect of on absolute growth rate is more remarkable; when is small, has minor promotion effect on absolute instability of the shear layer.
Now, we check the effect of obliquity of the disturbance wave on the absolute instability by varying the oblique angle, . Figure 5 shows the variation of absolute growth rate with oblique angle . It shows in Figure 5 that when the flow is two-dimensional (), the absolute growth rate is the smallest. When the oblique angle increases, the absolute growth rate increases. The increasing rate becomes larger when the oblique angle is greater. However, the absolute growth rate is always negative when varies from to . Hence, it can be concluded that the increase of oblique angle can increase the absolute growth rate, and the effect is more remarkable when the oblique angle is greater, or we can say that the three-dimensional disturbance waves are more dominant in inviscid shear layer than two-dimensional wave. But the increasing of would never change the flow into absolute instability.
The parameter is naturally supposed to be parameter influencing the absolute instability behavior of flow. Here we check the effect of and oblique angle comprehensively. Figure 6 gives the variation of absolute growth rate with and oblique angle. It shows that, for a constant M, the increase of oblique angle would increase the absolute growth rate, which can also be seen from Figure 5. For a given oblique angle, the effect of number on absolute growth rate differs when oblique angle is larger or relatively small. When the oblique angle is small, the increase of would decrease the absolute growth rate; while the oblique angle is greater (close to ), the effect of increasing is minor. That means the effect of number is opposite when the disturbances propagate upstream and downstream.
Figure 7 shows the critical value of Λ at which the flow transits into absolute instability at different and . It can be seen from Figure 7 that the increase of would cause the critical value of Λ decrease; that is, the increase of would expand the absolute instability region in -- space. However, the variation of has no influence on the absolute instability boundary, although the increase of can decrease the absolute instability growth rate.
In this work, we explore the transition of absolute instability and convective instability in a compressible inviscid shear layer, through a linear spatial-temporal instability analysis. From linearized governing equations of the shear layer and the ideal-gas equation of state, the dispersion relation for the pressure perturbation was obtained. The dispersion relation was solved using a shooting iterative procedure. The numerical results show that the increasing velocity ratio can make the inviscid shear layer transition to absolute instability. The increase of temperature ratio decreases the absolute growth rate when the temperature ratio is large; the effect of temperature ratio is opposite when the temperature ratio is relatively small. The obliquity of the perturbations would cause the increase of absolute growth rate. The effect of number is different when the oblique angle is large or small: when the oblique angle is relatively smaller, the increase of number would decrease the absolute growth rate; while the oblique angle is relatively large (close to ), the variation of number has no influence on the absolute growth rate. The absolute instability boundary has been given in the -- space. The increase of would expand the absolute instability region in the -- space, but the variation of has no influence on the absolute instability boundary.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
The financial support of National Nature Science Foundation of China (Support nos. 11672025 and 11302013) is gratefully acknowledged. The authors gratefully acknowledge financial support from China Scholarship Council (Support no. 201506025044).
M. Zhuang, T. Kubota, and P. E. Dimotakis, “On the instability of inviscid, compressible free shear layers,” AIAA Paper 1988-3538, AIAA, Washington, DC, USA, 1988.View at: Google Scholar
S. A. Ragab, “Instabilities in the wake mixing-layer region of a splitter plate separating two supersonic streams,” AIAA Paper 1988-3677, AIAA, Washington, DC, USA, 1988.View at: Google Scholar
R. J. Briggs, Electron Stream Interaction with Plasmas, MIT Press, Cambridge, Mass, USA, 1964.
M. P. Juniper, A. Hanifi, and V. Theofilis, “Modal stability theory,” in Applied Mechanics Reviews, vol. 66, 2014.View at: Google Scholar
F. Q. Hu, T. L. Jackson, D. G. Lasseigne, and C. E. Grosch, “Absolute-convective instabilities and their associated wave packets in a compressible reacting mixing layer,” Physics of Fluids A, vol. 5, no. 4, pp. 901–915, 1992.View at: Google Scholar
A. G. Kulikovskii and I. S. Shikina, “Influence of viscosity on stability of tangential discontinuity in incompressible liquid,” Moscow University Mechanics Bulletin, vol. 52, pp. 6–9, 1997.View at: Google Scholar
K. Manoharan and S. Hemchandra, “Absolute/convective instability transition in a backward facing step combustor: fundamental mechanism and influence of density gradient,” Journal of Engineering for Gas Turbines and Power, vol. 137, Article ID 021501, 2015.View at: Google Scholar
T. Liu, Shear-coaxial injection and mixing of cryogenic fluids under supercritical conditions [Ph.D. thesis], The Pennsylvania State University, 2007.
A. Bers, “Theory of absolute and convective instabilities,” in Proceedings of the International Congress on Waves and Instabilities in Plasma, Innsbruck, Austria, 1973.View at: Google Scholar