Transient Response Characteristics of Separated Flow and Heat Transfer in Enlarged Rectangular Channel ()
1. Introduction
Recently, prediction of separated flow and heat transfer in an enlarged channel is very important in relation to many types of fluid machinery and heat exchangers, and there have been so many experimental and numerical investigations [1-14]. It has been clarified in these studies that the flow in a two-dimensional symmetrically enlarged channel becomes asymmetric over the critical Reynolds number. The cause of becoming a steady asymmetric flow can be explained by the Coanda effect. The pressure decreases owing to an increase of the velocity near the lower or upper wall in the channel. The flow is biased to the lower pressure side when the pressure difference is once generated across the channel, and the asymmetry of flow is stably maintained. It seems that the instability of the separated shear layer disturbs the velocity resulting in the asymmetric flow. However, the details of these phenomena have not been clarified yet.
In engineering applications, the flow channel almost usually provides the side walls, which affect greatly the flow structure resulting in the three-dimensional flow [3, 13]. Accordingly, such the flow and heat transfer may be quite different from that obtained with the two-dimensional simulation. However, there has been little information about the three-dimensional structure of separated flow and heat transfer in an enlarged channel. Cherdron et al. [3] and Schreck and Schäfer [13] have shown the stabilizing effect of the side walls experimentally and numerically.
Several investigations regarding influences of the side walls upon the heat transfer have been performed. Yanaoka and Ota [15,16] ant Ota et al. [17] have investigated numerically a three-dimensional flow and heat transfer over a blunt flat plate in a channel and in a rectangular channel with a sudden expansion. It has been clarified that the flow in the separation bubble becomes extremely three-dimensional with an increase of the Reynolds number because of the side walls and the heat transfer is greatly enhanced near the side walls. Iwai et al. [18] have studied numerically a laminar flow and heat transfer over a backward-facing step in a rectangular duct and reported that the heat transfer is enhanced near the side walls. However, they have treated a steady asymmetric flow and there have been only a few works on an unsteady flow in the channel. The authors [19] have investigated effects of inlet flow imposed pulsating fluctuation [20] upon the separated flow and heat transfer in the enlarged rectangular channel. It is clarified that the pulsating fluctuation of St = 0.05 and 0.10 strongly affects on the flow in the channel and heat transfer on wall surfaces enhanced, especially on the wall covered with the long separation bubble.
However, transient flow and heat transfer with the inlet flow imposed fluctuation are not clarified in detail. The main objective of the present research is to clarify the transient response characteristics of the separated flow and heat transfer with only one cycle of a pulsating fluctuation in the enlarged rectangular channel.
2. Fundamental Equations and Numerical Procedures
In the present study, three-dimensional (3D) flow of an incompressible viscous fluid with constant properties is treated. Figure 1 illustrates a 3D enlarged channel along with the coordinates and several important geometrical parameters. The fully developed laminar flow enters into the inlet, but only one cycle of a pulsating fluctuation [20] is imposed initially. The flow separates at the step corner and reattaches to the downstream walls. The wall surfaces in the downstream of the step are heated. The present study simulates the forced convection and does not consider the influence of buoyancy. Fundamental equations are the continuity, momentum and energy ones.
(1)
(2)
(3)
In these equations, physical quantities are nondimen-
Figure 1. Flow configuration and coordinate system.
sionalized by, , , and, respectively. The superscript * indicates nondimensional quantities, and it is omitted in the above equations. The finite difference forms of these equations are obtained by the second-order Adams-Bashforth method for the time differentials. The fifth-order upwind difference is employed for the convection terms in order to make clear the details of flow structure, and the fourth-order central difference for the diffusion terms, and the second-order central difference for other terms. Resulting finite equations are solved using the SMAC method [21].
The numerical calculations are performed on the cases of the expansion ratio, ER = 2, the aspect ratio, AR = 8, the Reynolds number based on the mean flow velocity at the inlet and the step height, Re = 200, and the Prandtl number, Pr = 0.71, respectively. The computational region extends from 10 H upstream of the step to 100 H downstream of the step. The boundary conditions are as follows. As for the velocity, the no-slip condition is assumed on the walls and the zero gradient is used at the outlet. At the inlet, the fully developed duct flow velocity profile is presumed, though only one cycle of a pulsating fluctuation [20] is imposed initially. As for the temperature, it is assumed to be uniform at the inlet, and the lower and upper wall surfaces downstream of the step are heated under the condition of uniform temperature, and other walls are adiabatic. The second derivative is assumed to be zero at the outlet. Initial flow and temperature fields are steady state in the channel and the flow deflects downward, and after that the proper symmetric velocity profile with only one cycle of a pulsating fluctuation is imposed at the inlet.
The computational grids are generated using the Cartesian coordinates with non-uniform spacings and the grids of 251 × 155 × 61 are used.
3. Results and Discussion
In the preliminary calculations, in order to verify the present code, numerical streamwise velocity profiles on the center plane at several streamwise cross sections are compared with the previous measured ones by Fearn et al. [7] for ER = 3, AR = 24, Re = 186.7, and Durst et al. [9] for ER = 2, AR = 16, Re = 203.3 in Figure 2. In the case of Re = 186.7, the present results agree very well with the measured ones in the whole flow field. In both cases, the instability of shear layers separated from the step corner increases and the flow becomes asymmetric because of the pressure difference across the channel. The steady asymmetric flow is stably maintained. In the case of Re = 203.3, the present numerical results reasonably agree well with the measured ones, though small differences can be detected around the shear layers. It is necessary to confirm the grid dependency in the case of Re = 203.3 using a finer grid in order to clarify these small discrepancies. However, these comparisons show the good reliability of the present numerical scheme.
Figures 3 and 4 illustrate side and bottom views of instantaneous helicity isosurface. Two inclined vortices exist behind the step near the side walls symmetrically at t = 0. These vortices originated by a jet like flow in steady state [18] and almost keep the shape and location. These vortices are elongated toward spanwise center of the channel on the lower wall from t = 5. Corotating vortices appear downstream of the vortices at t = 17. At t = 20, the initial vortices shed two vortices. Counter rotating vortices are induced in front of the downstream vortices near the side walls. The initial shed vortices slightly incline to the upper wall and upper share layer sheds two vortices at t = 25. Downstream vortices successively induce counter rotating vortices from t = 25 to t = 80. The vortices gradually inclined from the walls from t = 47. The cloud of the vortices gradually decreases its intensity from t = 90 and are divided at and t = 95. These vortices monotonically decay and split and nearly disappear at t = 160.
Figures 5 and 6 present instantaneous local Nusselt number distribution on lower and upper walls. The distribution within the separation bubbles is little changed. On the lower wall, heat transfer is slightly enhanced in the spanwise center region from t = 25. Two, four and six areas of high Nusselt number exist near the side walls at t = 39, 45 and 55, respectively. Inclined vortices come close to the lower wall and a pair of counter rotating vortices entrains cold fluid to the wall. These regions are identical to area of high Nusselt number. At t = 62, high Nusselt number area appears in the spanwise center region and gradually increases the value. At t = 82, Nusselt number takes a peak in the spanwise channel center and
(a)(b)
Figure 2. Comparison of streamwise velocity profile. (a) ER = 3, AR = 24, Re = 186.7; (b) ER = 2, AR = 16, Re = 203.3.
Figure 3. Side view of helicity isosurface (ER = 2, AR = 8, Re = 200, St = 0.05, ■: h = −0.2, ■: h = +0.2).
Figure 4. Bottom view of helicity isosurface (ER = 2, AR = 8, Re = 200, St = 0.05, ■: h = −0.2, ■: h = +0.2)
Figure 5. Local Nusselt number distribution on lower wall (ER = 2, AR = 8, Re = 200, St = 0.05).
gradually decreases in accordance with the decay of the cloud vortices. On the other hand, effect of pulsation on the upper wall is observed from t = 40 and slower than the lower wall. Nusselt number on the upper wall takes a peak near the side walls similar to that on the lower wall. However the locations on the upper wall are slightly up-stream or downstream from them on the lower wall.
Time series of space averaged Nusselt number on the upper and lower walls is shown in Figure 7. It is clear that the transient response of heat transfer depends on
Figure 6. Local Nusselt number distribution on upper wall (ER = 2, AR = 8, Re = 200, St = 0.05).
Figure 7. Time series of space averaged Nusselt number (ER = 2, AR = 8, Re = 200, St = 0.05).
Strouhal number. Nusselt number of St = 0.05, 0.10 increases early and Numean of St = 0.025, 0.20 follows them. It seems that enhancement of heat transfer is almost independent on Strouhal number and the direction of wall.
4. Concluding Remarks
Numerical simulations of three-dimensional separated flow and heat transfer in the symmetric enlarged rectangular channel was conducted. Transient response of the flow and heat transfer in the channel with the fully developed inlet flow imposed one cycle of a pulsating fluctuation was investigated. Main results obtained are summarized as follows.
The vortices shed from behind of the step drift downstream along the side walls and sequentially induce some pairs of the streamwise vortices near the reattachment region. These cloud vortices incline the walls and gradually decay.
Nusselt number on the walls takes a peak near the side walls on the upper and lower walls similarly. However the locations on the upper wall are slightly upstream or downstream from them on the lower wall.
In the case of St = 0.05, 0.10, space averaged Nusselt number increases early. However, enhancement of heat transfer is almost independent on Strouhal number of the pulsating fluctuation and the direction of wall.
5. Acknowledgements
These results were partly obtained with a supercomputer SX-9 of Cyberscience Center, Tohoku University.
Nomenclature
AR Aspect ratio (= W3/H)
ER Expansion ratio (= W2/W1)
h Helicity
H Step height Nu
Nu sselt number (= aH/l)
Numean Space averaged Nusselt number
p Pressure
Pr Prandtl number (= n /a)
qw Heat flux
Re Reynolds number (= U H/n)
St Strouhal number (= f H/U)
t Time
u Velocity vector (u, v, w)
U Mean velocity at inlet
Umax Maximum velocity at inlet
x Coordinates (x, y, z)
Greek Letters
a Heat transfer coefficient
q Temperature
l Heat conductivity
v Kinematic viscosity
NOTES