Abstract
Numerical studies have been performed to visualize vortical flow structures emerged from jet cross-flow interactions. A single square jet issuing perpendicularly into a cross-flow was simulated first, followed by two additional scenarios, that is, inclined square jet at angles of 30° and 60° and round and elliptic jets at an angle of 90°, respectively. The simulation considers a jet to cross-flow velocity ratio of 2.5 and a Reynolds number of 225, based on the free-stream flow quantities and the jet exit width in case of square jet or minor axis length in case of elliptic jet. For the single square jet, the vortical flow structures simulated are in good qualitative agreement with the findings by other researchers. Further analysis reveals that the jet penetrates deeper into the cross-flow field for the normal jet, and the decrease of the jet inclination angle weakens the cross-flow entrainment in the near-wake region. For both noncircular and circular jet hole shapes, the flow field in the vicinity of the jet exit has been dominated by large-scale dynamic flow structures and it was found that the elliptic jet hole geometry has maximum “lifted-off” effect among three hole configurations studied. This finding is also in good qualitative agreement with existing experimental observations.
1. Introduction
The problem of injecting fluid through pipe/duct geometry into a mainstream cross-flow domain presents in many industrial and engineering applications, for example, turbine blade film cooling, fuel injection in IC engine, thrust, and noise control of S/VTOL aircraft, fuel-air mixing in gas turbine combustors, and pollutant dispersion from chimney stacks. Due to these wide range applications, the jet in cross-flow (JICF) configuration has been the subject for numerous experimental and theoretical studies. Since the observation of coherent structures in earlier 1970s, various experiments have been devoted to JICF research (see, e.g., Margason [1]) and it was believed, based on those studies, that the dynamic process of flow motions in the JICF is closely relevant to and, at some extent, predominated by the formation and evolution of large vortical structures, which are originated from the jet shear layer, evolved after the jet exit, and decayed further downstream. Fric and Roshko [2] identified four different vortical structures presented in a JICF system for the first time and pointed out that the counterrotating vortex pair (CRVP) is the dominant vortical structure in this kind of complex flow systems.
Despite the fact that the inclined jet and jet with various hole exit shape have been widely used for the blade film cooling by gas turbine and aeroengine industries, there are relatively few studies focusing on the understanding of the underlying flow physics. In case of the inclined jet, earlier experimental work by Lee et al. [3] revealed that significant flow structure changes happened at an inclined jet angle of 35°. Yuen and Martinez-Botas [4] carried out a systematic experimental study at three jet inclination angles of 30°, 60°, and 90° and identified the maximum cooling effectiveness, which was reached at a velocity ratio of 0.33 and jet angle of 30°. Similarly, the jet exit hole shape also plays an important role in determining the adiabatic cooling effectiveness and the cooling boundaries that are often used as a measure on how blade surface could be protected [5, 6]. Due to a “lifted-off” effect of counterrotating vortex pair (CRVP) in the near field of a hole exit, the “cooled” jet flow tends to be “separated or detached” from the blade surface and this will cause the significant reduction of the cooling effectiveness. However, some flow physics related to the underlying cooling mechanism is still not clear yet.
The present paper focuses on the investigation of the effect of different jet inclination angles and hole geometries on the dynamic evolution process of vortical structures associated with jet and cross-flow interactions using direct numerical simulation (DNS) approach. A parallel in-house DNS solver [7, 8] will be used to simulate the vortical flow development around square, round, and elliptic jets issuing normally into a cross-flow.
2. Numerical Method
In the present numerical paper, the compressible time-dependent three dimensional Navier-Stokes and energy equations are solved. By using reference values at free stream, the nondimensional form of these equations can be written as where represents shear stress tensor, Re the Reynolds number, the specific heats, Pr the Prandtl number, M the Mach number, and the dynamic viscosity. All other variables follow the standard notations as seen in Yao et al. [9]. The code was developed for compressible flow. In this application, we adopt low Mach number (~0.2) to reduce the compressibility effects.
The code has been parallelized using the MPI library and validated extensively for numerous configurations including laminar boundary layer and channel flows, turbulent boundary layer, turbulent channel flows, transonic flow over a bump geometry with shock/boundary layer, and jet in cross-flow. Recent code reengineering work has made it more robust and efficient [10].
3. Problem Definitions
The problem considered for simulations consists of a baseline configuration of a single square jet issuing perpendicularly into a main cross-flow field with two further scenarios, that is, (1) square jets issuing at an inclination angle of 60° and 30° against the horizontal walls and (2) round and elliptic jet hole shape issuing at a normal 90° angle against incoming cross-flow, respectively.
The computational box (baseline) has a cross-flow domain of in the streamwise, the wall normal and the lateral directions, respectively, and a jet domain of , where D is the jet width. In case of circular-type jets, the jet hole exit area is kept the same as the square jet case and then the radius of round jet and the major and minor axis length of elliptic jet (assuming a ratio of 2) can be derived accordingly. In case of square jet, it is located in 4D and 5D from the cross-flow inlet plane in the streamwise direction and 2.5D and 3.5D in the spanwise direction, respectively (see Figure 1). Following a previous numerical study of Sau et al. [11], a computational grid of 241 points in the streamwise (with 11 points in the jet domain), 81 points in the wall normal, and 61 points in the spanwise direction (with 11 points in the jet domain) is adopted. Both domain influence and grid refinement studies have been carried out [12]. For example, simulations were performed on domain sizes of and (using similar grid resolutions and distributions) and computational meshes of and (both applied to the baseline domain), respectively. The examination of streamwise and vertical velocity profiles at representative locations of the cross-flow domain has shown that very small differences were generated by the domain size increases and the grid refinements. It was therefore decided to retain both the baseline domain and the grids as described above for the flow analysis simulations.
The cross-flow velocity profile is initialized using a similarity solution of laminar boundary layer at Reynolds number of 225. A standard circular pipe flow Poiseuille-type velocity profile, that is, , is used at the inlet of a circular jet orifice and for 2D non-circular duct and elliptic pipe geometries, an interpolation procedure is applied. The characteristics boundary conditions are used at both outlet plane and upper surface and periodic conditions for the side planes. Simulations are performed for a jet to cross-flow velocity ratio and the Reynolds number , based on the free-stream quantities of the cross-flow and the jet width (D) and the Mach number , based on the free-stream quantities.
4. Results and Discussions
Similar to Sau et al. [11], results presented here were taken at an instantaneous dimensionless time unit of , as the main objective of the study is to visualize the vortex flow structures. Discussions will focus on understanding the underlying flow physics in JICF at different scenarios.
Figure 2(a) depicts contours of spanwise vorticity () at a midspan plane () and it reveals the presence of the shear layer vortex that emerges at the upstream side of the initial portion of the jet. The roll-up process of the shear layer upstream of the jet has been attributed to the Kelvin-Helmholtz (K-H) instabilities [2, 11, 13]. Figure 2(b) reveals the evolution of the three-dimensional streamlines participating in the formation of a K-H vortex system in the near-wake region that emerges at the upstream (wind side) of the initial portion of the jet flow. This type of roll-up process of the shear layer upstream of the jet has been attributed to the so-called Kelvin-Helmholtz instabilities as discussed in previous studies (e.g., Fric and Roshko [2], Smith and Mungal [14]). Due to the merging process whilst the vortex roller moving downstream, the scale of K-H vortices increases with the distance, similar to that in a free jet flow. Overall present observations are in good agreement with previous experimental and numerical predictions [12].
(a)
(b)
The influence of the streamline injection angle variations on the flow structures in jet cross-flow interaction has been assessed. Figure 3 presents the isosurfaces of spanwise vorticity () for 90°, 60°, and 30° injection angles (left). The CRVP is clearly present in all cases of injection angle with the shear layer (Kelvin-Helmholtz) vortices emerge at the upstream side of the initial portion of the jets. However, for the 30° injection case, as the jet centerline is much closer to the wall, the evidence of Kelvin-Helmholtz rollers is less clear. As injection angle increases, the jet penetration is higher and the spreading area is wider. The vorticity contours of are shown at location D from the wall in Figure 3 (right). Compared with the normal injection case at the same blowing ratio, the starting position of the CRVP is shifted downstream and its strength is reduced at both 60° and 30° injections. For the 30° injection, since the jet is much closer to the wall, the interaction between the jet and the boundary layer is more complex in the wake region. The evidence of wake vortices is not as clear as in the instantaneous plots in Figure 3(c) (right).
(a) Normal square jet at 90° angle
(b) Inclined square jet at 60° angle
(c) Inclined square jet at 30° angle
The study continues with the effect of jet exit hole shape variation with complicated nature of the JICF illustrated. Figure 4 depicts 3D isosurfaces of the spanwise vorticity and its 2D contours on the symmetric plane (, , ), that exhibit clear rollup of vortical structures, in addition to the tilting of the evolving vortex rings. The delayed rollup of vortical structures observed in the round and elliptic jets, in comparison to that of square jet, is probably due to the removal of the sharp corners, otherwise it will facilitate the earlier formation of vortex structures. Comparing to elliptic jet, the round jet has shown that the vortex rings form earlier downstream. This again could be attributed to the effect of jet orifice shape change. Note that the nature and the dynamics of these coherent structures play a significant role in the mixing and heat transfer process between the jet flow and the cross-flow. Figure 5 (left) presents the spanwise () vorticity contours in the (, D, ) plane for the three simulation cases. A different flow structure can be noted between the elliptic jet and that of circular and square jets. This is mainly due to the difference in vortex pairing process (see, e.g., Maidi and Yao [15]) which occurs further downstream from the jet orifice. It is important to recall that the presence of vortex pairing and the formation of kidney- and antikidney-shaped vortices depend on the aspect ratio of the elliptic jet [16]. Figure 5 (right) shows the lower-deck vortex pair (the large contrarotating vortex pair) appearing near the surface of the plate in all cases. The origin of these vortices is the lateral spillage and rollup of the jet sidewall boundary layer. On the top of the lower-deck vortex pair, an opposing upper-deck vortex pair (the small vortex pair) can be identified. The upper-deck vortex pair is the hole leading-edge boundary layer. Both lower- and upper-vortex pairs for a round hole were previously reported by other authors (e.g., Margason [1], Tyagi and Acharya [17]). These vortices have a sense of rotation that promotes both the jet “lifted-off” and “entrainment” of the cross-flow towards the wall, respectively. Thus they can seriously degrade the effectiveness of the protecting film layer near the surface. Based on Figure 5 (right), the jet “lifted-off” is more important in the elliptic jet hole case in which the upper-deck vortex pair, promoting the entrainment of the cross-flow toward the wall surface, is relatively weak compared to those identified in the other two hole geometry cases. The results reported here are in good qualitative agreement with the experimental studies of Haven and Kurosaka [16], New et al. [18], and DNS of Sau et al. [11].
(a) Round jet
(b) Elliptic jet
(a) Square jet
(b) Round jet
(c) Elliptic jet
5. Conclusions
The interactions between jet flow and cross-flow are examined by using a direct numerical simulation approach. Simulation considers two scenarios, including single square jet at three issuing angles of 30°, 60°, and 90° and single circular jet with round and elliptic hole exit geometries.
For single normal jet, the resulting vortex structures were in good agreement with both experiments and other numerical simulations in the literature. The computations successfully reproduced the important vortex patterns presented in this type of flow, including the counterrotating vortex pair (CRVP), horseshoe vortex, shear layer vortex, and wake vortex.
For single inclined jet at 30° and 60° angles, simulation results reveal the significant flow structure changes, with a local recirculation flow observed upstream of the jet orifice for a small injection angle of 30°. The vorticity contours also exhibit complicated flow patterns at two inclined jet cases, indicating strong 3D effects. The jet inclination decreases the lateral deflection of cross-flow.
For circular normal jet, it shows that the hole geometry has considerable influences on the near field characteristics of vortex structure. The proximity of these counterrotating vortices relative to each other is found to affect both the liftoff of the jet and the entrainment of cross-flow fluid toward the wall surface. The obtained results show that the jet liftoff as well as the fluid entrainment into the near wall region mechanisms depends on the hole geometry, in agreement with previous experimental observations in the literature. The maximum jet liftoff is found in the case of elliptic hole geometry.
Acknowledgment
The support of the UK Engineering and Physical Science Research Council (EPSRC) through the research Grant (EP/C014979/01) is gratefully acknowledged.