Study of flow and heat transfer in a 3D pin-fin channel

Size: px
Start display at page:

Download "Study of flow and heat transfer in a 3D pin-fin channel"

Transcription

1 Study of flow and heat transfer in a 3D pin-fin channel Jean-Eloi Lombard Swiss Institute of Technology Lausanne Supervisor : Dr. Mark Sawley Assistants : Laura Augello, Thierry Cornu December 27, 2011

2

3 Contents 1 Introduction 1 2 Problem definition 1 3 Techniques employed Mesh Turbulence model Numerical methods Results and discussion Convergence Conservation of heat Qualitative post-processing Contours Vector field Pressure and Temperature contours Velocity and Temperature profile Streamlines Velocity and Temperature plane-cut Quantitative post-processing Comparison with experimental extrapolation Total pressure drop Total heat flux Pin heat-flux and position Cooling Efficiency Proposed Improvements Numerical solution improvements Design improvement of the fin-pin channel geometry Conclusion 15

4

5 List of Figures 1 Rectangular pin-fin channel. The cross-flow enters through the front left side of the channel. The upper solid wall has been removed in the figure for visibility 1 2 Geometrical parameters. The array is centered on the origin, D = 12.5[mm], L x,y = 31.25[mm], X i = 125[mm], X c = 312.5[mm], Y c = 31.25, Z c = 12.5[mm] 1 3 Overview of the fine 59,964 hexahedral mesh from +Z View of the boundary layers around each pin in the fine 59,964 hexahedral mesh from +Z From left to right y+ values for each of the four full pins (Pin1, Pin2, Pin3 and Pin4) Convergence of the residuals of continuity, energy, k and ɛ as a function of iteration number for the three available wall functions : standard, unsteady and enhanced Convergence of the heat flux for all 10 pins as a function of iteration number for the three available wall functions : standard, unsteady and enhanced Relative pressure field distribution on the symmetry plane (z=0) Absolute temperature distribution on the symmetry plane (z=0) Velocity magnitude distribution on the symmetry plane (z=0) Velocity vector field in flow in near pin Enlarged view of the velocity vector field at the detachment zone on pin Relative pressure distribution on the pins and absolute temperature distribution on the symmetry plane Temperature and x component of velocity profile on line x = 0 z = 6.25[mm] Temperature and x component of velocity profile on line x = 0 y = A rake of 20 points (in z = 0[mm], x = 70[mm]) on the symmetry plane is defined from which the streamlines originate. They are computed both in the forward and backward directions x component of velocity in the plane (x = -20[mm]) The white interstice is Pin Absolute temperature in the plane (x = -20[mm]) The white interstice is Pin Aboslute pressure distribution on the symmetry plane. Each peak correspond to the point d arret of the pins Iso-surface of turbulence intensity. The down-wind pin keeps the high-turbulence flow near the upward pin. Note the difference in distance of the turbulence iso-surface for pins 3 and

6

7 1 Introduction The purpose of this study is to qualify and quantify the steady state flow and heat transfer capability of a 3D fin-pin array comprising of 32 pins arranged in a staggered arrangement located in a rectangular channel with solid walls. Figure 1 Rectangular pin-fin channel. The cross-flow enters through the front left side of the channel. The upper solid wall has been removed in the figure for visibility 2 Problem definition For the computations only a representative periodic section of the channel is considered the scheme of which is presented in (Fig. 2). The 3D model is obtained by extruding this 2D geometry along the z-axis. The cooling cross-flow is in the positive x direction. The pin-fin array consists of 7 stream-wise and 2 crosswise rows of equal-sized circular pins or radius (4 full pins and 6 half pins) in a uniformly spaced staggered arrangement. In this model, the two side surfaces (y = ±Y c /2) are considered to be periodic. The 3D geometry is obtained by extruding this 2D cross-section in the z direction. Furthermore since the flow is assumed symmetric in the z direction the flow is only computed within 0 z Z c /2. The cross-flow is air with the following properties : ρ [kg m 3 ] µ [P a s] C p [J kg 1 K 1 ] k [W m 1 K 1 ] Table 1 Properties of air Figure 2 Geometrical parameters. The array is centered on the origin, D = 12.5[mm], L x,y = 31.25[mm], X i = 125[mm], X c = 312.5[mm], Y c = 31.25, Z c = 12.5[mm] 1

8 Boundary conditions. At the inlet of the channel the flow velocity field U i = 5[m.s 1 ] is uniform and temperature is constant T i = 288[K]. The wall surface of each pin is at a fixed constant temperature T w = 398[K]. A uniform pressure is assumed at the outlet of the channel. The boundary at the side walls is defined as symmetric, the lower solid wall wall of the pin-fin array is adiabatic and has a no-slip boundary condition. Only half of the height in the z direction is computed so in z = 0[mm] correspond to a symmetric boundary condition. The effect of gravity is neglected. 3 Techniques employed 3.1 Mesh An automated mesh generator is used because it offers quick and descent meshes for computation use in simple flow cases. Both a coarse 33,006 and fine 59,964 hexahedral unstructured mesh are generated in two steps using : 1. first a 2D mesh is generated with a boundary layer around the pins (the radius of each pin is segmented into 51 uniform parts and the cells are formed with a uniform growth factor of 1.15 for 20 cells). For minimum computation error the regions before and after the array have been meshed by a cartesian grid. Hexahedral mapping is used for the region with the pins. Each pin has a boundary layer mesh. 2. the 2D mesh is extruded with the geometry along the z-direction. Figure 3 Overview of the fine 59,964 hexahedral mesh from +Z. Figure 4 View of the boundary layers around each pin in the fine 59,964 hexahedral mesh from +Z. 2

9 Mesh comparison. The coarse mesh is not able to correctly resolve the flow near the viscous sublayer even with a wall function (y+ max = 189 >> 30). For this reason the fine mesh is used to obtain all the significant results on flow and heat transfer in the remainder of this paper. Coarse Mesh Refined Mesh Cell # y+ max Cell quality % 80% Table 2 Mesh quality comparison Verification and validation For the non-equilibrium wall function to yield good results for the modeling of turbulence near the pin walls y+ should be smaller than 30 which is achieved with the fine mesh y+ [ ] x [m] Figure 5 From left to right y+ values for each of the four full pins (Pin1, Pin2, Pin3 and Pin4). 3.2 Turbulence model Although SST-k ω model has been shown to provide one of the most accurate results for heat and flow modeling amongst RANS methods in pin-fin configurations [6] the coarse cell limit imposed by the guidelines does not allow the construction of a mesh that guaranties y+ 1 at the walls. Due to this limitation it would not be feasible to use any 3 or more equation models or even the SST k ω which has been shown to provide most precise flow and heat transfer predictions. Hence the choice must be made between the k ɛ and 3

10 Spallart-Allmaras models. k ɛ being a two equation model is bound to better describe the flow. Even-though it requires solving an extra solution at every iteration the small grid used relative to computational power available renders this model the optimal choice. The choice of the k ɛ turbulence model is dictated by the size of the mesh and will not provide the most accurate prediction of flow since it doesn t predict flow detachment correctly [5] and since it is an important feature of cooling (an important fraction of cooling is done within the re-circulation zone) this model is not optimal. Cmu 0.09 C1-epsilon 1.44 C2-epsilon 1.92 TKE Prandtl Number 1 TDR Prandtl Number 1.3 Energy Prandtl Number 0.85 Wall Prandtl Number 0.85 Table 3 Parameters of k ɛ model Regarding wall treatment the non-equilibrium wall treatment is better suited for complex flows involving separation and strong pressure and temperature gradients (Fig. 6) for a mesh that does not guarantee a y+ 1 but still has y+ < Numerical methods A Pressure-Implicit with Splitting of Operator scheme 1 is used for the Pressure-Velocity coupling because it can take into account skewed cells and second-order spatial discretization is used for gradient, pressure, momentum, turbulent kinetic energy, turbulent dissipation rate and energy to minimize error. 1 the meshes are generated with ANSYS Gambit, the solutions computated with ANSYS Fluent and visualization with Paraview 4

11 4 Results and discussion 4.1 Convergence During the computations, monitor the residuals and the average heat flux at the surface of the pins. Provide detailed comments on the observed convergence history and whether it appears adequate to provide an accurate flow solution. A convergence of energy, continuity, k, and ɛ better than 5 orders of magnitude is obtained. Continuity Residual Std. k ε std. wall fct. Std. k ε enh. wall fct. Std. k ε unst. wall fct. Energy Residual Std. k ε std. Wall fct. Std. k ε Enhanced Wall fct. Std. k ε unst. wall fct Iterations Iterations Std. k ε std. Wall fct. Std. k ε Enhanced Wall fct. Std. k ε unst. wall fct Std. k ε std. Wall fct. Std. k ε Enhanced Wall fct. Std. k ε unst. wall fct. k Residual ε Residual Iterations Iterations Figure 6 Convergence of the residuals of continuity, energy, k and ɛ as a function of iteration number for the three available wall functions : standard, unsteady and enhanced. 5

12 Eventhough the k ɛ turbulence model with the unsteady wall function offers the best convergence properties in terms of residuals, there is no guarantee a priori it is converging to the correct physical solution. As can be seen in (Fig. 7) the three different wall functions result in completely different heat fluxes at convergence. There difference is of a factor 2 between unsteady wall function and the enhanced wall function while there is a 50% difference in predicted heat flux between the unsteady wall function and the standard one. The only to validate the correct wall function is by comparing the computed Nusselt number to the theoretical one using the Grimison correlation which give a 4% error for the unsteady wall function and more than 50% errors for both standard and enhanced wall functions Std. k ε std. Wall fct. Std. k ε Enhanced Wall fct. Std. k ε unst. wall fct. Total Pin Hit flux [W.m 2 ] Iterations Figure 7 Convergence of the heat flux for all 10 pins as a function of iteration number for the three available wall functions : standard, unsteady and enhanced. 4.2 Conservation of heat Thermal energy is conserved within in better than 5 orders of magnitude : Surface Heat transfer rate [W] Inlet Outlet Pin Pin Pin Pin Pin Pin Pin Pin Pin Pin Net Table 4 Energy conservation 6

13 5 Qualitative post-processing As expected the flow around a cylinder at Re 7000 exhibits a turbulent recirculation zone but the use of RANS does not allow for an appreciation of the vortex shedding taking place behind each pin. 5.1 Contours The relative pressure distribution betrays the recirculation zones behind the pins. Figure 8 Relative pressure field distribution on the symmetry plane (z=0). The absolute temperature distribution shows how the recirculation zones behind the pins works to advect the heat from the pins. Figure 9 Absolute temperature distribution on the symmetry plane (z=0). The velocity magnitude distribution shows a higher velocity of the flow between the rows of pins, with vanishing velocity at the boundary of the pins. Also the recirculation zones behind the pins have low velocity magnitude relative to the main flow. Figure 10 Velocity magnitude distribution on the symmetry plane (z=0). 7

14 5.2 Vector field Four main features of the flow can be seen in (Fig. 11), a broad view of the velocity vector field around pin 1. Firstly the impinging point in the bottom left-corner, secondly, the Figure 11 Velocity vector field in flow in near pin 1. recirculation zone in the down-flow of pin 1, thirdly the high-velocity main flow with a higher density of longer arrows and finally the detachment zone an enlarged view of which is provided in (Fig. 12). Figure 12 Enlarged view of the velocity vector field at the detachment zone on pin 1. 8

15 5.3 Pressure and Temperature contours As can be seen in (Fig. 13) the two components of heat dissipation in the flow are the horse shoe vortex on the leading edge of the pins with high temperature gradient and high pressure and the turbulent mixing in the recirculation behind the trailing-edge. Note the ripples in Figure 13 Relative pressure distribution on the pins and absolute temperature distribution on the symmetry plane. the back flow of Pin1 in the pressure distribution on the symmetry plane 5.4 Velocity and Temperature profile In (Fig. 14) it is interesting to note that temperature reaches the boundary condition temperature of the pins of 398[K] at the boundaries of pins 6 and 9 and also the recirculation zone downwind of pin 2 that betrays the a higher temperature relative to the high velocity region where it is minimal. The correlation is evident when looking at the x component of velocity on the same line. The velocity vanishes at the wall of pins 6 and 9 and even though the x component of velocity is barely negative one can still see the recirculation zone. Temperature [K] v x [K] y [mm] y [mm] Figure 14 Temperature and x component of velocity profile on line x = 0 z = 6.25[mm] In (Fig. 15) it is interesting to note that the temperature of the flow increases from the bottom to the symmetry plane but that the x component of the velocity profile is more 9

16 complex. The no-slip boundary condition is verified on the bottom plane the negative x component of the velocity betrays the recirculation zone down-wind of pin 2. Temperature [K] v x [K] z [mm] x z [mm] x 10 3 Figure 15 Temperature and x component of velocity profile on line x = 0 y = Streamlines (Fig. 16) presents the streamlines on the symmetry plane and allows for an appreciation of the size of the recirculation zones behind the full pins (1,2,3, and 4) and the relative proximity of adjacent streamlines. The biggest recirculation zone is the last pin (pin 4) and pins 1 through 3 have roughly the same recirculation geometry. The streamlines closest one to another near pins 2 and 3 yielding high velocity flow near the walls of pins 2 and 3 which allows for a delayed separation and better heat advection. Figure 16 A rake of 20 points (in z = 0[mm], x = 70[mm]) on the symmetry plane is defined from which the streamlines originate. They are computed both in the forward and backward directions. It is also interesting to note the recirculation that is clearly put in evidence behind pins 2 and 3 and that because of the turbulent flow it is not perfectly symmetric to the X Z plane. 10

17 5.6 Velocity and Temperature plane-cut The distribution of the x component of velocity in the plane (x = -20[mm] i.e. intersecting Pin2 ) in (Fig. 17) shows the vanishing velocity (in white) near the wall of the pin and the recirculation zone behind pins 5 and 8 with negative and zero x component of velocity (in black) and the high velocity flow between these regions (in light-grey). Figure 17 x component of velocity in the plane (x = -20[mm]) The white interstice is Pin 2. The temperature distribution somewhat follows the velocity distribution in the plane x = 20[mm] with the high velocity regions remaining at a temperature near the inlet conditions 288[K] because heat is advected by the thin region in red (temperature of 398[K] at the pin wall) near the pin wall were the velocity is much lower. Note in light blue the temperature of the recirculation zones behind pins 5 and 8. Figure 18 Absolute temperature in the plane (x = -20[mm]) The white interstice is Pin 2. 11

18 6 Quantitative post-processing 6.1 Comparison with experimental extrapolation For the fin-pin channel flow the Reynolds number is defined as : Re D,max = ρu maxd µ with U max = and the average Nusselt number 2 Nu D = qd AkT in U il y L y D = where q = 34.33[W ] it the total heat flow normal to the pin walls, k = [W m 1 K 1 ] si the thermal conductivity, A = π D/2 10 Z c /2 = π / [m 2 ] is the total surface area of the pins and T in = T o T ( i ) log Tw Ti T w T o with T i,t o and T w the inlet, outlet and pin-wall temperature respectively. The average Nusselt number for the fin-pin array can be approximated using the Grimison correlation : Nu D = C 1 Re m D,max with c 1 = 0.521, m = Hence the error in computed Nusselt number relative to the Grimison correlation is 4% : Nu err D = 3.94% One possible source of error is the k ɛ model that cannot (because of the isotropic assumption) accurately predict separation and suffers from over prediction of turbulent energy at the impinging point. 2 is the ratio of convective to conductive heat transfer (i.e. a small Nusselt number indicate heat is mainly transferred through conduction whereas a large Nusselt number indicates heat is mostly dissipated by convection). 12

19 6.2 Total pressure drop Pressure [Pa] x [m] Figure 19 Aboslute pressure distribution on the symmetry plane. Each peak correspond to the point d arret of the pins. The absolute pressure difference between inlet and outlet is computed to be : p = = 140 ± 1[P a] hence a relative low pressure drop of 0.14% ± Hence the energy required to push the air through the pin-fin array is small. 6.3 Total heat flux Pin Heat transfer rate [W] Pin Pin Pin Pin Pin Pin Pin Pin Pin Pin Total Pin Heat transfer rate per unit area [mw m 2 ] Pin Pin Pin Pin Pin Pin Pin Pin Pin Pin Table 5 Total heat flux computed from the heat transfer rate of each pin and heat transfer rate per unit area. The pin-fin array section is computed to have a heat transfer rate of 48.80±0.01[W ] with ± 0.01 [W] from the row of 4 pins and ± 0.01[W ] for the row with 3 pins. Hence the global design with 5 rows of 4 pins and 3 rows of 3 pins has total power dissipation of ± 0.01[W ]. 13

20 6.4 Pin heat-flux and position Pin1 is exposed to natural flow (i.e. identical to cylinder in free flow), whereas second, third, fourth pins are exposed to highly turbulent flow similarly to pins in the rows of 3. The extend and efficiency of the turbulence is due to two factors : upwind flow they are in the recirculation zones, and secondly the downwind pin limits the size of the recirculation zone. Hence the most efficient pins are Pins 2 and 3. Full Pin Heat transfer rate [W] Pin Pin Pin Pin Table 6 Heat transfer rate for the 4 pins in the row of 4. Figure 20 Iso-surface of turbulence intensity. The down-wind pin keeps the highturbulence flow near the upward pin. Note the difference in distance of the turbulence iso-surface for pins 3 and 4. 14

21 6.5 Cooling Efficiency Provide insights into the heat transfer between the cross-flow fluid and the pins. Comment on the cooling efficiency of the pin-fin channel. For efficient cooling a pin-fin array must maximize heat transfer augmentation while minimizing pressure drop (i.e. drop in velocity magnitude). Cooling efficiency can be computed as the ratio between heat transferred from the pins and the kinetic energy of the inlet air. v in = 5[m s 1 ] and v out = 4.87[m s 1 ] and T in = 288[K] whereas T out = [K]. The coefficient of performance (COP) [10] of the fin-pin array is : COP = Q flow W where Q flow is the difference between heat of flow and outlet and inlet and W the amount of work required for insuring the inlet conditions, hence : COP = cṁr T 2M ṁ in v 2 in 2 = 5R T 2Mv 2 in 1.17 where R = 8.314[J K 1 mol 1 ] is the ideal gas constant, M = 28.97[g mol 1 ] is molecular mass of air, v in = 5[m s 1 ] the inflow velocity. Hence the pin-fin array with the given boundary conditions evacuates 17% more power than is required cost to produce the flow. 7 Proposed Improvements 7.1 Numerical solution improvements Even-though a high order LES with a much finer grid would yield more accurate results, if the purpose of this work is to determine the efficiency of the proposed pin-fin array for engineering purposes, then limit in size of the grid is not a limiting factor for the quality of the flow and heat transfer properties required. 7.2 Design improvement of the fin-pin channel geometry. Since we have shown that the staggered array is efficient insofar that produces turbulent flow throughout the array it might be interesting to add the number of pins in each row. Also the velocity of the inlet flow might influence the overall efficiency by developing more turbulent flow in the array. 8 Conclusion The 32 pin-fin array is computed to dissipate ± 0.01[W ] with a COP of The high average Nusselt number indicates a high ratio of convective to conductive heat transfer. Correlation to experimental data by means of the Grimison correlation shows good agreement with 4% error. A parametric study with different inlet velocities and pin temperatures is required to optimize the operating efficiency. 15

22 References [1] W. L. Oberkampf and T. G. Trucano. Verification and validation in computational fluid dynamics. Progress in Aerospace Sciences, 38(3): , [2] S. Krajnović. Large eddy simulation of the flow over a three-dimensional hill. Flow, Turbulence and Combustion, [3] R. Webb and E. Eckert. Application of rough surfaces to heat exchanger design. International Journal of Heat and Mass Transfer, 15(9): , [4] M. E. Lyall, A. A. Thrift, K. A. Thole, and A. Kohli. Heat transfer from low aspect ratio pin fins. Journal of Turbomachinery, 133(1):011001, [5] R. Bahadur and A. Bar-Cohen. Orthotropic thermal conductivity effect on cylindrical pin fin heat transfer. International Journal of Heat and Mass Transfer, 50(5-6): , [6] W. Khan, J. Culham, and M. Yovanovich. Convection heat transfer from tube banks in crossflow: Analytical approach. International Journal of Heat and Mass Transfer, 49(25-26): , [7] Stephen B. Pope Turbulent Flows Cambridge University Press, 2000 [8] Incropera F.P., DeWitt D.P., Bergman T.L. & Lavine A.S. Fundamentals of Heat and Mass Transfer Wiley, 6th edition, [9] Lienhard J.H. & Lienhard J.H. A Heat Transfer Handbook 4th edition, 2011 [10] Hans U. Fuchs The Dynamics of Heat Springer-Verlag, 1996, p

ENERGY PERFORMANCE IMPROVEMENT, FLOW BEHAVIOR AND HEAT TRANSFER INVESTIGATION IN A CIRCULAR TUBE WITH V-DOWNSTREAM DISCRETE BAFFLES

ENERGY PERFORMANCE IMPROVEMENT, FLOW BEHAVIOR AND HEAT TRANSFER INVESTIGATION IN A CIRCULAR TUBE WITH V-DOWNSTREAM DISCRETE BAFFLES Journal of Mathematics and Statistics 9 (4): 339-348, 2013 ISSN: 1549-3644 2013 doi:10.3844/jmssp.2013.339.348 Published Online 9 (4) 2013 (http://www.thescipub.com/jmss.toc) ENERGY PERFORMANCE IMPROVEMENT,

More information

Analysis of the Cooling Design in Electrical Transformer

Analysis of the Cooling Design in Electrical Transformer Analysis of the Cooling Design in Electrical Transformer Joel de Almeida Mendes E-mail: joeldealmeidamendes@hotmail.com Abstract This work presents the application of a CFD code Fluent to simulate the

More information

Turbulent Boundary Layers & Turbulence Models. Lecture 09

Turbulent Boundary Layers & Turbulence Models. Lecture 09 Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects

More information

International Conference on Energy Efficient Technologies For Automobiles (EETA 15) Journal of Chemical and Pharmaceutical Sciences ISSN:

International Conference on Energy Efficient Technologies For Automobiles (EETA 15) Journal of Chemical and Pharmaceutical Sciences ISSN: HEAT TRANSFER ENHANCEMENT WITH PRESSURE LOSS REDUCTION IN COMPACT HEAT EXCHANGERS USING VORTEX GENERATORS Viswajith M V*, Gireesh Kumaran Thampi, James Varghese Department of Mechanical Engineering, School

More information

Validation 3. Laminar Flow Around a Circular Cylinder

Validation 3. Laminar Flow Around a Circular Cylinder Validation 3. Laminar Flow Around a Circular Cylinder 3.1 Introduction Steady and unsteady laminar flow behind a circular cylinder, representing flow around bluff bodies, has been subjected to numerous

More information

Calculations on a heated cylinder case

Calculations on a heated cylinder case Calculations on a heated cylinder case J. C. Uribe and D. Laurence 1 Introduction In order to evaluate the wall functions in version 1.3 of Code Saturne, a heated cylinder case has been chosen. The case

More information

3D Numerical Study on Laminar Forced Convection in V-Baffled Square Channel

3D Numerical Study on Laminar Forced Convection in V-Baffled Square Channel American Journal of Applied Sciences 10 (10): 1287-1297, 2013 ISSN: 1546-9239 2013 Boonloi and Jedsadaratanachai, This open access article is distributed under a Creative Commons Attribution (CC-BY) 3.0

More information

CFD Analysis of Forced Convection Flow and Heat Transfer in Semi-Circular Cross-Sectioned Micro-Channel

CFD Analysis of Forced Convection Flow and Heat Transfer in Semi-Circular Cross-Sectioned Micro-Channel CFD Analysis of Forced Convection Flow and Heat Transfer in Semi-Circular Cross-Sectioned Micro-Channel *1 Hüseyin Kaya, 2 Kamil Arslan 1 Bartın University, Mechanical Engineering Department, Bartın, Turkey

More information

Numerical simulation of fluid flow in a monolithic exchanger related to high temperature and high pressure operating conditions

Numerical simulation of fluid flow in a monolithic exchanger related to high temperature and high pressure operating conditions Advanced Computational Methods in Heat Transfer X 25 Numerical simulation of fluid flow in a monolithic exchanger related to high temperature and high pressure operating conditions F. Selimovic & B. Sundén

More information

Performance of Elliptical Pin Fin Heat Exchanger with Three Elliptical Perforations

Performance of Elliptical Pin Fin Heat Exchanger with Three Elliptical Perforations www.cfdl.issres.net Vol. 3 (2) June 2011 Performance of Elliptical Pin Fin Heat Exchanger with Three Elliptical Perforations Monoj Baruah 1, Anupam Dewan 2c and P. Mahanta 1 1 Department of Mechanical

More information

On the transient modelling of impinging jets heat transfer. A practical approach

On the transient modelling of impinging jets heat transfer. A practical approach Turbulence, Heat and Mass Transfer 7 2012 Begell House, Inc. On the transient modelling of impinging jets heat transfer. A practical approach M. Bovo 1,2 and L. Davidson 1 1 Dept. of Applied Mechanics,

More information

External Forced Convection :

External Forced Convection : External Forced Convection : Flow over Bluff Objects (Cylinders, Spheres, Packed Beds) and Impinging Jets Chapter 7 Sections 7.4 through 7.8 7.4 The Cylinder in Cross Flow Conditions depend on special

More information

International Journal of Scientific & Engineering Research, Volume 6, Issue 5, May ISSN

International Journal of Scientific & Engineering Research, Volume 6, Issue 5, May ISSN International Journal of Scientific & Engineering Research, Volume 6, Issue 5, May-2015 28 CFD BASED HEAT TRANSFER ANALYSIS OF SOLAR AIR HEATER DUCT PROVIDED WITH ARTIFICIAL ROUGHNESS Vivek Rao, Dr. Ajay

More information

Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders

Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders A. Jugal M. Panchal, B. A M Lakdawala 2 A. M. Tech student, Mechanical Engineering Department, Institute

More information

Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza

Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza REPRESENTING PRESENCE OF SUBSURFACE CURRENT TURBINES IN OCEAN MODELS Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza 1 Momentum Equations 2 Effect of inclusion of Coriolis force

More information

A Numerical Study of Forced Convection Heat Transfer for Staggered Tube Banks in Cross-Flow

A Numerical Study of Forced Convection Heat Transfer for Staggered Tube Banks in Cross-Flow A Numerical Study of Forced Convection Heat Transfer for Staggered Tube Banks in Cross-Flow T. A. Tahseen 1, M. Ishak 1,2 and M. M. Rahman 1,2 1 Faculty of Mechanical Engineering, University Malaysia Pahang

More information

CFD Time Evolution of Heat Transfer Around A Bundle of Tubes In Staggered Configuration. G.S.T.A. Bangga 1*, W.A. Widodo 2

CFD Time Evolution of Heat Transfer Around A Bundle of Tubes In Staggered Configuration. G.S.T.A. Bangga 1*, W.A. Widodo 2 CFD Time Evolution of Heat Transfer Around A Bundle of Tubes In Staggered Configuration G.S.T.A. Bangga 1*, W.A. Widodo 2 1,2 Department of mechanical engineering Field of study energy conversion Institut

More information

COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF A V-RIB WITH GAP ROUGHENED SOLAR AIR HEATER

COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF A V-RIB WITH GAP ROUGHENED SOLAR AIR HEATER THERMAL SCIENCE: Year 2018, Vol. 22, No. 2, pp. 963-972 963 COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF A V-RIB WITH GAP ROUGHENED SOLAR AIR HEATER by Jitesh RANA, Anshuman SILORI, Rajesh MAITHANI *, and

More information

Numerical Heat Transfer Study of Turbulent Square Duct Flow through W-Type Turbulators

Numerical Heat Transfer Study of Turbulent Square Duct Flow through W-Type Turbulators search Article International Journal of Thermal Technologies E-ISSN 2277 4114 214 INPRESSCO, All Rights served Available at http://inpressco.com/category/ijtt/ merical Heat Transfer Study of Turbulent

More information

CFD Analysis for Thermal Behavior of Turbulent Channel Flow of Different Geometry of Bottom Plate

CFD Analysis for Thermal Behavior of Turbulent Channel Flow of Different Geometry of Bottom Plate International Journal Of Engineering Research And Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 13, Issue 9 (September 2017), PP.12-19 CFD Analysis for Thermal Behavior of Turbulent

More information

Natural Convection from a Long Horizontal Cylinder

Natural Convection from a Long Horizontal Cylinder Natural Convection from a Long Horizontal Cylinder Hussein Awad Kurdi Saad Engineering Technical College of Al Najaf, Al-Furat Al-Awsat Technical University, Iraq ABSTRACT: Natural convection from a Long

More information

Active Control of Separated Cascade Flow

Active Control of Separated Cascade Flow Chapter 5 Active Control of Separated Cascade Flow In this chapter, the possibility of active control using a synthetic jet applied to an unconventional axial stator-rotor arrangement is investigated.

More information

NUMERICAL SIMULATION OF FLOW THROUGH TURBINE BLADE INTERNAL COOLING CHANNEL USING COMSOL MULTIPHYSICS

NUMERICAL SIMULATION OF FLOW THROUGH TURBINE BLADE INTERNAL COOLING CHANNEL USING COMSOL MULTIPHYSICS International Journal of Emerging Technology and Innovative Engineering Volume 1, Issue 12, December 2015 (ISSN: 2394 6598) NUMERICAL SIMULATION OF FLOW THROUGH TURBINE BLADE INTERNAL COOLING CHANNEL USING

More information

A CFD Analysis Of A Solar Air Heater Having Triangular Rib Roughness On The Absorber Plate

A CFD Analysis Of A Solar Air Heater Having Triangular Rib Roughness On The Absorber Plate International Journal of ChemTech Research CODEN( USA): IJCRGG ISSN : 0974-4290 Vol.5, No.2, pp 964-971, April-June 2013 ICGSEE-2013[14th 16th March 2013] International Conference on Global Scenario in

More information

MAE 598 Project #1 Jeremiah Dwight

MAE 598 Project #1 Jeremiah Dwight MAE 598 Project #1 Jeremiah Dwight OVERVIEW A simple hot water tank, illustrated in Figures 1 through 3 below, consists of a main cylindrical tank and two small side pipes for the inlet and outlet. All

More information

CHAPTER 4 OPTIMIZATION OF COEFFICIENT OF LIFT, DRAG AND POWER - AN ITERATIVE APPROACH

CHAPTER 4 OPTIMIZATION OF COEFFICIENT OF LIFT, DRAG AND POWER - AN ITERATIVE APPROACH 82 CHAPTER 4 OPTIMIZATION OF COEFFICIENT OF LIFT, DRAG AND POWER - AN ITERATIVE APPROACH The coefficient of lift, drag and power for wind turbine rotor is optimized using an iterative approach. The coefficient

More information

Project #1 Internal flow with thermal convection

Project #1 Internal flow with thermal convection Project #1 Internal flow with thermal convection MAE 494/598, Fall 2017, Project 1 (20 points) Hard copy of report is due at the start of class on the due date. The rules on collaboration will be released

More information

Numerical Study of Convective Heat Transfer for Flat Unglazed Transpired Solar Collectors

Numerical Study of Convective Heat Transfer for Flat Unglazed Transpired Solar Collectors Purdue University Purdue e-pubs International High Performance Buildings Conference School of Mechanical Engineering 2012 Numerical Study of Convective Heat Transfer for Flat Unglazed Transpired Solar

More information

ABSTRACT I. INTRODUCTION

ABSTRACT I. INTRODUCTION 2016 IJSRSET Volume 2 Issue 4 Print ISSN : 2395-1990 Online ISSN : 2394-4099 Themed Section: Engineering and Technology Analysis of Compressible Effect in the Flow Metering By Orifice Plate Using Prasanna

More information

Investigation of Jet Impingement on Flat Plate Using Triangular and Trapezoid Vortex Generators

Investigation of Jet Impingement on Flat Plate Using Triangular and Trapezoid Vortex Generators ISSN 2395-1621 Investigation of Jet Impingement on Flat Plate Using Triangular and Trapezoid Vortex Generators #1 Sonali S Nagawade, #2 Prof. S Y Bhosale, #3 Prof. N K Chougule 1 Sonalinagawade1@gmail.com

More information

Numerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel

Numerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel Numerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel Arunanshu Chakravarty 1* 1 CTU in Prague, Faculty of Mechanical Engineering, Department of Process Engineering,Technická

More information

Laminar flow heat transfer studies in a twisted square duct for constant wall heat flux boundary condition

Laminar flow heat transfer studies in a twisted square duct for constant wall heat flux boundary condition Sādhanā Vol. 40, Part 2, April 2015, pp. 467 485. c Indian Academy of Sciences Laminar flow heat transfer studies in a twisted square duct for constant wall heat flux boundary condition RAMBIR BHADOURIYA,

More information

Numerical Investigation on Turbulent Forced Convection in Heating Channel Inserted with Discrete V-Shaped Baffles

Numerical Investigation on Turbulent Forced Convection in Heating Channel Inserted with Discrete V-Shaped Baffles Journal of Mathematics and Statistics Original Research Paper Numerical Investigation on Turbulent Forced Convection in Heating Channel Inserted with Discrete V-Shaped Baffles 1 Amnart Boonloi and 2 Withada

More information

External Forced Convection. Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.

External Forced Convection. Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display. External Forced Convection Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display. Drag and Heat Transfer in External flow Fluid flow over solid bodies is responsible

More information

Effective Heat Transfer Enhancement in Finned Tube Heat Exchanger with Different Fin Profiles

Effective Heat Transfer Enhancement in Finned Tube Heat Exchanger with Different Fin Profiles International Journal of Engineering Research (ISSN : 2319-6890) Volume No.2, Issue No.2, pp : 83-87 01 April 2013 Effective Heat Transfer Enhancement in Finned Tube Heat Exchanger with Different Fin Profiles

More information

FLOW SIMULATION AND HEAT TRANSFER ANALYSIS USING COMPUTATIONAL FLUID DYNAMICS (CFD)

FLOW SIMULATION AND HEAT TRANSFER ANALYSIS USING COMPUTATIONAL FLUID DYNAMICS (CFD) PERIODIC FLOW SIMULATION AND HEAT TRANSFER ANALYSIS USING COMPUTATIONAL FLUID DYNAMICS (CFD) *M.SNEHA PRIYA, **G. JAMUNA RANI * (Department of Mechanical Engineering, koneru lakshmai University, Guntur)

More information

CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE

CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral

More information

University of Maiduguri Faculty of Engineering Seminar Series Volume 6, december Seminar Series Volume 6, 2015 Page 58

University of Maiduguri Faculty of Engineering Seminar Series Volume 6, december Seminar Series Volume 6, 2015 Page 58 University of Maiduguri Faculty of Engineering Seminar Series Volume 6, december 2015 IMPINGEMENT JET COOLING OF GAS TURBINE COMBUSTOR WALL OF HEAT FLUX IMPOSED HOT - SIDE: CONJUGATE HEAT TRANSFER INVESTIGATIONS

More information

Convection. forced convection when the flow is caused by external means, such as by a fan, a pump, or atmospheric winds.

Convection. forced convection when the flow is caused by external means, such as by a fan, a pump, or atmospheric winds. Convection The convection heat transfer mode is comprised of two mechanisms. In addition to energy transfer due to random molecular motion (diffusion), energy is also transferred by the bulk, or macroscopic,

More information

Lab 5: Post Processing and Solving Conduction Problems. Objective:

Lab 5: Post Processing and Solving Conduction Problems. Objective: Lab 5: Post Processing and Solving Conduction Problems Objective: The objective of this lab is to use the tools we have developed in MATLAB and SolidWorks to solve conduction heat transfer problems that

More information

Simulation and improvement of the ventilation of a welding workshop using a Finite volume scheme code

Simulation and improvement of the ventilation of a welding workshop using a Finite volume scheme code 1 st. Annual (National) Conference on Industrial Ventilation-IVC2010 Feb 24-25, 2010, Sharif University of Technology, Tehran, Iran IVC2010 Simulation and improvement of the ventilation of a welding workshop

More information

THERMO-FLOW CHARACTERISTICS OF A PIN-FIN RADIAL HEAT SINKS ACCORDING TO THEIR FIN HEIGHT PROFILE

THERMO-FLOW CHARACTERISTICS OF A PIN-FIN RADIAL HEAT SINKS ACCORDING TO THEIR FIN HEIGHT PROFILE HEFAT2012 9 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 16 18 July 2012 Malta THERMO-FLOW CHARACTERISTICS OF A PIN-FIN RADIAL HEAT SINKS ACCORDING TO THEIR FIN HEIGHT

More information

AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract

AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract nd Workshop on CFD Uncertainty Analysis - Lisbon, 19th and 0th October 006 AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION Alfredo Iranzo 1, Jesús Valle, Ignacio Trejo 3, Jerónimo

More information

The effect of geometric parameters on the head loss factor in headers

The effect of geometric parameters on the head loss factor in headers Fluid Structure Interaction V 355 The effect of geometric parameters on the head loss factor in headers A. Mansourpour & S. Shayamehr Mechanical Engineering Department, Azad University of Karaj, Iran Abstract

More information

A numerical study of heat transfer and fluid flow over an in-line tube bank

A numerical study of heat transfer and fluid flow over an in-line tube bank Fluid Structure Interaction VI 295 A numerical study of heat transfer and fluid flow over an in-line tube bank Z. S. Abdel-Rehim Mechanical Engineering Department, National Research Center, Egypt Abstract

More information

Numerical Investigation of Convective Heat Transfer in Pin Fin Type Heat Sink used for Led Application by using CFD

Numerical Investigation of Convective Heat Transfer in Pin Fin Type Heat Sink used for Led Application by using CFD GRD Journals- Global Research and Development Journal for Engineering Volume 1 Issue 8 July 2016 ISSN: 2455-5703 Numerical Investigation of Convective Heat Transfer in Pin Fin Type Heat Sink used for Led

More information

NUMERICAL HEAT TRANSFER ENHANCEMENT IN SQUARE DUCT WITH INTERNAL RIB

NUMERICAL HEAT TRANSFER ENHANCEMENT IN SQUARE DUCT WITH INTERNAL RIB NUMERICAL HEAT TRANSFER ENHANCEMENT IN SQUARE DUCT WITH INTERNAL RIB University of Technology Department Mechanical engineering Baghdad, Iraq ABSTRACT - This paper presents numerical investigation of heat

More information

CFD Study of Flow Over Parallel Ridges with Varying Height and Spacing

CFD Study of Flow Over Parallel Ridges with Varying Height and Spacing Proceedings of the World Congress on Engineering 21 Vol II WCE 21, June 3 - July 2, 21, London, U.K. CFD Study of Flow Over Parallel Ridges with Varying Height and Spacing Lee Chin Yik, Salim Mohamed Salim,

More information

CFD in Heat Transfer Equipment Professor Bengt Sunden Division of Heat Transfer Department of Energy Sciences Lund University

CFD in Heat Transfer Equipment Professor Bengt Sunden Division of Heat Transfer Department of Energy Sciences Lund University CFD in Heat Transfer Equipment Professor Bengt Sunden Division of Heat Transfer Department of Energy Sciences Lund University email: bengt.sunden@energy.lth.se CFD? CFD = Computational Fluid Dynamics;

More information

ENGR Heat Transfer II

ENGR Heat Transfer II ENGR 7901 - Heat Transfer II External Flows 1 Introduction In this chapter we will consider several fundamental flows, namely: the flat plate, the cylinder, the sphere, several other body shapes, and banks

More information

LES ANALYSIS ON CYLINDER CASCADE FLOW BASED ON ENERGY RATIO COEFFICIENT

LES ANALYSIS ON CYLINDER CASCADE FLOW BASED ON ENERGY RATIO COEFFICIENT 2th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics ANALYSIS ON CYLINDER CASCADE FLOW BASED ON ENERGY RATIO COEFFICIENT Wang T.*, Gao S.F., Liu Y.W., Lu Z.H. and Hu H.P. *Author

More information

Enhancement of Heat Transfer Effectiveness of Plate-pin fin heat sinks With Central hole and Staggered positioning of Pin fins

Enhancement of Heat Transfer Effectiveness of Plate-pin fin heat sinks With Central hole and Staggered positioning of Pin fins Enhancement of Heat Transfer Effectiveness of Plate-pin fin heat sinks With Central hole and Staggered positioning of Pin fins Jubin Jose 1, Reji Mathew 2 1Student, Dept. of Mechanical Engineering, M A

More information

PHYSICAL MECHANISM OF CONVECTION

PHYSICAL MECHANISM OF CONVECTION Tue 8:54:24 AM Slide Nr. 0 of 33 Slides PHYSICAL MECHANISM OF CONVECTION Heat transfer through a fluid is by convection in the presence of bulk fluid motion and by conduction in the absence of it. Chapter

More information

UNIT II CONVECTION HEAT TRANSFER

UNIT II CONVECTION HEAT TRANSFER UNIT II CONVECTION HEAT TRANSFER Convection is the mode of heat transfer between a surface and a fluid moving over it. The energy transfer in convection is predominately due to the bulk motion of the fluid

More information

Heat Transfer Simulation by CFD from Fins of an Air Cooled Motorcycle Engine under Varying Climatic Conditions

Heat Transfer Simulation by CFD from Fins of an Air Cooled Motorcycle Engine under Varying Climatic Conditions , July 6-8, 2011, London, U.K. Heat Transfer Simulation by CFD from Fins of an Air Cooled Motorcycle Engine under Varying Climatic Conditions Pulkit Agarwal, Mayur Shrikhande and P. Srinivasan Abstract

More information

Due Tuesday, November 23 nd, 12:00 midnight

Due Tuesday, November 23 nd, 12:00 midnight Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem

More information

Numerical Study of the Three-dimensional Flow in a Flat Plate Solar Collector with Baffles

Numerical Study of the Three-dimensional Flow in a Flat Plate Solar Collector with Baffles Numerical Study of the Three-dimensional Flow in a Flat Plate Solar Collector with Baffles M.A. AMRAOUI a, K. ALIANE b a. Department of Mechanical Engineering, Faculty of Technology, University of Tlemcen,

More information

HEFAT th International Conference on Heat Transfer, Fluid Mechanics, and Thermodynamics September 2005, Cairo, Egypt AA10

HEFAT th International Conference on Heat Transfer, Fluid Mechanics, and Thermodynamics September 2005, Cairo, Egypt AA10 HEFAT5 4 th International Conference on Heat Transfer, Fluid Mechanics, and Thermodynamics 9- September 5, Cairo, Egypt AA Numerical Study of Natural Convection Heat Transfer in Enclosures with Conducting

More information

Chapter 7: External Forced Convection. Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University

Chapter 7: External Forced Convection. Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University Chapter 7: External Forced Convection Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University Objectives When you finish studying this chapter, you should be able to: Distinguish between

More information

PERFORMANCE SCREENING OF A LOUVERED FIN AND VORTEX GENERATOR COMBINATION

PERFORMANCE SCREENING OF A LOUVERED FIN AND VORTEX GENERATOR COMBINATION HEFAT2014 10 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 14 26 July 2014 Orlando, Florida PERFORMANCE SCREENING OF A LOUVERED FIN AND VORTEX GENERATOR COMBINATION Bernd

More information

Explicit algebraic Reynolds stress models for internal flows

Explicit algebraic Reynolds stress models for internal flows 5. Double Circular Arc (DCA) cascade blade flow, problem statement The second test case deals with a DCA compressor cascade, which is considered a severe challenge for the CFD codes, due to the presence

More information

SIMULATION OF THERMAL CHARACTERISTICS OF RADIATORS USING A POROUS MODEL. YETSAN Auto Radiator Co. Inc Çorum, Turkey NOMENCLATURE

SIMULATION OF THERMAL CHARACTERISTICS OF RADIATORS USING A POROUS MODEL. YETSAN Auto Radiator Co. Inc Çorum, Turkey NOMENCLATURE Proceedings of CONV-14: Int. Symp. on ConvectiveHeatandMass Transfer June8 13, 2014, Turkey CONV-14 176 SIMULATION OF THERMAL CHARACTERISTICS OF RADIATORS USING A POROUS MODEL Kadir G. Güler 1,2 and BarbarosÇetin

More information

COMPUTATIONAL ANALYSIS OF LAMINAR FORCED CONVECTION IN RECTANGULAR ENCLOSURES OF DIFFERENT ASPECT RATIOS

COMPUTATIONAL ANALYSIS OF LAMINAR FORCED CONVECTION IN RECTANGULAR ENCLOSURES OF DIFFERENT ASPECT RATIOS HEFAT214 1 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 14 16 July 214 Orlando, Florida COMPUTATIONAL ANALYSIS OF LAMINAR FORCED CONVECTION IN RECTANGULAR ENCLOSURES

More information

Numerical Investigation of Multijet Air Impingement on Pin Fin Heat Sink with Effusion Slots

Numerical Investigation of Multijet Air Impingement on Pin Fin Heat Sink with Effusion Slots , 23-25 October, 2013, San Francisco, USA Numerical Investigation of Multijet Air Impingement on Pin Fin Heat Sink with Effusion Slots N. K. Chougule G. V. Parishwad A. R. Nadgire Abstract The work reported

More information

Numerical Simulation of the Hagemann Entrainment Experiments

Numerical Simulation of the Hagemann Entrainment Experiments CCC Annual Report UIUC, August 14, 2013 Numerical Simulation of the Hagemann Entrainment Experiments Kenneth Swartz (BSME Student) Lance C. Hibbeler (Ph.D. Student) Department of Mechanical Science & Engineering

More information

INTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 2, No 1, 2011

INTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 2, No 1, 2011 Experimental and Numerical comparison between the performance of Helical cone coils and ordinary helical coils used as dehumidifier for humidification dehumidification in desalination units Abo Elazm M.M.

More information

Heat Transfer Enhancement in Turbines using Dimples- A CFD Study Nitish Kumar Jain 1 Dr.U.S. Mallikarjun 2

Heat Transfer Enhancement in Turbines using Dimples- A CFD Study Nitish Kumar Jain 1 Dr.U.S. Mallikarjun 2 IJSRD - International Journal for Scientific Research & Development Vol. 3, Issue 05, 2015 ISSN (online): 2321-0613 Heat Transfer Enhancement in Turbines using Dimples- A CFD Study Nitish Kumar Jain 1

More information

FLOW DISTRIBUTION IN A HYDRAULIC HEADBOX

FLOW DISTRIBUTION IN A HYDRAULIC HEADBOX FLOW DISTRIBUTION IN A HYDRAULIC HEADBOX Lu Hua, Eric Bibeau _, Pingfan He, Martha Salcudean and Ian Gartshore Department of Mechanical Engineering, The University of British Columbia, Vancouver, BC V6T

More information

A Computational Fluid Dynamics Investigation of Solar Air Heater Duct Provided with Inclined Circular Ribs as Artificial Roughness

A Computational Fluid Dynamics Investigation of Solar Air Heater Duct Provided with Inclined Circular Ribs as Artificial Roughness Bonfring International Journal of Industrial Engineering and Management Science, Vol. 4, No. 3, August 2014 115 A Computational Fluid Dynamics Investigation of Solar Air Heater Duct Provided with Inclined

More information

ME 144: Heat Transfer Introduction to Convection. J. M. Meyers

ME 144: Heat Transfer Introduction to Convection. J. M. Meyers ME 144: Heat Transfer Introduction to Convection Introductory Remarks Convection heat transfer differs from diffusion heat transfer in that a bulk fluid motion is present which augments the overall heat

More information

SHELL SIDE NUMERICAL ANALYSIS OF A SHELL AND TUBE HEAT EXCHANGER CONSIDERING THE EFFECTS OF BAFFLE INCLINATION ANGLE ON FLUID FLOW

SHELL SIDE NUMERICAL ANALYSIS OF A SHELL AND TUBE HEAT EXCHANGER CONSIDERING THE EFFECTS OF BAFFLE INCLINATION ANGLE ON FLUID FLOW THERMAL SCIENCE: Year 2012, Vol. 16, No. 4, pp. 1165-1174 1165 SHELL SIDE NUMERICAL ANALYSIS OF A SHELL AND TUBE HEAT EXCHANGER CONSIDERING THE EFFECTS OF BAFFLE INCLINATION ANGLE ON FLUID FLOW by Rajagapal

More information

Heat Transfer F12-ENG Lab #4 Forced convection School of Engineering, UC Merced.

Heat Transfer F12-ENG Lab #4 Forced convection School of Engineering, UC Merced. 1 Heat Transfer F12-ENG-135 - Lab #4 Forced convection School of Engineering, UC Merced. October 23, 2012 1 General purpose of the Laboratory To gain a physical understanding of the behavior of the average

More information

Numerical Analysis of Tube-Fin Heat Exchanger using Fluent

Numerical Analysis of Tube-Fin Heat Exchanger using Fluent Numerical Analysis of Tube-Fin Heat Exchanger using Fluent M. V. Ghori & R. K. Kirar Patel college of Science and Technology, Rajiv Gandhi Proudyogiki Vishwavidhyalaya, Ratibad, Bhopal- 462036, India E-mail

More information

Applied CFD Project 1. Christopher Light MAE 598

Applied CFD Project 1. Christopher Light MAE 598 Applied CFD Project 1 Christopher Light MAE 598 October 5, 2017 Task 1 The hot water tank shown in Fig 1 is used for analysis of cool water flow with the heat from a hot plate at the bottom. For all tasks,

More information

HEAT TRANSFER CAPABILITY OF A THERMOSYPHON HEAT TRANSPORT DEVICE WITH EXPERIMENTAL AND CFD STUDIES

HEAT TRANSFER CAPABILITY OF A THERMOSYPHON HEAT TRANSPORT DEVICE WITH EXPERIMENTAL AND CFD STUDIES HEAT TRANSFER CAPABILITY OF A THERMOSYPHON HEAT TRANSPORT DEVICE WITH EXPERIMENTAL AND CFD STUDIES B.M. Lingade a*, Elizabeth Raju b, A Borgohain a, N.K. Maheshwari a, P.K.Vijayan a a Reactor Engineering

More information

PERIODICALLY FULLY DEVELOPED LAMINAR FLOW AND HEAT TRANSFER IN A 2-D HORIZONTAL CHANNEL WITH STAGGERED FINS

PERIODICALLY FULLY DEVELOPED LAMINAR FLOW AND HEAT TRANSFER IN A 2-D HORIZONTAL CHANNEL WITH STAGGERED FINS THERMAL SCIENCE: Year 2017, Vol. 21, No. 6A, pp. 2443-2455 2443 PERIODICALLY FULLY DEVELOPED LAMINAR FLOW AND HEAT TRANSFER IN A 2-D HORIZONTAL CHANNEL WITH STAGGERED FINS by Oguz TURGUT a* and Kamil ARSLAN

More information

DNS of Buoyancy Driven Flow Inside a Horizontal Coaxial Cylinder

DNS of Buoyancy Driven Flow Inside a Horizontal Coaxial Cylinder DNS of Buoyancy Driven Flow Inside a Horizontal Coaxial Cylinder Imama Zaidi 1, Yacine Addad 2, Dominique Laurence 1 1 The University of Manchester, School of Mechanical, Aerospace and Civil Eng., M60

More information

IJSER. PRESENT study considers numerical investigation for 1 INTRODUCTION. Vivek Shrivastava, Pavan Badami, Saravanan V, K N Seetharamu

IJSER. PRESENT study considers numerical investigation for 1 INTRODUCTION. Vivek Shrivastava, Pavan Badami, Saravanan V, K N Seetharamu International Journal of Scientific & Engineering Research, Volume 5, Issue 4, April-2014 554 Effect of Triangular Wake Splitter on Flow and Heat Transfer over a Circular Cylinder for Various Chord Lengths

More information

Computational and Experimental Studies of Fluid flow and Heat Transfer in a Calandria Based Reactor

Computational and Experimental Studies of Fluid flow and Heat Transfer in a Calandria Based Reactor Computational and Experimental Studies of Fluid flow and Heat Transfer in a Calandria Based Reactor SD Ravi 1, NKS Rajan 2 and PS Kulkarni 3 1 Dept. of Aerospace Engg., IISc, Bangalore, India. ravi@cgpl.iisc.ernet.in

More information

Thermo-Hydraulic performance of Internal finned tube Automobile Radiator

Thermo-Hydraulic performance of Internal finned tube Automobile Radiator Thermo-Hydraulic performance of Internal finned tube Automobile Radiator Dr.Kailash Mohapatra 1, Deepiarani Swain 2 1 Department of Mechanical Engineering, Raajdhani Engineering College, Bhubaneswar, 751017,

More information

Heat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay

Heat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay Heat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay Lecture No. 18 Forced Convection-1 Welcome. We now begin our study of forced convection

More information

JJMIE Jordan Journal of Mechanical and Industrial Engineering

JJMIE Jordan Journal of Mechanical and Industrial Engineering JJMIE Jordan Journal of Mechanical and Industrial Engineering Volume Number, June.6 ISSN 995-6665 Pages 99-4 Computational Fluid Dynamics of Plate Fin and Circular Pin Fin Heat Sins Mohammad Saraireh *

More information

EFFECT OF ROTOR S ANGULAR SPEED AND ASPECT RATIO ON ENTROPY GENERATION IN A ROTOR-CASING ASSEMBLY

EFFECT OF ROTOR S ANGULAR SPEED AND ASPECT RATIO ON ENTROPY GENERATION IN A ROTOR-CASING ASSEMBLY Third International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 1-1 December 3 EFFECT OF ROTOR S ANGUL SPEED AND ASPECT RATIO ON ENTROPY GENERATION IN A ROTOR-CASING

More information

Modeling of Fluid Flow and Heat Transfer for Optimization of Pin-Fin Heat Sinks

Modeling of Fluid Flow and Heat Transfer for Optimization of Pin-Fin Heat Sinks Modeling of Fluid Flow and Heat Transfer for Optimization of Pin-Fin Heat Sinks by Waqar Ahmed Khan Athesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree

More information

CFD ANALYSIS OF CD NOZZLE AND EFFECT OF NOZZLE PRESSURE RATIO ON PRESSURE AND VELOCITY FOR SUDDENLY EXPANDED FLOWS. Kuala Lumpur, Malaysia

CFD ANALYSIS OF CD NOZZLE AND EFFECT OF NOZZLE PRESSURE RATIO ON PRESSURE AND VELOCITY FOR SUDDENLY EXPANDED FLOWS. Kuala Lumpur, Malaysia International Journal of Mechanical and Production Engineering Research and Development (IJMPERD) ISSN(P): 2249-6890; ISSN(E): 2249-8001 Vol. 8, Issue 3, Jun 2018, 1147-1158 TJPRC Pvt. Ltd. CFD ANALYSIS

More information

Numerical Analysis of a Multi-Row Multi-Column Compact Heat Exchanger

Numerical Analysis of a Multi-Row Multi-Column Compact Heat Exchanger Numerical Analysis of a Multi-Row Multi-Column Compact Heat Exchanger Ashkan Motamedi 1, Arturo Pacheco-Vega 1,, J. Rafael Pacheco 2 1 Department of Mechanical Engineering, California State University,

More information

CFD SIMULATION OF A SINGLE PHASE FLOW IN A PIPE SEPARATOR USING REYNOLDS STRESS METHOD

CFD SIMULATION OF A SINGLE PHASE FLOW IN A PIPE SEPARATOR USING REYNOLDS STRESS METHOD CFD SIMULATION OF A SINGLE PHASE FLOW IN A PIPE SEPARATOR USING REYNOLDS STRESS METHOD Eyitayo A. Afolabi 1 and J. G. M. Lee 2 1 Department of Chemical Engineering, Federal University of Technology, Minna,

More information

SELF-SUSTAINED OSCILLATIONS AND BIFURCATIONS OF MIXED CONVECTION IN A MULTIPLE VENTILATED ENCLOSURE

SELF-SUSTAINED OSCILLATIONS AND BIFURCATIONS OF MIXED CONVECTION IN A MULTIPLE VENTILATED ENCLOSURE Computational Thermal Sciences, 3 (1): 63 72 (2011) SELF-SUSTAINED OSCILLATIONS AND BIFURCATIONS OF MIXED CONVECTION IN A MULTIPLE VENTILATED ENCLOSURE M. Zhao, 1, M. Yang, 1 M. Lu, 1 & Y. W. Zhang 2 1

More information

INTERNATIONAL JOURNAL OF SCIENTIFIC & TECHNOLOGY RESEARCH VOLUME 5, ISSUE 09, SEPTEMBER 2016 ISSN

INTERNATIONAL JOURNAL OF SCIENTIFIC & TECHNOLOGY RESEARCH VOLUME 5, ISSUE 09, SEPTEMBER 2016 ISSN Numerical Analysis Of Heat Transfer And Fluid Flow Characteristics In Different V-Shaped Roughness Elements On The Absorber Plate Of Solar Air Heater Duct Jitesh Rana, Anshuman Silori, Rohan Ramola Abstract:

More information

Open boundary conditions in numerical simulations of unsteady incompressible flow

Open boundary conditions in numerical simulations of unsteady incompressible flow Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can

More information

This chapter focuses on the study of the numerical approximation of threedimensional

This chapter focuses on the study of the numerical approximation of threedimensional 6 CHAPTER 6: NUMERICAL OPTIMISATION OF CONJUGATE HEAT TRANSFER IN COOLING CHANNELS WITH DIFFERENT CROSS-SECTIONAL SHAPES 3, 4 6.1. INTRODUCTION This chapter focuses on the study of the numerical approximation

More information

Drag Coefficient and Nusselt Number for Laminar Pulsating Flow in Porous Media

Drag Coefficient and Nusselt Number for Laminar Pulsating Flow in Porous Media Drag Coefficient and Nusselt Number for Laminar Pulsating Flow in Porous Media T.I. Mulcahey, M.G. Pathak, S.M. Ghiaasiaan Cryo Lab, School of Mechanical Engineering Georgia Institute of Technology, Atlanta,

More information

This section develops numerically and analytically the geometric optimisation of

This section develops numerically and analytically the geometric optimisation of 7 CHAPTER 7: MATHEMATICAL OPTIMISATION OF LAMINAR-FORCED CONVECTION HEAT TRANSFER THROUGH A VASCULARISED SOLID WITH COOLING CHANNELS 5 7.1. INTRODUCTION This section develops numerically and analytically

More information

On the influence of tube row number for mixed convection around micro tubes

On the influence of tube row number for mixed convection around micro tubes Thessaloniki, Greece, 22-24 August 211 On the influence of tube row number for mixed convection around micro tubes Chuanshan DAI, Qiuxiang WANG, Biao LI * Corresponding author: Tel.: +86-22274183; Fax:

More information

Is My CFD Mesh Adequate? A Quantitative Answer

Is My CFD Mesh Adequate? A Quantitative Answer Is My CFD Mesh Adequate? A Quantitative Answer Krzysztof J. Fidkowski Gas Dynamics Research Colloqium Aerospace Engineering Department University of Michigan January 26, 2011 K.J. Fidkowski (UM) GDRC 2011

More information

Department of Mechanical Engineering

Department of Mechanical Engineering Department of Mechanical Engineering AMEE401 / AUTO400 Aerodynamics Instructor: Marios M. Fyrillas Email: eng.fm@fit.ac.cy HOMEWORK ASSIGNMENT #2 QUESTION 1 Clearly there are two mechanisms responsible

More information

Theoretical Gas Flow through Gaps in Screw-type Machines

Theoretical Gas Flow through Gaps in Screw-type Machines Theoretical Gas Flow through Gaps in Screw-type Machines Prof. Dr.-Ing. K. Kauder, Dipl.-Ing. D. Stratmann University of Dortmund, Fachgebiet Fluidenergiemaschinen (The experimental part of these studies

More information

Documentation of the Solutions to the SFPE Heat Transfer Verification Cases

Documentation of the Solutions to the SFPE Heat Transfer Verification Cases Documentation of the Solutions to the SFPE Heat Transfer Verification Cases Prepared by a Task Group of the SFPE Standards Making Committee on Predicting the Thermal Performance of Fire Resistive Assemblies

More information

Comparison Of Square-hole And Round-hole Film Cooling: A Computational Study

Comparison Of Square-hole And Round-hole Film Cooling: A Computational Study University of Central Florida Electronic Theses and Dissertations Masters Thesis (Open Access) Comparison Of Square-hole And Round-hole Film Cooling: A Computational Study 2004 Michael Glenn Durham University

More information

TURBULENT FLOW ACROSS A ROTATING CYLINDER WITH SURFACE ROUGHNESS

TURBULENT FLOW ACROSS A ROTATING CYLINDER WITH SURFACE ROUGHNESS HEFAT2014 10 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 14 16 July 2014 Orlando, Florida TURBULENT FLOW ACROSS A ROTATING CYLINDER WITH SURFACE ROUGHNESS Everts, M.,

More information