optimal thermohydraulic performance in three sides

0 downloads 0 Views 3MB Size Report
laminar flow and heat transfer characteristics in a three dimensional .... winglet type vortex generators (DWs) on forced convection heat transfer and friction loss ...... and the balance length of 630 mm served as the entry length for flow ...
OPTIMAL THERMOHYDRAULIC PERFORMANCE IN THREE SIDES ARTIFICIALLY ROUGHENED SOLAR AIR HEATER

A THESIS Submitted in fulfillment of the requirements of the award of the degree of DOCTOR OF PHILOSPHY In MECHANICAL ENGINEERING By ASHWINI KUMAR (2013PGPHDME08)

Under the supervision of Prof. Bishwa Nath Prasad & Dr. Krishna Deo Prasad Singh Department of Mechanical Engineering National Institute of Technology, Jamshedpur Jamshedpur-831014 (India) January 2017

NATIONAL INSTITUTE OF TECHNOLOGY JAMSHEDPUR CERTIFICATE I hereby certify that the work presented in the thesis entitled " Optimal Thermohydraulic Performance In Three Sides Artificially Roughened Solar Air Heater" in partial fulfillment of the requirements for the award of the degree of Doctor of Philosophy in Mechanical Engineering, submitted in the Department of Mechanical Engineering, National Institute of Technology, Jamshedpur is an authentic record of my own work carried out under the supervision of Prof. B.N. Prasad and Dr. K.D.P. Singh, Department of Mechanical Engineering, National Institute of Technology, Jamshedpur. The matter presented in the thesis has not been submitted by me for the award of any other degree of this or any other institute.

Date:

(Ashwini Kumar)

This is to certify that the above statement made by the candidate is correct to the best of our knowledge.

(Dr. B.N. Prasad)

(Dr. K.D.P. Singh)

Professor

Associate Professor

Department of Mechanical Engineering

Department of Mechanical Engineering

NIT, Jamshedpur, India 831014

NIT, Jamshedpur, India 831014

The Ph.D. thesis Viva-voce Examination of Mr. Ashwini Kumar, Research Scholar was held on ______________.

Signature of Internal Examiner

Signature of External Examiner i

ACKNOWLEDGEMENTS The author expresses his deep sense of gratitude and indebtedness to Dr.B.N.Prasad, Professor and Dr.K.D.P.Singh, Associate Professor both from Mechanical Engineering Department, NIT Jamshedpur, who provided their valuable guidance, unfailing inspiration, whole hearted cooperation and help in carrying out this work. Their encouragement at each step of the work and painstaking efforts in providing in final shape to the text are gratefully acknowledged. The author gratefully acknowledges the facility provide by the institute and the department. Thanks are due to staff of Mechanical Engineering Department for their co-operation at the every step of the work. Last but not the least, the moral support and patience of my family members especially my parents Mr.Yogendra Roy & Mrs. Kiran Yadav, my brother Amit Kumar, my wife Mrs. Bandana Singh and motivation given by my fellow friends especially Ravi Kumar, Abhishek Priyam and Aruna Kumar Behura is highly acknowledged.

ASHWINI KUMAR

ii

My Work Is Dedicated To

My Parents And My Little Angel

Kashvi Kavya iii

ABSTRACT Solar air heaters are one of the simplest and cost effective solar energy utilization systems for various drying applications. The thermal efficiency of solar air heaters is found to be low due to low heat transfer coefficient on the air side. Attempts have been made to enhance the heat transfer rate from the absorber plate to air by using different roughness elements. Enhancement of heat transfer coefficient is achieved by providing artificial roughness, always goes with an increment of pressure drop and more pumping power is required. So, there is a need to optimize the system parameters to maximize heat transfer, while keeping the friction losses as low as possible. The present work aimed at obtaining the optimal thermo hydraulic performance condition in three sides roughened and glass covered solar air heaters to get the maximum heat transfer for minimum pumping power. An extensive experimental investigation has been carried out with the objective of obtaining and validating the results on heat transfer, friction factor, and thermal performance of three sides roughened and glass covered solar air heaters along with those on smooth and three sides glass covered solar air heaters. Roughened absorber plates were prepared by soldering small diameter wires on them. The operating parameters covered the range of relative roughness height, e/D, 8.35 × 10-3 – 3.75 × 10-2, relative roughness pitch, p/e, 10 – 30, flow Reynolds number, Re, 4000 – 20000 and Prandtl number, Pr, 0.7 for air. Data were collected on a similar sized three sides glass covered smooth solar air heater simultaneously. Calibrated copper constantan thermocouples were soldered on the absorber plates to measure local plate temperature, whereas, the temperature of air was measured by digital thermometer inserted midway into the air duct. Air was sucked through the collectors using a blower with electric power input controlled by an auto-variac to obtain variable flow rates, the flow rate being measured by a U-tube manometer fitted in a calibrated flange-tap orifice-meter. The intensity of global radiations was measured by digital pyranometer. The reduced data were further utilized to obtain the thermo hydraulic results. The results have been found to compare well with the available ones. Correlations have been developed for heat transfer and friction coefficient in three sides roughened and glass covered as well as three sides glass covered smooth solar air heaters, written under as:-

iv

𝑁𝑢3𝑟 = 0.1415 𝑝/𝑒

−0.0713

𝑒/𝐷

0.0871

𝑅𝑒

0.746

for 𝑒 + ≤ 23

𝑁𝑢3𝑟 = 0.0422 𝑝/𝑒

−0.0461

𝑒/𝐷

0.0375

𝑅𝑒

1.208

for 𝑒 + > 23

𝑓3𝑟 = 0.372 𝑝/𝑒

−0.319

𝑁𝑢3𝑔𝑠 = 0.0782 𝑃𝑟

0.5

𝑒/𝐷

0.354

0.7

𝑓𝑆

𝑅𝑒

𝑅𝑒

−1.76

0.932

A good agreement has been found between predicted and experimental values with a deviation of ±12% and±11% respectively, for Nusselt number and friction factor. The effects of roughness and flow parameters (p/e, e/D, Re) on heat transfer, friction factor and thermo hydraulic performance have been determined. Results show that the values of Stanton number are higher in three sides roughened and glass covered solar air heater and they are in the range of 36% to 65% more as compared to three sides glass covered smooth ones and the values of friction factor in three sides roughened and glass covered collector have been found to increase in the range of 28% to 38%, as compared to three sides glass covered smooth ones. The value of heat transfer enhancement factor of such solar air heaters have been found to lie in the range of0.378 to 0.487,resulting in an enhancement of about 50 to 59% in thermal performance over those of three sides glass covered smooth solar air heaters, in the range of the parameters investigated. The values of performance parameters, FR and F’ have been found to be more than those of three sides glass covered smooth solar air heaters by an amount of 51-60% and 41-49% respectively. Thermal performance equations in terms of the performance parameters have been derived. The roughness Reynolds numbers, 𝑒 + = 𝑒/𝐷

𝑓3𝑟 2

𝑅𝑒, have

been considered for optimizing the thermo hydraulic performance of three sides roughened and glass covered solar air heaters. The optimal thermo hydraulic performance of such solar + air heaters has been achieved at an optimal value of the optimization parameter, 𝑒𝑜𝑝𝑡 = 23.

The experimental results have been utilized to develop the correlations for thermo hydraulic efficiency written under as:Ƞ𝑡𝑕𝑒𝑟𝑚𝑜 = 0.102 𝑒 +

0.658

for 𝑒 + ≤ 23

Ƞ𝑡𝑕𝑒𝑟𝑚𝑜 = 11.91 𝑒 +

−0.913

for 𝑒 + > 23

The deviation between the predicted and the experimental values of thermo hydraulic efficiency was found to be ±7.23%. v

An uncertainty of +1.89% to +7.06% in the values of thermo hydraulic performance has + been found. The value of optimal thermo hydraulic performance corresponding to 𝑒𝑜𝑝𝑡 = 23

have been found to be 78.8% within the range of the parameters investigated. It can therefore be concluded that a solar air heater with three sides artificially roughened absorber plates and glass covers perform considerably better with respect to energy collection as compared to three sides glass covered smooth solar air heaters and the optimal thermohydraulic performance of three sides roughened and glass covered solar air heaters is + obtained when 𝑒𝑜𝑝𝑡 = 23.

Keywords: Relative roughness pitch, (p/e); Relative roughness height, (e/D); Flow Reynolds number, (Re); Stanton number, (St); Collector roughness and flow parameters; Collector performance parameters; Thermal performance; Roughness Reynolds number; and Thermo hydraulic performance.

vi

CONTENTS CERTIFICATE

i

ACKNOWLEDGEMENT

ii

ABSTRACT

iii

CONTENTS

vi

LIST OF FIGURES

ix

LIST OF TABLES

xii

NOMENCLATURE

xiii

Chapter No. 1

Title INTRODUCTION AND LITERATURE SURVEY

1-22

1.1 Solar thermal collectors

1

1.2 Solar air heaters

1

1.3 Solar air heater performance

2

1.4 Effect of roughness and flow parameters in solar air heaters

7

1.4.1

20

Optimization of the system parameters

1.5 Objectives of the present work

2

Page No.

OPTIMAL

THERMOHYDRAULIC

22

PERFORMANCE 23-40

ANALYSIS 2.1

Introduction

23

2.2

Optimization methodology

23

2.3

Effect of efficiency parameters

26

2.4

Thermohydraulic performance

36

vii

3

4

SYSTEM DEVELOPMENT AND INVESTIGATIONS

41-59

3.1

Introduction

41

3.2

Solar air heater ducts models

41

3.3

Experimental Set-up

42

3.3.1 Absorber plates and glass covers

43

3.3.2 Test section and entry section of the duct

43

3.4

48

Instrumentation

3.4.1 Blower and electric motor

48

3.4.2 Flow control arrangement

48

3.4.3 Air flow measurement

48

3.4.4 Temperature and Pressure drop measurement

51

3.4.5 Solar radiation intensity measurement

54

3.5

54

Experimental procedure

RESULTS AND DISCUSSIONS

60-110

4.1

Introduction

60

4.2

Data presentation, reduction and validation

60

Sample calculation

61

4.3

Variation of intensity and ambient temperature

62

4.4

Validation of the experimental data

66

4.5

Development of correlations for Nusselt number and 71

4.2.1

friction factor 4.5.1 Correlation for Nusselt number of three sides glass covered 78 smooth solar air heater 4.5.2 Correlation for friction factor of three sides roughened and 82 glass covered solar air heater 4.6 Effect of roughness and flow parameters on heat transfer and 91 friction factor 4.7 Heat transfer enhancement factor and friction loss factor

96

4.8 Thermal performance of solar air heater

96

viii

5

4.9 Collector heat removal factor and collector efficiency factor

101

4.10 Thermal performance curves

104

4.11 Optimal thermo hydraulic performance

106

CONCLUSIONS AND FUTURE SCOPE

111-112

5.1 Conclusions

111

5.2 Future scope of the work

112

REFERENCES

113-121

UNCERTAINTY ANALYSIS

122-124

LIST OF PUBLICATIONS FROM THIS WORK

125

ix

LIST OF FIGURES Figure No.

Title of Figures

Page No.

2.1(a)

Three sides roughened solar air heater duct

24

2.1(b)

Absorber plate with artificial roughness

24

2.2

Effect of e+ on B-1 for varying p/e

27

2.3

Effect of e+ on C-1 for varying p/e

29

+

-1

2.4

Effect of e on L for varying p/e

30

2.5

The efficiency parameter L-1 as a function of p/e and e+

31

2.6.(1-3)

Optimal thermo hydraulic performance curve for three sides 33-35 artificially roughened solar air heater

2.7

Optimum thermo hydraulic performance curves

37

2.8

Optimal range of e/D

39

3.1 (a,b)

Solar air heater duct models

42

3.2(a)

Schematic diagram of experimental set-up

44

3.2 (b)

Top view of the experimental set-up

45

3.3

Photograph of experimental set-up with two ducts

45

3.4

Photograph of experimental set-up with measuring instruments

46

3.5

Photograph of three sides roughened absorber plate

46

3.6

Typical three sides artificially roughened absorber plate

47

3.7

Typical three side smooth absorber plate

47

3.8

Photograph of digital pyranometer system

48

3.9

Orifice plate location and specification

49

3.10

Details of flange-taps orifice meter construction

50

3.11

Orifice meter calibration

52

3.12

Thermocouple arrangement on top and side of the absorber plates

53

3.13

Arrangement of air temperature and pressure taps along the duct 53 length

4.1

Typical variation of solar radiation and ambient temperature on a 64 day

4.2

Typical air and plate temperature in solar air heaters

65

4.3

Comparison of heat transfer data in smooth solar air heater

67

x

4.4

Comparison of friction factor data in smooth solar air heater

68

4.5

Comparison of heat transfer data in three sides roughened solar air 69 heater

4.6

Comparison of friction factor data in three sides roughened solar air 70 heater

4.7

Variation of Nusselt number with Reynolds number (𝑒 + ≤ 23)

4.8

3𝑟 Plot of (𝑅𝑒)0.746 versus p/e

𝑁𝑢

72 73

𝑁𝑢 3𝑟 (𝑅𝑒)0.746

4.9

Plot of

4.10

Variation of Nusselt number with Reynolds number (𝑒 + > 23)

4.11

3𝑟 Plot of (𝑅𝑒)1.208 versus p/e

(𝑝/𝑒)−0.0713 versus e /D

𝑁𝑢

74

75 76

𝑁𝑢 3𝑟 (𝑅𝑒)1.208

4.12

Plot of

4.13

Nusselt number versus flow Reynolds number for smooth solar air 79

(𝑝/𝑒)−0.0461 versus e/D

77

heater 4.14

𝑁𝑢 3𝑔𝑠

80

Plot of (𝑅𝑒)0.932 versus Pr 𝑁𝑢 3𝑔𝑠 0.932

4.15

Plot of (𝑅𝑒)

(𝑃𝑟 )0.5 versus fs

81

4.16

Variation of friction factor with Reynolds number

83

4.17

Plot of

𝑓3𝑟 (𝑅𝑒)−1.76

versus p/e

84

(𝑝/𝑒)−0.319 versus e/D

85

𝑓3𝑟 (𝑅𝑒)−1.76

4.18

Plot of

4.19

Comparison of experimental Nusselt number values to the predicted 86 ones

4.20

Comparison of experimental Nusselt number values to the predicted 87 ones

4.21

Predicted and experimental heat transfer data in solar air heater xi

89

4.22

Predicted and experimental values of friction factor data in solar air 90 heater

4.23

Effect of p/e on Stanton number in three sides artificially roughened 92 solar air heater

4.24

Effect of e/D on Stanton number in three sides artificially 93 roughened solar air heater

4.25

Effect of p/e on friction factor in three sides artificially roughened 94 solar air heater

4.26

Effect of e/D on friction factor in three sides artificially roughened 95 solar air heater

4.27

Heat transfer enhancement factor and friction loss factor versus 97 Reynolds number

4.28

Effect of roughened sides and glass covers on Stanton number in 99 solar air heaters

4.29

Thermal performance characteristics of solar air heaters

100

4.30

Collector heat removal factor in solar air heaters

102

4.31

Collector efficiency factor in solar air heaters

103

4.32

Thermal performance curves for solar air heaters

105

4.33

Effect of p/e and e+ on thermo hydraulic performance

107

+

4.34

Effect of e/D and e on thermo hydraulic performance

108

4.35

Thermohydraulic efficiency versus roughness Reynolds number

110

xii

LIST OF TABLES Table No.

Title of Tables

Page No.

2.1

Range of efficiency parameter, L-1, as a function of e+ and p/e

2.2

+ Values of 𝑒𝑜𝑝𝑡 , for different roughness geometries for fully 36

32

developed turbulent flow 2.3

Values of Re, at varying value of e/D and p/e=10, for the optimal 40 condition

3.1

Test data collected during a typical run (Run no. 56)

55

3.2

Range of parameters investigated

56

4.1

Values of I, Ta, h, St and Nu as a function of time during the day 63 (Run no. 56)

4.2

Values of slope, intercept, collector heat removal factor and collector 104 efficiency factor

xiii

NOMENCLATURE

𝑨𝒄

Collector area, m2

𝑨𝟐

Throat area of orifice meter, m2

𝑨𝒇

Friction similarity function, m2

𝑩

Solar air heater duct height, m

𝑩−𝟏

Stanton number roughness parameter

𝑪𝑭

Dimensionless friction factor coefficient

𝑪𝑯

Dimensionless heat transfer coefficient

𝑪−𝟏

Efficiency roughness parameter

Cp

Specific heat at constant pressure of air, kJ/kg.K

𝑪𝒅

Discharge coefficient

𝑫, 𝑫𝒉

Hydraulic diameter of solar air heater duct, m

d1

Flow pipe diameter, m

d2

Throat diameter of orifice meter, m

d/W

Relative gap position

d/D

Relative print diameter

Dj/Dh Ratio of jet diameter to hydraulic diameter e 𝒆+ 𝒆+ 𝒆+ 𝒐𝒑𝒕

Artificial roughness height, m Roughness Reynolds number = 𝑒/𝐷 Roughness Reynolds number = 𝑒/𝐷

𝑓3𝑟 2 𝑓 2

𝑅𝑒 (Present case) 𝑅𝑒 (Prasad and Saini, 1991)

Optimal roughness Reynolds number 𝒆/𝑫 𝒆/𝑫𝒐𝒑𝒕

Relative roughness height Optimal relative roughness height

𝒇𝒔

Friction factor in four sided smooth duct

𝒇𝒓

Average friction factor in three sides roughened duct (Prasad et al., 2014)

𝒇𝟑𝒓

Friction factor for in three sides roughened duct (Present case)

𝒇

Average friction factor=(𝑓𝑠 + 𝑓𝑟 )/2, (Prasad and Saini, 1991)

𝒇

Friction factor

𝑭′

Collector efficiency factor xiv

𝑭"

Collector flow factor

𝑭𝑹

Collector heat removal factor

𝑭𝟎

Collector heat removal factor based on outlet temperature, 𝑡0

𝑭𝟎 𝑼 𝑳 𝑭𝟎 (𝝉𝜶)

Slope of the thermal performance curve Intercept of the thermal performance curve

𝑭𝒂𝒗

Average collector efficiency factor

𝒇𝑭

Friction loss factor

𝒇

Moderate pressure drop increase 𝒇𝟎 𝑮

Mass flow rate per unit collector area, kg/s.m2

g

Acceleration due to gravity, m/s2

g’

Heat transfer function

g/e

Relative gap width

g/p

Groove position to pitch ratio

GH

Heat transfer roughness function

𝒉

Convective heat transfer coefficient, W/m2.K

𝑰

Intensity of solar radiation, W/m2

k

Thermal conductivity of collector material, W/m.K

𝑲′

Geometric coefficient of collector, 𝐾 ′ = 𝐿/ 𝐵𝐷

L

Collector length, m

L-1

Efficiency parameter, L-1 = C-1 - B-1

L/e

Relative long way length

l/s

Relative length of metal grid

𝒎

Mass flow rate of air, kg/s

0.25

𝑵𝒖𝒓 , 𝑵𝒖 Nusselt number for top side roughened duct 𝑵𝒖𝒓

Average Nusselt number for three sided roughened duct (Prasad et al., 2014)

𝑵𝒖𝟑𝒓

Nusselt number for three sides roughened duct (Present case)

𝑵𝒖𝟑𝒈𝒔

Nusselt number for three sides glass covered smooth duct (Present case)

𝑵𝒖

Nusselt number

xv

𝑵𝒖𝑭

Heat transfer enhancement factor

𝑵𝒖𝟎

Nusselt number for fully developed flow smooth channel

𝑵𝒖 𝑵𝒖𝟎

Heat transfer augmentation

𝒑

Roughness pitch, m

P

Motor power

P1, P2 Pressure in two limbs of the manometer, N-m2 p/e

Relative roughness pitch

𝑷𝒓

Prandtl number

𝑷𝒓𝒕

Turbulent Prandtl number

QU

Useful heat gain, W

𝑹

Roughness function

𝑹𝒆

Reynolds number

𝑹𝒅

Throat Reynolds number

𝑺𝒕𝟑𝒓

Average Stanton number for three sides roughened duct (Present case)

𝑺𝒕 , 𝑺𝒕

Stanton number for top side roughened duct

𝑺𝒕𝒓 𝑺𝒕𝑺 𝑺𝒕𝟑𝒈𝒔

Average Stanton number for three sides roughened duct (Prasad et al., 2014) Stanton number for smooth solar air heaters Stanton number for three sides glass covered smooth solar air heaters

SWG Standard Wire Gauge 𝑻𝒂

Ambient air temperature, K

𝑻𝒊

Inlet temperature of air, K

𝑻𝟎

Outlet temperature of air

𝑻𝒇

Fluid temperature at different location of the collector

𝑻𝒑

Plate temperature at different location of the absorber plate

𝐓𝐟

Average temperature of fluid

𝑻𝒑

Average temperature of plate

𝜶

Rib angle of attack

𝝉𝜶

Transmittance-absorptance product

ƞ𝒕𝒉

Collector thermal efficiency

xvi

Ƞ

Optimization parameter

ƞ𝒕𝒉𝒆𝒓𝒎𝒐

Thermo hydraulic efficiency

𝑼𝑳

Overall heat transfer coefficient, W/m2.K

V2

Fluid velocity in the throat, m-s-1

Va

Volume flow rate of air, m3/s

𝑾

Solar air heater duct width

w

Width of rib, m

W/w Relative roughness width W/H Duct aspect ratio X

Distance in x-direction, m

Y

Distance in y-direction, m

X/Dh Ratio of streamwise pitch to hydraulic diameter Y/Dh Ratio of spanwise pitch to hydraulic diameter 𝝆

Air density, kg/m3

∆𝐩 Pressure drop, N/m2 ∆𝐓 Temperature rise 𝝑 kinematic viscosity 𝜷𝟏

Diameter ratio, 𝛽1 =

𝑑2 𝑑1

1r,1G One side roughened collector with only top side glass cover 3r,3G Three sides roughened collector with three sides glass covers (top and sides) S,1G Smooth collector with only top side glass cover S, 3G Smooth collector with three sides glass covers (top and sides)

xvii

CHAPTER- 1 INTRODUCTION AND LITERATURE SURVEY

Since the beginning of time, people have been fascinated by the power of sun. Ancient civilizations personified the sun, worshipping as a God. Throughout history, farming and agriculture efforts have relied upon the sun’s rays to grow crops and sustain population. However, we have developed the ability to harness the sun’s awesome power. The resulting technologies have promising implications for the future of renewable energy and sustainability. 1.1 SOLAR THERMAL COLLECTORS A solar thermal collector collects heat by absorbing sunlight. A collector is a device for capturing solar radiation. Solar radiation is energy in the form of electromagnetic radiation from the infrared (long) to the ultraviolet (short) wavelengths. The quantity of solar energy striking the earth's surface (solar constant) averages about 1000 watts per square meter under clear skies, depending upon weather conditions, location and orientation. The term ‘solar collector’ commonly refers to solar hot water panels, but may refer to installations such as solar parabolic troughs and solar or basic installations such as solar air heaters. Concentrated solar power plants usually use the more complex collectors to generate electricity by heating a fluid to drive a turbine connected to an electrical generator. Simple collectors are typically used in residential and commercial buildings for space heating. Solar collectors are either non-concentrating or concentrating. In the nonconcentrating type, the collector area (i.e., the area that intercepts the solar radiation) is the same as the absorber area (i.e., the area absorbing the radiation). In these types the whole solar panel absorbs light. Concentrating collectors have a bigger interceptor than absorber. 1.2 SOLAR AIR HEATERS Solar air heating is a solar thermal technology in which the energy from the sun is captured by an absorber and used to heat air flowing over it. Solar air heating is a renewable

1|Page

energy heating technology used to heat or condition air for buildings or process heat applications. It is typically the most cost-effective out of all the solar technologies. Solar air heaters are a system that collect solar energy and transfers the heat to passing air, which is either stored or used for space heating. The collectors are often black to absorb more of the sun's energy and a conductive material, often metal, acts as a heat exchanger. Solar air heaters can compliment traditional indoor heating systems by providing a free and clean source of heat. 1.3 SOLAR AIR HEATER PERFORMANCE The useful thermal energy for quasi-steady state condition, delivered by a solar collector is obtained from the energy balance between the energy absorption rate of the absorber & thermal energy losses directly and indirectly to the surroundings. It can be written as: Qu Ac

= I τα − UL (Tp − Ta )

(1.1)

The difficulty is to utilize the above equation, which lies in the calculation of mean absorber plate temperature, Tp , which is a function of the collector configuration, the incident solar radiation and entering fluid condition. Therefore, Eq. (1.1) is usually represented in a more practical form as: Qu Ac

= FR [I τα − UL Ti − Ta ]

(1.2)

where, the heat removal factor FR is defined as the ratio of rate of actual useful energy collected to the rate if the entire collector surface were at fluid inlet temperature and can be written as: FR = A

m C P (T o −T i ) c [I

(1.3)

τα −U L T i −T a ]

A very widely used basic equation (1.2), known as (Hottel-Whillier-Bliss equation) can be modified to include thermal efficiency as: Ƞth = FR

τα − UL

T i −T a I

(1.4)

2|Page

This equation expresses the thermal efficiency of the collector in terms of two major operating parameters, incidence solar radiation normal to the collector plate and temperature difference Ti − Ta . Eqs. (1.2) and (1.4) have been in wide use to compare the performance characteristics of different types of collectors and to investigate the effect of various design parameters. Heat removal factor FR is a function of plate efficiency factor F ′ defined by Bliss (1959), written under as: h

F ′ = h+U

(1.5)

L

The value of F ′ is a strong function of configuration parameters given by (Duffie et al. 1980). A third factor known as collector flow factor F ′′ has been included (Duffie et al. 1980) to present equations given below: m CP

FR = A F ′′ =

c UL

FR F′

=

1 − exp − m CP A c U L F′

A c U L F′

(1.6)

m CP

1 − exp −

A c U L F′

(1.7)

m CP

The collector flow factor F ′′ is a function of dimensionless collector capacitance rate,

m CP A c UL F′

.

The collector thermal efficiency, Ƞth , given by Eq. (1.4) can be written as: Ƞth = FR τα − FR UL

(T i −T a ) I

(1.8)

This can be seen to have the form of the equation of straight line, (y = mx + c), as the term FR and UL are usually nearly constant over the operating range of the collector. Thus, FR τα , is the intercept and −FR UL is the slope of the straight line drawn on a graph with abscissa units of

(T i −T a ) I

and ordinate of Ƞth .

The terms FR τα and FR UL are used as thermal performance parameters of the flat plate solar collectors. A solar air heater working on an open cycle which is drawing ambient air outside, the inlet temperature always coincides with the ambient temperature i.e. Ti = Ta. Under such operating conditions, the conventional equation (1.4) will not be meaningful as this reduces 3|Page

toȠth = FR τα . Therefore, for a solar air heater, it was considered useful to utilize the following equation given by (Duffie and Beckman, 1980; Reddy and Gupta, 1980; Biondi et al., 1988): Ƞth = F0

(T 0 −T a )

τα − UL

(1.9)

I

where,F0 is the heat removal factor referred to the outlet temperature, also expressed by Biondi et al (1988) as: F0 = GCP

e

U L F ′ /GC P

−1

(1.10)

UL

Further performance can also be expressed as: Ƞth = GCP

(T 0 −T i )

(1.11)

I

For Ti = Ta, the variable reduced temperatures if Eqns. (1.9) and (1.11) coincide. Therefore, Eqns. (1.9) and (1.11) could be represented on a single diagram having the same quantity in the abscissa (Reddy and Gupta, 1980; Biondi et al, 1988). This was considered the most suitable way to describe the characteristic behaviors of solar air collectors (Reddy and Gupta, 1908; Biondi et al, 1988; Khalid et al, 1987). According to Duffie and Beckman (1980), an approximate analytical correction can be made to FR τα and FR UL for flow rates for which actual test results are not available. It was proposed that the data be plotted as Ƞth versus Tav =

T i +T 0 2

T av −T a I

instead of usual

(T 0 −T i ) I

taking

. Then the performance equation could be in the form of:

Ƞth = Fav τα − Fav UL

T av −T a I

(1.12)

where, Fav approximated to F ′ and corrections for FR τα and FR UL could be obtained from the equations (Duffie and Beckman, 1980) given below: m CP

FR τα = Fav τα

m CP

AC F av U L 2 AC +

(1.13)

m CP

FR UL = Fav UL

m CP

AC F av U L 2 AC +

(1.14) 4|Page

However, if Ƞth is plotted against

(T 0 −T a ) I

, as is often done for air heater test data (Reddy

and Gupta, 1980; Biondi et al, 1988), the equation for Ƞth could be written as: Ƞth = F0 τα − F0 UL

(T 0 −T a ) I

(1.15)

and then the correction for FR τα and FR UL would be written as: m CP

FR τα = F0 τα

m CP m CP

F R UL = F 0 UL

m CP

AC

A C +F 0 U L

AC

A C +F 0 U L

(1.16)

(1.17)

(Gillet et al. 1983; Bernier and Plett, 1988) declared that the thermal performance representation, testing procedures and mathematical models associated with the solar air collectors had not been developed to the same extent as for liquid collector. The two basic differences between air and liquid collectors are (i) relatively poor thermal conductivity of the air resulting in lower absorber plate to air heat transfer coefficient as compared to that for liquid collectors, which implies that the former have a strong dependence on flow rate, (ii) most air collectors leak since it is not economical to build a perfectly sealed solar air collector. An important aspect of the analytical methods of performance prediction of solar air heaters is the heat transfer coefficient in between the absorber plate and air, and the pressure drop in the collector. A number of empirical relations have been predicted for different configurations and working conditions of the collectors. The correlation given by Kays (1966) for heat transfer for flow between parallel plates, with the upper being heated and lower plate insulated has been utilized by Duffie and Beckman (1980) and Khalid et al (1987) for fully developed turbulent flow as: NuS = 0.0158 Re0.8

(1.18)

For solar air heaters operating in the range of 2000 < Re < 10000, Jansen (1985) utilized the relation recommended by Shewen and Hollands (1980 &1978) based on the data of Kays and London (1964) as: NuS = 0.00269 Re

(1.19) 5|Page

The corresponding friction factor was given by Shewen and Hollands (1978) as: fS = 0.1335 Re−0.317

(1.20)

The equation suggested by Altemani and Sparrow (1980) for flow in a triangular duct with heated inclined surfaces and an insulated third surface utilized by Khalid et al (1987) for the mass flow rate between 0.01 – 0.05 kg/s.m2 was given as: NuS = 0.019 Re0.781

(1.21)

Biondi et al (1988) expressed Kays’ Eq. (1.18) by introducing into it, the geometric coefficient of collector given by: L

K ′ = BD 0.25

(1.22)

NuS = 0.0158

GK ′

0.8

D

μ

(1.23)

For fully developed turbulent flow inside tubes, the correlation for heat transfer and friction factor as recommended by Duffie and Beckman (1980) are: fs

Nus = 1.07+12.7√f8

Re P r s √8(p r

2/3 −1)

fS = 0.79 ln Re − 1.64

−2

(1.24)

(1.25)

Prasad and Mullick (1983) recommended for bare cover solar air heater duct in Reynolds number range of 5000 – 25000: Nus = 0.014 pr 0.5 fs 0.5 Re

(1.26)

Equation of Malik and Buelow (1975) which agreed within ±10 % for Reynolds number ranging from 10000 – 20000 in solar air heater could be written as: 0.01344 Re 0.75

Nus = 1−1.586 Re −0.125

(1.27)

The well known equation of Blasius utilized for the calculation of pressure drop for smooth surfaces in solar air heaters (Sukhatme and Nayak, 2008) for 5000 < Re < 30000 is written as:

6|Page

fs = 0.079Re−0.25

(1.28)

Another equation in this range (Kays, 1966) is given by: 1 √4f s

= −0.8 + 0.87 ln Re√4fs

(1.29)

and in the range of 30000 < Re < 100000 Eq. (1.29) closely can be approximated as: fs = 0.046 Re−0.2

(1.30)

Moody chart could also be utilized for friction data in fully developed flow in smooth ducts. However, a number of researchers have focused their study towards the enhancement of thermal performance of conventional solar air heaters, but are limited to only one side glass cover. As the glass covers provided in solar air heaters affect the thermal performance, the increment in the number of glass covers would be the reason behind the better thermal performance. As such, smooth solar air heaters if provided with glass covers on three absorber plate sides might result in better performance. 1.4 EFFECT OF ROUGHNESS AND FLOW PARAMETERS IN SOLAR AIR HEATERS Use of artificial roughness of different configurations has been used in plenty to enhance the heat transfer rate in solar air heaters during the last decades. Varying magnitudes of roughness result in varying values of heat transfer and friction factor enhancement. The wide range of relative roughness height covered was 0.001-0.0333. The data was correlated according to the law of wall similarity, by the friction similarity function (Webb et.al., 1971, a, b). + u+ = e e

2

2e D + 3.75 f + 2.5 ln

(1.31)

Dipprey and Sabersky, (1963) developed a heat momentum transfer analogy relation for flow in sand-grain roughened tubes. Three rough tubes were tested, having effective roughness ratios equal to 0.0488, 0.0138 and 0.0024 in the Reynolds number range of 60000 to 500000 in while the Prandtl number ranged from 1.2 to 5.95. The correlation developed was:-

7|Page

CF

2C H −1 CF

+ Af = g e+, Pr

(1.32)

2

Based on the approach considered by Han, (1984), analysis for the effect of artificial roughness was made by Prasad and Saini, (1988), for heat transfer and friction factor in a solar air heater using small diameter wire artificial roughness on the absorber surface, having relative roughness pitch of 10, 15, 20 and relative roughness height of 0.020, 0.027 and 0.033 to predict for a correlation for the average Nusselt number written under as: f/2

Nu = 1+

Re Pr

p e

(1.33)

f/2 4.5(e + )0.28 Pr 0.57 −0.95( )0.53

(Han et al., 1991) studied the effect of parallel and V-shaped staggered discrete ribs and found out that 600 staggered discrete V-shaped ribs provide higher heat transfer than that for the parallel discrete ribs. Experiments were conducted by (Lau et al., 1991 a, b, c) to study the effect of discrete V-shaped ribs on turbulent heat transfer and friction for fully developed flow of air in a square channel. They found that the average Stanton number for the inclined 450 and 600 ribs was 20-35% higher than in the 900 full rib cases. (Gupta et al, 1993) used transverse wire roughness on the top surface in a solar air heater to investigate for the effect of solar air heater duct aspect ratio and relative roughness height for a relative roughness pitch of 10 and flow Reynolds number of 3000-18000 to arrive at the following correlations: e

W

for e+ < 35, Nur = 0.000824(D )−0.178 ( H )0.288 (Re)1.62 e

W

for e+ ≥ 35, Nur = 0.00307(D )0.469 ( H )0.245 (Re)0.812

(1.34) (1.35)

Studies conducted by Saini and Saini, (1997) deals with the effect of expanded wire mesh geometry parameter of relative long way length of mesh, relative short way length of mesh and relative roughness height of mesh on heat transfer. Comparison of thermal performance of roughened absorber plate fixed with staggered discrete V-apex (up and down) was done by (Muluwork et al., 1998, 2000). The increase of relative roughness length ratio in the range of 3-7, increases Stanton number. Boosting in Stanton number ratio was found to be of order of 1.32-2.47. (Karwa et. al., 1999) developed heat transfer coefficient and friction factor correlations in rib-roughened solar air heater duct for transitional flow. 8|Page

Verma and Prasad, (2000) developed the correlation for heat transfer in a top side artificially roughened solar air heater for fully developed turbulent flow as under: p

e

p

e

Nur = 0.08596(e )−0.054 (D )0.072 (Re)0.728 , for

e+ ≤ 24

(1.36)

Nur = 0.02954(e )−0.016 (D )0.021 (Re)0.802 , for e+ > 24

(1.37)

(Bhagoria et al., 2002) used wedge shaped transverse repeated rib roughness on one broad heated wall of solar air heater duct and generated data pertinent to friction and heat transfer. They analyzed that the presence of wedge shape ribs yield maximum enhancement in Nusselt number, about 2.4 times as compared to smooth duct. Nusselt number increases and attains maximum value at a wedge angle of about 100 and then sharply decreases with increasing wedge angle beyond 100. The correlations developed were:−4

e

1.21

Nu = 1.89 × 10 Re exp

−1.5 ln

φ

0.426

D

p

2.94

e

× exp

−0.71 ln

p

2

e

φ

−0.018

10

×

2

10

f = 12.44 e D

(1.38) 0.99

φ

0.49

10

Re−0.18

p

−0.52

(1.39)

e

(Momin et al., 2002) experimented on flow through duct roughened with V-shape ribs attached to the underside of one broad wall of the duct, to collect data on heat transfer and fluid flow characteristics. They observed that the Nusselt number increases with an increase of Reynolds number. It was found that for relative roughness height of 0.034 and for angle of attack of 600, the V-shaped ribs enhance the values of Nusselt number by 1.14 and 2.3 times respectively over inclined ribs and smooth plate case at Reynolds number of 17034 and arrived the correlations as:Nu = 0.067Re0.888 e D f = 6.266 e D

0.565

α

0.424

α

−0.093

60

−0.077

60

× exp

2 −0.782 ln α 60

2 Re−0.425 exp[−0.719ln α 60 ]

(1.40) (1.41)

Karwa, (2003) had conducted experimental study on heat transfer and friction characteristics in a high aspect ratio duct with transverse, inclined, V-up continuous and V-down 9|Page

continuous, V-up discrete ribs and V-down discrete ribs. He investigated that enhancement in Stanton number over smooth duct was found to be 65-90%, 87-112%, 102-137%, 110147%, 93-134% and 102-142% respectively. The friction factor ratio for V-up continuous and V-down continuous, V-up discrete ribs and V-down discrete ribs was found to be 3.92, 3.65, 2.47 and 2.58 respectively and the correlations arrived are: G = 32.26e−0.006 W/H

0.5

p/e

2.56

−0.4

p/e

exp 0.7343 ln p/e

2

e+

−0.08

(1.42)

for 7 ≤ e+ ≤ 20 R = 1.66 e+

−0.0078

W/H

2.695

exp −0.762 ln p/e

2

e+

−0.075

(1.43)

for 20 ≤ e+ ≤ 60 R = 1.325 e+

−0.0078

W/H

−0.4

p/e

2.695

exp −0.762 ln p/e

2

(1.44)

(Cho et al., 2003) experimentally investigated the effect of gap in the inclined ribs on heat transfer in square duct with rib to pitch height ratio of 8 and angle of attack of 600.Experiments were carried out by maintaining the gap width same as the rib width and by varying the gap position over the duct width for parallel and cross rib arrangement on two opposite walls, they investigated that, the inclined rib with a downstream gap shows significant enhancement in heat transfer as compared to that of continuous inclined rib arrangement. Sahu and Bhagoria, (2005) investigated effect of 900 broken transverse ribs on heat and fluid flow characteristics using roughness height of 1.5 mm, duct aspect ratio value of 8, pitch in the range of 10-30 mm and Reynolds number in the range of 3000-12000. Heat transfer enhancement was reported to be 1.25-1.4 times over smooth duct. From the experiment maximum thermal efficiency was found in the order of 83.5%. The effect of relative roughness pitch, relative roughness height and relative groove position on the heat transfer coefficient and friction factor has been studied experimentally by (Jaurkur et al., 2006) using rib-grooved roughness and the correlations arrived are:

10 | P a g e

Nu = 0.00206Re0.936 e/D p

2

0.349

2

3.318

exp −0.868 ln p/e

2

g/p

1.108

exp 2.486 ln g/

3

+ 1.405 ln g/p

(1.45)

f = 0.0012Re−0.199 e/D p

p/e

0.585

p/e

7.19

exp −1.854 ln p/e

3

+ 0.8662 ln g/p

2

g/p

0.645

exp 1.513 ln g/ (1.46)

CFD investigation using ten different ribs namely rectangular, square, chamfered, triangular etc. is done by (Chaube et al., 2006). FLUENT CFD code is used for analysis using the SST k-ἐ turbulence model. The heat flux of 1100 W/m2 is provided on the absorber plate only and rib surface is kept adiabatic. Their study reported that highest heat transfer is achieved with chamfered ribs and the best performance index is found with rectangular ribs of size 3×5. (Karmare et al., 2007) experimentally study the effect of heat transfer performance of a rectangular duct with metal grit ribs as roughness element, employed on one broad wall, transferring heat to the air flowing through it. The effect of relative length, height and pitch of the metal grid ribs on the heat transfer and friction factor has been studied for the flow range of Reynolds numbers 4000-17000. It has been found that maximum heat transfer rate was reported for the set of parameters and maximum friction factor was observed for the set of parameters. Correlations for Nusselt number and friction factor were developed based on the experimental results and were given by: Nu = 2.4Re1/3 e/D

0.42

f = 15.55Re−0.26 e/D

l/s

0.94

−0.146

l/s

p/e

−0.27

−0.27

p/e

−0.51

(1.47) (1.48)

(Varun et al., 2008) experimentally used a combination of transverse and inclined ribs as roughness geometry and examined the thermal performance for the range of Reynolds number (Re), 2000-14000, pitch of ribs (p), 5-13mm, roughness height (e), 1.6mm and aspect ratio (W/H), of 10. Results show that the collector roughened with this type of roughness provides best performance at relative roughness pitch (p/e) of 8 and the correlations developed were:-

11 | P a g e

p f = 1.0858(Re)−0.3685 ( e)0.0114

(1.49)

p Nu = 0.0006(Re)1.213 ( e)0.0104

(1.50)

(Ahrwal et al., 2008) studied the effect of width and position of gap in inclined split-ribs having square cross section on heat transfer and friction characteristics of a rectangular air heater duct. The increase in Nusselt number and friction factor was in the range of 1.48-2.59 times and 2.26-2.9 times of the smooth duct respectively for the range of Reynolds numbers from 3000-18000 and the correlations arrived are: Nu = 0.0102Re1.148 e/D f = 0.5Re−0.0836 e/D

0.51

1 − 0.25 − d/w 2 0.01 1 − g/e

2

0.72

(1.51) (1.52)

Experimental study was performed by (Huang et al., 2008) for a measurement of heat transfer coefficients in a square channel with a perforation baffle on Nusselt numbers. Local heat transfer enhancement at the downstream of step baffle was found to be more significant when Reynolds number was large and when baffle height was higher. Heat transfer coefficient off centre was found to be better because of secondary flows that appeared off centre after the air flow passes through the baffle. The heat transfer enhancement in case of a baffle with holes was reported greater than that without holes. Karwa and Maheshwari, (2009) made an experimental study of heat transfer and friction in a rectangular section duct with fully perforated baffles or half perforated baffles at relative roughness pitch of 7.2-28.8 affixed to one of the broader walls. They investigated that there was enhancement of 79-169% in Nusselt number over the smooth duct for the fully perforated baffles and 133-274% for the half perforated baffles while the friction factor for fully perforated baffles was found to be 2.98-8.02 times of that for the smooth duct and 4.42-17.5 times for the half perforated baffles. The correlations arrived are: Nu = 0.0893 Re0.7608

(1.53)

f = 0.1673 Re−0.0213

(1.54)

12 | P a g e

An experimental investigation had been carried out by Bopche and Tandale, (2009) to study the heat transfer coefficient and friction factor by using artificially roughness in the form of inverted U-shaped turbulators, on the absorber surface of an air heater duct. As compared to smooth one, the enhancement in heat transfer and friction factor was reported of the order of 2.82 and 3.72 times, respectively and the correlations developed are: Nu = 0.5429Re0.7054 e/D f = 1.2134Re−0.2076 e/D

0.3619

0.3285

p/e p/e

−0.1592

−0.4259

(1.55) (1.56)

A 3D solar air heater duct provided with artificial roughness in the form of thin circular wire in arc shaped has been analyzed by (Kumar et al., 2009), using CFD. Renormalization k-ἐ turbulence model is used for analysis in FLUENT CFD code. Their study reveals the maximum value of overall enhancement ratio as 1.7 for the range of parameters investigated. A parametric study of artificial roughness geometry of expanded metal mesh type in the absorber plate of solar air heater duct has been carried out by (Gupta et al. 2009). The performance evaluation in terms of energy augmentation ratio (EAR), effective energy augmentation ratio (EEAR) and exergy augmentation ratio (EXAR) has been carried out for various values of Reynolds number and roughness parameters of expanded metal mesh roughness geometry in the absorber plate of solar air heater duct. It was investigated that the augmentation ratios decrease at faster rate with Reynolds number in order of EAR, EEAR and EXAR. It was also studied that augmentation ratios increase with increase in duct depth and intensity of solar radiation. Thermal performance of wire-screen metal mesh was experimentally studied by (Paswan et al., 2009) as artificial roughness geometry on the underside of absorber plate. The thermal performance enhancement of a solar air heater was depended strongly upon the diameter of the wire and pitch, mass flow rate, insulation and inlet temperature. It was investigated that decreasing wire pitch, thermal performance of a solar air heater increases. Sriromrein and Promvong, (2010) studied that heat transfer and friction characteristic of a rectangular duct roughened artificially with Z-shaped ribs. The enhancement in heat transfer rate and best thermal performance was reported for Z-rib inclined at 45º. (Hans et al., 2010) developed multiple V-rib roughness and performed extensive experimentation to collect data 13 | P a g e

on heat transfer and fluid flow characteristics of a roughened duct. Maximum enhancement in Nusselt number and friction factor due to presence of multi v-rib roughness has been found to be 6 and 5 times, respectively, in comparison to the smooth duct. The maximum enhancement in heat transfer has been achieved corresponding to relative roughness width (W/w) of 6. It has been observed that Nusselt number attain maximum value at 600angle of attack. The correlations proposed were:Nu = 3.35 × 10−5 Re0.92 e D exp

2 −0.61 ln α 90

p

8.54

e

0.77

× exp

W

0.43

w

−2.0407 ln

α p

−0.49

90

× exp −0.1177

2 ln W w

×

2

e

(1.57)

f= 4.47 × 10−4 Re−0.3188 e D exp −0.52 ln

α 90

2

p

0.73 8.9

e

W

0.22

w

α

−0.39

90

exp −2.133 ln

p

× 2

e

(1.58)

(Champookham et al., 2010) carried out an experimental investigation to study the effect of combined wedge ribs and winglet type vortex generators on heat transfer and friction loss behaviors for turbulent air flow through a constant heat flux channel. Experimental investigations of heat transfer, flow friction and thermal performance factor characteristics in a tube fitted with delta winglet twisted tape, using water as working fluid was done by Eiamsa-ard et al., (2010). In the study, they described the influences of the oblique deltawinglet twisted tape and straight delta-winglet twisted tape arrangements and found that the oblique delta-winglet twisted tape is more effective turbulators giving higher heat transfer coefficient than the straight delta-winglet twisted tape. A numerical investigation had been carried out by Kwankaomeng et al., (2010) to study laminar flow and heat transfer characteristics in a 3D isothermal wall square-channel with 300 staggered angled-baffles. The studied the effect of different baffle heights at a single pitch ratio of 3 on heat transfer and pressure loss in the channel. The vortex flow was created by using the 300 angled baffle cavities leading to drastic increase in heat transfer in the channel. The order of enhancement is about 100-650% for using the angled baffles with blockage ratio of 0.1-0.3. The heat transfer augmentation was associated with enlarged pressure loss ranging from 1 to 17 times above the smooth channel. A thermal enhancement factor for the inclined baffles at 14 | P a g e

blockage ratio of 0.15 was found to be highest about 2.9. (Karmare et al., 2010) had carried out analysis using CFD on solar air heater with a absorber plate roughened with metal ribs of circular, square and triangular cross-section having 540, 560, 580,600 and 620 inclinations to the air flow. The system and operating parameters studied were: e/D = 0.044, p/e = 17.5 and l/ s = 1.72. Range of Reynolds number was 3600-17000. Their study reported that maximum heat transfer is obtained with the square cross-section ribs 580 angle of attack and there is 30% enhancement in the heat transfer for square plate over smooth surface. A numerical investigation had been carried out by Promvong et al., (2010) to examine the laminar flow and heat transfer characteristics in a three dimensional isothermal wall square channel with 450 angled baffles. In order to generate a pair of mainstream wise vortex flows through the tested section, baffles with an attack angle of 450 were mounted in tandem and in-line arrangement on the lower and upper walls of the channel. Effects of different baffle heights on heat transfer and pressure loss in the channel were studied and results of the 45 0 in-line baffle were also compared with those of the 900 transverse baffle and the 450 staggered baffle. The effect of wedge-shaped transverse ribs roughness on flow through a rectangular duct of a solar air heater is studied using CFD by Gandhi et al., (2010). A 2D CFD simulation is carried out using FLUENT on a duct having relative roughness pitch of 4.5, relative roughness height of 0.022 and rib wedge angle of 150. The results of numerical analysis are in good agreement with the experimental results in terms of velocity profile and turbulence intensity. To study the effect on the thermal performance of the ribbed channel, the parallel ribs have been installed either onto one wall of the channel or, in-line, onto two opposite walls with a rib pitch-to-height ratio ranging from 6.66-20.0. (Kumar et al., 2011) have investigated for the enhancement of heat transfer coefficient of a solar air heater having roughened air duct with a provision of artificial roughness in the form of 600 inclined discrete ribs. For the relative roughness pitch of 12, relative gap position of 0.35 and relative roughness height of 0.0498, the maximum heat transfer enhancement takes place. Statistical correlations for friction factor and Nusselt number have been derived as a function of gap position, rib depth, pitch and Reynolds number. The correlations have been established to predict the values of friction factor and Nusselt number with an average absolute standard deviation of 3.4% and 3.8% respectively and have been given by:15 | P a g e

Nu = 3 × 10−5 Re0.947 e D f = 0.0014Re−0.23 e D

p

0.29

0.804

p

5.885

e 4.516

e

d

d

0.115

W 0.097

W

× exp −1.237 ln

× exp −0.944 ln

p

p

2

(1.59)

e 2

(1.60)

e

An experimental investigation for heat transfer and friction factor was carried out by (Lanjewar et al., 2011), having the characteristics of a rectangular duct roughened W-shaped ribs, that is arranged at an inclination with respect to the flow direction. The duct used, had a width to height ratio of 8, relative roughness pitch of 10, relative roughness height of 0.03375 and angle of attack of flow of 300-750 and the correlations proposed were:R=

2

2e D + 3.75 f + 2.5ln

(1.61)

e+ = √ f 2 Re e D g′ =

f

(1.62) 2

2St − 1

f +R

(1.63)

An experimental investigation was conducted by (Bhushan et al., 2011) for a range of system and operating parameters in order to analyze the effect of artificial roughness on heat transfer and friction in solar air heater duct having protrusions as roughness geometry. Maximum enhancement of Nusselt number and friction factor has been found 3.8 and 2.2 times respectively in comparison to smooth duct for the investigated range of parameters. Maximum enhancement in heat transfer coefficient has been found to occur for relative short way length (S/e) of 31.25, relative long way length (L/e) of 31.25 and relative print diameter (d/D) of 0.294. The correlations developed are: Nu = 2.1 × 10−88 Re1.452 S/e D

−3.9

12.94

L/e

99.2

d/

exp −10.4{log S/e }2 exp −77.2{log

f = 2.32Re−0.201 S/e

−0.383

L/e

−0.484

d/D

L e

}2 exp −7.83{ln

0.133

d D

}2

(1.64) (1.65)

16 | P a g e

(Pomovonge et al., 2011) experimentally studied the effects of combined ribs and deltawinglet type vortex generators (DWs) on forced convection heat transfer and friction loss behaviors for turbulent air flow through a solar air heater channel. Rectangular channel with aspect ratio of 10, height of 30 mm with Reynolds number in the range of 5000-20000 was used. They reported that combined rib and the DW provides considerable heat transfer 𝑁𝑢

augmentations, 𝑁𝑢 = 2.3 − 2.6 and also causes a moderate pressure drop increase, 𝑓/𝑓0 = 0

4.7 − 10.1, depending on the attack angle and Re values. An experimental study has been carried out by (Yadav et al., 2012) to study the effect of heat transfer and friction characteristics of turbulent flow of air passing through rectangular duct which is roughened by circular protrusions arranged in angular arc fashion. The investigation reported the maximum enhancement in heat transfer and friction factor as 2.89 and 2.93 times as compared to smooth duct. The correlations arrived are: Nu = 0.154Re1.017 p/e f = 7.207Re−0.56 p/e

−0.38

−0.18

e/D

e/D

0.521

0.176

α/60

α/60

−0.213

0.038

exp −2.023 ln α/60

exp 1.412 ln α/60

2

2

(1.66) (1.67)

CFD analysis on solar air heater is done by (Sharma et al., 2012) using V-shaped ribs roughness on the underside of the absorber plate. Finite volume method with semi-implicit method for pressure linked equation algorithm is used for computations. Reynolds number range was from 5000-15000. The results obtained show that the combined effect of swirling motion, detachment and reattachment of the fluid were responsible for the increase of heat transfer rate during CFD analysis. Nusselt number increases and friction factor decreases with increase in Reynolds number for all combination of relative roughness height (e/D) and relative roughness pitch (p/e). An average percentage deviation predicted between CFD and exact solution was found less than ±3% V-shaped rib roughness found to give high rate of heat transfer. A numerical investigation on laminar flow and heat transfer characteristics in a 3D isothermal wall square-channel fitted with in-line 450 V-shaped baffles on two opposite wall was done by (Promvong et al., 2012). The in-line V-baffles with its V-tip pointing downstream and the angle of attack of 450 relative to the flow direction were mounted repeatedly on the lower and upper walls. The V-baffled channel flow was found to be fully developed periodic flow and heat transfer profiles at about x/D=8 downstream of the tested 17 | P a g e

channel inlet. It was reported that the P-vortex flow caused by the V-baffle induced impingement flows on the channel walls leading to drastic increase in the heat transfer rate. The heat transfer in the V-baffled channel was found to be about 1-21 times higher than the smooth channel with no baffle. Experimental investigation has been carried out by (Chauhan et al., 2013) to study heat transfer and friction factor characteristics using impinging jets in solar air heaters duct. The study shows that there is considerable enhancement in heat transfer and friction factor by 2.67 and 3.5 times respectively and the correlations arrived were: Nu = 1.658 × 10−3 Re0.8512 X/Dh

0.1761

Y/Dh

0.141

Dj /Dh

−1.9854

exp⁡−0.3498{ln⁡Dj /

Dh }2

f=

(1.68)

0.3475Re−0.5244

X

0.4169

Dh

Y

D j −1.4848

0.5321

Dh

Dh

exp −0.22210{ln Dj /Dh }

(1.69)

2

Prasad, (2013) represented the experimental results on heat transfer and thereby thermal performance of artificially roughened solar air heaters for fully developed turbulent flow. Such solar air heaters have been found to give considerably high value of collector heat removal factor (FR ), collector efficiency factor (F ′ ) and thermal efficiency (ƞth ) as compared to the corresponding values of those of smooth collectors. In the range of the operating parameters investigated, the ratio of the respective values of the parametersFR , F ′ and ƞth for the roughened collectors to the smooth collectors have been found to be 1.786, 1.806 and 1.842 respectively and the correlations arrived were:-

Nu = f 2 RePr/ 1 +

f

2

2

f= [0.95

p

0.53

e

+2.5 ln D 2e −3.75]2

4.5 e+

0.28

Pr 0.57 − 0.95

p

0.53

e

(1.70)

(1.71)

A CFD based investigation of turbulent flow through a solar air heater is carried out by Yadav and Bhagoria, (2013, a) with square sectioned transverse rib roughness. The 2D 18 | P a g e

analysis of rectangular duct is performed using ANSYS FLUENT 12.1 software and using (RNG k-ἐ) turbulence model. Analysis using four different configurations of rib roughness keeping relative roughness pitch constant at p/e = 14.29 and six different values of Reynolds number, ranging from 3800-18000, revels that the relative roughness height is a vital factor and mainly affects the rate of heat transfer rate and the increase in flow friction. An investigation has been conducted by Yadav and Bhagoria, (2013, b) by using FLUENT 12.1 and (RNG k-ἐ) turbulence model on solar air heater, using square sectioned transverse rib roughness, with different values of relative roughness pitch (7.14 ≤ p/e ≤ 17.86). 2D CFD analysis to predict the effect of small diameter of transverse wire rib roughness on heat transfer enhancement in solar air heater is done by Yadav and Bhagoria, (2013, c). In this analysis also FLUENT 12.1 software is used with (RNG k-ἐ) turbulence model. (Prasad et al., 2014) have analyzed with respect to fluid flow and heat transfer in a novel solar air heater having artificial roughness on three sides (the two side walls and the top side) of the rectangular solar air heater duct, with three sides glass covers. Equations for friction factor and heat transfer parameter have been developed. The analytical values of friction factor and heat transfer parameter have been found to be 2 to 40% more and 20 to 75% more than those of the respective values of (Prasad and Saini, 1988) for the same range of the values of operating parameters p/e, e/D and Re and fixed values of W and B. The correlations developed are:f r /2

Nur = 1+

p f r /2 4.5(e + )0.28 Pr 0.57 −0.95( e )0.53

2

W +2B 0.95

fr =

p

0.53 e

2 +2.5 ln D 2e −3.75

2 W +B

Re Pr

(1.72)

+W f s

(1.73)

Behura et al., (2016), investigated experimentally the heat transfer, friction factor and thermal performance of a novel type of three sides artificially roughened and three sides glass covered solar air heaters under fully developed turbulent flow conditions and compare well with analytical values of Prasad et al., (2014). Theoretical investigations on thermal and thermohydraulic performance of wavy finned absorber plate solar air heaters have been 19 | P a g e

reported recently by Priyam and Chand (2016).Chouksey and Sharma, (2016), investigated for the thermal performance characteristics of solar air heater having its duct with blackened wire screen matrices. However, in all the above cases, provision of artificial roughness and glass cover has remained limited to only one side (top side) of the solar air heater duct except those of the recent ones (Prasad et al., 2014 and Behura et al., 2016), wherein it has been concluded that three sides roughened and glass covered solar air heaters perform even better than those of one side roughened and glass covered solar air heaters, but friction factor also increases. 1.4.1

Optimization of the System Parameters As the enhancement of heat transfer coefficient is attempted by providing artificially roughness, it always goes with an increment of pressure drop and the requirement of pumping power is increased. So, there is a need to optimize the system parameters to maximize heat transfer while keeping friction losses as low as possible. Despite plenty of works dealing with the effect of artificial roughness on heat transfer and friction factor in solar air heaters, very few dealt with the thermo hydraulic performance in 𝑒

them. The thermo hydraulic optimization parameter, 𝑒 + = 𝐷

𝑓 2

𝑅𝑒, known as roughness

Reynolds number, has been considered to arrive at the optimal thermo hydraulic performance condition. Covering a wide range of the values of heat transfer surface area, overall heat conductance and flow friction power, (Webb and Eckert, 1971), arrived at the conclusion that the value of the parameter, 𝑒 + = 20, gives the optimal thermo hydraulic performance. For optimal thermo hydraulic performance of circular tube roughened with ribs, Lewis, (1975), introduced new efficiency parameters (𝐿−1 , 𝐵 −1 , 𝐶 −1 , 𝐺𝐻 𝑎𝑛𝑑 𝑅𝑀 ) and arrived at the conclusion that the value of the roughness Reynolds number, 𝑕+ = 𝑒 + = 20, corresponds to the optimal thermo hydraulic condition. Sheriff and Gumbley, (1966), has studied for annulus with wire type roughness and found the value of roughness Reynolds number, 𝑒 + = 35, for the optimal thermo hydraulic condition. Optimal thermo hydraulic performance analysis in one side artificially roughened solar air heater of (Prasad and Saini, 1991), has obtained a particular set of values of roughness and flow parameters to give the value of roughness Reynolds number, 𝑒 + = 24, for optimal thermo hydraulic condition. 20 | P a g e

Verma et al., (1991) have given a technical note on ‘optimization of solar air heaters of different designs’, to obtain the optimum flow channel depth and mass flow rate in 10 different designs of solar air heaters and have found that a single glazing solar air heater operating under double flow configuration gives the best performance. Hegazy, (1995) has developed an analytical criterion to determine the optimal flow channel depth of conventional flat plate solar air heaters. Gupta et al., (1997) have shown the effect of roughness parameters on the thermal as well as the thermo hydraulic performance of roughened solar air heaters to design the optimal operating conditions. A technical note on ‘optimal thermo hydraulic performance of flat plate solar air heaters, operating with fixed or variable pumping power’, has been given by Hegazy, (1999). Optimal thermo hydraulic performance of solar air heaters has been investigated (Verma and Prasad, 2000), for the maximum heat transfer and minimum pressure drop to arrive at the conclusion that the value of 𝑒 + = 24,corresponds to the optimal thermo hydraulic performance. Saha and Mahanta, (2001) have investigated for the performance of second-law modeling of flat plate collectors with the objective of describing collector performance when entropy generation is minimum. Karwa et al., (2001) have investigated for thermo hydraulic performance of solar air heaters with chamfered repeated rib-roughness on the air flow side of the absorber plates for the range of flow Reynolds number, 3750-16350. Mittal and Varshney, (2005), has worked on optimal thermo hydraulic performance of wire mesh packed solar heater. Second law optimization of solar air heater having chamfered rib groove as a roughness element has been analyzed by Layek, et al., (2007). Karmare and Tikekar, (2008), have optimized the thermo hydraulic performance of solar air heater integrated with metal rib as roughness element. Karmare and Tikekar, (2009) have investigated for the effect of metal rib grifts on absorber plate at 600 of angle of attack to the flow direction of air on the thermo hydraulic performance of solar air heater. Siddhartha et al., (2012) have undertaken the work with an objective to develop and implement a trained particle swarm optimization algorithm for prediction of an optimized set of design and operating parameters for a smooth flat plate solar air heater. Yadav and Bhagoria, (2014), investigated numerically a CFD based thermo hydraulic performance analysis of an artificially roughened solar air heater having equilateral triangular sectioned rib roughness on the absorber plate. A numerical investigation is conducted (Yadav and Bhagoria, 2014) to analyze the two-dimensional 21 | P a g e

incompressible Navier–Stokes flows through the artificially roughened solar air heater for relevant Reynolds number range of 3800 to 18,000and twelve different configurations of equilateral

triangular

sectioned

rib

roughness (𝑝/𝑒 = 7.14– 35.71 and 𝑒/𝐷 =

0.021– 0.042). Sharma and Kalamkar, (2015), reviewed the thermo hydraulic performance of solar air heater having artificial roughness elements and reported the various artificial roughness geometries used by investigators to enhance the heat transfer, which is attained by increasing the pressure drop across the surfaces of the absorber plate of the solar air heater and its effect on heat transfer and friction factor through experiments.

1.5 OBJECTIVES OF THE PRESENT WORK All the above mentioned works on thermo hydraulic and optimal thermo hydraulic performance deal with one side roughened solar air heaters. Three sides roughened solar air heater of Prasad et al, (2014) and Behura et al, (2016), have been reported to be even more efficient than those of one side roughened solar air heater. Therefore, optimal thermohydraulic performance analysis and investigation of three sides roughened solar air heaters is required and the present work has been taken up with the following objectives: 1) To analyze for the optimal thermo hydraulic performance of the three sides artificially roughened solar air heater for maximizing heat transfer while keeping friction loss minimum. 2) To develop the experimental set-up to perform experiments under actual outdoor conditions for data collection with a view towards validation of the optimal thermo hydraulic performance analysis for such solar air heater. 3) To obtain the results with respect to Nusselt number, Stanton number and friction factor from the data collected in such solar air heater. 4) To derive the correlations for heat transfer and friction factor in such solar air heaters. 5) To represent and discuss the effect of roughness and flow parameters on heat transfer, friction factor, collector performance parameters and thermo hydraulic performance in such solar air heater. 6) To obtain and represent the thermal efficiency curves for such solar air heaters. 7) To obtain the optimal thermo hydraulic performance results and represent the optimal thermo hydraulic performance curves for such solar air heater. 22 | P a g e

CHAPTER - 2 OPTIMAL THERMOHYDRAULIC PERFORMANCE ANALYSIS

2.1 INTRODUCTION It has already been pointed out in the previous chapter that the provision of artificial roughness on the underside of the absorber plate of solar air heaters leads to the enhancement of heat transfer coefficient and friction losses. In view of this, there is a need to critically examine the quantitative influence of the roughness and flow parameters to attain the optimal thermo hydraulic performance condition, while keeping the pumping power as low as possible. The fundamental parameters which control the momentum and heat transfer for a roughened surface are momentum transfer function (RM) and heat transfer roughness function (GH). The present chapter has been devoted to the analysis of optimization aspects of three sides roughened and glass covered duct system that would possibly lead to the quantitative evaluation of various physical and flow parameters which would help in designing such solar air heaters. The effect of various efficiency parameters have been compared to those of the referred ones. 2.2 OPTIMIZATION METHODOLGY The solar air heater duct model of width, W, height, B and length, L, as shown in Fig 2.1a with provision of roughness elements at pitch, p of height, e on the absorber plates has been considered for analysis. The methodology of optimization adopted is similar to those of (Lewis, 1975 a, b; Prasad and Saini, 1991), that a dimensionless parameter, e+ (roughness Reynolds number), could plausibly form the basis to analyze the optimal thermo hydraulic performance of three sides artificially roughened solar air heaters. The correlations developed (Prasad et al., 2014), for friction factor and heat transfer written here under as Eqs. (2.1) and (2.2) respectively has been utilized. Fig. 2.1(b) shows typically the pitch and height of the roughened elements.

23 | P a g e

f3r = St 3r =

W +2B f r +W f s 2(W +B)

(2.1) f 3r /2

f 1+√( 3r )[4.5 2

p 0.53 e + 0.28 Pr 0.57 −0.95 ]

(2.2)

e

24 | P a g e

where, W/B>>1, and fr is as that of (Webb et al., 1971), considered by (Prasad and Saini, 1991), written under as Eq. (2.3) 2

fr = 0.95

p 0.53 D +2.5 ln e 2e

(2.3)

2

−3.75

Consequently, substituting forfr , Eq. (2.1) could be written to be Eq. (2.4) as under: 2 2 p 0.53 D 0.95 +2.5 ln −3.75 e 2e

W +2B

f3r =

+W f s

(2.4)

2 W +B

Now,

substituting

for

the

parameters,

heat

transfer

roughness

function, GH = 4.5(e+)0.28 pr0.57 and momentum transfer roughness function, R M = 0.95 P e St 3r =

0.53

, Eq. (2.2) could be written as Eq. (2.5) under:

f 3r /2

1

f 1+√( 3r ) 2

G H −R M

where, e+ = e D

f 3r 2

(2.5)

Re

(2.6)

The optimization parameter, Ƞ,for the present case of three sides roughened solar air heater is defined as (Lewis, 1975; Prasad and Saini, 1991) S t 3r

Ƞ=

3

S t 3gS f 3r

= f 3r , GH , R M

(2.7)

fs

which can further be written as: Ƞ = Ƞ f 3r , B−1 , C −1

(2.8)

and could further also be represented as: Ƞ = Ƞ f 3r , L−1

(2.9)

25 | P a g e

where, B−1 = GH − Pr t R M = 4.5(e+)0.28 Pr 0.57 − Pr t × 0.95(p/e)0.53 C−1 = 2.5ln⁡ (e+) + 5.5 − R M = 2.5ln⁡ (e+)5.5 − 0.95

p

(2.10)

0.53

(2.11)

e

L−1 = C−1 − B−1 = 2.5ln⁡ (e+) + 5.5 − GH − R M 1 − Pr t

(2.12)

Also, that, the following equations for smooth solar air heater are used in Eq. (2.7) fS = 0.079Re−0.25

(2.13)

St s = 0.023Re0.2 pr −2/3

(2.14)

The values of fr , fs and St s have been worked out by substituting the values of geometrical and operating parameters W, B, p/e, e/D and Re in Eqs. (2.3), (2.13) and (2.14) respectively, which have further been utilized to work out for the values of f 3r , St r and Ƞ, by the respective equations. 2.3 EFFECT OF EFFICIENCY PARAMETERS The values of the parameters, e+, B-1, C-1, and L-1 have also been worked out by the respective equations, taking the value of the turbulent Prandtl number, Pr t , equal to 0.90. The values of the parameters as above have been worked out, assuming that W/B >>1, p/e ≥ 10 and e/D, so selected that the roughness height is of the order of or slightly greater than the viscous sub-layer thickness of the hydrodynamic boundary layer over a flat plate. The values of flow Reynolds number vary from 4000 to 20,000. Since, the parameter, e+, has remained the optimization parameter in literature, Figs. 2.2, 2.3 & 2.4 have been drawn to see the effect of e+ on the parameters C-1, B-1, and L-1, for the present analysis and that of (Prasad and Saini, 1991). It could be seen from Figs. 2.2 & 2.3 that the values of the parameters B-1 and C-1 increase with increasing values of e+.

. 26 | P a g e

(1) p/e=10(Preasent case) (2) p/e=10(Prasad and Saini,1991) (3) p/e=20(Present case) (4) p/e=20(Prasad and Saini,1991) (5) p/e=30(Present case) (6) p/e=30(Prasad and Saini,1991)

Re  4000  20000 e/D=0.0225 (Preseant Case) e/D=0.0205 (Prasad and Saini,1991)

(1)

10 (3)

8

(2)

(5)

(4)

6

-1

B 4

(6)

2

0 5

10 15 20 25 30 35 40 45 50 55 60 65 70

e   e / D f 3r / 2 Re +

-1

Fig. 2.2. Effect of e on B for varying p/e

27 | P a g e

Figs. 2.2 & 2.3 also show the curves for one side artificially roughened solar air heater, based on the heat transfer as friction factor correlations developed (Prasad and Saini, 1988), and written under as Eqs. (2.15) and (2.16): f=

St =

W +2B f S +W f r

(2.15)

2(W +B) f/2

(2.16)

p 0.53 1+ f 2[4.5 e + 0.28 Pr 0.57 −0.95 ] e

It is clear from Figs. 2.2 & 2.3 that values of both the parameters, B-1 and C-1 are higher in the case of three sides artificially roughened solar air heater than those of only one side artificially roughened solar air heater (Prasad and Saini, 1988). Fig. 2.4 shows that the value of the efficiency parameter L-1,increases with the increasing value of e+, up to certain maximum value when after its value decreases with increasing value of e +, for every curve drawn for varying values of the relative roughness pitch, p/e, equal to 10, 20 and 30. It appears that the maximum values of the parameter, L-1, could be expected between 20-30. To ascertain the particular value of e+, Table-2.1 and Fig. 2.5 have been prepared, which show that e+ = 23 gives the maximum value of the efficiency parameter, L-1, for all values of p/e. Therefore, the optimal thermo hydraulic performance equation for three sides artificially roughened solar air heater could be written as: e e+ opt = D

f 3r 2

Re = 23

Hence, the roughness Reynolds number, e+ = e D

(2.17) f 3r 2

Re constituting a specific set of

values of roughness and flow parameters, p/e, e/D and Re to obtaine+ opt = 23, could result in the maximum heat transfer for minimum pressure drop. The set of values of p/e, e/D and Re, vary with variance in either of the values of p/e, e/D and Re, individually or combined, which may form the basis for preparing design curves to give the optimal thermo hydraulic performance of three sided artificially roughened solar air heaters in a similar way as those of (Prasad and Saini, 1991).

28 | P a g e

Re  4000  20000

(1) p/e=10(Present Case) (2) p/e=10(Prasad and Saini,1991) (3) p/e=20(Present Case) (4) p/e=20(Prasad and Saini,1991) (5) p/e=30(Present Case) (6) p/e=30(Prasad and Saini,1991)

e/D=0.0225 (Present Case) e/D=0.0205 (Prasad and Saini,1991)

14

(1)

12

(3)

(2)

10

(4)

(5)

8 C

(6)

-1

6

4

2 5 10 15 20 25 30 35 40 45 50 55 60 65 70

e   e / D f 3r / 2 Re +

-1

Fig. 2.3. Effect of e on C for varying p/e 29 | P a g e

e/D = 0.0225 (constant for p/e= 10,20 &30) (1) p/e=10 Re  4000  20000 (2) p/e=20 4.8 (3) p/e=30

4.6 4.4 4.2 4.0 3.8 3.6 -1

L

(1)

(2)

(3)

3.4 3.2 3.0 2.8 2.6 2.4 2.2 2.0 5

10

15

20

25

30

35

40

45

50

55

60

65

e   e / D f 3r / 2 Re +

-1

Fig.2.4. Effect of e on L for varying p/e

30 | P a g e

4.2

eopt  23

e

+

23 25 22 26 24 28 30 17 34

4.1

4.0 -1

L

3.9

3.8

10

15

20

25

30

35

40

p/e -1

+

Fig.2.5. The efficiency parameter L as a function of p/e & e

31 | P a g e

Table-2.1 Range of efficiency parameter L-1 as a function of e+ and p/e p/e = 10 e+

p/e =20 L-1

e+

p/e =30 L-1

e+

L-1

17

4.143224

17

4.000322

17

3.793982

22

4.180075

22

4.037173

22

3.830833

23

4.181923

23

4.039021

23

3.832681

24

4.177088

24

4.034186

24

3.827846

25

4.181683

25

4.038781

25

3.832441

26

4.179871

26

4.036969

26

3.830629

28

4.17343

28

4.030528

28

3.824188

30

4.163822

30

4.02092

30

3.81458

34

4.137277

34

3.994375

34

3.788035

Therefore, the optimal thermo hydraulic performance equation given under as, e+ opt = e

D

f 3r 2

Re = 23, has formed the basis for preparing optimal performance curves for three

sides artificially roughened solar air heater as shown in Figs. 2.6.1, 2.6.2 and 2.6.3.

32 | P a g e

0.09

0.08

0.07

0.06

e/D

4000

0.05

5000 6000

0.04

7000 8000

0.03

0.02 10

15

20

25

30

35

40

p/e Fig.2.6.1 Optimal thermohydraulic performance curves for three sides artificially roughened solar air heater

33 | P a g e

0.030 0.028 Re=9000

0.026

10000 0.024

11000 0.022

12000 13000

0.020

e/D

14000

0.018 0.016 0.014 0.012 10

15

20

25

30

35

40

p/e Fig.2.6.2 Optimum thermohydraulic performance curves for three sides artificially roughened solar air heater

34 | P a g e

0.018

Re=15000

0.016

16000 17000 e/D

0.014

18000 19000 20000

0.012

0.010

10

15

20

25

30

35

40

p/e Fig.2.6.3 Optimum thermohydraulic performance curves for three sides artificially roughened solar air heater

35 | P a g e

Table – 2.2 shows the optimum values of roughness Reynolds number e+ opt found in literature for different roughness geometry with range of values of roughness and flow parameters for fully developed turbulent flow, with the present value of e+ opt = 23. Table-2.2 Values of 𝒆+ 𝒐𝒑𝒕 for different roughness geometries for fully developed turbulent flow Sl.No

References

Roughness and flow p/e

e/D

Re

0.018-

1×104–

0.035

2×103

0.01-

6×103–

0.04

10×103

0.02-

1×104–

0.10

1×107

Prasad and Saini, Rectangular with top 10-40

0.020-

3×103–

(1991)

0.033

20×103

0.010-

5×103– 20×103

. 1.

duct geometry Sheriff

and Annulus with wires

10

Gumley, (1966) 2.

Webb and Eckert, Rectangular with ribs

10-40

(1971) 3.

Lewis, (1975 b)

Circular

tubes

with 2-60

ribs 4.

5.

7.

+ 𝑒𝑜𝑝𝑡

side wires

Verma

and Rectangular with top 10-40

Prasad, (2000)

side wires

0.030

Present case

Rectangular with three 10-30

0.013 – 4×103–

sides wires

0.025

35

20

20

24

24

23

20×103

2.4 THERMO HYDRAULIC PERFORMANCE Fig. 2.7 compares the optimal thermo hydraulic performance curves, corresponding to the e equation e+ opt = D

f 3r 2

Re = 23 for the present solar air heater and those of (Prasad and

e Saini, 1991), corresponding to the equation, e+ opt = D

f 2

Re = 24 for given values of

Reynolds number. This figure shows that for the same value of flow Reynolds number and relative roughness pitch, p/e, values of relative roughness height,(e/D) is low for the present solar air heater.

36 | P a g e

(1) Re=4000

0.10 0.09

(2) Re=4000

0.08 0.07 +

0.06

e =24 (Prasad and Saini,1991)- 1,3,5 +

e =23 (Present case)- 2,4,6

e/D 0.05 0.04 (3) Re=10000 (4) Re=10000

0.03 0.02

(5) Re=20000 (6) Re=20000

0.01 10

15

20

25

30

35

40

p/e Fig. 2.7 optimal thermohydraulic performance curves

37 | P a g e

Fig. 2.8 shows the optimal range of relative roughness height, (e/D), for 10 ≤ p/e ≤ 30, at varying values of flow Reynolds number from 4000-20000, for the present case and that of (Prasad and Saini, 1991), which shows that the range of optimal values of (e/D),for the present case is lower than that of (Prasad and Saini, 1991), and that too decreases with increasing values of flow Reynolds number. It could also be seen that the range of the value of (e/D) , goes on decreasing for increasing values of flow Reynolds number. Table -2.3 shows the values of flow Reynolds number for varying values of e/D at given values of relative roughness pitch, p/e, equal to 10 for the present case of e+ opt = 23 and that of (Prasad and Saini, 1991), e+ opt = 24.It could be seen that for the same values of p/e and e/D, lower values of flow Reynolds number for the present case, correspond to the optimal thermo hydraulic performance. Lower value of flow Reynolds number yields in higher grade energy collection, since mass flow rate involved will be less. For three sides roughened solar air heater (Prasad et al., 2014) the value of the average friction factor, f3r will be more than that of the one side roughened solar air heater, f, of (Prasad and Saini, 1988) for comparison, the present optimal thermo hydraulic performance Eq. (2.18) is written under together with the optimal thermo hydraulic performance equation of (Prasad and Saini, 1991), as Eq. (2.19): e e+ opt = D e e+ opt = D

f 3r 2

f 2

Re = 23

Re = 24, (Prasad and Saini, 1991)

(2.18)

(2.19)

Now, it can be attributed that smaller value of e+ opt , for the present case, will lead to smaller value of flow Reynolds number to result in higher grade energy collection. Therefore, it could be concluded from the present analysis that three sides roughened solar air heaters are quantitatively even better than those of one side roughened solar air heaters of (Prasad and Saini, 1991), w.r.t. optimal thermo hydraulic performance conditions.

38 | P a g e

(e / D)opt for10  p / e  30

0.10

e/D

0.08

1

0.06

2 1

+

(e/D)optRange for e =24 (Prasad and Saini,1991) [1-1-1-1] +

(e/D)optRange for e =23 (Present case) [2-2-2-2]

2 0.04

0.02 1 2 1

2

0.00 0

Re2Re1

4000

8000

12000

16000

20000

Re Fig.2.8 .Optimal range of e/D

39 | P a g e

Table-2.3 Values of Re at varying value of e/D and p/e=10, for the optimal thermo hydraulic condition

e/D

Re1

(e+=24),

Prasad Re2 (e+=23), present case

&Saini,1991

0.01

24413.65

19901.10

0.02

12206.83

11450.553

0.03

8137.885

7300.369

0.04

6103.413

5725.277

0.05

4882.731

4780.221

0.06

4368.942

4074.184

0.07

4187.665

3904.158

0.08

3951.707

3862.638

From Fig. 2.8, it can be concluded that the range of the values of the parameter, (e/D) at lower values of flow Reynolds number is more than those at higher values of flow Reynolds number, providing more options for selecting the values of (e/D) . The optimal range of (e/D),(Prasad and Saini, 1991), represented above that of the present case, in this figure, would attribute that for the same range of the value of(e/D) , one side roughened and glass covered solar air heaters (Prasad and Saini, 1991), will have higher value of flow Reynolds number, resulting in lower grade of energy collection under optimal thermo hydraulic condition than that of three sides roughened and glass covered solar air heaters. For example, at a value of 𝑒/𝐷 = 0.06, the value of flow Reynolds number for the present case of Re2 lower than that of Re1 (Prasad and Saini, 1991). Therefore, for the similar operating condition three sides roughened and glass covered solar air heaters will perform thermo hydraulically better than those of one side roughened and glass covered solar air heaters.

40 | P a g e

CHAPTER – 3 SYSTEM DEVELOPMENT AND INVESTIGATIONS

3.1 INTRODUCTION To verify the expressions developed for heat transfer and friction factor parameters (Prasad et al, 2014) in three sides artificially roughened solar air heater, experimental results have been reported by the authors (Behura et al, 2016) with respect to the effect of roughness and flow parameters on heat transfer and friction factor. Similarly, an experimental set up has been developed to perform experiment and collect wide range of data, to validate the analyzed expressions developed for optimal thermo hydraulic performance of three sides roughened and glass covered solar air heater. Similar size of glass covered smooth solar air heater has been developed for experimentation. 3.2 SOLAR AIR HEATER DUCT MODELS Figs. 3.1 (a) and (b) show the two solar air heater duct models of similar size having high aspect ratio: (i) smooth with three sides glass covers and (ii) three sides roughened with three sides glass covers. Circular wire of different diameters has been provided on the absorber plate at varying pitches to serve as an artificial roughness element. The two ducts have been made similar in dimensions under similar conditions. Ambient air has been induced by two separate inlets but, exhausted by means of a single outlet. Wooden bottom sheet of 25 mm thickness has been used for back insulation in both solar air heater ducts and two sides insulation in smooth solar air heater duct. 4mm thick glass covers have been provided. An air gap of 25 mm has been provided for each duct resulting in a duct cross section of 200 mm × 25 mm each.

41 | P a g e

Figs. 3.1 (a) and (b) Solar air heater duct models

Out of the total duct length of 2130 mm, 1500 mm portion was utilized for instrumentation and the balance length of 630 mm served as the entry length for flow stabilization. In three sides roughened and glass covered solar air heater duct, the roughness elements of small diameter wires have been provided on the flow side of the absorber plate normal to the ambient air flow direction at varying pitch.

3.3 EXPERIMENTAL SET-UP Figs. 3.2 (a) and (b) show the schematic diagram and top view of experimental set-up respectively, combining the two solar air heater ducts (three sides roughened and smooth) with a single blower to run simultaneously. Directions denoted by arrows X and Y represent the ambient air flow in the roughened and smooth ducts respectively. Both the ducts are having three sides glass covers. Thermocouples were provided on the absorber plates to measure temperature, while digital thermometers were used to measure the air temperature in the air duct. The photographs of experimental set- up with instruments used have been shown in Figs. 3.3, while Fig 3.4 shows the experimental set up with instrumentation.

42 | P a g e

3.3.1

Absorber Plates and glass covers Black painted 22 SWG, G.I. sheets of 2000 mm length, 200 mm width and 25 mm height has been used as absorber plates in both the ducts. Fig.3.5 shows a typical artificially roughened absorber plate. Thin circular wires were soldered to the underside of the absorber plate to provide artificial roughness. Figs. 3.6 and 3.7 show the respective three sides roughened and smooth absorber plates.

3.3.2

Test Section and Entry Section of the Duct As mentioned in the previous section, an absorber plate length of 2000 mm has been instrumented to form test section, in which the entry length of 500 mm which also included about 100 mm long bell-mounted-converging section, that was 10 times the hydraulic diameter. Test sections and entry sections for both the roughened and smooth ducts were similar in dimensions. Entry and test length segments of both the ducts were provided with wooden covers and glass covers of 4 mm thickness. The entry section has the roughness geometry similar to the test section in the roughened duct.

43 | P a g e

All dimensions are in mm X-Three sides roughened and glass covered solar air heater duct inlet Y- Smooth and three sides glass covered solar air heater duct inlet P- Unheated wooden covered entry section Q- Glass covered test section R- Glass covers S- Flange tape orifice-meters

Fig. 3.2(a) Schematic diagram of the experimental set-up

44 | P a g e

Fig. 3.2 (b) Top view of the experimental set-up

Fig. 3.3 Photograph of experimental set-up with two ducts

45 | P a g e

Fig. 3.4 Photograph of experimental set-up with measuring instruments

Fig.3.5 Photograph of the roughness layout on top side of three sides roughened solar air heater

46 | P a g e

Fig. 3.6 Typical three sides artificially roughened absorber plates

Fig. 3.7 Typical three sides smooth absorber plates

47 | P a g e

Fig. 3.8 Photograph of digital pyranometer system 3.4

INSTRUMENTATION

3.4.1

Blower and Electric Motor The ambient air was sucked through the solar air heater duct system by means of a blower powered by a 3-phase, 5 H.P. motor. The power of the motor has to be such that it is capable of sucking air simultaneously through both the duct systems at desired mass flow rates to provide the maximum Reynolds number of 20000.

3.4.2

Flow Control Arrangement Air flow rate through the smooth and three sides roughened duct system was controlled by means of 3- phase auto-variac connected in series with the blower motor. Voltage regulation of the variac provided controlled supply of power to motor resulting in desired blower speeds.

3.4.3

Air Flow Measurements Air flow measurement was accomplished by a calibrated flange-taps type orifice-meter which was designed, fabricated and finally fitted in the 76.2 mm diameter pipe that carried the air flow from the duct system. Figs.3.9 and 3.10 show the details, specification and location of orifice-meter. The mass flow rate through the orifice-meter is expressed in the

48 | P a g e

terms of physical parameters of the orifice-meter and the pressure difference across it, by Eq. (3.1) written under: m = Cd A2

2ρ P 1 −P 2 1/2 1−β 41

(3.1)

where, Cd is the coefficient of discharge.

Fig. 3.9 Orifice plate location and specification 49 | P a g e

Fig. 3.10 Details of the flange-taps orifice meter construction 50 | P a g e

Coefficient of discharge of the orifice-meter has been determined by the observation of flow rate values through a pitot-tube. The flow rate has been determined by measuring average flow velocity using a pitot-tube, using the relationship: V = √2∆P/ρ

(3.2)

where, ∆P is pressure difference by pitot-tube. The mass flow rate has been determined by using the relationship: m = ρAC V

(3.3)

where, AC is the cross-sectional area of the pipe. The experimental values of Cd in Eqn. (3.1) have been plotted as a function of throat Reynolds number in Fig. 3.11, by using Murdock equation (Ozisik, 1985) for an orificemeter with flange-taps and are given by: Cd = C0 + C

10 4

a

Rd

(3.4)

where, the values of C0 and C for 76.2 mm diameter pipe, can be seen from Table no. 21.5 of the reference (Ozisik, 1985) as, C0 = 0.60394, C = 0.03270 and the exponent a=1, for flange-taps. Mass flow rate in the orifice-meter pipe bas been determined from Eq. (3.1) using experimental values of pressure p1 and p2. 3.4.4

Temperature and Pressure Drop Measurement Calibrated copper-constantan thermocouples of 28 SWG have been used to measure the plate temperature and the air temperature was measured by a digital thermometer. Thermocouples were soldered on the absorber plates at six sections. The thermocouple wires passing through connectors were connected to selector switch through copper connecting wires of 28 SWG. Plate and air temperature readings were noted from the digital display. Thermocouple arrangement for the plate temperature measurement has been shown in Fig. 3.12.The duct pressure drop was measured by means of a multi-tube manometer. The pressure signals were fed to the manometer through PVC piping attached to pressure taps of 6 mm diameter PU-hose connector inserted along the centre line of the bottom wall of the 51 | P a g e

duct. Pressure taps arrangement for pressure drop measurement and digital thermometer for air temperature measurement has been shown in Fig. 3.13.

1.0

 10 4  C d  0.60394  0.03270   Rd 

a

Experimental values of C

d

0.8

0.6

Cd

0.4

0.2

0.0 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200

(Throat Reynolds Number) Rd  103 Fig. 3.11 Orifice-meter Calibration

52 | P a g e

Dimensions are in mm Fig. 3.12 Thermocouple arrangement on top and side of the absorber plates

Fig. 3.13 Arrangement of air temperature sensors and pressure taps along the duct length 53 | P a g e

3.4.5

Solar Radiation Intensity Measurement Intensity of solar radiation was measured by means of a standard digital display pyranometer system in terms of W/m2.

3.5

EXPERIMENTAL PROCEDURE Test runs were conducted under controlled conditions of mass flow rate (i.e. flow Reynolds number, 𝑅𝑒) and the roughness parameters (𝑝/𝑒, 𝑒/𝐷), under actual outdoor conditions with varying values of intensity of solar radiation at quarter hourly intervals each day between 11 AM to 2 PM in the months of February to May 2015. Parameters measured were: 1) Pressure difference across orifice-meter 2) Pressure drop across inclined tube manometer 3) Temperature of the absorber plates 4) Temperature of air flowing in the duct 5) Ambient temperature 6) Intensity of solar radiation In 105 numbers of runs, 15 sets of three sides roughened absorber plate, along with the smooth absorber plate, each under 7 values of mass flow rate were tested. Typical test data collected during a test run are shown in Table 3.1 and range of parameter investigated is shown in table 3.2.

54 | P a g e

Table-3.1 Test Data Collected During a Typical Run (Run No. 56)

55 | P a g e

Table-3.2 Range of Parameters Investigated Set No.

1

2

3

4

Run no. 1

Date

𝑚,(Kg/s)

𝑅𝑒

15.02.2015

0.0375

19498

2

16.02.2015

0.0282

13572

3

17.02.2015

0.0271

14127

4

18.02.2015

0.0187

9749

5

19.02.2015

0.0175

6894

6

20.02.2015

0.0118

6166

7

21.02.2015

0.0085

4356

8

22.02.2015

0.0354

18576

9

23.02.2015

0.0227

12575

10

24.02.2015

0.0251

13125

11

25.02.2015

0.0243

9347

12

26.02.2015

0.0169

6263

13

27.02.2015

0.0108

6145

14

28.02.2015

0.0076

4243

15

01.03.2015

0.0361

18778

16

02.03.2015

0.0262

13451

17

03.03.2015

0.0284

14124

18

04.03.2015

0.0194

9554

19

05.03.2015

0.0161

6374

20

06.03.2015

0.0112

6151

21

07.03.2015

0.0077

4254

22

08.03.2015

0.0345

18463

23

09.03.2015

0.0232

13472

24

10.03.2015

0.0231

14117

25

11.03.2015

0.0157

9249

26

12.03.2015

0.0155

6194

27

13.03.2015

0.0098

6106

28

14.03.2015

0.0069

4156

𝑝,(mm) 𝑒,(mm)

𝑝/𝑒

𝑒/𝐷

11.0

1.1

10

0.0250

16.5

1.1

15

0.0250

22.0

1.1

20

0.0250

27.5

1.1

25

0.0250

56 | P a g e

5

6

7

8

29

15.03.2015

0.0376

19698

30

16.03.2015

0.0296

13772

31

17.03.2015

0.0262

14227

32

18.03.2015

0.0191

9049

33

19.03.2015

0.0189

6994

34

20.03.2015

0.0121

6366

35

21.03.2015

0.0085

4368

36

22.03.2015

0.0344

19198

37

23.03.2015

0.0252

13272

38

24.03.2015

0.0241

13827

39

25.03.2015

0.0157

9449

40

26.03.2015

0.0135

7079

41

27.03.2015

0.0109

5966

42

28.03.2015

0.0063

4056

43

29.03.2015

0.0364

19398

44

30.03.2015

0.0272

13472

45

31.03.2015

0.0261

14027

46

01.04.2015

0.0177

9649

47

02.04.2015

0.0165

6794

48

03.04.2015

0.0108

6076

49

04.04.2015

0.0073

4256

50

05.04.2015

0.0304

18498

51

06.04.2015

0.0242

12572

52

07.04.2015

0.0231

13127

53

08.04.2015

0.0117

9649

54

09.04.2015

0.0135

6794

55

10.04.2015

0.0101

5865

56

11.04.2015

0.0094

4886

33.0

1.1

30

0.0250

10.0

1.0

10

0.0225

15.0

1.0

15

0.0225

20.0

1.0

20

0.0225

57 | P a g e

9

10

11

12

57

12.04.2015

0.0354

19476

58

13.04.2015

0.0262

13549

59

14.04.2015

0.0252

14105

60

15.04.2015

0.0167

9727

61

16.04.2015

0.0155

6872

62

17.04.2015

0.0102

6143

63

18.04.2015

0.0079

4334

64

19.04.2015

0.0363

19483

65

20.04.2015

0.0271

13557

66

21.04.2015

0.0259

14112

67

22.04.2015

0.0176

9734

68

23.04.2015

0.0164

6879

69

24.04.2015

0.0107

6151

70

25.04.2015

0.0072

4342

71

26.04.2015

0.0361

19208

72

27.04.2015

0.0269

13282

73

28.04.2015

0.0258

13837

74

29.04.2015

0.0174

9459

75

30.04.2015

0.0162

6604

76

01.05.2015

0.0052

5876

77

02.05.2015

0.0071

4066

78

03.05.2015

0.0375

19512

79

04.05.2015

0.0251

11450

80

05.05.2015

0.0268

14124

81

06.05.2015

0.0191

9781

82

07.05.2015

0.0181

7301

83

08.05.2015

0.0123

6503

84

09.05.2015

0.0079

4101

25.0

1.0

25

0.0225

30.0

1.0

30

0.0225

6.0

0.6

10

0.0130

9.0

0.6

15

0.0130

58 | P a g e

13

14

15

85

10.05.2015

0.0352

19228

86

11.05.2015

0.0259

13302

87

12.05.2015

0.0248

13857

88

13.05.2015

0.0164

9479

89

14.05.2015

0.0152

6624

90

15.05.2015

0.0095

5896

91

16.05.2015

0.0062

4086

92

17.05.2015

0.0377

19499

93

18.05.2015

0.0283

13573

94

19.05.2015

0.0272

14128

95

20.05.2015

0.0188

9750

96

21.05.2015

0.0176

6895

97

22.05.2015

0.0119

6167

98

23.05.2015

0.0086

4357

99

24.05.2015

0.0378

19504

100

25.05.2015

0.0284

13577

101

26.05.2015

0.0273

14132

102

27.05.2015

0.0189

9754

103

28.05.2015

0.0177

6899

104

29.05.2015

0.0121

6172

105

30.05.2015

0.0087

4361

12.0

0.6

20

0.0130

15.0

0.6

25

0.0130

18.0

0.6

30

0.0130

59 | P a g e

CHAPTER – 4 RESULTS AND DISCUSSIONS 4.1 INTRODUCTION The present chapter reports the data presentation, reduction and validation. Based on the experimental results, correlations of Nusselt number and friction factor have been developed as a function of flow and roughness parameters. The effect of roughness and flow parameters on heat transfer, friction factor, collector performance parameters and thermo hydraulic performance has been presented. Thermo hydraulic performance results have been obtained. Thermal and thermo hydraulic performance curves have been presented to design such solar air heaters

4.2 DATA PRESENTATION, REDUCTION AND VALIDATION The experimental data with respect to the relative roughness height, relative roughness pitch, flow Reynolds number, intensity of solar radiation, plate and air temperatures, pressure drops along the duct length and orifice-meters have been reduced to obtain the results, using the relevant expressions. The experimental values of heat transfer coefficient have been obtained using the following Eq. (4.1):

𝑚𝐶𝑃 𝑇0 − 𝑇𝑖 = 𝑕𝐴𝐶 𝑇𝑃 − 𝑇𝑓

(4.1)

where, the values of mass flow rate have been calculated by using Eq. (3.4). The values of heat transfer coefficient, h, have been further used to calculate the values of Nusselt number by the following Eq. (4.2):

𝑁𝑢 =

𝑕𝐷 𝐾

(4.2)

The values of Nusselt number have been further used to obtain the values of Stanton number by Eq. (4.3) written under: 𝑁𝑢

𝑆𝑡 = 𝑅𝑒 𝑃𝑟

(4.3) 60 | P a g e

The respective values of Stanton numbers 𝑆𝑡3𝑟 and 𝑆𝑡3𝑔𝑠 , for the respective solar air heaters have been worked out using Eq. (4.3). However, to obtain the analytical values of Stanton number, the following Eq. (4.4) of Prasad et al. (2014) has been utilized: 𝑁𝑢3𝑟 =

𝑓3𝑟 /2 𝑓 1+√( 3𝑟 )[4.5 2

𝑝 0.53 𝑒 + 0.28 𝑃𝑟 0.57 −0.95 ]

𝑅𝑒 𝑃𝑟

(4.4)

𝑒

where, 𝑓3𝑟 , has been taken from Eq.(2.4). For smooth solar air heater, Eq. (2.14) has been used to obtain the value of 𝑆𝑡3𝑔𝑠 . The experimental values of friction factor, 𝑓3𝑟 , have been worked out from Eq. (3.13) and the values of 𝑓𝑆 have been found out from Eq.(2.13). 4.2.1

Sample Calculation Experimental results have been presented in the form of Nusselt Reynolds number, friction factor, temperature rise of air, intensity of solar radiation and the roughness parameters (𝑝/𝑒 and 𝑒/𝐷). Sample calculations for Run no. 56 (Table 3.1 and 3.2) have been illustrated as under:

(a) Average Temperature Average fluid temperature and average plate temperature for three sides roughened and the smooth ones have been calculated by using the following Eqs. (4.5) and (4.6) to be 35.290C & 50.140C respectively for three sides roughened duct and 33.60C & 51.430C for the smooth ones. 𝑇𝑓 = 𝑇𝑃 =

𝑇𝑓1 +𝑇𝑓2 +𝑇𝑓3 +𝑇𝑓4 +𝑇𝑓5 +𝑇𝑓6 +𝑇𝑓7 7 𝑇𝑃 1 +𝑇𝑃 2 +𝑇𝑃 3 +𝑇𝑓𝑃 4 +𝑇𝑃 5 +𝑇𝑃 6 +𝑇𝑃 7 7

(4.5) (4.6)

(b) Reynolds Number For a pressure drop of 25 mm of water across orifice-meter, throat velocity of air (V2) can be calculated using Eq. (7) under: 𝑉2 =

1/2

2∆𝑝 𝜌 1−𝛽 1

4

(4.7)

which works out to be 22.95 m/s for, 𝛽1 = 0.5. And hence, throat Reynolds number (Rd) as calculated using Eq. (8) under: 61 | P a g e

𝑅𝑑 =

𝑉2 𝑑 2 𝜗

(4.8)

is 45105. Discharge coefficient from Eq. (3.4) using the value of throat Reynolds number is worked out to be 0.6112. Mass flow rate given by Eq. (3.1), works out to be 0.00941kg/s, which gives a value of air flow velocity of 1.47 m/s through each duct, corresponding to a Reynolds number values of 4886 in each of the duct (three sides roughened and smooth) which remains invariant during a day. (c) Heat Transfer Coefficient Values of heat transfer coefficient from Eq. (4.1) work out to be 47.73 W/m2.K and 27.94 W/m2.K for three sides roughened and smooth respectively, giving the value of respective Nusselt numbers equal to 78.27& 45.81which yield in the value of Stanton numbers equal to 0.0228 & 0.0134 for the respective cases. (d) Friction Factor Values of friction factor have been calculated using Eq. (4.9) under: 2∆𝑃𝐷

𝑓 = 4𝐿𝜌 𝑣 2

(4.9)

are worked out be 0.0353 and 0.0251 for three sides roughened and smooth one respectively. (e) Temperature Rise of Air The difference of outlet and inlet temperatures of air through Eq. (4.10) in three sides roughened and smooth solar air heaters. ∆𝑇 = 𝑇0 − 𝑇𝑖

(4.10)

Table 4.1 shows the values of different resulting parameters for test Run no. 56. 4.3 VARIATION OF INTENSITY AND AMBIENT TEMPERATURE Fig. 4.1 shows typically the variation of intensity of solar radiation and ambient temperature. Fig. 4.2 shows the instantaneous plate and fluid (air) temperature in three sides roughened and smooth solar air heaters, across the collector length. Fig. 4.2 shows that the plate temperature in smooth solar air heater is more than that in three sides roughened one, but air temperature in three sides roughened solar air heater is more than that in smooth one. Higher values of plate temperature and lower values of air temperature in three sides glass covered 62 | P a g e

smooth solar air heater than those in three sides roughened and glass covered solar air heater indicate higher rate of heat transfer in three sides roughened collector than that in smooth collector. Table- 4.1Values of I, Ta, h, St and Nu as a Function of Time During the Day (Run No. 56)

Time

𝐼

, 𝑇𝑎 , 0C

𝑆𝑡3𝑟

𝑆𝑡3𝑔𝑠

(hours) w/m2

𝑕3𝑟 ,

𝑁𝑢3𝑟

W/m2K

𝑕3𝑔𝑠

,

𝑁𝑢3𝑔𝑠

W/m2K

11.00

733

29.0

0.02286

0.01337

46.73

78.18

26.37

45.73

11.15

746

29.2

0.02287

0.01341

47.84

78.23

27.48

45.84

11.30

756

29.6

0.02257

0.01308

46.72

77.18

26.27

44.72

11.45

766

30.4

0.02288

0.01339

47.73

78.27

27.94

45.81

12.00

811

30.9

0.02259

0.01292

46.62

77.27

26.27

44.17

12.15

844

31.3

0.02253

0.01289

46.09

77.06

26.91

44.06

12.30

877

31.8

0.02283

0.01319

47.13

78.09

27.31

45.11

12.45

899

32.3

0.02284

0.01318

47.12

78.11

27.21

45.09

13.00

842

33.1

0.02289

0.01327

47.39

78.29

27.93

45.37

13.15

816

33.5

0.02262

0.01295

46.29

77.37

26.92

44.29

13.30

786

34.1

0.02299

0.01317

47.02

78.02

27.21

45.02

13.45

777

34.7

0.02274

0.01311

46.81

77.79

26.18

44.81

14.00

748

35.1

0.02304

0.01340

47.78

78.81

27.87

45.79

63 | P a g e

Fig. 4.1 Typical variation of solar radiation and ambient temperature on a day

64 | P a g e

.

m  0.00835  0.0375Kg / s

56 54

Legends TP Tf S, 3G 3r, 3G

52 50 48 46 44

36

0

Tf ( C)

40

TP ( C),

42

0

38

34 32 30 28 0

150 300 450 600 750 900 1050 1200 1350 1500

Collector length, L (mm) Fig. 4.2 Typical air and plate temperature in solar air heaters

65 | P a g e

4.4 VALIDATION OF THE EXPERIMENTAL DATA The experimental values of Nusselt number and friction factor obtained for smooth solar air heater have been compared with those calculated by using the relevant Eqs. (4.11) and (4.12) respectively and have been shown in Figs. 4.3 and 4.4. It could be seen that the experimental values of Nusselt number and friction factor compare well with a deviation of ±2.7% and ±2.5% respectively. 𝑁𝑢𝑆 = 0.023𝑅𝑒 0.8 𝑃𝑟 0.4 (Dittus-Boelter equation)

(4.11)

𝑓𝑆 = 0.085𝑅𝑒 −0.25 (Modified Blasius equation)

(4.12)

Similarly, the present experimental values of heat transfer and friction factor data in three sides roughened solar air heaters have been compared with those of the experimental values of Behura et al., (2016) and are shown in Figs. 4.5 and 4.6 respectively. It could be seen from Figs. 4.5 and 4.6 that the present values of heat transfer and friction factor data compare well with those of the experimental values of Behura et al., (2016), with a maximum variation of ±4.5% and ±4.1% respectively. Since, the present results on heat transfer and friction factor have been validated, correlations for heat transfer and friction factor have been developed utilizing the present data value results.

66 | P a g e

Present experimental values Dittus -Boelter Eq.(4.11)

70

60

50

40

Nu 30

20

10 3000

6000

9000

12000

15000

18000

21000

Re Fig.4.3 Comparison of heat transfer data in smooth solar air heaters

67 | P a g e

Present experimental values Modified Blasius Eq.(4.12) 0.036 0.032 0.028

f 0.024 0.020 0.016 3000

6000

9000

12000

15000

18000

21000

Re

Fig. 4.4 Comparison of friction factor data in smooth solar air heaters

68 | P a g e

Present experimental values Behura et al. (2016) 80 70

60 50

Nu3r , Nur

40 30

20 10 4000 6000 8000 10000 12000 14000 16000 18000 20000 Re Fig. 4.5 Comparison of heat transfer data in three sides roughened solar air heater

69 | P a g e

Behura et al. (2016) Present experimental values 0.052 0.048 0.044 0.040

f3r , fr 0.036 0.032 0.028 0.024 4000 6000 8000 10000 12000 14000 16000 18000 20000 Re

Fig. 4.6 Comparison of friction factor data in three sides roughened solar air heaters

70 | P a g e

4.5 DEVELOPMENT

OF

CORRELATIONS

FOR

NUSSELT

NUMBER

AND

FRICTION FACTOR On the basis of experimental data of heat transfer and friction factor, the correlations for three sides roughened and glass covered solar air heaters, have been developed by using the method of regression analysis, as a function of roughness and flow parameters. A correlation for heat transfer for three sides glass covered smooth solar air heater has also been developed. Since, Nusselt number and friction factor both are the strong function of roughness and flow parameters, the following functional relationship for 𝑁𝑢3𝑟 and 𝑓3𝑟 can be written for three sides roughened solar air heater as under: 𝑝 𝑒

𝑁𝑢3𝑟 = 𝑓1 (𝑅𝑒, 𝑒 , 𝐷 )

(4.13)

𝑝 𝑒

𝑓3𝑟 = 𝑓3 (𝑅𝑒, 𝑒 , 𝐷 )

(4.14)

In order to determine a functional relationship between Nu and Re, all the experimental values of Nusselt number, for 𝑒 + ≤ 23, were plotted against entire range of Reynolds number, as shown in Fig. 4.7. The well known regression analysis method (curve fitting method by using origin pro 9.0, 32 Bit) for forced convective heat transfer has been utilized to determine the functional equation, written under as: 𝑁𝑢3𝑟 = 𝐴0 𝑅𝑒

0.746

(4.15)

where, 𝐴0 is the function of relative roughness pitch (𝑝/𝑒). 𝐴0 =

𝑁𝑢 3𝑟 𝑅𝑒 0.746

, is plotted

against all the experimental values of 𝑝/𝑒, as shown in Fig. 4.8. The regression analysis to fit a straight line through these points has been given by: 𝑁𝑢 3𝑟 𝑅𝑒 0.746

= 𝐵0 (𝑝/𝑒)−0.0713

(4.16) 𝑁𝑢 3𝑟

where, 𝐵0 is the function of relative roughness height (𝑒/𝐷). 𝐵0 =

𝑅𝑒 0.746

(𝑝/𝑒)−0.0713 ,

is plotted against entire experimental values of e/D, as shown in Fig. 4.9. The regression

71 | P a g e

analysis to fit a straight line through these points, yield the following correlation for Nusselt number of three sides roughened and glass covered solar air heaters for 𝑒 + ≤ 23 : 𝑁𝑢3𝑟 = 0.1415 𝑝/𝑒

−0.0713

𝑒/𝐷

0.0871

𝑅𝑒

0.746

(4.17)

Similarly, the correlation for Nusselt number corresponding to (𝑒 + > 23), were developed by using the method, same as that of the above. Figs. (4.10-4.12), show the plots between the obtained results, whenafter, applying regression analysis method, the following correlation for Nusselt number for three sides roughened and glass covered solar air heaters corresponding to 𝑒 + > 23, were developed:

𝑁𝑢3𝑟 = 0.0422 𝑝/𝑒

−0.0461

𝑒/𝐷

0.0375

𝑅𝑒

1.208

Nu3r=A0(Re)

for

𝑒 + > 23

12000

16000

(4.18)

0.746

Nusselt number

90

60

30

0 0

4000

8000

20000

Reynolds number

Fig.4.7 Variation of Nusselt number with Reynolds number

72 | P a g e

0.00045 0.00040 0.746

0.00035

Nu3r=B0(Re) (p/e)

-0.0713

Nu3r/(Re)

0.746

0.00030 0.00025 0.00020 0.00015 0.00010 0.00005 8

12

16

20

24

28

32

p/e

Fig. 4.8 Plot of

𝑁𝑢 3𝑟 𝑅𝑒 0.746

versus p/e

73 | P a g e

0.0016

0.746

-0.0713

0.0871

0.030

0.035

Nu3r=C0(Re) (p/e)

0.0012

[{Nu3r/(Re)

0.746

}/(p/e)

-0.0713

]

0.0014

(e/D)

0.0010

0.0008

0.0006 0.015

0.020

0.025

0.040

e/D

Fig.4.9 Plot of

𝑁𝑢 3𝑟 𝑅𝑒 0.746

(𝑝/𝑒)−0.0713 versus e/D

74 | P a g e

100

Nu3r=A0(Re)

1.208

Nusselt number

80

60

40

20

0 0

4000

8000

12000

16000

20000

Reynolds number Fig.4.10 Variation of Nusselt number with Reynolds number

75 | P a g e

0.00045 0.00040 0.00035

1.208

Nu3r=B0(Re) (p/e)

-0.0461

Nu3r/(Re)

1.208

0.00030 0.00025 0.00020 0.00015 0.00010 0.00005 8

12

16

20

24

28

p/e Fig. 4.11 Plot of

𝑁𝑢 3𝑟 versus 𝑅𝑒 1.208

p/e

76 | P a g e

0.0018

0.0016 ] -0.0461

}/(p/e) 1.208

[{Nu3r/(Re)

1.208

-0.0461

0.030

0.035

Nu3r=C0(Re) (p/e) 0.0014

0.0012

0.0010

0.0008

0.0006 0.015

0.020

0.025

0.040

0.045

e/D

Fig.4.12 Plot of

𝑁𝑢 3𝑟 𝑅𝑒 1.208

(𝑝/𝑒)−0.0461 versus e/D

77 | P a g e

4.5.1

Correlation for Nusselt Number for Three Sides Glass Covered Smooth Solar Air Heater Nusselt number of smooth solar air heaters is the function of Prandtl number (Pr) friction factor (fS) and flow Reynolds number (Re), expressed as Eq. (4.19) under: 𝑁𝑢3𝑔𝑠 = 𝑓(𝑅𝑒, 𝑃𝑟, 𝑓𝑆 )

(4.19)

and hence, to determine the functional relationship between the flow Reynolds number and Nusselt number for three sides smooth solar air heater, all the experimental values of flow Reynolds number were plotted against all the experimental value of Nusselt number, as shown in Fig.4.13. The regression analysis method has been utilized to determine the functional equation, written by Eq. (4.20) under: 𝑁𝑢3𝑔𝑠 = 𝐴1 𝑅𝑒

0.932

(4.20)

where, 𝐴1 is the function of 𝑃𝑟. 𝐴1 =

𝑁𝑢 3𝑔𝑠 𝑅𝑒 0.932

, is plotted against all the experimental values

of 𝑃𝑟 , as shown in Fig. 4.14. The regression analysis to fit a straight line through these points has been given by: 𝑁𝑢 3𝑔𝑠 𝑅𝑒 0.932

= 𝐵1 (𝑃𝑟)0.5

(4.21) 𝑁𝑢 3𝑔𝑠

where, 𝐵1 is the function of 𝑓𝑆 . 𝐵1 =

𝑅𝑒 0.932

(𝑃𝑟)0.5 , is plotted against entire

experimental values of 𝑓𝑆 , as shown in Fig.4.15. The regression analysis to fit a straight line through these points, yields the following correlation for Nusselt number in three sides glass covered smooth solar air heaters:

𝑁𝑢3𝑔𝑠 = 0.0782 𝑃𝑟

0.5

𝑓𝑆

0.7

𝑅𝑒

0.932

(4.22)

78 | P a g e

70

60

0.932

Nu3gS=A1(Re)

50

40

Nu3sg 30

20

10 3000

6000

9000

12000

15000

18000

21000

Re Fig.4.13Nusselt number versus flow Reynolds number for smooth solar air heater

79 | P a g e

0.0018 0.0016 0.0014 0.932

Nu3gS=B1(Re) (Pr)

Nu3gS/(Re)

0.932

0.0012

0.5

0.0010 0.0008 0.0006 0.0004 0.0002 0.690

0.692

0.694

0.696

0.698

0.700

Pr Fig.4.14 Plot of

𝑁𝑢 3𝑔𝑠 𝑅𝑒 0.932

versus 𝑃𝑟

80 | P a g e

0.0037 0.0036

[{Nu3gS/(Re)

0.932

}/(Pr)

0.5

]

0.0035

0.932

0.5

Nu3gS=C1(Re) (Pr) (fS)

0.7

0.0034 0.0033 0.0032 0.0031 0.0030 0.025 0.026 0.027 0.028 0.029 0.030 0.031 0.032

fS 𝑁𝑢 3𝑔𝑠

Fig.4.15 Plot of

𝑅𝑒 0.932

(𝑃𝑟 )0.5 versus 𝑓𝑆

81 | P a g e

4.5.2

Correlation for Friction Factor of Three Sides Roughened and Glass Covered Solar Air Heater For the determination of functional relationship between the flow Reynolds number and friction factor, all the experimental values of flow Reynolds number were plotted against all the experimental value of friction factor, as shown in Fig. 4.16. The regression analysis method has been utilized to determine the functional equation, written under as: 𝑓3𝑟 = 𝐴2 𝑅𝑒

−1.76

(4.23)

Where, 𝐴2 is the function of /𝑒. 𝐴2 =

𝑓3𝑟 𝑅𝑒 −1.76

, is plotted against all the experimental

values of 𝑝/𝑒, as shown in Fig. 4.17. The regression analysis to fit a straight line through these points has been given by: 𝑓3𝑟 𝑅𝑒 −1.76

= 𝐵2 (𝑝/𝑒)−0.319

(4.24)

Where, 𝐵2 is the function of 𝑒/𝐷 . 𝐵2 =

𝑓3𝑟 𝑅𝑒 −1.76

(𝑝/𝑒)−0.319 , is plotted against entire

experimental values of 𝑒/𝐷, as shown in Fig. 4.16. The regression analysis to fit a straight line through these points, yield the following correlation for Nusselt number of three sides glass covered smooth solar air heaters : 𝑓3𝑟 = 0.372 𝑝/𝑒

−0.319

𝑒/𝐷

0.354

𝑅𝑒

−1.76

(4.25)

82 | P a g e

0.044 -1.76

f3r=A2(Re) 0.040

0.036

f 3r 0.032

0.028

0.024 4000 6000 8000 10000 12000 14000 16000 18000 20000

Re

Fig. 4.16 Variation of friction factor with Reynolds number

83 | P a g e

0.5

-1.76

f3r=B2(Re) (p/e)

-0.319

f3r/(Re)

-1.76

0.4

0.3

0.2

8

12

16

20

24

28

p/e Fig. 4.17 Plot of

𝑓3𝑟 𝑅𝑒 −1.76

versus 𝑝/𝑒

84 | P a g e

0.18 -1.76

-0.319

0.354

0.016

0.018

0.020

f3r=C2(Re) (p/e)

(e/D)

0.14

[{f3r/(Re)

-1.76

}/(p/e)

-0.319

]

0.16

0.12

0.10

0.012

0.014

0.022

0.024

0.026

e/D

Fig.4.18 Plot of

𝑓3𝑟 𝑅𝑒 −1.76

(𝑝/𝑒)−0.319 versus 𝑒/𝐷

85 | P a g e

100

Predicted Nusselt number

80

60

+12% 40

-12% 20

0 0

20

40

60

80

100

Experimental Nusselt number Fig.4.19 Comparison of experimental Nusselt number values to the Predicted ones

86 | P a g e

0.05

Predicted friction factor

0.04

0.03

+11% 0.02

-11%

0.01

0.00 0.00

0.01

0.02

0.03

0.04

0.05

Experimental friction factor Fig. 4.20 Comparison of experimental friction factor values to the Predicted ones

87 | P a g e

It can seen from the Figs. 4.19 and 4.20, that the experimental values of the Nusselt number and friction factor compare well with that of the predicted ones with the average deviations of ±12% and ±11%, respectively, which are within the acceptable limit. Hence, the correlations developed are reasonably satisfactory for the prediction of Nusselt number and friction factor in three sides artificially roughened and glass covered solar air heaters. Fig. 4.21 shows the comparison of heat transfer data for three sides roughened and glass covered solar air heater and three sides glass covered smooth one with those of predicted Eq. (4.17) and (4.22) respectively. The experimental values of Nusselt number for both three sides roughened and three sides glass covered smooth solar air heaters are seen to compare well with those of the predicted Eqs. (4.17) and (4.22) respectively. For the entire range of experimental values of heat transfer for three sides roughened and three sides glass covered smooth solar air heaters, the maximum deviation was found to be about ±12% and ±10%, respectively, which are reasonably satisfactory for the predicted correlations. Fig. 4.22 shows the comparison of friction factor values for three sides roughened and glass covered solar air heater with those of predicted Eq. (4.25).The experimental values of friction factor are seen to compare well with those of the predicted Eq. (4.25).The maximum deviation was found to be about ±11%.

88 | P a g e

(1)

Predicted Eq. (4.17) Present experimental (Nu3r)

(2)

Predicted Eq. (4.22) Present experimental (Nu3gs)

80

(1) 70

60

50

Nu3gs Nu3r

(2)

40

30

20

10 4000 6000 8000 10000 12000 14000 16000 18000 20000

Re Fig. 4.21Predicted and experimental heat transfer data in solar air heater

89 | P a g e

Predicted Eq. (25) Present experimental (f3r) 0.052 0.048 0.044 0.040

f 3r

0.036 0.032 0.028 0.024 4000 6000 8000 10000 12000 14000 16000 18000 20000 Re

Fig. 4.22 Predicted and experimental values of friction factor data in solar air heater

90 | P a g e

4.6 EFFECT OF ROUGHNESS AND FLOW PARAMETERS ON HEAT TRANSFER AND FRICTION FACTOR The experimental values of Nusselt number have been used to get the values of Stanton number by the following conventional Eq. (4.26), which have been further taken as the heat transfer data to see the effect of roughness and flow parameters on it. 𝑆𝑡 = 𝑁𝑢/ 𝑅𝑒 𝑃𝑟

(4.26)

Figs.4.23 and 4.24 show the effect of the roughness and flow parameters, 𝑝/𝑒, 𝑒/𝐷 and 𝑅𝑒 on Stanton number in three sides roughened solar air heater. It could be seen from these figures that the values of Stanton number in three sides roughened collector decrease with increasing values of relative roughness pitch, 𝑝/𝑒, for a given value of relative roughness height, 𝑒/𝐷 and decreasing values of relative roughness height, 𝑒/𝐷, for a given value of relative roughness pitch, 𝑝/𝑒. The values of Stanton number are found to decrease with increasing values of flow Reynolds number. Stanton number for three sides smooth solar air heater for varying values of Reynolds number is also shown in these Figs. 4.23 and 4.24. It could be worked out from Figs. 4.23 and 4.24 that the values of Stanton number are more in three sides roughened and glass covered collector in the range of 36% to 65%, as compared to the three sides glass covered smooth ones. Fig. 4.25 shows the effect of 𝑝/𝑒 on friction factor for a given value of 𝑒/𝐷, equal to 0.0225. It could be seen from the figure that the values of friction factor increase with decrease in the value of the relative roughness pitch, 𝑝/𝑒, and decrease with increasing values of the flow Reynolds number, 𝑅𝑒. Fig. 4.26 shows the effect of 𝑒/𝐷 on friction factor for a given value of 𝑝/𝑒, equal to 20. It is quite clear from the figure that the values of friction factor increase with the increase in the value of relative roughness height, 𝑒/𝐷 and decrease with increasing values of flow Reynolds number, 𝑅𝑒. The values of friction factor in three sides roughened and glass covered collector have been found to increase in the range of 28% to 38%, as compared to the three sides glass covered smooth ones.

91 | P a g e

p/e=10 p/e=20 p/e=30 smooth

e/D= 0.0225

0.008

0.007

0.006

St 3 r , St3 gs

0.005

0.004

0.003

0.002 4000

8000

12000

16000

20000

Re Fig. 4.23 Effect of p/e on Stanton number in three sides artificially roughened solar air heater

92 | P a g e

e/D=0.0250 e/D=0.0225 e/D=0.0130 smooth

p/e=20 0.008

0.007

0.006

St3r , St3gs

0.005

0.004

0.003

0.002 4000

8000

12000

16000

20000

Re Fig. 4.24 Effect of e/D on Stanton number in three sides artificially roughened solar air heater

93 | P a g e

p/e= 10 p/e= 20 p/e= 30 smooth

0.044

e/D=0.0225

0.040

0.036

f 3r , fS

0.032

0.028

0.024 4000

6000 8000 10000 12000 14000 16000 18000 20000

Re

Fig. 4.25 Effect of p/e on friction factor in three sides artificially roughened solar air heater

94 | P a g e

0.044

p/e=20

e/D= 0.0250 e/D= 0.0225 e/D= 0.0130 smooth

0.040

0.036

f 3r , fS

0.032

0.028

0.024

0.020 4000

6000

8000 10000 12000 14000 16000 18000 20000

Re

Fig. 4.26 Effect of e/D on friction factor in three sides artificially roughened solar air heater

95 | P a g e

4.7 HEAT TRANSFER ENHANCEMENT FACTOR AND FRICTION LOSS FACTOR Fig. 4.27 shows the results of heat transfer enhancement factor, 𝑁𝑢𝐹 and friction loss factor, 𝑓𝐹 with increasing values of Reynolds number for a given value of 𝑒/𝐷, equal to 0.0225 at varying values of 𝑝/𝑒. The values of 𝑁𝑢𝐹 and 𝑓𝐹 have been given by Eqs. (4.27) and (4.28) respectively:

𝑁𝑢𝐹 = 𝑓𝐹 =

𝑁𝑢 𝑟 −𝑁𝑢 𝑆 𝑁𝑢 𝑆

𝑓𝑟 −𝑓 𝑆 𝑓𝑆

(4.27) (4.28)

It could be seen from the figure that the values of both heat transfer enhancement factor and friction loss factor increase with the increasing values of Reynolds number. Fig. 4.27 also shows that for a given value of e/D the rate of increment of 𝑓𝐹 is slight higher than that in 𝑁𝑢𝐹 at varying values of p/e. Therefore, the three sides roughened solar air heaters have been found to perform better both quantitatively and qualitatively as compared to the smooth ones. 4.8 THERMAL PERFORMANCE OF SOLAR AIR HEATER Three sides artificially roughened and glass covered solar air heaters have been found more efficient than those of only top side roughened ones (Prasad et al., 2014; Behura et al., 2016). This is due to the fact that three sides roughened and glass covered solar air heaters have more absorber plate area and enhanced heat transfer coefficient than those of only one side roughened and glass covered solar air heaters. The absorber plate area of only top side glass covered roughened as well as smooth (Prasad, 2013) solar air heater is equal to: 𝐴𝐶 = 𝑊𝐿

(4.29)

whereas, the absorber plate area of three sides roughened and glass covered solar air heater and the three sides glass covered smooth one is equal to: 𝐴𝐶 = 𝑊𝐿 + 2(𝐵𝐿)

(4.30)

96 | P a g e

p/e=10(NuF)

0.50

p/e=10(fF) p/e=20(NuF)

e/D=0.0225

p/e=20(fF)

0.48

p/e=30(NuF) p/e=30(fF)

0.46

NuF fF

0.44

0.42

0.40

0.38 2000 4000 6000 8000 10000 12000 14000 16000 18000 20000

Re Fig. 4.27 Heat transfer enhancement factor and friction loss factor versus Reynolds number

97 | P a g e

The absorber plate area for three sides roughened and glass covered solar air heater and three sides glass covered smooth solar air heater is more by a magnitude of, 2(𝐵𝐿), than those of one side glass covered roughened as well as smooth solar air heater, which results in more heat collection for the same value of mass flow rate. Fig. 4.28 shows the comparative values of Stanton number for the respective solar air heaters. It could be observed from this figure that the solar air heater having three roughened sides and glass covers has the higher value of Stanton number, followed by the solar air heaters: one roughened side and glass cover; smooth and three sides glass covers; smooth and one side glass cover respectively. Thermal performance for three sides roughened & glass covered and three sides glass covered smooth solar air heaters, operating without recycling of air, has been presented on the basis of outlet fluid temperature. Based on the conventional performance equation (1.8) of (Bliss, 1959), replaced by Eq. (1.9), for solar air heaters without air recycling, thermal performance has been presented in Fig. 4.29. This figure represents typically the performance characteristics of the roughened solar air heaters for p/e =10, e/D = 0.0247: (i) three sides roughened and glass covered, (ii) top side roughened and glass covered, as well as of the smooth solar air heater: (i) smooth and three sides glass covered and (ii) smooth and only top side glass cover, all at varying mass flow rates. The efficiency data values given by the following equation: ƞ𝑡𝑕 = 𝑚𝐶𝑝 (𝑇0 − 𝑇𝑖 )/𝐼𝐴𝐶

(4.31)

shown by the respective lines are represented for each mass flow rate from the origin, resulting in the respective lines for the respective mass flow rates. Utilizing the efficiency data values, the respective curves A−𝐵3𝑟,3𝐺 , 𝐶 − 𝐷1𝑟,1𝐺 , 𝐸 − 𝐹𝑆,3𝐺 and 𝐺 − 𝐻𝑆,1𝐺 , have been obtained by the least square fit method for three sides roughened and glass covered solar air heater (present case), one side roughened and glass covered solar air heater (Prasad, 2013), smooth and three sides glass covered solar air heater (present case) and smooth and only top side glass covered solar air heater (Prasad, 2013), respectively. As could be seen from this figure that three sides roughened and glass covered solar air heater has higher value of thermal efficiency than those of only one side roughened and glass covered solar air

98 | P a g e

heater, followed by smooth and three sides glass covered solar air heater having slightly more thermal efficiency than that of smooth and only top side glass covered solar air heater.

p/e=20 e/D=0.0225

0.008

.

m  0.00835  0.0375Kg / s

S, 3G Present case 1r, 1G Prasad (2013) S, 1G Prasad (2013) 3r, 3G Present case

0.007

0.006

S t3r , St3 gs

0.005

0.004

0.003

0.002 4000

8000

12000

16000

20000

Re Fig. 4.28 Effect of roughened sides and glass covers on Stanton number in solar air heaters

99 | P a g e

m

kg/s 3r,3G(Present case) 1r,1G(Prasad, 2013) S,3G(Present case) S,1G(Prasad, 2013)

m1 0.0085 p/e=10 e/D=0.0247

m2 0.0118 m3 0.0175 m4 0.0187 m5 0.0271 m6 0.0282 m7 0.0375

90

m7

80

m6

70

m5

A

m4

60

m3

th%

50

m2

C

40 30

m1

E B3r,3G D1r,1G(Prasad, 2013)

G

20

FS,3G HS,1G(Prasad, 2013)

10 0 0.000

0.004

0.008

0.012 0

0.016

0.020

2

(T0-Ti)/I, C-m /W Fig 4.29 Thermal performance characteristics of solar air heaters

100 | P a g e

4.9 COLLECTOR HEAT REMOVAL FACTOR (𝐅𝐑 ) AND COLLECTOR EFFICIENCY FACTOR (𝐅 ′ ) Using the values of the slope,F0 UL and intercept, F0 τα ,of the respective performance characteristic curves, obtained from Fig. 4.29, values of the performance parameters FR UL and FR τα have been obtained using equations (1.10, 1.16 and 1.17) of Duffie and Beckman, (1980). In three sides roughened and glass covered solar air heater, for mass flow rate, 𝑚, equal to 0.0375 Kg/s and Cp equal to 1005 J/Kg K, the value of 𝐹𝑅 𝑈𝐿 and 𝐹𝑅 𝜏𝛼 , worked out from Eqs. (1.16) and (1.17) comes out equal to 8.881 W/m2K and 0.697 respectively. For 𝜏𝛼 , equal to 0.84, the value of 𝐹𝑅 coming to 0.823, gives the value of 𝑈𝐿 equal to 10.701 W/m2K and that of 𝐹 ′ equal to 0.891. Figs. 4.30 and 4.31 show the values of collector heat removal factor, (𝐹𝑅 ) and collector efficiency factor,(𝐹 ′ ), respectively for three sides roughened & glass covered and three sides glass covered smooth ones.It could be seen that the values of collector heat removal factor, 𝐹𝑅 , and collector efficiency factor, 𝐹 ′ , increases in three sides roughened and glass covered solar air heater as compared to smooth and three sides glass covered solar air heater. The values of the collector heat removal factor,𝐹𝑅 and plate efficiency factor, 𝐹 ′ , have been found to enhance in the range of 51 to 60% and 41 to 49% respectively, as compared to three sides glass covered smooth solar air heater, within the range of the parameters investigated.

101 | P a g e

1.1

p/e= 20,e/D=0.0225

3r,3G S, 3G

.

1.0

m  0.00835  0.0375Kg / s

0.9 0.8 0.7

FR 0.6 0.5 0.4 0.3 0.2 0.1 4000 6000 8000 10000 12000 14000 16000 18000 20000

Re Fig. 4.30 Collector heat removal factor in solar air heaters

102 | P a g e

e/D=0.0225 .

1.1

m  0.00835  0.0375Kg / s

p/e=20

3r,3G S, 3G

1.0

0.9

0.8

0.7

F' 0.6

0.5

0.4

4000

6000

8000 10000 12000 14000 16000 18000 20000

Re

Fig. 4.31 Collector efficiency factor in solar air heaters

103 | P a g e

Table- 4.2 represents the comparative view of the values of the slope, intercept, collector heat removal factor and collector efficiency factor obtained presently and those of Prasad, (2013), respectively. Higher values of the parameters, 𝐹𝑅 and 𝐹 ′ , indicate higher value of thermal efficiency. Table- 4.2 Values of Slope, Intercept, Collector Heat Removal Factor and Collector Efficiency Factor Slope,

𝐹0 𝑈𝐿 Intercept,

2

𝐹0 𝜏𝛼

(w/m .K)

3𝑟, 3𝐺

Collector heat Collector removal

efficiency

factor, 𝐹𝑅

factor, 𝐹 ′

(Present 8.56

0.725

0.822

0.918

(Prasad, 6.50

0.536

0.545

0.591

(Present 3.94

0.321

0.410

0.504

(Prasad, 3.60

0.291

0.391

0.496

case) 1𝑟, 1𝐺 2013) 𝑆, 3𝐺 case) 𝑆, 1𝐺 2013)

4.10 THERMAL PERFORMANCE CURVES Thermal performance equations (4.32) and (4.33) have been derived based on the conventional Eq. (1.8) and presented by thermal performance curves respectively in Fig. 4.32. It has been found that three sides roughened & glass covered solar air heater is 50 to 59% more efficient than that of three sides glass covered smooth one. 𝑇𝑖 −𝑇𝑎

ƞ𝑡𝑕(3𝑟,3𝐺) = 0.822 𝜏𝛼 − 0.822𝑈𝐿

ƞ𝑡𝑕(𝑆,3𝐺) = 0.41 𝜏𝛼 − 0.41𝑈𝐿

𝐼

𝑇𝑖 −𝑇𝑎 𝐼

(4.32)

(4.33)

104 | P a g e

th 

90 85 80 75 70 65 60 55 50 45 40 35 30 25 20 15 10 5 0

3r, 3G (Present work)

S, 3G (Present work)

0.002 0.004 0.006 0.008 0.010 0.012 0.014 0.016 0.018 0

2

(Ti-Ta)/I, C-m /w

Fig. 4.32 Thermal performance curves for solar air heaters

105 | P a g e

4.11 OPTIMAL THERMOHYDRAULIC PERFORMANCE Thermo hydraulic performance of solar air heater is expressed by Eq. (4.34) written under as: 𝑆𝑡 3𝑟

Ƞ𝑡𝑕𝑒𝑟𝑚𝑜 =

𝑆𝑡 3𝑔𝑠

𝑓3𝑟

(4.34)

𝑓𝑆

The values of thermo hydraulic performance worked out by utilizing Eq.(4.34) have been represented in Fig. 4.33, as a function of the parameter, e+, given by Eq.(2.6) to see the effect of the parameter, e+, on thermo hydraulic performance for varying values of p/e and given value of e/D, equal to 0.0225. Similarly, Fig. 4.34 represents the data for varying values of e/d and given value of p/e, equal to 20. It could be seen from these figures that for every value of p/e and e/D, the value of thermo hydraulic performance increases with increasing values of e+, and attains the maximum value at 𝑒 + = 23, whenafter, its value decreases for increasing values of e+. Since, e+, has remained the optimization parameter to arrive at the optimal thermo hydraulic performance + condition in literature, the value of 𝑒𝑜𝑝𝑡 = 23 , has the optimal thermo hydraulic

performance in three sides roughened and glass covered solar air heaters which validates the thermo hydraulic optimization analysis presently.

106 | P a g e

p/e=30 p/e=20 p/e=10

90 85

Re=4000-20000

e/D=0.0225

80 75

thermo %

70 65 60 55 50 45 40 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40

e   e / D f 3r / 2 Re

Fig. 4.33 Effect of p/e and e+ on thermohydraulic performance

107 | P a g e

e/D=0.0130 e/D=0.0225 e/D=0.0250

90 85

Re=4000-20000

p/e=20

80 75

thermo %

70 65 60 55 50 45 40 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40

e   e / D f 3r / 2 Re

Fig. 4.34 Effect of e/D and e+ on thermohydraulic performance

108 | P a g e

Covering the total range of experimental data, Fig. 4.35 has been drawn to see the effect of the optimization parameter, e+, on thermo hydraulic performance. It could be seen from this figure that the optimal value of thermo hydraulic performance is achieved at the value of + 𝑒𝑜𝑝𝑡 = 23 , which validates the present analysis. Based on the experimental results,

correlations for thermo hydraulic efficiency have been obtained by using regression analysis method. The correlations are given below as: Ƞ𝑡𝑕𝑒𝑟𝑚𝑜 = 0.102 𝑒 +

0.658

for

𝑒 + ≤ 23

(4.34)

Ƞ𝑡𝑕𝑒𝑟𝑚𝑜 = 11.91 𝑒 +

−0.913

for

𝑒 + > 23

(4.35)

The present optimal value of thermo hydraulic performance in such solar air heaters could + be read from the Fig.4.35 to be about 78.8% for optimal value of 𝑒𝑜𝑝𝑡 = 23.

The value of the optimization parameter, 𝑒 +, combines in itself the values of the roughness and flow parameters, under varying conditions of the roughness and flow parameters. As also, the effect of either of the parameters on it may affect the values of other parameters to result in the same value of it. An increase in the value of one parameter, might be compensated by an equivalent decrease in the values of other parameters, to result in the + similar value of 𝑒𝑜𝑝𝑡 . Since, performance of solar air heaters is very much mass flow rate

dependent, flow Reynolds number greatly affects the heat transfer coefficient as well as the friction factor. Results show that heat transfer coefficient increases with increasing values of flow Reynolds number but friction factor decreases with increasing values of it. The value of + 𝑒𝑜𝑝𝑡 = 24, for one side roughened solar air heater (Prasad and saini, 1991; Verma and + Prasad, 2000), and 𝑒𝑜𝑝𝑡 = 23, for three sides roughened solar air heaters, presently, differ

marginally, which might be due to the two additional roughened sides of the present solar air heater duct . Three sides roughened and glass covered solar air heater has been found to have more value of the heat transfer coefficient than that of one side glass covered + roughened solar air heater (Prasad et al., 2014; Behura et al., 2016). The value of 𝑒𝑜𝑝𝑡 , is

mainly affected by the values of p/e, e/D and Re. As the heat transfer and pressure drop results have been found to compare well within the parameters investigated with those of the available ones, it gives the authentication of the 109 | P a g e

+ present investigation and therefore, 𝑒𝑜𝑝𝑡 = 23, is the only condition established for three

sides roughened and glass covered solar air heaters. The correlations developed for heat transfer, friction factor and thermo hydraulic efficiency are also valid within the parameters investigated.

85

Re=4000-20000 p/e=10-30 e/D=0.0130-0.0250

80 75

thermo(%)

70 65 60 55 50 6

8 10 12 14 16 18 20 22 24 26 28 30 32 34 36

e

+

Fig. 4.35 Thermohydraulic efficiency versus roughness Reynolds number

110 | P a g e

CHAPTER 5 CONCLUSIONS AND FUTURE SCOPE 5.1 CONCLUSIONS The following conclusions have been derived based on the present results:1. The optimal value of the optimization parameter, 𝑒 + , for optimal thermo hydraulic performance in three sides roughened and glass covered solar air heaters has been obtained + to be 𝑒𝑜𝑝𝑡 = 23 and validated.

2. The heat transfer and fluid flow characteristics of three sides artificially roughened and glass covered solar air heaters as also three sides glass covered smooth solar air heater have been investigated under actual outdoor conditions. 3. Experimental data on heat transfer and friction factor compare well those of the analytical data within the range of the parameters investigated. 4. Correlations for Nusselt number and friction factor have been developed in terms of roughness and flow parameters as under: 𝑁𝑢3𝑟 = 0.1415 𝑝/𝑒

−0.0713

𝑒/𝐷

0.0871

𝑅𝑒

0.746

for

𝑒 + ≤ 23

𝑁𝑢3𝑟 = 0.0422 𝑝/𝑒

−0.0461

𝑒/𝐷

0.0375

𝑅𝑒

1.208

for

𝑒 + > 23

𝑓3𝑟 = 0.372 𝑝/𝑒

−0.319

𝑁𝑢3𝑔𝑠 = 0.0782 𝑃𝑟

0.5

𝑒/𝐷 𝑓𝑆

0.7

0.354

𝑅𝑒

𝑅𝑒

0.932

−1.76

5. The values of Stanton number worked out using the experimental value of Nusselt number in three sides roughened and glass covered solar air heater have been utilized to see the effect of roughness and flow parameter on it and represent the thermo hydraulic performance results. 6. The values of Stanton number in three sides artificially roughened and glass covered solar air heaters enhance by 36% to 65% as compared to smooth and three sides glass covered solar air heaters for the same range of values of parameters investigated. 7. The value of heat transfer enhancement factor of three sides roughened and glass covered solar air heater have been found to be 0.378 to 0.487 times over and above the three sides glass covered smooth solar air heater in the range of parameters investigated. 111 | P a g e

8. Collector heat removal factor, 𝐹𝑅 , and collector efficiency factor, 𝐹 ′ , in three sides artificially roughened and glass covered solar air heaters have been found to enhance by about 51 to 60% and 41 to 49% respectively, over those of three sides glass covered smooth solar air heaters. 9. Three sides roughened and glass covered solar air heaters are thermally 53-57%, more efficient than those of three sides glass covered smooth solar air heaters. 10. Correlations for thermo hydraulic efficiency have been developed for three sides glass covered and roughened solar air heater as under: Ƞ𝑡𝑕𝑒𝑟𝑚𝑜 = 0.102 𝑒 +

0.658

for

𝑒 + ≤ 23

Ƞ𝑡𝑕𝑒𝑟𝑚𝑜 = 11.91 𝑒 +

−0.913

for 𝑒 + > 23

+ 11. The value of optimal thermo hydraulic performance corresponding to 𝑒𝑜𝑝𝑡 = 23, have been

found to be 78.8%, within the range of the parameters investigated. 12. Thermal performance equation have been derived for three sides roughened and smooth solar air heaters with three sides glass covers as under: 𝑇𝑖 −𝑇𝑎

ƞ𝑡𝑕(3𝑟,3𝐺) = 0.822 𝜏𝛼 − 0.822𝑈𝐿 ƞ𝑡𝑕(𝑆,3𝐺) = 0.41 𝜏𝛼 − 0.41𝑈𝐿

𝐼

𝑇𝑖 −𝑇𝑎 𝐼

Therefore, three sides glass covered and roughened solar air heaters are thermally more efficient than those of three sides glass covered smooth solar air heaters, and there exists an + optimal thermo hydraulic performance condition at 𝑒𝑜𝑝𝑡 = 23, in three sides roughened

solar air heaters results in the maximum heat transfer for minimum friction factor. 5.2 FUTURE SCOPE OF THE WORK The present work has been investigated for the optimal thermo hydraulic performance of three sides roughened and glass covered solar air heaters by providing transverse wire roughness. The work may be extended to: 1. Use of three sides inclined wire roughness in such solar air heaters may enhance the heat transfer rate and hence may give better thermal and thermo hydraulic performance. 2. Use of booster mirror in such solar air heaters may be the reason behind better performance. 112 | P a g e

REFRENCES Aharwal, K. R., Gandhi, B. k. and Saini, J. S., 2008. Experimental investigation on heat transfer enhancement due to a gap in an inclined continuous rib arrangement in a rectangular duct of solar air heater, Renew. Energy, 33, 585-596. Altermani, C.A.C. and Sparow, E.M., 1980. Turbulent heat transfer and fluid flow in an unsymmetrically heated triangular duct, Trans. ASME, J. Heat Mass Transfer, 102, 590. Behura, K.Arun, Prasad, B. N., Prasad, L., 2016. Heat transfer, friction factor and thermal performance of three sides artificially roughened solar air heaters. Sol. Energy. 130, 46-59. Bernier, M.A. and Plett, E.G., 1988. Thermal performance representation and testing of air solar collectors, Trans. ASME, 110, 74-81. Bhagoria, J. L., Saini, J. S. and Solanki, S. C., 2002. Heat transfer coefficient and friction factor correlations for rectangular solar air heater duct having transverse wedge shaped rib roughness on the absorber plate, Renew. Energy, 25, 341-369. Bhushan, B. and Singh, R., 2011. Nusselt number and friction factor correlations for solar air heater duct having artificially roughened absorber plate, Sol. Energy, 133, 1277-1287. Biondi, P., Cicala, L. and Farina, G., 1988. Performance analysis of solar air heaters of conventional design, Sol. Energy, 41 (1), 101-107. Bliss, R.W., 1959. The derivation of several plate efficiency factors useful in the design of flat plate solar heat collectors, Sol. Energy, 3, 55-64. Bopche, S. B. and Tandale, M. S., 2009. Experimental investigation on heat transfer and frictional characteristics of turbulator roughened solar air heater duct, Int. J. Heat Mass Transfer, 52, 2834-2848. Chaube, A., Sahoo, P. K. and Solanki, S. C., 2006.Analysis of heat transfer augmentation and flow characteristics due to rib roughness over absorber plate of a solar air heater, Renew. Energy, 31 (3), 317-331.

113 | P a g e

Cho, H. H., Kim, Y. Y., Rhee, D. H. and Lee, S. Y., 2003.The effect of gap position in discrete ribs on local heat/mass transfer in a square duct, J. Enhanced Heat Trans., 10, 287300. Champookham, T., Thianpong, C., Kwankaomeng, S. and Promvonge, P., 2010. Heat transfer augmentation in a wedge ribbed channel using winglet vortex generators, Int. Commun Heat Mass Transfer, 37, 163-169. Chauhan, R. and Thakur, N. S., 2013. Heat transfer and friction factor correlations for imping jet solar air heater, Exp. Therm. Fluid Sci., 14, 760-767. Chouksey, V.K. and Sharma, S.P., 2016. Investigations on thermal performance characteristics of wire screen packed bed solar air heater, Sol. Energy, 132, 511-605. Dipprey, D. F. and Sabersky, R. H., 1963. Heat and momentum transfer in smooth and rough tubes at various Prandtl numbers, Int. J. Heat Trans. 6, 329-353. Duffie, J.A. and Beckman, W.A., 1980. Solar energy thermal processes, Willey Interscience, New York. Eiamsa-ard, S., Wongcharee, K., Eiamsa-ard, P. and Thianpong, C., 2010. Heat transfer enhancement in a tube using delta winglet twisted tapes inserts, Appl. Therm. Engg., 30, 310-318. Gandhi, B. K. and Singh, K. M., J-Mc, 2010. Experimental and numerical investigations on flow through wedge shape rib roughened duct, Inst. Engg. (India), 90, 13-18. Gillet, W.B., Aranviteh, E. and Moon, J.E., 1983. Solar collector testing in the European Community, Int. J. Sol. Energy, 1, 317-341. Gupta, D. and Solanki, S. C., 1993. Heat and fluid flow in rectangular solar air heater ducts having transverse rib roughness on absorber plates, Sol. Energy, 51 (1), 31-37. Gupta, M. K. and Solanki, S. K., 2009.Performance evaluation of solar air heater having expanded metal mesh as artificial roughness on absorber plate, Int. J. Therm. Sci., 48, 10071016. 114 | P a g e

Gupta, D., Solanki, S. C. and Saini, J. S., 1997. Thermo-hydraulic performance of solar air heaters with roughened absorber plates, Sol. Energy. 61, 33-42. Han, J.C., 1984. Heat transfer and friction channels with two opposite rib-roughened walls, Trans. ASME J. of Heat Trans. 106, 774-781. Hans V. S., Saini, R. P. and Saini, J. S., 2010. Heat transfer and friction factor correlations for a solar air heater duct roughened artificially with multiple V-ribs, Sol. Energy, 84, 898911. Han, J. C., Zhang, Y. M. and Lee, C. P., 1991.Augmented heat transfer in square channels with parallel, crossed and V-shaped angled ribs, Trans., ASME, J. Heat Trans., 113, 590596. Hegazy, A.A., 1995. Optimization of flow channel depth for conventional flat plate solar air heaters, Renew. Energy, 1, 15-21. Hegazy, A.A., 1999. Optimizing the thermo hydraulic performance of flat plate solar air heaters operating with fixed/variable pumping power, Renew. Energy, 18, 283-304. Huang, K. D., Tzeng, S.C., Jeng, T. M., Wang, J. R. Cheng, S. Y. and Tseng, K. T., 2008. Experimental study of fluid flow and heat transfer characteristics in the square channel with a perforation baffle, Int. Commun Heat Mass Transfer, 36, 264-268. Jaurkur, A. R., Saini, J. S. and Gandhi, B. K., 2006. Heat transfer and friction factor characteristics of rectangular solar air heater duct using rib-grooved artificially roughness, Sol. Energy, 80 (8), 895-907. Karwa, R., Solanki, S.C., Saini, J.S., 1999. Heat transfer coefficient and friction factor correlation for the transitional flow regime in rib roughened rectangular ducts. Int. J. Heat Mass Trans. 42, 1597–1615. Karwa, R., Solanki, S.C., Saini, J.S., 2001. Thermo hydraulic performance of solar air heaters having integral chamfered rib roughness on absorber plates, Energy, 26,161-176.

115 | P a g e

Karwa, R., 2003. Experimental studies of augmented heat transfer and friction in asymmetrical heated rectangular ducts with ribs on the heated wall in transverse inclined, Vcontinuous and V-discrete pattern, Int. Commun Heat Mass Transfer, 30, 241-250. Karwa, R. and Maheshwari, B. K., 2009. Heat transfer and friction in an asymmetrically heated rectangular duct with half and fully perforated baffles at different pitches, Int. Commun Heat Mass Transfer, 36, 264-268. Karmare, S. V. and Tikekar, A. N., 2007. Heat transfer and friction factor correlation for artificially roughened duct with metal grift ribs, Int. J. Heat Mass Transfer, 50, 4342-4351. Karmare, S. V. and Tikekar, A. N., 2009. Experimental investigation of optimum thermo hydraulic performance of solar air heaters with metal rib grits roughness, Sol. Energy. 83, 6-13. Karmare, S. V. and Tikekar, A. N., 2010. Analysis of fluid flow and heat transfer in a rib grid roughened solar air heater using CFD, Sol. Energy, 84 (3), 409-417. Kays, W.M., 1966. Convective Heat and Mass Transfer, McGraw Hill, New York. Kays, W.M. and London, A.L., 1964. Compact Heat Exchangers, McGraw Hill, New York. Khalid, A.J. and Shaker, M.M., 1987. An experimental investigation into a solar assisted desiccant-evaporative air conditioning system, Sol. Energy, 39 (2), 97-107. Kumar, S. and Saini, R. P., 2009. CFD based performance analysis of a solar air heater duct provided with artificial roughness, Renew. Energy, 34, 1285-1291. Kumar, S. T., Mittal, V., Thakur, N. S. and Kumar, A., 2011. Heat transfer and friction factor correlations for rectangular solar air heater duct having 600 inclined continuous discrete rib arrangements, Br. J. Appl. Sci. Technol., 3, 67-93. Kwankaomeng, S., Jed Sadaratanachri, W. and Promvong, P., 2-4 June, 2010. Laminar periodic flow and heat transfer in square channel with 300 inclined baffles, PEA-AIT, International conference on energy and sustainable development: Issues and strategies (ESD 2010), The Empress Hotel, Chiang Mai, Thailand. 116 | P a g e

Lanjewar, A., Bhagoria, J. L. and Sarviya, R. M., 2011.Experimental study of augmented heat transfer and friction in solar air heater with different orientations of W-rib roughness, EXP Therm. Fluid Sci., 35 986-995. Layek, A., Saini, J. S, and Solanki, S. C., 2007.Second law optimization of a solar air heater having chamfered rib–groove roughness on absorber plate, Renew. Energy, 32, 1967-1980. Lau, S. C., Kakraja, R. T. and McMillin, R. D., 1991 (a).Effect of V-shaped rib arrays on turbulent heat transfer and friction of fully developed flow in a square channel, Int. L. Heat Mass Transfer, 34, 1605-1616. Lau, S. C., McMillin, R. D. and Han, J. C., 1991 (b).Heat transfer characteristics of turbulent flow in a square channel with angle discrete ribs, ASME, J. Turbo-mech, 113, 367-374. Lau, S. C., McMillin, R. D. and Han, J. C., 1991 (c). Turbulent heat transfer and friction in a square channel with discrete and turbulators, ASME, J. Turbo-mech, 113, 360-366. Lewis, M. J., 1975 (a).An elementary analysis for predicting the momentum and heat transfer characteristics of hydraulically rough surface, Trans. ASME J. of Heat Trans. 97, 249-254. Lewis, M. J., 1975 (b).Optimizing the thermo hydraulic performance of rough surfaces, Int. J. Heat Mass Trans. 18, 1243-1248. Malik, M.A.S. and Buelow, F.H., 1975. Hydro-dynamic and heat transfer characteristics of a heated duct, Helio-technique and Development, 2, 3. Mittal, M. k. and Varshney, L., 2005. Optimal thermo hydraulic performance of a wire mesh Packed solar air heater, Sol. Energy, 80, 1112-1120. Momin, A. M. E., Saini, J. S. and Solanki, S. C., 2002. Heat transfer and friction in solar air heater duct with V-shaped rib roughness on absorber plate, Int. J. Heat Mass Transfer, 45, 3383-3396. Muluwork, K. B., Saini, J. S. and Solanki, S. C., 1998. Studies on discrete rib roughened solar air heaters, In: Proceedings of national solar energy convention, 75-84. 117 | P a g e

Muluwork, K. B., 2000. Investigations on fluid flow and heat transfer in roughened absorber solar heaters, [PhD. Thesis] IIT Roorkee. Ozisik, M. N., 1985. Heat transfer: A basic approach, McGraw Hill, Book Company, Singapore. Paswan, M. K. and Sharma, S. P., 2009. Thermal performance of wire- mesh roughened solar air heaters, Arab Res. Inst. Sci. Engg., 1, 31-40. Prasad, B. N. and Saini, J. S., 1988.Effect of artificial roughness on heat transfer and friction factor in a solar air heater, Sol. Energy, 41, 555-560. Prasad, B. N. and Saini, J.S., 1991.Optimal thermo-hydraulic performance of artificially roughened solar air heaters, Sol. Energy, 47, 91-96. Prasad, B. N., 2013.Thermal performance of artificially roughened solar air heaters, Sol. Energy, 91, 59-67. Prasad, B. N., Behura, K. Arun. and Prasad, L., 2014. Fluid flow and heat transfer analysis for heat transfer enhancement in three sided artificially roughened solar air heater, Sol. Energy, 105, 27-35. Prasad, K. and Mullick, S.C., 1985. Heat transfer characteristics of a solar air heater used for drying purposes, Appl. Energy 13, 83–93. Priyam, A. and Chand, P., 2016. Thermal and thermo hydraulic performance of wavy finned absorber solar air heater, Sol. Energy, 130, 250-259. Promvong, P., Kwankaomeng, S. and Thianpong, C., 2011. Thermal behavior in solar air heater channel fitted with combined rib and delta winglet, Int. Commun Heat Mass Transfer, 38, 749-756. Promvong, P., Jedsadaratanachai, W., Kwankaomeng, S. and Thianpong, C., 2012.3D simulation of laminar flow and heat transfer in V-baffled square channel, Int. Commun Heat Mass Transfer, 39, 85-93.

118 | P a g e

Promvong, P., Stripattanapipat, S. and Kwankaomeng, S., 2010. Laminar periodic flow and heat transfer in square channel with 450 in-line baffles on two opposite walls, Int. J. Therm. Sci., 49, 963-975. Reddy, T.A. and Gupta, C.L., 1980.Generating application design data for solar air heating system, Sol. Energy, 25, 527-530. Saha, S.K. and Mahanta, D.K., 2001. Thermodynamic optimization of solar flat-plate collector, Renew. Energy, 23, 181-193. Sahu, M. M. and Bhagoria, J. L., 2005. Augmentation of heat transfer coefficient by using 900 broken transverse ribs on absorber plate of solar air heater, Renew. Energy, 30, 20632075. Saini, R.P. and Saini, J.S., 1997.Heat transfer and friction factor correlations for artificially roughened duct with expanded metal mesh as roughness element, Int. journal of Heat Mass Trans. 40, 973-986. Sheriff, N and Gumley, P., 1966. Heat transfer and friction properties of surfaces with discrete roughness, Int. J. Heat Mass Trans. 9, 1297-1320. Sharma, A. K. and Thakur, N. S., 2012. CFD based fluid flow and heat transfer analysis of a V-shaped roughened surface solar air heater, Int. J. Engg. Sci. Technol., 4 (5), 2115-2121. Sharma, S. K. and Kalamkar, V. R., 2015.Thermo-hydraulic performance analysis of solar air heaters having artificial roughness–A review, Renewable and sustainable energy reviews., 41, 413-435. Shewen, E.C. and Hollands, K.G.T., March-1980.Development of a fixed position, fixed configuration air heating solar collector, National Research council of Canada. Shewen, E.C. and Hollands, K.G.T., 1978.Optimization of flow passage geometry in solar air heaters, Paper 1-1-2, Proc. of Solar Energy Society of Canada, Inc. London. Siddhartha., Sharma, N. and Varun., 2012. A practical swarm optimization algorithm for optimization of thermal performance of a smooth flat plate solar air heater, 38, 406-413. 119 | P a g e

Sriromrein, P. and Promvong, P., 2-4 June, 2010. Augmented heat transfer in rectangular duct with angled Z-shaped ribs, In: Proceedings of the International Conference on energy and sustainable development, Thailand. Sukhatme, S.P. and Nayak, J.K., 2008. Solar Energy, Principles of Thermal Collection and Storage, Tata McGraw Hill, New Delhi Ted, J, Jansen, 1985. Solar Engineering Technology, Prentice Hall, New Jersey. Varun, Saini, R. P. and Singal, S. K., 2008.Investigation of thermal performance of solar air heater having roughness elements as a combination of inclined and transverse ribs on the absorber plate, Renew. Energy, 33, 1398-1405. Verma, R., Chandra, R. and Garg, H. P., 1991.A technical note on optimization of solar air heaters of different designs, Renew. Energy, 2, (4/5), 521-531. Verma, S.K. and Prasad, B.N., 2000. Investigation for the optimal thermo hydraulic performance of artificially roughened solar air heaters, Renew. Energy, 20, 19-36. Webb, R.L., Eckert, E.R.G and Goldstein, R.J., 1971. Heat transfer and friction in tubes with repeated-rib roughness, Int. J. of Heat and Mass Trans. 14, 601-617. Webb, R. L. and Eckert, E. R. G., 1971.Application of rough surfaces to heat exchanger design, Int. J. Heat Mass Trans. 15, 1647-1658. Yadav, S., ManeeshKausal, Varun and Siddhartha., 2012. Nusselt number and friction factor correlations for solar air heater duct having protrusions as roughness elements on absorber. Yadav, A. S. and Bhagoria, J. L., 2013a. Modeling and simulation of turbulent flows through a solar air heater having square-sectioned transverse rib roughness on the absorber plate, Sci. world, J. Yadav, A. S. and Bhagoria, J. L., 2013b.Numerical investigation of flow through an artificially roughened solar heater, Int. J. Ambient Energy.

120 | P a g e

Yadav, A. S. and Bhagoria, J. L., 2013c. A CFD Based heat transfer and fluid flow analysis of a solar air heater provided with circular transverse wire rib roughness on the absorber plate, Energy, 55, 1127-1142. Yadav, A. S. and Bhagoria, J.L., 2014. A CFD based thermo hydraulic performance analysis of an artificially roughened solar air heater having equilateral triangular sectioned rib roughness on the absorber plate, Int. J. Heat Mass Trans. 70, 1016-1039.

121 | P a g e

UNCERTAINTY ANALYSIS Uncertainty analysis has been done by using the method [J.P.Holman], given by the following form of equation: 𝑈𝑥 =

𝜕𝑥 𝜕𝑦 1

𝑈𝑦1

2

+

𝜕𝑥 𝜕𝑦 2

𝑈𝑦2

2

+ ……………+

𝜕𝑥 𝜕𝑦 𝑛

2

𝑈𝑦𝑛

(A-1)

where, 𝑦1 , 𝑦2 … . 𝑦𝑛 are the variables affecting the parameters X and U stands for the uncertainty. One of the most important parameter in this work is the thermal efficiency calculation, given by the following equation: ƞ𝑡𝑕 = 𝑚𝐶𝑝 (𝑇0 − 𝑇𝑖 )/𝐼𝐴𝐶

(A-2)

From Eqs. (A-1) and (A-2), the uncertainty in the thermal efficiency of solar air heaters, U is given by: 𝑈Ƞ = −

𝐶𝑝 𝑇0 −𝑇𝑖 𝐼𝐴𝐶

𝑚 𝐶𝑃 𝑇0 −𝑇𝑖 𝐴𝐶 𝐼 2

𝑈𝐼

2

𝑈𝑚

+

𝑚 𝑇0 −𝑇𝑖 𝐼𝐴𝐶

1 2

𝑈𝐶𝑝

2

+

𝑚 𝐼𝐴𝐶

𝑈 𝑇0 −𝑇𝑖

2

+ −

𝑚 𝐶𝑃 𝑇0 −𝑇𝑖 𝐼𝐴𝐶 2

𝑈𝐴𝐶 + (A-3)

Where, 𝑈𝑚 , 𝑈𝐶𝑝 , 𝑈(𝑇0 −𝑇𝑖 ) , 𝑈𝐴𝐶 and 𝑈𝐼 are the uncertainties in mass flow rate of air, specific heat of air, air temperature rise, area of collector and intensity of solar radiation respectively. A.1 EXPERIMENTAL DATA For a particular test run (Test Run no. 56) typical data have been used to illustrate and workout for uncertainties, which are length of test section of solar air heater, L = 1500mm, width of test length of solar air heater, W = 200mm, internal diameter of orifice pipe, d1 = 76mm, throat diameter of orifice plate, d2 = 38mm, pressure drop across orifice meter, ∆p = 22mm of water, air temperature rise, (T0 – Ti), intensity of solar radiation, I = 899 W/m2. Dimensions of orifice pipe and throat diameters are measured with the help of vernier caliper for which uncertainties are: 𝑈𝑑1 = 𝑈𝑑2 = 0.02 𝑚𝑚

(A-4) 122 | P a g e

Length and width of solar air heater duct are measured with standard steel scale for least count equal to 1mm for which uncertainties are: 𝑈𝐿 = 𝑈𝑤 = 1 𝑚𝑚

(A-5)

Uncertainties in pressure drop and temperature are: 𝑈∆𝑝 = 1𝑚𝑚

(A-6)

𝑈𝑇0 −𝑇𝑖 = 0.10 𝐶

(A-7)

The value of intensity of solar radiation has been taken from standard digital pyranometer and therefore, the uncertainty in intensity of solar radiation is: 𝑈𝐼 = 0

(A-8)

A.2 UNCERTAINTY IN MASS FLOW RATE MEASUREMENT Mass flow rate through the duct has been calculated by the Eq. (3.4) and rewritten as: m = Cd A2

2ρ P 1 −P 2 1/2

(A-9)

1−β 41

Where, ∆p is pressure drop orifice meter. An uncertainty in the mass flow rate is given by: 2

𝑈𝑚 =

Cd

Cd A2 2ρ

1 A 2 2∆p 2 1 1 ρ 2 1−β 41 2

1

U∆p

1

2 1−β 41 ∆p 2

+

1−β 41

1 2

2

UA 2

2

+ 2Cd A2 2ρ∆p

1 2

β 31

3

2 1−β 41 2

Uβ 1

+

1

2



2ρ∆p

+

A 2 2ρ∆p 1

1−β 41 2

1 2

2 2

UC d

(A-10)

Where, U∆p , UA 2 , Uβ 1 , Uρ and UC d are uncertainties in pressure drop across orifice meter, area of throat, diameter ratio, density of air and coefficient of discharge of orifice meter respectively. The values of density of air and coefficient of discharge are taken from standard tables and therefore, the uncertainties are: Uρ = 0

(A-11)

123 | P a g e

UC d = 0

(A-12)

Uncertainties in UA 2 and Uβ 1 have been worked out to be 0.0003×10-5 and 17.04×10-7 respectively. Therefore, uncertainties in the mass flow rate work out to be: 𝑈𝑚 = 2.02287 × 10−9

(A-13)

A.3 UNCERTAINTY IN COLLECTOR AREA 𝑈𝐴𝐶 =

𝐿𝑈𝑤

2

+ 𝑊𝑈𝐿

2 1/2

(A-14)

= [(1500×1)2 + (200×1)2]1/2 = 1.5133 × 10-3 A.4 UNCERTAINTY IN SPECIFIC HEAT The value of specific heat of air has been taken from standard table, therefore, the uncertainty in specific heat is: 𝑈𝐶𝑝 = 0

(A-15)

A.5 UNCERTAINTY INTHERMAL EFFICIENCY When substituted for the values of various quantities in Eq. (A-3), the maximum uncertainty in the value of thermal efficiency is worked out to be 1.76%.

124 | P a g e

LIST OF PUBLICATIONS FROM THIS WORK 1. Prasad B.N., Kumar Ashwini and Singh K.D.P., 2015. Optimization Of Thermo Hydraulic Performance In Three Sides Artificially Roughened Solar Air Heaters, Sol. Energy, 111, 313-319. 2. Kumar Ashwini, Prasad B.N., and Singh K.D.P., 2016. Performance Characteristics Of Three Sides Glass Covered Smooth Solar Air Heater, Transylvanian Review, 24, (11), 32473256. 3. Kumar Ashwini, Prasad B.N. and Singh K.D.P., 2018. Development Of Correlations Of Heat Transfer And Fluid Flow Characteristics For Three Sides Artificially Roughened Solar Air Heaters, J P Journal of Heat and Mass Transfer, 15, (1), 93-123. 4. Kumar Ashwini, Prasad B.N., and Singh K.D.P., 2017. Thermal And Thermo Hydraulic Performance Of Three Sides Artificially Roughened Solar Air Heaters, Int. Res. J. of Advanced Engineering and Science, 02, (2), 215-231. 5. Kumar Ashwini, Prasad B.N., and Singh K.D.P., 2017. Analysis Of Collector Performance Parameters In Three Sides Artificially Roughened And Glass Covered Solar Air Heater, Int. Res. J. of Advanced Engineering and Science, 02, (3), 224-229. 6. Kumar Ashwini, Prasad B.N., and Singh K.D.P., 2017. Heat Transfer And Fluid Flow Characteristics Of Three Sides Artificially Roughened And Glass Covered Solar Air Heaters, Int. Res. J. of Advanced Engineering and Science, 02, (3), 230-244. 7. Kumar Ashwini, Prasad B.N and Singh K.D.P., 2016. Performance Characteristics Of Three Sides Artificially Roughened Solar Air Heaters, proceedings, Global Conference on Renewable Energy, National Institute of Technology, Patna, India, GCRE-2016, 14, 132137. 8. Kumar Ashwini, Singh K.D.P. and Prasad B.N., 2016. Enhancement Of Collector Performance Parameters In Three Sides Artificially Roughened Solar Air Heater, proceedings, International Conference on Advances in Steel, Power and Construction Technology, O P Jindal University, Raigarh, India. 9. Kumar Ashwini, Prasad B.N., and Singh K.D.P., 2017. Investigation For Optimal Thermo Hydraulic Performance Of Three Sides Roughened And Glass Covered Solar Air Heater, Int. J. of Green Energy, Taylor and Francis. (Under Review).

125 | P a g e

10. Kumar Ashwini, Prasad B.N. and Singh, K.D.P., 2018. Optimal Thermo Hydraulic Performance Of Three Sides Roughened And Glass Covered Solar Air Heaters, Int. J. of Heat and Fluid Flow, Elsevier. (Under Review). 11. Kumar, Ashwini, Prasad, B.N., Singh, K.D.P., 2018. Thermo Hydraulic Performance Of Three Sides Roughened And Glass Covered Solar Air Heater, International Energy Journal, Thailand. (Under Review).

126 | P a g e