EFFECT OF THE SIZE OF GRADED BAFFLES ON THE PERFORMANCE OF CHANNEL HEAT EXCHANGERS

Similar documents
A Numerical Study of Heat Transfer and Fluid Flow past Single Tube

Numerical Analysis of Heat Transfer and Pressure Drop in a Channel Equipped with Triangular Bodies in Side-By-Side Arrangement

FORCED CONVECTION HEAT TRANSFER FROM A RECTANGULAR CYLINDER: EFFECT OF ASPECT RATIO

Numerical Heat and Mass Transfer

A Comparative Investigation into Aerodynamic Performances of Two Set Finned Bodies with Circular and Non Circular Cross Sections

UNSTEADY LAMINAR CONVECTION FLOW OVER PERIODIC GROOVES BY USING SIO 2 -WATER NANOFLUID

OPTIMIZATION ANALYSIS FOR LOUVER-FINNED HEAT EXCHANGERS

Research & Reviews: Journal of Engineering and Technology

THE NEAR-WALL INFLUENCE ON THE FLOW AROUND A SINGLE SQUARE CYLINDER.

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

TURBULENT WALL JET OVER A FORWARD-BACKWARD FACING STEP PAIR

Turbulent Flow. Turbulent Flow

modeling of equilibrium and dynamic multi-component adsorption in a two-layered fixed bed for purification of hydrogen from methane reforming products

Heat Transfer and Turbulent Nanofluid Flow over a Double Forward- Facing Step

Simulation and experiment of the effect of clearance of impeller wear-rings on the performance of centrifugal pump

NUMERICAL SIMULATION OF FLOW OVER STEPPED SPILLWAYS

Turbulent Flow in Curved Square Duct: Prediction of Fluid flow and Heat transfer Characteristics

PERMEABILITY-POROSITY RELATIONSHIP ASSESSMENT BY 2-D NUMERICAL SIMULATIONS

A large scale tsunami run-up simulation and numerical evaluation of fluid force during tsunami by using a particle method

Calculation of Aerodynamic Characteristics of NACA 2415, 23012, Airfoils Using Computational Fluid Dynamics (CFD)

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests

A Hybrid Variational Iteration Method for Blasius Equation

A NUMERICAL STUDY OF HEAT TRANSFER AND FLUID FLOW IN A BANK OF TUBES WITH INTEGRAL WAKE SPLITTER

Simulation of Flow Pattern in Open Channels with Sudden Expansions

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

A numerical study on the shell-side turbulent heat transfer enhancement of shell-and-tube heat exchanger with trefoilhole

HEAT TRANSFER THROUGH ANNULAR COMPOSITE FINS

2 MODELS A typcal hgh-sded artculated lorry, whch was nvestgated extensvely by Baker and hs colleagues n wnd tunnels and later used n the dynamc analy

Numerical Analysis and Optimization on Vortex Shedder by the Linearity of Strouhal Number against Reynold Number

EXAMPLES of THEORETICAL PROBLEMS in the COURSE MMV031 HEAT TRANSFER, version 2017

Energy configuration optimization of submerged propeller in oxidation ditch based on CFD

Annexes. EC.1. Cycle-base move illustration. EC.2. Problem Instances

Flow equations To simulate the flow, the Navier-Stokes system that includes continuity and momentum equations is solved

Principles of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer

DUE: WEDS FEB 21ST 2018

2 Finite difference basics

INTERROGATING THE FLOW BEHAVIOUR IN A NOVEL MAGNETIC DESICCANT VENTILATION SYSTEM USING COMPUTATIONAL FLUID DYNAMICS (CFD)

The Finite Element Method

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

One-sided finite-difference approximations suitable for use with Richardson extrapolation

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

Global Sensitivity. Tuesday 20 th February, 2018

Design optimization of a high specific speed Francis turbine runner

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

GeoSteamNet: 2. STEAM FLOW SIMULATION IN A PIPELINE

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

Computational investigation of the external excitation frequency effect on liquid sloshing phenomenon

NUMERICAL ANALYSIS OF TURBULENT FLOW WITH HEAT TRANSFER IN A SQUARE DUCT WITH 45 DEGREE RIBS

Numerical simulation of water flow in an axial flow pump with adjustable guide vanes

Lattice Boltzmann simulation of nucleate boiling in micro-pillar structured surface

Numerical Study of Propane-Air Mixture Combustion in a Burner Element

An Experimental and Numerical Study on Pressure Drop Coefficient of Ball Valves

NUMERICAL SIMULATION OF HEAT TRANSFER CHARACTERISTICS IN THE ABSORBER TUBE OF PARABOLIC TROUGH COLLECTOR WITH INTERNAL FLOW OBSTRUCTIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018

The unsteady flow characteristic research on the initial period flow of micro channel

CFD VALIDATION OF STRATIFIED TWO-PHASE FLOWS IN A HORIZONTAL CHANNEL

Investigation of a New Monte Carlo Method for the Transitional Gas Flow

Nomenclature. I. Introduction

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

Numerical simulation of flow past circular duct

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Constitutive Modelling of Superplastic AA-5083

FEATURES OF TURBULENT TRANSPORT OF MOMENTUM AND HEAT IN STABLY STRATIFIED BOUNDARY LAYERS AND THEIR REPRODUCTION IN ATMOSPHERIC MESOSCALE MODELS

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)

Basic concept of reactive flows. Basic concept of reactive flows Combustion Mixing and reaction in high viscous fluid Application of Chaos

Introduction to Computational Fluid Dynamics

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

Inductance Calculation for Conductors of Arbitrary Shape

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient

829. An adaptive method for inertia force identification in cantilever under moving mass

A Cartesian-grid integrated-rbf method for viscoelastic flows

Polynomial Regression Models

OFF-AXIS MECHANICAL PROPERTIES OF FRP COMPOSITES

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

Handout: Large Eddy Simulation I. Introduction to Subgrid-Scale (SGS) Models

STUDY OF A THREE-AXIS PIEZORESISTIVE ACCELEROMETER WITH UNIFORM AXIAL SENSITIVITIES

A NEW FILTERED DYNAMIC SUBGRID-SCALE MODEL FOR LARGE EDDY SIMULATION OF INDOOR AIRFLOW

Developing of laminar fluid flow in rectangular microchannels

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

3-D NUMERICAL SIMULATION OF FLOW AND CLEAR WATER SCOUR BY INTERACTION BETWEEN BRIDGE PIERS

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)

CHAPTER 9 CONCLUSIONS

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

Simulation of Turbulent Flow Using FEM

Numerical and Experimental Analysis of Elliptic Finned-Tube Heat Exchangers under Misted Conditions

Statistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA

IC Engine Flow Simulation using KIVA code and A Modified Reynolds Stress Turbulence Model

Numerical Studies on Flow Features past a Backward Facing Sharp Edge Step Introducing Hybrid RANS-LES

2010 Black Engineering Building, Department of Mechanical Engineering. Iowa State University, Ames, IA, 50011

Conduction Shape Factor Models for Three-Dimensional Enclosures

Application of the Adjoint Method for Vehicle Aerodynamic Optimization. Dr. Thomas Blacha, Audi AG

Study on Effect of Different Porosity on Thermally Driven Heat Transfer in a Centrifugal Force Field

Survey of applications of discrete vortex method in civil engineering

Analysis of Unsteady Aerodynamics of a Car Model with Radiator in Dynamic Pitching Motion using LS-DYNA

Transcription:

EFFECT OF THE SIZE OF GRADED BAFFLES ON THE PERFORMANCE OF CHANNEL HEAT EXCHANGERS Djamel SAHEL 1, Houar AMEUR 2,*, and Bak TOUHAMI 3 1 Techncal Scences Department, Unversty of Amar Teldj, Laghouat, Algera *2 Insttute of Scence and Technology, Unversty Center of Naama (Ctr Unv Naama), 45000, Algera 3 Mechancal Engneerng Department, Unversty of scence and the technology, Oran, Algera * Correspondng author; E-mal: houar_ameur@yahoo.fr; Tel: +213770343722 The bafflng technque s well-known for ts effcency n terms of enhancement of heat transfer rates throught channels. However, the baffles nsert s accompaned by an ncrease n the frcton factor. Ths ssue remans a great challenge for the desgners of heat exchangers. To overcome ths ssue, we suggest n the present paper a new desgn of baffles whch s here called graded baffle-desgn. The baffles have an up- or downgraded heght along the channel length. Ths geometry s characterzed by two ratos: up-graded baffle rato and down-graded baffle rato whch are vared from 0 to 0.08. For a range of Reynolds number varyng from 10 4 to 2x10 4, the turbulent flow and heat transfer characterstcs of a heat exchanger channel are numercally studed by the computer code Fluent. The obtaned results revealed an enhancement n the thermo-hydraulc performance offered by the new suggested desgn. For the channel wth a down-graded baffle rato equal to 0.08, the frcton factors decreased by 4-8%. Key words: Heat exchanger channel; Graded baffles; Turbulent flow; Frcton factor; CFD. 1. Introducton Several technques are developed to reduce the energy cost and to enhance the thermal performance of varous devces used n many nductral applcatons. In ths framework, rbs,baffles or obstacles play an mportant role n for the mprovement of the thermal transfer rate n varous thermal systems lke the heat exchangers, the gas turbnes blades, the solar collectors, the car radators and other applcatons of fluds mechancs such as the water desalnaton. Some expermental and numercal works have been acheved on the effect of dfferent geometrcal parameters of baffles as the shape, sze and spacng between baffles, the attack angle, the blockage rato and the porous space. The shape of baffles s a necessary parameter for the generaton of vortces [1-3]. Therefore, varous shapes have been developed, for example, the V-baffle shape [4-1

9], the W-baffle shape [10], the Z-baffle shape [11], mult V-type perforated baffles [12-17] and damond shape proposed by Srpattanappat et al. [18]. All these shapes permt an enhancement n the heat transfer rates, but wth catastrophc pressure drops. Dutta et al. [19] studed by experments the effct of baffle attack angle on the characterstcs of heat transfer and pressure drop n a rectangular channel wth nclned plate and perforated baffles. They showed the man dependency of the thermal transfer rates on the geometry, the orentaton and the poston of the second baffle. For a square channel Promvonge et al. [20] explored the effect of baffle nclnaton wth varous values of blockage ratos (BR). Wth the 45 nclned baffles and BR = 0-0.05, they found an ncrease n the thermal transfer rate from 150 to 850% accompaned wth pressure losses varyng from 2 to 70 tmes compared wth the smooth channel. In other works [21, 22], these authors found that the combnaton between the attack angle and the V-baffle shape allow an enhancement n heat transfer rates assocated wth losses n pressure. In the same framwork, Promvonge et al. [23] studed by experments the effects of the delta baffles shapes on the thermal performance of solar ar heater. They reported that the combnaton between the delta wnglet and the rb ensures consderable heat transfer mprovements (Nu/Nu 0 = 2.3 2.6) and also yelds a moderate pressure drop ncrease (f/f 0 = 4.7 10.1). Tandroglu et al. [24] explored n ther expermental study three cases of nclned baffles (45, 90 and 180 ) wth three values of dameter rato H/D = 1, 2 and 3 and they obtaned the best performance wth the 90 case for the forced convecton. In another work, Tandroglu [25] developed emprcal formulas for the calculaton of Nusselt numbers and the frcton coeffcents. The same author, tested n another artcle the entropy generaton for dfferent types of baffles [26]. The V-shaped baffle or rb s another technque to enhance the thermal performance due to the presence of two hgh thermal regons. Some researcher [27, 28] reported that the staggered arrangement of rbs gves less hgher heat transfer coeffcent than the symmetrc arrangement rbs. Further studes on the enhancement of heat transfer by the bafflng technque can be found elsewhere [29-32]. The effcency of porous baffles s confrmed by many researchers because the perforated space enables the agtaton of the flow stagnated behnd of baffles [33, 34], t ncreases the thermal transfer ntensty and elmnates the lowers heat transfer areas (LHTA) [35-38]. All studes n lterature confrm that the bafflng technque ncreases sgnfcantly the heat transfer rates. However, these works dd not correct the major problems of ths technque: the formaton of LHTA and the pressure losses. Amng to reduce the frcton factor n baffled channels, we propose here new desgned baffles wth up-graded and down-graded heght along the channel length. The performance of the new desgn has been nvestgated for dfferent flow rates. 2. Geometry of the problem studed The purpose of the present work s to study the flow and heat transfer characterstcs n a horzontal channel wth a baffle seres of up- and down-graded heght along the channel length (L). The baffles are nserted n a staggered array on the upper and lower walls of the channel (Fgure 1). Varatons n the baffle heght are defned by DBR (Fgure 1a) and GBR (Fgure 1b) for the decreased (.e. down-graded) and ncreased (.e. up-graded) baffle ratos, respectvely. The obtaned results from the baffled cases are compared wth those for a smooth channel (wthout baffles) havng the followng geometrcal parameters: H s the channel heght fxed to 0.1 m, 2

b s the baffle heght, e s the baffle thckness (e = 0.02H) and b/h s known as the blockage rato (BR). The spacng rato between baffles s S = s/h = 1. Effects of DBR and GBR where nvestgated by changng ther values for 0 to 0.08. Fgure 1. Geometry of the problem studed 3. Mathematcal formulaton The numercal model for heat transfer and flud flow n the two-dmensonal channel s based on the above assumptons.the flow through the duct s governed by the Reynolds averaged Naver Stokes (RANS) equatons [39] and the energy equaton: Contnuty Equaton u 0 (1) x Momentum equaton u ' ' P u u j u u j x x x j x Where u' s a fluctuatng component of velocty. (2) Energy equaton T ut Γ Γ x x x t j velocty gradents as seen n the equaton below: 3 (3) Where: Γ s the molecular thermal dffusvty gven by: Pr Γ t s the turbulent thermal dffusvty gven by: t t Prt. The Reynolds-averaged approach to turbulence model requres the modelng of Reynolds ' ' stresses uu n Eq. (2). The Boussnesq hypothess relates the Reynolds stresses to the mean j

' ' u u 2 u u u j t k t j x j x 3 x j (4) 1 ' ' Where k s the turbulent knetc energy defned by, k u u and δ j s the Kronecker delta. 2 j The RNG-based k-ε turbulence model s derved from the nstantaneous Naver-Stokes equatons [39], usng a mathematcal technque called renormalzaton group (RNG) methods. The steady state transport equatons are expressed as: t k k u Gk x j x j k x j t u C1 C 2 x x j x j k k 2 (5) (6) Where G k s the rate of generaton of the k-ε model C 1ε and C 2ε are constants. μ t s the turbulent vscosty defned by 2 t C ( k / ),C μ s a constant set to 0.0845 and derved usng the RNG theory. The frcton factor fs calculated by the followng equaton: 2 P f L D U 2 / h where Ps the pressure drop across the length of the channel L. To determne the heat transfer characterstcs, we defne the local Nusselt number as: hx Dx Nu( x ) (8) k f and the average Nusselt number as: 1 Nu ( x ) Nu( x ) A A (9) The thermal enhancement factor η s defned by: Nu / Nu / f / f 1/3 (10) 0 0 Where Nu 0 and f 0 stand for Nusselt number and frcton factor for the smooth channel, respectvely. 4. Numercal smulaton 4.1. Boundary Condtons Ar was used as a workng flud n all cases wth a unform nlet temperature (300 K, Pr = 0.7). Impermeable boundary and no-slp wall condtons have been mplemented over the channel wall as well as the baffle. The temperature of the bottom and upper plates s mantaned constant at 330 K whle the baffle s assumed as adabatc wall condton. The Reynolds number s changed from 10000 to 20000. (7) 4

4.2. Grd Generaton Geometry and mesh of the computatonal doman were created wth the computer tool Gambt. The mesh generated s quadrlateral wth regular grd elements. Mesh tests were realzed by a seres of grd selected n the unbaffled channel as follows: 25608, 36100, 40200, 59015 and 71112 elements. From the mesh wth 40200 elements, the devaton of the Nusselt number s not mportant (do not exceed 2%) wth the rase of the grd densty. Thus and for the next calculatons, the fnal selected grd number had 40200 elements. Usng the same approach for the other cases wth baffles, the selected grds were vared between 37000 and 39500 elements. The necessary detals on mesh test are resumed on Table 1. Computatons were acheved wth the CFD FLUENT, whch s based on the method of fnte volume to solve the governng equatons. The flow regme s consdered as turbulent and steady state. The QUICK numercal scheme coupled wth the SIMPLE (Sem-Implct Pressure Lnked Equaton) algorthm were used n the fnte volume method to solve all equatons. The RNGk-ε model was employed for the closure of the equatons. Default under-relaxaton factors of the solver are used to control the update of computed varables at each teraton. These factors are: 0.3, 1, 0.7 and 1 for pressure, densty, momentum and energy, respectvely. Solutons were consdered converged when the normalzed resdual values were below 10-5 for all varables and below 10-7 only for the energy equaton. Table 1. Tests of grd sensvlty for all confgurtatons Test 1 Test 2 Test 3 Test 4 Test 5 Grd 1 Nu Grd 2 Nu Grd 3 Nu Grd 4 Nu Grd 5 Nu Smooth channel 25608 29.17 36100 19.28 40200 48.8 59015 49.10 71112 48.16 DBR=0.02 25788 162.35 37254 212.36 39554 305.49 60002 308.19 73254 307.88 DBR=0.04 24772 190.00 36992 119.24 39490 281.09 59322 278.01 70786 278.92 DBR=0.06 26710 122.89 35998 252.10 38112 273.77 58088 274.14 71250 272.63 DBR=0.08 25718 198.58 34048 212.12 38586 267.42 61002 269.78 73254 268.55 GBR=0.02 23711 50.11 36254 201.33 39554 314.76 60997 313.44 69011 313.78 GBR=0.04 24111 141.88 37982 252.02 38999 310.32 52049 308.19 65114 307.93 GBR=0.06 24352 197.28 37449 199.78 39500 314.75 58021 314.77 66194 313.58 GBR=0.08 24228 98.39 36229 177.99 38511 310.48 48049 313.13 62004 312.82 5. Results and dscusson 5.1. Verfcaton of results for a smooth channel Our predcted results for the Nusselt number and frcton factor through a smooth channel are respectvely compared wth Dttus-Boelter and Blasus correlatons. These correlatons, whch are used by for the characterzaton of flud flow and heat transfer n ducts [24], are avalable n the lterature [40] as follows: 0.4 0.8 Nu 0.023 ReD h Pr (13) 0.25 Re (14) f 0 0.0791 5

1000. f Nu Fgure 2a and 2b show, respectvely, a comparson of Nusselt number and frcton factor obtaned from the present study wth those from equatons (11) and (12). In these fgures, the presented results agree reasonably well wthn ±2.9% and ±5.5% for both frcton factor correlaton of Blasus and Nusselt number correlaton of Dttus-Boelter, respectvely. 55 Dttus-Boelter Smooth channel 50 45 40 35 30 10000 12000 14000 16000 18000 20000 Re (a) 8,0 7,8 Smooth channel of present work Blasus correlaton 7,6 7,4 7,2 7,0 6,8 6,6 10000 12000 14000 16000 18000 20000 Re (b) Fgure 2.Verfcaton of (a) Nusselt number and (b) frcton factor for a smooth channel 5.2. Flow Structure For Re = 18000 and S = 1, the axal velocty s presented on Fgure 3 for DBR = 0.08, DBR =0.04, DBR/GBR = 0, GBR = 0.04 and GBR = 0.08. Ths fgure shows the aerodynamc nature of the flow n the presence of baffles wth down- and up-graded heght. At the baffle, the flow s separated n three zones: a zone of recrculated flow behnd the baffles where the vortex appears clearly. The second zone s located n front of the frst zone where the axal velocty s very ntense, and the thrd zone s located between the frst and the second zones where the flow passes wth a mean velocty. Ths aerodynamc phenomenon depends on the sze of baffles and the drecton of changes n ther heght (up- or down-graded heght). Consequently, ths fgure shows that the rase of DBR/GBR decreases the formaton of the vortex structures. The effect of down- and up-graded sze of baffles on the formaton of vortex s llustrated n Fgure 3. The formaton of the vortex s clearly appearng n the case of channel wth down-graded sze of baffle DBR = 0.08 and 0.04. On the other hand, the vortex sze n the case of channel wth upgraded sze of baffle GBR = 0.08 and 0.04 s small compared wth DBR = 0.08 and 0.04, respectvely. 6

Ths phenomenon s due to the development of the flow along the channel n the presence of baffles. Consequently, the down-graded sze of baffle ensures further ntensfed vortces despte of the reducton of baffles sze (DBR = 0.08). Therefore, the weaker vortex s observed wth the case GBR = 0.08, because that the small baffles located at the entrance of channel do not allow an effcent developed flow. (a) (b) (c) (d) (e) Fgure 3. Velocty contours for dfferent baffled channels (a) DBR = 0.08; (b) DBR = 0.04; (c) DBR, GBR = 0; (d) GBR = 0.04; (e) GBR = 0.08 7

5.3. Pressure losses 3,1 3,0 2,9 2,8 2,7 f / f 0 2,6 2,5 2,4 2,3 2,2 2,1 DBR, GBR = 0 (Smple baffle) DBR = 0.02 DBR = 0.04 DBR = 0.06 DBR = 0.08 GBR = 0.02 GBR = 0.04 GBR = 0.06 GBR = 0.08 10000 12000 14000 16000 18000 20000 Re Fgure 4.Varaton of f/f 0 wth Reynolds number for varous baffled channels Fgure 4 presents the varaton of the normalzed frcton factor, f/f 0 wth Reynolds number values for dfferent GBR and DBR. As shown on ths fgure, the use of down-/up-graded heght of baffle leads to a consderable decrease n the frcton factor compared wth the smple baffles (GBR, DBR = 0). The DBR values present further ncreases n the frcton factor compared wth GBR values. Dependng on Reynolds number, the lowest value of frcton factor s obtaned wth DBR = 0.08, whch s consdered to be smaller by about 4-8% thant that of a smple baffle channel. 5.4. Heat transfer 1200 1000 DBR, GBR = 0 (Smple baffle) DBR = 0.02 DBR = 0.04 DBR = 0.06 DBR = 0.08 GBR = 0.02 GBR = 0.04 GBR = 0.06 GBR = 0.08 800 Re = 18000 Nu x 600 400 200 0 0,0 0,2 0,4 X/L 0,6 0,8 1,0 Fgure 5. Axal varaton of Nu X along lower channel wall for varous baffled channels, Re = 18000 At Re = 18000, the axal varaton of Nu x along the lower channel wall s presented on Fgure 5 for dfferent values of DBR/GBR. The followng phenomena can be observed: The weak thermal transfer rate remarked on the bass of baffle, and partcularly behnd t, s caused by the formaton of hot pockets n ths regon. The hghest thermal transfer rates are observed n the zones where the flow s reattached (.e. whch between two successve baffles). 8

The smple baffle allows the hghest thermal transfer rates. However, the case GBR = 0.08 gves the lowest heat transfer rates. 7,5 7,0 6,5 6,0 DBR = 0.02 DBR = 0.04 DBR = 0.06 DBR = 0.08 DBR, GBR = 0 (Smple baffle) GBR = 0.02 GBR = 0.04 GBR = 0.06 GBR = 0.08 Nu / Nu 0 5,5 5,0 4,5 4,0 3,5 3,0 10000 12000 14000 16000 18000 20000 Re Fgure 6. Varaton of Nu/Nu 0 wth Reynolds number for varous baffled channels 3,8 3,6 3,4 3,2 3,0 2,8 2,6 2,4 2,2 DBR, GBR = 0 (Smple baffle) DBR = 0.02 DBR = 0.04 DBR = 0.06 DBR = 0.08 GBR = 0.02 GBR = 0.04 GBR = 0.06 GBR = 0.08 10000 12000 14000 16000 18000 20000 Re Fgure 7. Thermal enhancement factor for varous baffled channels For a smple baffle, S = 1 and BR = 0.5, varatons of the average Nu/Nu 0 rato wth Reynolds number for dfferent down-/up-graded heghts of baffles are presented on Fgure 6. Ths fgure shows that the Nu/Nu 0 value ncreases wth the rase of Reynolds number values for all GBR and DBR cases. The use of smple baffle ensures a good thermal transfer rate and the baffles wth GBR = 0.02 reduces the heat transfer rate by 5% compared wth the smple baffle. The lowest thermal transfer rate whch s gven by the DBR = 0.08 tends to 16% compared wth the smple baffle. Fgure 7 shows the varaton of the thermal enhancement factor (η) wth Reynolds number for dfferent down-/up-graded heghts of baffles. The enhancement factor tends to ncrease wth the rase of Re for all cases. The cases of GBR = 0.02 and smple baffle (GBR = 0) present the same best enhancement factors (η) up to 3.6 tmes at the hghest value of Reynolds number. 6. Conclusons The turbulent flow and heat transfer characterstcs n a two-dmensonal channel equpped wth down-/up-graded heght of baffles have been explored numercally. The predcted results confrm that 9

the new desgn of baffles provdes an adequate reducton n frcton factors. Dependng on Reynolds number, the case of DBR = 0.08 presents the mnmum value of frcton factor whch s about 4 8% less than that of the smple baffled channel (DBR, GBR = 0). However, the frcton factor reducton s assocated wth a reducton n heat transfer rates tends to 5% as a maxmum value compared wth the smple baffled channel. Nomenclature b Smple baffle heght, [m] D Hydraulc dameter, [m] f Frcton factor, [-] h Heat transfer coeffcent, [Wm -2 K -1 ] H Channel heght, [m] G k Turbulent knetc energy producton k Turbulent knetc energy, [m 2 s -2 ] k f Thermal conductvty, [Wm -1 K -1 ] L Channel length, [m] Nu Nusselt Number Re Reynolds number,(= ρud/μ) S Spacng rato (= s/d) T Temperature,[K] u Mean velocty at the channel,[ms -1 ] Greek Letter ε Turbulent energy dsspaton μ Dynamc vscosty, [kgm -1 s -1 ] η Thermal enhancement factor, [-] ρ Densty,[kgm -3 ] Subscrpts j x drecton y drecton Abbrevatons BR Blockage rato (=b/d) DBR Down graded baffle rato (= Dhb/H) Dhb Down graded heght of baffle [m] GBR Up-Graded Baffle Rato (= Ghb/H) Ghb Up-Graded heght of baffle, [m] 10

References [1] Turgut, O., Arslan, K., Perodcally Fully Developed Lamnar Flow and Heat Transfer n a Two- Dmensonal Horzontal Channel Wth Staggered Fns, Thermal Scence, 21 (2017), 6, pp. 2443-2455 [2] Turgut, O., Arslan, K., Perodcally Fully Developed Lamnar Flow and Heat Transfer n a 2-D Horzontal Channel wth Staggered Fns, Thermal Scence, 21 (2017), 6a, pp. 2443-2455 [3] Kumar, R., et al., Expermental Investgaton on Overall Thermal Performance of Flud Flow n a Rectangular Channel wth Dscrete V-Pattern Baffle, Thermal Scence, 22 (2018), 1, pp. 183-191 [4] Lee, D. H., et al., Detaled Measurement of Heat/Mass Transfer Wth Contnuous and Multple V-Shaped Rbs n Rectangular Channel, Energy, 34 (2009), 11, pp. 1770-1778 [5] Peng, W., et al., Expermental and Numercal Investgaton of Convecton Heat Transfer n Channels wth Dfferent Types of Rbs, Appled Thermal Engneerng, 31 (2011), 14-15, pp.2702-2708 [6] Promvonge, P., et al., Numercal Heat Transfer Study of Turbulent Square-Duct Flow Through Inlne V-Shaped Dscrete Rbs, Internatonal Communcatons n Heat and Mass Transfer, 38 (2011), 10, pp. 1392-1399 [7] Promvonge, P., et al., 3D Smulaton of Lamnar Flow and Heat Transfer n V-Baffled Square Channel, Internatonal Communcatons n Heat and Mass Transfer, 39 (2012), 1, pp. 85-93 [8] Hans, V. S., et al., Heat Transfer and Frcton Factor Correlatons for a Solar Ar Heater Duct Roughened Artfcally wth Multple V-Rbs, Solar Energy,84 (2010), 6, pp. 898-911 [9] Sngh, S., et al., Exergy Based Analyss of Solar Ar Heater Havng Dscrete V-Down Rb Roughness on Absorber Plate, Energy, 37 (2012), 1, pp. 749-758 [10] Lanjewar, A., et al., Heat Transfer and Frcton n Solar Ar Heater Duct wth W-Shaped Rb Roughness on Absorber Plate, Energy, 36 (2011), 7, pp. 4531-4541 [11] Srromreun, P., et al., Expermental and Numercal Study on Heat Transfer Enhancement n a Channel wth Z-Shaped Baffles, Internatonal Communcatons n Heat and Mass Transfer, 39 (2012), 7, pp. 945-952 [12] Kumar, A., Km, M. H., Thermal Hydraulc Performance n a Solar Ar Heater Channel wth Mult V-Type Perforated Baffles, Energes 9 (2016), 7, pp. 564 [13] Kumar, R., et al., Expermental Study and Correlaton Development for Nusselt Number and Frcton Factor for Dscretzed Broken V-Pattern Baffle Solar Ar Channel, Expermental Thermal and Flud Scence, 81 (2016), 2, pp. 56-75 [14] Kumar, R., et al., Expermental Study of Heat Transfer Enhancement n a Rectangular Duct Dstrbuted by Mult V Perforated Baffle of Dfferent Relatve Baffle Wdth, Heat and Mass Transfer, 53 (2017), 4, pp. 1289-1304 [15] Kumar, R., et al., Expermental Study of Enhancement of Heat Transfer and Pressure Drop n a Solar Ar Channel wth Dscretzed Broken V-Pattern Baffle, Renewable Energy, 101 (2017), 2, pp. 856-872 [16] Kumar, A., et al., Developng Heat Transfer and Pressure Loss n an Ar Passage wth Mult Dscrete V-Blockages, Expermental Thermal and Flud Scence, 84 (2017), 6, pp. 266-278 11

[17] Sharma, A., et al., Optmzng Dscrete V Obstacle Parameters Usng a Novel Entropy-VIKOR Approach n a Solar Ar Flow Channel, Renewable Energy, 106 (2017), 6, pp. 310-320 [18] Srpattanappat, S., Promvonge, P., Numercal Analyss of Lamnar Heat Transfer n a Channel wth Damond-Shaped Baffles, Internatonal Communcatons n Heat and Mass Transfer, 36 (2009), 1, pp. 32-38 [19] Dutta, P., Hossan, A., Internal Coolng Augmentaton n Rectangular Channel Usng Two Inclned Baffles, Internatonal Journal of Heat and Flud Flow, 26 (2005), 2, pp. 223-232 [20] Promvonge, P., Thanpong, C., Thermal Performance Assessment of Turbulent Channel Flows over Dfferent Shaped Rbs, Internatonal Communcatons n Heat and Mass Transfer, 35 (2008), 10, pp.1327-1334 [21] Promvonge, P., Kwankaomeng, S., Perodc Lamnar Flow and Heat Transfer n a Channel wth 45 Staggered V-Baffles, Internatonal Communcatons n Heat and Mass Transfer, 37 (2010), 7, pp. 841-849 [22] Promvonge, P., Heat Transfer and Pressure Drop n a Channel wth Multple 60 V-Baffles, Internatonal Communcatons n Heat and Mass Transfer, 37 (2010), 7, pp. 835-840 [23] Promvonge, P., et al., Thermal Behavor n Solar Ar Heater Channel Ftted wth Combned Rb and Delta-Wnglet, Internatonal Communcatons n Heat and Mass Transfer, 38 (2011), 6, pp. 749-756 [24] Tandroglu, A., Ayhan, T., Energy Dsspaton Analyss of Transent Heat Transfer for Turbulent Flow n a Crcular Tube wth Baffle Inserts, Appled Thermal Engneerng, 26 (2006), 2, pp. 178-185 [25] Tandroglu, A., Effect of Flow Geometry Parameters on Transent Heat Transfer for Turbulent Flow n a Crcular Tube wth Baffle Inserts, Internatonal Journal of Heat and Mass Transfer, 49 (2006), 9-10, pp. 1559-1567 [26] Tandroglu, A., Effect of Flow Geometry Parameters on Transent Entropy Generaton for Turbulent Flow n Crcular Tube wth Baffle Inserts, Energy Converson and Management, 48 (2007), 3, pp. 898-906 [27] Sngh, S., et al., CFD (Computatonal Flud Dynamcs) Investgaton on Nusselt Number and Frcton Factor of Solar Ar Heater Duct Roughened wth Non-Unform Cross-Secton Transverse Rb, Energy 84 (2015), 5, pp. 509-517 [28] Yang, W., Expermental Study on The Heat Transfer Characterstcs of Hgh Blockage Rbs Channel. Expermental Thermal and Flud Scence 83 (2017), 5, pp. 248-259 [29] Skullong, S., et al., Effects of Rb Sze and Arrangement on Forced Convectve Heat Transfer n a Solar Ar Heater Channel, Heat and Mass Transfer 51 (2015), 10, pp. 1475-1485 [30] Sahel, D., et al., Enhancement of Heat Transfer n a Rectangular Channel wth Perforated Baffles, Appled Thermal Engneerng, 101 (2016), pp. 156-164 [31] Sahel, D., et al., A Numercal Study of Flud Flow and Heat Transfer over a Fn and Flat Tube Heat Exchangers wth Complex Vortex Generators, The European Physcal Journal Appled Physcs, 78 (2017), pp. 34805. [32] Kumar, R., et al., Comparatve Study of Effect of Varous Blockage Arrangements on Thermal Hydraulc Performance n a Roughened Ar Passage, Renewable and Sustanable Energy Revews 81 (2018), 1, pp. 447-463 12

[33] Davd Huang, K., et al., Expermental Study of Flud Flow and Heat Transfer Characterstcs n the Square Channel wth a Perforaton Baffle, Internatonal Communcatons n Heat and Mass Transfer, 35 (2008), 9, pp. 1106-1112 [34] Karwa, R., Maheshwar, B.K., Heat Transfer and Frcton n an Asymmetrcally Heated Rectangular Duct wth Half and Fully Perforated Baffles at Dfferent Ptches, Internatonal Communcatons n Heat and Mass Transfer, 36 (2009), 3, pp. 264-268 [35] Kumar, A., et al., Effect of Roughness wdth Rato n Dscrete Mult V-Shaped Rb Roughness on Thermo-Hydraulc Performance of Solar Ar Heater, Heat Mass Transfer, 51 (2015), 2, pp. 209-220 [36] Boukhada, K., et al., Effect of the Perforaton Desgn on the Flud Flow and Heat Transfer Characterstcs of a Plate Fn Heat Exchanger, Internatonal Journal of Thermal Scences, 126 (2018), 4, pp. 172-180 [37] Alem, K., et al., CFD Investgatons of Thermal and Dynamc Behavors n a Tubular Heat Exchanger wth Butterfly Baffles, Fronters n Heat and Mass Transfer (FHMT), 10, (2018), pp. 27 [38] Mellal, M., et al., Hydro-Thermal Shell-Sde Performance Evaluaton of a Shell and Tube Heat Exchanger under dfferent Baffle Arrangement and Orentaton, Internatonal Journal of Thermal Scences, 121 (2017), 11, pp. 138-149 [39] Fluent Inc,. (2005), User s Gude 6.2, Centerra Resource Park 10 Cavendsh Court Lebanon, NH 03766. [40] Incropera, F., Dewt, P. D., Introducton to Heat Transfer, Ffth edton John Wley & Sons Inc (2006) 13