Ninth International Conference on Computational Fluid Dynamics (ICCFD9), Istanbul, Turkey, July -, ICCFD9-xxxx Inclined slot et impinging on a moving wall Benmouhoub Dahbia & Mataoui Amina USTHB, Laboratoire de mécanique des fluides théorique et appliquée, Faculté de physique, Université des sciences et de la technologie Houari Boumediene, B.P. Bab Ezzouar Al Alia, Alger, Algérie. Corresponding author: benmouhoub_d@yahoo.com & mataoui_amina@yahoo.fr Abstract: In this work, heat transfer from moving wall due an inclined plane turbulent et has been numerical investigated. The inclination of nozzle is required for the control of stagnation point location. The numerical predictions are based on RSM second order turbulence model, coupled with enhanced wall treatment. For a given impinging distance of nozzle thickness (e) and a Reynolds number of,the main parameter of this study is the surface-to-et velocity ratio Rs varying from to.7, requiring a et inclination between to.the calculations are in good agreement with available data. The numerical results show that heat transfer is greatly influenced by moving wall effect. The local Nusselt number decreases with increasing surface-to-et velocity ratios (R s =). However, the accurate inclination of the et enhances heat transfer and recovers the stagnation point location, of the interaction on a stationary wall. Keywords: Inclined slot et, Moving wall, CFD, Heat transfer, Nusselt number, stagnation point location. Introduction Heat transfer by impinging et are encountered in many industrial applications such as tempering and shaping glass, draying textiles and paper, cooling of electronic components and turbine blades, and annealing of continuous sheets of materials, etc. Particular interests were devoted to impinging ets due to their characteristics of heat and mass transfer. Impinging et on a moving wall were little bit investigated compared to the case of immobile wall. Among work implying a plane turbulent air et impinging a moving wall, we quote the experimental studies of Subba Rau and Schlünder [] of Van Heiningen [] and Senter [] or the numerical studies of Chattopadhay et al [, ], Huang et al [], Sharif et al [], Aghahani et al [9] and Benmouhoub and Mataoui []. All these works have reported the effect of the impingement surface motion on the flow field and heat transfer. The movement of the impacted wall causes the displacement of stagnation point and reduces its heat transfer. Senter [], from measurements of the case of perpendicular et impingement on a mobile wall, has examined the inclined case of in order to recover the characteristics of perpendicular impingement.
Free Boundary Free Boundary Indeed, for Reynolds number of and surface-to-et velocity ratio of, Senter [] confirmed an improvement of the Nusselt average number of % for a et inclination of. In this paper, we present the results of preliminary study in order to find the accurate inclination giving the stagnation point location of orthogonal et impinging an immobile wall Problem Statement. Geometry of problem The geometry and boundary condition are summarized in figure (). The flow is steady in average and fully turbulent. The fluid (air) is incompressible with constant thermo-physical properties. The et exit velocity is V. The moving wall is moving with velocity U corresponding to the surface-to-et velocity U ratio R s defined by: R s. The et at the exit is maintained at ambient temperature T V C, while the moving wall is heated at T H (T C =9K and T H =K). Furthermore due to the et inclination, the velocity components are: U V.cos( ) ; V V.sin( ) and Kinetic energy and dissipation rate are deduced from et parameters as follows: k k I. V where I. and.e Confinement plate (U=, V=, T= Tc) Inclined et U=V cos() V=V sin() T=Tc Confinement plate (U=, V=, T= Tc) Impinging immobile plate ( U=, V=, T= Tc) y O e Impinging moving hot plate ( U=U, V=, T= TH) x H Impinging immobile plate (U=, V=, T= Tc) L L Figure: Configuration and boundary conditions.. Governing equations Mass, momentum, and energy conservation equations for incompressible flow, constant fluid properties and steady state conditions are as follows: Mass conservative equation: Ui () x Momentum conservation equation: Ui P U i U uiu x xi x x Energy conservation equation: T T Ui ui () xi xi xi. Turbulence modeling i ()
In equations () and (), Reynolds stress component uu i and turbulent heat flux ui require modeling. The closure of these equations is achieved by means of linear strain pressure - Reynolds stress second order model. This model does not require eddy viscosity hypothesis. The differential equation Reynolds stress transport consists of standard Reynolds stress model based on dissipation equation ε. There are three versions of the standard Reynolds stress models. They are called: LRR-IP, LRRQI and SSG. Reece and Rodi [] developed the LRR-IP and LRR-QI models. The pressurestrain expression consists of an anisotropy tensor, mean strain rate tensor and vorticity tensor. The production due to the buoyancy is neglected in this study. By analogy with molecular transport, for all models (first or second order models), the Simple Gradient Diffusion Hypothesis (SGDH) is used. The following algebraic constitutive law is allows deducing the velocity-temperature correlation: T ui t x () t k where t and t= C Prt Near the wall, specific treatment is required. After several tests, enhanced wall treatment predicts the flow fields with a best accuracy [].. Numerical procedure The numerical predictions based on finite volume method are performed by ANSYS FLYENT. CFD code. The discritisation of each equation is achieved on collocated meshes. The algorithm SIMPLEC is applied for pressure-velocity coupling. For the interpolation of convection-diffusion terms: second order scheme for pressure and POWER LAW scheme for the other variable, are applied Patankar []... Grid arrangement A two dimensional structural non-uniform grid is generated. The enhanced wall treatment is used, requiring a fine grid size in the vicinity to each wall (Fig. ) - - - - - - Figure : Typical grid arrangement Several grids were tested in order to examine their independence on the solution. Fig. shows the effect of cells size on the evolution of local Nusselt number. It is noticed that the independence of the grid on the solution is obtained for 9 x 9 nodes beyond which no further signification change is observed. This grid distribution is therefore used in all calculations knowing that all geometrical parameters are not modified.
Figure : Effect of grid refinement on local Nusselt number along the moving wall. Re=, R s = H/e=, α=... Validation The validation of numerical technique is checked on the dynamic and thermal fields. The normalized velocities (U/V ) and (V/V ) for two cross sections (y=. and ) for two velocity ratios: Rs=, α= and Rs=, α=9 and Re=, were compared to PIV experimental data of Senter []. As shows figure, a good agreement was obtained. Figure also shows a good agreement with the experimental data of Gardon and Akfirat [] for the case of an immobile wall corresponding to R s =, α=, H/e= and Re=. Figure : Validation; mean velocity components U/V and V/V, Re=, for horizontal section =. and.
Figure : Comparison of local Nusselt number prediction with the experiment data of Gardon and Akfirat [] Results and Discussion Figure () shows the streamlines contours for Re=, for the cases of surface-to-et velocities ratios such R s.7 and their corresponding inclined cases which allow to get the case of a et impinging perpendicularly on a stationary wall ( α ). All results are compared to those of available experimental measurements (PIV) []. A good agreement is obtained with PIV measurement [] excepting the cases where R s. For R s =, α=, relating to stationary plate, the flow topology is perfectly symmetric relative to the et axis. Two identical eddies are observed, they are generated by the flow driving and by the upper wall confinement. When the wall moves from left to right < R s, a third contra-rotating eddy appears at the left part of the moving wall. This new recirculation reduces the size of the left principal vortex and thus the et is deviated towards the right part of the mobile wall. The location of the stagnation point is therefore shifted. For R s >, the et is completely deviated by the movement of the wall. The et does not attain the moving wall and its streamlines become progressively parallel to the wall when the surface-to-et velocities ratio increases. The shear drives the flow by the moving wall is much larger than that of the left region. The inclination of the nozzle is required in order to recover the case of et impinging perpendicularly to stationary wall. These inclinations are obtained following several tests. When the velocity ratio increases, the suitable inclination of the et augments. Figure (7) illustrates a comparison for each components of the average velocity U/V and U/ V, between the inclined and orthogonal cases, through the cross section =.. From this figure, for < R s : we observe a good similarity of the value of orthogonal et impinging on a fixed surface and those of the inclined et interaction on a moving surface in the stagnation points region. However, for < R s.7, the profiles are different. For this range of velocity ratio, the flow is parallel to the wall even with the inclination of the nozzle; the et does not hit the moving wall.
7 - - - - - - PIV [] RSM (a) Rs=, α= PIV [] : Rs=., α= PIV [] : Rs=, α= PIV [] : Rs=., α= 7 - - - - - - - - - Rs=., α= Rs=., α= - - - Rs=, α= 7 - - - - - - - - Rs=., α= - - Rs=., α= - - - - - Rs=., α= Rs=.7, α= - - - Rs=, α=9 - - - - Rs=., α=. - Rs=., α=.7 - - - - Rs=., α=7.7 - - - Rs=.7, α= (b) Rs, α Figure : Streamlines contours: effect of plate velocity ratio and inclination angle Re=; Rs.7; α
Figure 7: Effect of the inclination of the et on transverse profiles of average velocity U/V and U/ V Re=, =.. The dimensionless Nusselt number defines the heat transfer of fluid along the heated impinging wall, which characterizes the ratio of convective to conductive heat transfer across the corresponding wall. The local Nusselt number distribution over the impingement wall is defined in Eq.: e T Nu( x) T T n H C ywall Where n is the perpendicular direction to the corresponding wall. Figure (a) shows the evolution of local Nusselt number Nu(x) along the moving plate for Re= and for various values of R s. It is observed from this figure that, for stationary plate (R s =), the curve is symmetrical, it is characterized by a maximum at the stagnation points which are located on the et axis. On either side of this axis the local Nusselt number decreases gradually. While for the case of moving wall R s, we note a peak appears at the left end of the impinging wall (=-), which augments when R s increases. The development of this maximum can be explained by the velocity discontinuity between the moving part and the left immobile wall part []. For the cases:. R s The maximum value decreases and its position is deflected in the same direction as the moving wall, this is due to the shear generated by plate movement. A drop occurs in the left part of the plate. This reduction is related with the apparition of the new recirculation []. Local Nusselt number augments on the right part of the plate when R s increases by the effect of the progressive et deviation (Benmouhoub and Mataoui[9]). For R s > For these cases, the plate velocity exceeds that of the et. The maximum value of local Nusselt number progressively disappears at the left side of the wall. The Nu(x) increases significantly, this is due to the shear generated by the movement of the wall. Figure (b) illustrates the distribution of the local Nusselt number for each surface-to-et velocity ratio R s, the inclined et is in the opposite direction to the movement of the wall with an suitable angle to reach the maximum at the same point as the orthogonal impact et on an immobile wall. It was found that the inclination of the nozzle increases with increasing velocity ratio R s. The value of the local Nusselt number at the stagnation point increases compared to the case of the vertical et (Fig. 9). () 7
Table shows the increasing percentage of local Nusselt at the stagnation points between the two cases (inclined and non-inclined et). (a) (b) Figure : Distribution of local Nusselt number along the hot moving wall for Re=: (a) perpendicular et; (b) inclined et [9] Figure 9: Nusselt number at stagnation point Re=. Table. Gain of local Nusselt number at stagnation points R s...7...7 δ(%).. 9.7..77.79.79 Conclusion The global structure and heat transfer due to confined slot et impinging on a moving wall was investigated numerically for a given range of plate velocity ratios and et inclination. The RSM turbulence model with enhanced wall treatment for near wall turbulence modeling is used. Results are in good agreement with previous experimental measurement. For Rs, the deviation of the stagnation point from the et exit axis generates an important variation of the heat transfer. With increasing velocity plate ratio, the magnitude of maximum local Nusselt number decreases and its position is deviated in the direction of movement of the wall. However, the accurate inclination of the et in the opposite direction to the movement of the wall considerably improves local heat transfer in the region of the stagnation point.
References [] Subba Rau, K., & Schlünder, E. U. (977). Heat transfer between an impinging et and a continuously moving surface. Wärme-stoffübertr,, -. [] Van Heiningen, A. R. P., Muumdar, A. S., & Douglas, W. J. M. (977). Flow and heat transfer characteristics of a turbulent slot et impinging on a moving wall. In Symposium on Turbulent Shear Flows (Vol., pp. -9). [] Senter, J,. Senter. Analyse experimentale et numérique des écoulements et des transferts de chaleur convectifs produits par un et plan impactant une plaque plane mobile. Phd, thesis Université de Nanates, France,. [] Chattopadhyay, H., Biswas, G., & Mitra, N. K. (). Heat transfer from a moving surface due to impinging slot ets. Journal of heat transfer, (), -. [] Chattopadhyay, H., & Saha, S. K. (). Turbulent flow and heat transfer from a slot et impinging on a moving plate. International Journal of Heat and Fluid Flow, (), -97. [] Huang, P. G., Muumdar, A. S., & Douglas, W. J. M. (9). Numerical prediction of fluid flow and heat transfer under a turbulent impinging slot et with surface motion and crossflow. ASME paper, (-WA), -. [7] Sharif, M. A. R., & Baneree, A. (9). Numerical analysis of heat transfer due to confined slotet impingement on a moving plate. Applied Thermal Engineering, 9(), -. [] Aghahani, M., Eslami, G., & Hadidi, A. (). Heat transfer in a turbulent et impinging on a moving plate considering high plate-to-et velocity ratios. Journal of Mechanical Science and Technology, (), 9-. [9] Benmouhoub, D., & Mataoui, A. (). Inclination of an impinging et on a moving wall to control the stagnation point location. International Journal of Thermal Sciences, 9, 9-. [] Mobtil, M., Bougeard, D., & Solliec, C. (). Inverse determination of convective heat transfer between an impinging et and a continuously moving flat surface. International Journal of Heat and Fluid Flow,, -9. [] Launder, B. E., Reece, G. J., & Rodi, W. (97). Progress in the development of a Reynoldsstress turbulence closure. Journal of fluid mechanics, (), 7-. [] Kader, B. A. (9). Temperature and concentration profiles in fully turbulent boundary layers. International Journal of Heat and Mass Transfer, (9), -. [] Patankar, S. V. (9). Numerical heat transfer and fluid flow: Computational methods in mechanics and thermal science. [] Gardon, R., & Akfirat, J. C. (9). Heat transfer characteristics of impinging two-dimensional air ets. Journal of Heat Transfer, (), -. 9