Jtam.dvi JOURNAL OF THEORETICAL AND APPLIED MECHANICS 45, 3, pp. 479-488, Warsaw 2007 ABSOLUTE INSTABILITY OF A DOUBLE RING JET – NUMERICAL STUDY Jarosław Bijak Andrzej Bogusławski Czestochowa University of Technology, Institute of Thermal Machinery, Poland e-mail: abogus@imc.pcz.czest.pl; jbijak@imc.pcz.czest.pl The paper is aimed at analyzing the influence of external flow on abso- lutely unstable modes of a ring jet with recirculation zone. The double ring jet is also considered. The investigation is carried out by means of the linearized, inviscid spatio-temporal stability theory. Calculations are based on the shootingmethod for asymmetric azimuthal mode. During the numerical experiment, absolutely unstable modes were iden- tified. No major influence of the external stream on the stability of the flowwas also proved. Key words: absolute instability, double jet, shootingmethod 1. Introduction The linear spatio-temporal theory is applied inorder to showthemostunstable modes of axisymmetric jets and to analyze the influence of an external co- flow on the convective and absolute instabilities. The applied spatio-temporal analysis follows the formulation of Briggs (1964) and Bers (1983) for plasma physics. The complex wave number and the complex frequency are used in the perturbation equation. Using this approach, an absolutely unstable flow can be defined as the one which has a mode with zero group velocity in the upper half of the complex frequency plane. If the zero group velocity lies in the lower half plane, the flow is called convectively unstable. In the field of fluid mechanics, the concept was comprehensively described by Huerre and Monkewitz (1990). The analysis is based mainly on the numerical work of Michalke (1999) concerning a single ring jet with recirculation zone. Michalke has proved, on the basis of the inviscid linearized theory, that the flow is only convectively unstable for the axisymmetric mode and becomes absolutely unstable for the 480 J. Bijak, A. Bogusławski first azimuthal mode with back flow greater than −0.3 of themaximum axial velocity. It seems to be interesting to analyze whether an additional outer ring jet stabilizes or destabilizes the flow and what is the effect of the second recirculation zone between both co-axial jets. The starting point in the research is the configuration considered by Mi- chalke, i.e. a single jet with recirculation zone. The known solution of the flow is used to speed up the saddle point searching for further slightly mo- dified configurations. The iterative procedure is applied, i.e. solution of one configuration serves as an estimation for the next one. The base flow is first modified by increasing the outer co-flow up to 0.9 of the maximum velocity. The changes are sufficiently small to keep the solution close to the previous one. Next, the second recirculation zone is created by decreasing velocity in the region between the jets to the level of −0.3 of themaximum velocity. The final velocity profile is the double annular jet. Contents of the paper is as follows: Section 2 contains a short introduction to absolute/convective instability concepts, in Sections 3 and 4 the scope of the work and numerical procedure are described, results are in Section 5, the paper is summarized in the last Section 6. 2. Spatio-temporal stability analysis The classical linearized stability analysis involves a mean parallel flow and an infinitesimal perturbation superimposed on it. The fluctuations are decompo- sed into elementary instability waves p̂(r)exp[i(αx−ωt)] (2.1) of complex wave number α and complex frequency ω; in the spatial stability theory the frequency is real whilst the wave number is complex, similarly the temporal theory considers complex frequency and real wave number. The x variable is the spatial coordinate in flow direction, r is the radial direction and t is time. The p̂(r) is an amplitude which satisfies the Orr-Sommerfeld equation. In these circumstances the differential equation leads to eigenvalue and the eigenfunction problem, where ω/α is the eigenvalue and p̂(r) the eigenfunction. Next, according to the stability theory, the flow is stable if the Green function of the differential equation is zero along all rays x/t= const lim t→∞ G(x,t) = 0 for all x t = const (2.2) Absolute instability of a double ring jet... 481 and the flow is unstable if the Green function along at least one ray x/t= const is infinite lim t→∞ G(x,t) =∞ for at least one x t = const (2.3) Contrary to the predecessors, the spatio-temporal theory distinguishes two types of instability: convective and absolute. The first one means that initial perturbation is amplified and convected away from its origin. In the second type, the perturbation spreads in all directions and eventually contaminates all spatial positions. It follows that the unstable flow is convectively unstable if lim t→∞ G(x,t) = 0 for x t =0 (2.4) and is absolutely unstable if lim t→∞ G(x,t) =∞ for x t =0 (2.5) The above condition for absolutely unstable flow is necessary but not suffi- cient. The additional constraint was formulated in plasma physics and is often referred to Briggs (1964) or Fainberg et al. (1961) condition. It states that the absolutely unstable flow has a saddle point in the space (α,ω) and the point is formed by coalescence of an upstream and a downstream α branches – it is also called a pinching requirement. To find such a saddle point, a map of α(ω) is constructed where the α branch coalescence can be easily identified. Additionally, as it was shown by Chomaz et al. (1988), a flow is absolutely unstable if the saddle point lies within the region 0< Im(ω)¬ Im(ω(αm)) (2.6) The αm is the wave number of the most amplified wave in time. Exhaustive description of the spatio-temporal theory can be found in the papers of Bers (1983) or Huerre andMonkewitz (1990). 3. Scope of the work The starting point of the numerical experiment is theMichalke velocity profile Michalke (1999) U(r)= 4BF(r)[1−BF(r)] F(r)= 1 1+[exp(ar2)−1]N (3.1) B= 1 2 (1+ √ 1−U0) 482 J. Bijak, A. Bogusławski With the parameters U0 = −0.3, N = 1, the profile represents a ring jet with a back flow on the jet axis. Themaximum jet velocity defines the radius of the jet R, so that Umax = U(R). The nondimensional Michalke profile – normalized by length R and velocity Umax – is shown in Figure 1. The nondimensional form will be used in the calculation. For convenience, the introducednotation is preserved and further in the text all symbols are treated as nondimensional. Fig. 1. Velocity profiles: for the base, Michalke profile 3, flowwith small external co-flow (0.3) + and final profile with co-flow (0.9) and back-flow (−0.3) 2 For the velocity shape mentioned above, Michalke had found one saddle point which satisfied the pinching requirement. The influence of an external co-flow and secondary back flowon the saddle point locuswill be investigated. The modification of the base profile is done by addition of an outer velocity profile of the formMichalke and Hermann (1982) Uouter = U2 2 − U2 2 tanh [ (|r−R2|−R2w) θ 2 ] (3.2) The profile has a single ring jet shape with the velocity magnitude U2. R2 is the distance of the ring from the axis, R2w is thewidth of the ring and θ con- trols the mixing layer thickness. The profile is superimposed on the base one (proposedbyMichalke (1982)) in order to keep theMichalke profileunaffected. Superposition of two velocity profiles created by the above equationwith posi- tive and negative U2 andwith the base profile, leads to the final co-axial flow with two recirculation zones – one on the axis and the second one between two inflow streams. During the numerical experiment the velocity magnitu- de U2 varied from 0.1 to 0.9, in the case of co-flow, and from 0 to −0.3 for the back-flow. The recirculation zone is taken into account only for flows with co-flow greater than 0.3. The remaining parameters of the external flow are summarized in Table 1. The flows with andwithout recirculation zone are shown in Figure 1. Absolute instability of a double ring jet... 483 Table 1.Parameters of the external flow U2 R2 R2w θ co-flow 0.1 to 0.9 4 0.4 15 back flow 0 to−0.3 2.8 0.5 10 4. Numerical procedure The inviscid, linear stability equation (Boguslawski, 2002; Jendoubi and Stry- kowski, 1994) d2p̂(r) dr2 + (1 r − 2 U− c dU dr )dp̂(r) dr − ( α2+ m2 r2 ) p̂(r)= 0 (4.1) where: c= ω/α – the complex phase velocity, has an asymptotic solution in the area of zero velocity gradient dU/dr of the form p̂(r)=C1Im(αr)+C2Km(αr) (4.2) where: C1, C2 – arbitrary constants; Im, Km – modified Bessel functions of order m. The second order differential equation (??) can be transformed to the first order using the substitution χ=−iα p̂ ν̂ (4.3) and dp̂(r) dr =−iα(U− c)ν̂ (4.4) The resulting relation dχ dr =−α2(U − c)+χ [ 1 U −c (α2+ ( m r )2 α2 χ− dU dr ) + 1 r ] (4.5) is solved numerically. The boundary conditions for the above equation is derived from equation (4.2) and the presented substitution. The equations (4.3) and (4.4) lead to the relation for χ χ=−α2(U −c) p̂ dp̂(r) dr (4.6) 484 J. Bijak, A. Bogusławski At the boundaries r=0 and r=∞ the pressure and its derivatives are p̂0 =C1Im p̂∞ =C2Km dp̂ dr0 =C1I ′ m dp̂ dr∞ =C2K ′ m (4.7) Substituting theboundaryvalues to the equation (4.6) andneglecting constant terms, the boundary conditions have the form χ0 =−α 2(U − c) Im I′m χ∞ =−α 2(U −c) Km K′m (4.8) For a value of ω, a guessed value of α is chosen, then Eq. (4.5) is integrated bymeans of the sixth orderRunge-Kutta scheme fromboth sides χ0, the axis, and χ∞, a point sufficiently far from the axis. The two solutions are compared at r=1and if the curves donotmeet, anewvalueof α is calculated according to the Newton-Ralphsonmethod α∗ =α−J−1(χL−χR) (4.9) Value α∗ is only anestimation of thepropervalue for ω, so theproceduremust be repeated until both sides appear with the same left (χL) and right (χR) value at r = 1, including a small error of the order of 10−5. The J in Eq. (4.9) is the difference J = dχL dα − dχR dα (4.10) Finally, ω or c is an eigenvalue and the continuous curve χ(r) is an eigen- function of Equation (4.5). The above procedure is performrd for complex ω in a wide range in order to obtain a α(ω) map in the α complex space. Figu- re 2 presents themap for configuration with the back-flow (−0.3) and co-flow (0.3). The pressure distribution p̂ can be easily obtained from χ using the rela- tion p̂(r)= p̂(r∗)exp ( − r∫ r∗ α2(U − c) χ(r) dr ) (4.11) p̂(r∗) is a boundaryvalue at a certain position r= r∗ and is unknown.Because of that, the pressure is normalized to obtain p̂max = 1. The authors have integrated the integral of Eq. (4.11) from the right side (r∗ = r∞) towards zero, which seems to bemore accurate. Absolute instability of a double ring jet... 485 Fig. 2. Map of α(ω) for configuration with back-flow (−0.3) and co-flow (0.3); a line is the image of α function for constant Im(ω), values denoted near the line 5. Results 5.1. Localization of absolutely unstable modes The experiment revealed numerous saddle points for the considered configura- tions. Nevertheless, only one for each case satisfied the pinching requirement. The remaining points have been eliminated because they were located in a region above Im(ω(αm)) of Eq. (2.6) or ωi was negative. For example, in the double jet with recirculation zone (−0.3) and inflow (0.3), three saddle points were found with: (+) ω=0.5014+0.0510i, (3) ω=1.2742+0.7146i and (2) ω = 1.4514+0.7323i. The map α(ω) of the configuration, Fig. 2, shows that the line of constant ωi = 0.19 has the highest ω imaginery part among the lines crossing the real α axis. This means that perturbation with the ωi = 0.19 is the wave most amplified in time. Accordingly, only the first saddle point is absolutely unstable. As it was noted, no influence of the external co-flow on the stability (in fact on the locus of the single saddle point) was observed. Figure 3 presents the amplification rate ωi as a function of the co-flow velocity. Similarly, no change of the saddle point was observed in the presence of back-flow, neither in α nor in ω complex spaces. In the result, all configura- tions have the same absolutely unstable complex frequency ω=0.501+0.051i and complex wave number α=1.696−0.782i. 5.2. Pressure distribution of the absolutely unstable mode Surprisingly, the external flow with and without recirculation zone did not change the shape of the amplitude pressure disturbance. The pressure ampli- 486 J. Bijak, A. Bogusławski Fig. 3. Amplification rate ωi as a function of co-flow velocity Fig. 4. Pressure distribution along the radial direction. The shape is the same for all configurations tude, shown in Figure 4, is confined to the region of the first ring jet. It is zero on the axis and does not act on the secondary jet. The peak value of the eigenfunction is located in the vicinity of the first inflection point. Fig. 5. Imaginary part of the pressure amplitude for flows without back-flow and with co-flow 0.0 (3), 0.1 (+), 0.2 (2), 0.3 (×), 0.6 (△), 0.9 (⋆) The presence of the additional flow is manifested only in the perturbation phase. For the configuration without the back-flow the phase is increasing with external stream,which is shown inFigure 5.The creation of recirculation zone has a similar effect. The phase shift is in the negative direction with the Absolute instability of a double ring jet... 487 increase of back-flow – Figure 6. However, the phase variation does not have any significant effect on physical behavior of the flow. Fig. 6. Influence of the recirculation zone on the pressure amplitude phase. The co-flow is 0.3 and back-flow 0.0 (3), −0.1 (+), −0.2 (2), −0.3 (×) 6. Summary The linear spatio-temporal stability analysis allowed to identify complex para- meters (α,ω) of the absolutely unstable modes in the presence of outer flow. Amongmany saddle points only one for each of the corresponding configura- tions was proved to be absolutely unstable. No influence of the external flow on the growth rate and wave length of the absolutely unstable mode was ob- served. Analysis of eigenfunctins showed also no significant change. Only a phase variation was noted, but it had no major physical consequences. The investigation was focused on asymmetricmode (m=1)without swirl. Recent measurements of double ring jets Frania (2006), Frania andHirsch (2005) sho- wed that a flow without inlet swirl had a nonzero annular velocity close to the nozzle. It permits to expect additional unstable modes with the presence of a swirl. However, it is only an assumption of authors of the text based on themean velocity data. In future the work will be extended to the analysis of double annular jet with inner and outer swirl. Acknowledgment The authors would like to acknowledge the State Committee for Scientific Rese- arch for financial support. The work was financed from the Statutory Research (BS) funds. References 1. Bers A., 1983, Space-time evolution of plasma instabilities – absolute and co- nvective,Handbook of Plasma Physics, 1, 451-517,Amsterdam, North-Holland 488 J. Bijak, A. Bogusławski 2. BogusławskiA., 2002,Niestabilnosc absolutna i konwekcyjna swobodnej osio- wosymetrycznej strugi plynu o niejednorodnej gestosci, Monografia nr 85, Poli- technika Częstochowska 3. Briggs R.J., 1964, Electron-Stream Interaction With Plasmas, Cambridge, Mass: MIT Press 4. Chomaz J.M., Huerre P., Redekopp L.G., 1988, Bifurcations to local and global modes in spatially-developing flows,Phys. Rev. Lett., 60, 25 5. FainbergYa.,KurilkoB., ShapiroV., 1961, Instabilities in the interaction of charged particle beams with plasmas, Sov. Phys. Tech. Phys., 6, 6. Frania T., 2006,Experimental and Numerical Analysis of Annular Jet Flows, PhD. Thesis, Vrije Universiteit Brussel 7. Frania T., Hirsch Ch., 2005,Measurement of the 3D turbulent flow field of confined double annular jet, Abstract submitted for the 4th AIAA Theoretical Fluid Mechanics Conference, Toronto 8. HuerreP.,MonkewitzP.A., 1990,Local andglobal instabilities in spatially developing flows,Ann. Rev. Fluid Mech., 22, 473-537 9. Jendoubi J., Strykowski P.J., 1994, Absolute and convective instability of axisymetric jets with external flow,Phys. Fluids, 6, 9 10. Michalke A., 1999, Absolute inviscid instability of a ring jet with back flow and swirl,Eur. J. Mech. B/Fluids, 18 11. MichalkeA., HermannG., 1982,On the inviscid instability of a circular jet with external flow, J. Fluid Mech., 114, 343 Niestabilność absolutna podwójnej strugi pierścieniowej – analiza numeryczna Streszczenie W pracy analizuje się wpływ zewnętrznego przepływu na mody niestabilności absolutnej strugi pierścieniowej ze strefą recyrkulacji oraz podwójnej strugi pierście- niowej. Badania przeprowadzono przy użyciu liniowej, czasowo-przestrzennej teorii stabilności dla przepływównielepkich. Do rozwiązania zagadnienia brzegowegowyko- rzystanometodynumerycznebazującenaalgorytmie strzałów.Wpracyuwzględniono tylko asymetryczne obwodowemody niestabilności. Wtrakcie eksperymentunumerycznegozidentyfikowanomodyniestabilności abso- lutnej, nie stwierdzającwpływu zewnętrznego przepływu formującego zarówno strefę recyrkulacji, jak i drugą strugę pierścieniową, na stabilność przepływu. Manuscript redceived February 12, 2007; accepted for print March 21, 2007