Optimum Fin Profile under Dry and Wet Surface Conditions

Similar documents
OPTIMIZATION OF THE LONGITUDINAL FINS WITH DIFFERENT GEOMETRIES FOR INCREASING THE HEAT TRANSFER

Fluid flow consideration in fin-tube heat exchanger optimization

TOPICAL PROBLEMS OF FLUID MECHANICS 17 ONE-DIMENSIONAL TEMPERATURE DISTRIBUTION OF CONDENSING ANNULAR FINS OF DIFFERENT PROFILES

Modelling a Spiral Porous Fin with Temperature Dependent and Independent Thermal Conductivity

Available online at ScienceDirect. International Conference On DESIGN AND MANUFACTURING, IConDM 2013

HEAT TRANSFER AND TEMPERATURE DISTRIBUTION OF DIFFERENT FIN GEOMETRY USING NUMERICAL METHOD

DIFFERENTIAL TRANSFORMATION METHOD TO DETERMINE TEMPERATURE DISTRIBUTION OF HEAT RADIATING FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY

INVESTIGATIONS ON THE INTERNAL SHAPE OF CONSTRUCTAL CAVITIES INTRUDING A HEAT GENERATING BODY

Safa Bozkurt Coşkun and Mehmet Tarik Atay. Received 13 December 2006; Revised 11 April 2007; Accepted 22 September 2007 Recommended by Josef Malek

Introduction to Heat and Mass Transfer. Week 5

Convective Mass Transfer

Maximal Thermo-geometric Parameter in a Nonlinear Heat Conduction Equation

Performance of Annular Fin with Different Profile Subjected to Heat Transfer Coefficient by using FEA

CONVECTION WITH SIMPLIFIED CONSTRAINTS 1.1 INTRODUCTION Three quarters of a century ago, a paper by Harper and Brown (1922) appeared as an NACA

Estimation of Heat Dissipation from Plate with Multiple Tapered and Rectangular Fins

Applied Thermodynamics for Marine Systems Prof. P. K. Das Department of Mechanical Engineering Indian Institute of Technology, Kharagpur

Experiment 1. Measurement of Thermal Conductivity of a Metal (Brass) Bar

Study of Temperature Distribution Along the Fin Length

Chapter 9 NATURAL CONVECTION

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

Level 7 Post Graduate Diploma in Engineering Heat and mass transfer

ANALYSIS OF A ONE-DIMENSIONAL FIN USING THE ANALYTIC METHOD AND THE FINITE DIFFERENCE METHOD

EXPERIMENTAL AND NUMERICAL STUDIES OF A SPIRAL PLATE HEAT EXCHANGER

Numerical Study of Laminar Free Convection About a Horizontal Cylinder with Longitudinal Fins of Finite Thickness

Minimum volume of the longitudinal fin with rectangular and triangular profile by a modified Newton-Raphson method Nguyen Quan¹,

ME6301- ENGINEERING THERMODYNAMICS UNIT I BASIC CONCEPT AND FIRST LAW PART-A

A NUMERICAL APPROACH FOR ESTIMATING THE ENTROPY GENERATION IN FLAT HEAT PIPES

Maximum Heat Transfer Density From Finned Tubes Cooled By Natural Convection

Thermal Analysis of Convective-Radiative Circular Porous Fins with Sinusoidal Cross-Section Geometries using Galerkin Method

CHAPTER 4 MATHEMATICAL MODELING AND SIMULATION

Experimental Analysis of Wire Sandwiched Micro Heat Pipes

A REPORT ON PERFORMANCE OF ANNULAR FINS HAVING VARYING THICKNESS

Documentation of the Solutions to the SFPE Heat Transfer Verification Cases

PREDICTION OF MASS FLOW RATE AND PRESSURE DROP IN THE COOLANT CHANNEL OF THE TRIGA 2000 REACTOR CORE

The Numerical Psychrometric Analysis

[Yadav*, 4.(5): May, 2015] ISSN: (I2OR), Publication Impact Factor: (ISRA), Journal Impact Factor: 2.114

Characteristics of forced convection heat transfer in vertical internally finned tube B

HEAT TRANSFER FROM FINNED SURFACES

4. Analysis of heat conduction

Mathematical Modelling for Refrigerant Flow in Diabatic Capillary Tube

Optimization of Peripheral Finned-Tube Evaporators Using Entropy Generation Minimization

Minhhung Doan, Thanhtrung Dang

Chapter 10: Steady Heat Conduction

Chapter 10: Boiling and Condensation 1. Based on lecture by Yoav Peles, Mech. Aero. Nuc. Eng., RPI.

Chapter 4 Water Vapor

Analyzing Mass and Heat Transfer Equipment

The Choice of Supply Design Conditions

Keywords: Spiral plate heat exchanger, Heat transfer, Nusselt number

C ONTENTS CHAPTER TWO HEAT CONDUCTION EQUATION 61 CHAPTER ONE BASICS OF HEAT TRANSFER 1 CHAPTER THREE STEADY HEAT CONDUCTION 127

S.E. (Chemical) (Second Semester) EXAMINATION, 2012 HEAT TRANSFER (2008 PATTERN) Time : Three Hours Maximum Marks : 100

10 minutes reading time is allowed for this paper.

Chapter 3: Steady Heat Conduction. Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University

1. Water Vapor in Air

Heat and Mass Transfer Unit-1 Conduction

INTRODUCTION: Shell and tube heat exchangers are one of the most common equipment found in all plants. How it works?

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

Optimization of flue gas turbulent heat transfer with condensation in a tube

7LAMINAR CONVECTIVE HEAT TRANSFER FOR IN-PLANE SPIRAL COILS OF NON-CIRCULAR CROSS-SECTIONS DUCTS A Computational Fluid Dynamics Study

heat transfer process where a liquid undergoes a phase change into a vapor (gas)

Introduction to Heat and Mass Transfer. Week 12

Numerical Heat and Mass Transfer

Simulation of Heat and Mass Transfer in the Corrugated Packing of the Counter Flow Cooling Tower

Convective Heat and Mass Transfer Prof. A.W. Date Department of Mechanical Engineering Indian Institute of Technology, Bombay

Chapter 6 Fundamental Concepts of Convection

ENTROPY GENERATION DUE TO EXTERNAL FLUID FLOW AND HEAT TRANSFER FROM A CYLINDER BETWEEN PARALLEL PLANES

The Effect Of MHD On Laminar Mixed Convection Of Newtonian Fluid Between Vertical Parallel Plates Channel

A SIMILARITY SOLUTION OF FIN EQUATION WITH VARIABLE THERMAL CONDUCTIVITY AND HEAT TRANSFER COEFFICIENT

Chapter 5 Time-Dependent Conduction

Pin Fin Lab Report Example. Names. ME331 Lab

COOLING TOWER MODEL IN MNR

Film condensation on horizontal tube with wall suction effects

Axial profiles of heat transfer coefficients in a liquid film evaporator

NUMERICAL INVESTIGATION AND STATISTICAL ANALYSIS OF LASER BENDING OF TITANIUM SHEETS

FLOW BOILING HEAT-TRANSFER IN PLATE MICRO- CHANNEL HEAT SINK

Subject: Principles of Refrigeration and Air Conditioning Lecturer: Assistant Professor Dr. Waheed Shaty Mohammed

Theory. Humidity h of an air-vapor mixture is defined as the mass ratio of water vapor and dry air,

Heat Transfer Characteristics and Performance of a Spirally Coiled Heat Exchanger under Sensible Cooling Conditions

R13 SET - 1 '' ''' '' ' '''' Code No RT21033

Meteorology 6150 Cloud System Modeling

Pressure Distribution of Refrigerant Flow in an Adiabatic Capillary Tube

Free Convective Heat Transfer From A Vertical Surface For The Case Of Linearly Varying Thermal Potential

Thermal Systems. What and How? Physical Mechanisms and Rate Equations Conservation of Energy Requirement Control Volume Surface Energy Balance

We are IntechOpen, the first native scientific publisher of Open Access books. International authors and editors. Our authors are among the TOP 1%

Basic Thermodynamics Prof. S. K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur. Lecture No 16

NUMERICAL STUDY OF HEAT AND MASS TRANSFER DURING EVAPORATION OF A THIN LIQUID FILM

Diffusional Growth of Liquid Phase Hydrometeros.

Heat processes. Heat exchange

IHTC DRAFT MEASUREMENT OF LIQUID FILM THICKNESS IN MICRO TUBE ANNULAR FLOW

(Refer Slide Time: 00:00:43 min) Welcome back in the last few lectures we discussed compression refrigeration systems.

FINITE ELEMENT ANALYSIS OF THERMAL CHARACTERISTICS OF ANNULAR FINS WITH DIFFERENT PROFILES

A COMPUTATIONAL ALGORITHM FOR THE ANALYSIS OF FREEZE DRYING PROCESSES WITH SPECIAL REFERENCE TO FOOD PRODUCTS / MATERIALS

Experimental Analysis on Pin Fin Heat Exchanger Using Different Shapes of the Same Material

FINITE ELEMENT ANALYSIS OF MIXED CONVECTION HEAT TRANSFER ENHANCEMENT OF A HEATED SQUARE HOLLOW CYLINDER IN A LID-DRIVEN RECTANGULAR ENCLOSURE

Arctice Engineering Module 3a Page 1 of 32

Chapter 3 NATURAL CONVECTION

DETERMINATION OF R134A S CONVECTIVE HEAT TRANSFER COEFFICIENT IN HORIZONTAL EVAPORATORS HAVING SMOOTH AND CORRUGATED TUBES

Performance evaluation of heat transfer enhancement for internal flow based on exergy analysis. S.A. Abdel-Moneim and R.K. Ali*

PROBLEM (a) Long duct (L): By inspection, F12. By reciprocity, (b) Small sphere, A 1, under concentric hemisphere, A 2, where A 2 = 2A

1 Introduction. EPJ Web of Conferences 110,

GRAVITY EFFECT ON THE DISTRIBUTION OF REFRIGERANT FLOW IN A MULTI-CIRCUITED CONDENSER

Transcription:

Optimum Fin Profile under Dry Wet Surface Conditions Balaram Kundu Somchai Wongwises Department of Mechanical Engineering Jadavpur University, Kolkata 7 3 Fluid Mechanics, Thermal Engineering Multiphase Flow Research ab (FUTURE), Department of Mechanical Engineering King Mongkut s University of Technology Thonburi (KMUTT) Bangmod,Bangkok 4 India Thail. Introduction Fins or extended surfaces are frequently employed in heat exchangers for effectively improving the overall heat transfer performance. The simple design of fins their stability in different surface conditions have created them a popular augmentation device. The different fin shapes are available in the literature. The geometry of the fin may be dependent upon the primary surface also. For circular primary surface, the attachment of circumferential fins is a common choice. The longitudinal pin fins are generally used to the flat primary surface. However, due to attachment of fins with the primary surface, the heat transfer augments but the volume, weight, cost of the heat exchanger equipments increase as well. Hence, it is a challenge to the designer to minimize the cost for the attachment of fins. This can be done by determining the optimum shape of a fin satisfying the maximization of heat transfer rate for a given fin volume. In general, two different approaches are considered for the optimization of any fin design problem. Through a rigorous technique, the profile of a fin for a particular geometry (flat or curved primary surface) may be obtained such that the criteria of the maximum heat transfer for a given fin volume or equivalently minimum fin volume for a given heat transfer duty is satisfied. In a parallel activity, the optimum dimensions of a fin of given profile (rectangular, triangular etc.) are determined from the solution of the optimality criteria. The resulting profile obtained from the first case of optimum design is superior in respect to heat transfer rate per unit volume thus it is very much important in fin design problems. However, it may be limited to use in actual practice because the resulting profile shape would be slightly difficult to manufacture fabricate. Alternatively, such theoretical shape would first be calculated then a triangular profile approximating the base two thirds of the fin would be used. Such a triangular fin transfers heat per unit weight, which is closer to that of the analytical optimum value. Under a convective environmental condition, Schmidt (96) was the first researcher to forward a systematic approach for the optimum design of fins. He proposed heuristically that for an optimum shape of a cooling fin, the fin temperature must be a linear function

4 Heat Transfer - Mathematical Modelling, Numerical Methods Information Technology with the fin length. ater, through the calculus of variation, Duffin (959) exhibited a rigorous proof on the optimality criteria of Schmidt. iu (96) extended the variational principle to find out the optimum profile of fins with internal heat generation. iu (96) Wilkins (96) addressed for the optimization of radiating fins. Solov ev (968) determined the optimum radiator fin profile. The performance parameter of annular fins of different profiles subject to locally variable heat transfer coefficient had been investigated by Mokheimer (). From the above literature works, it can be indicated that the above works were formulated based on the length of arc idealization (AI). Maday (974) was the first researcher to eliminate AI obtained the optimum profile through a numerical integration. It is interesting to note that an optimum convecting fin neither has a linear temperature profile nor possesses a concave parabolic shape suggested by Maday. The profile shape contains a number of ripples denoted as a wavy fin. The same exercise was carried out for radial fins by Guceri Maday (975). ater Razelos Imre (983) applied Pontryagin's minimum principle to find out the minimum mass of convective fins with variable heat transfer coefficient. Zubair et al. (996) determined the optimum dimensions of circular fins with variable profiles temperature dependent thermal conductivity. They found an increasing heat transfer rate through the optimum profile fin by % as compare to the constant thickness fin. A variational method was adopted by Kundu Das (998) to determine the optimum shape of three types of fins namely the longitudinal fin, spine disc fin. A generalized approach of analysis based on a common form of differential equations a set of boundary conditions had been described. For all the fin geometries, it was shown that the temperature gradient is constant the excess temperature at the tip vanishes. By taking into account the AI, Hanin Campo (3) forecasted a shape of a straight cooling fin for the minimum envelop. From the result, they have highlighted that the volume of the optimum circular fin with consideration of AI found is 6.-8 times smaller than the volume of the corresponding Schmidt s parabolic optimum fin. A new methodological determination for the optimum design of thin fins with uniform volumetric heat generation had been done by Kundu Das (5). There are ample of practical applications in which extended surface heat transfer is involved in two-phase flow conditions. For example, when humid air encounter into a cold surface of cooling coils whose temperature is maintained below the dew point temperature, condensation of moisture will take place, mass heat transfer occur simultaneously. The fin--tube heat exchangers are widely used in conventional air conditioning systems for air cooling dehumidifying. In the evaporator of air conditioning equipment, the fin surface becomes dry, partially or fully wet depending upon the thermogeometric psychrometric conditions involved in the design process. If the temperature of the entire fin surface is lower than the dew point of the surrounding air, there may occur both sensible latent heat transferred from the air to the fin so the fin is fully wet. The fin is partially wet if the finbase temperature is below the dew point while fin-tip temperature is above the dew point of the surrounding air. If the temperature of the entire fin surface is higher than the dew point, only sensible heat is transferred so the fin is fully dry. For wet surface, the moisture is condensed on the fin surface, latent heat evolves mass transfer occurs simultaneously with the heat transfer. Thermal performance of different surface conditions of a fin depends on the fin shape, thermophysical psychrometric properties of air. Many investigations have been devoted to analyze the effect of condensation on the performance of different geometric fins. It is noteworthy to mention that for each instance, a

Optimum Fin Profile under Dry Wet Surface Conditions 5 suitable fin geometry has been selected a priory to make the analysis. For the combined heat mass transfer, the mathematical formulation becomes complex to determine the overall performance analysis of a wet fin. Based on the dry fin formula, Threlkeld (97) McQuiston (975) determined the one-dimensional fin efficiency of a rectangular longitudinal fin for a fully wet surface condition. An analytical solution for the efficiency of a longitudinal straight fin under dry, fully wet partially wet surface conditions was introduced elaborately by Wu Bong (994) first with considering temperature humidity ratio differences as the driving forces for heat mass transfer. For the establishment of an analytical solution, a linear relationship between humidity ratio the corresponding saturation temperature of air was taken. ater an extensive analytical works on the performance optimization analysis of wet fins was carried out by applying this linear relationship. A technique to determine the performance optimization of straight tapered longitudinal fins subject to simultaneous heat mass transfer has been established analytically by Kundu () Kundu Das (4). The performance optimum dimensions of a new fin, namely, SRC profile subject to simultaneous heat mass transfer have been investigated by Kundu (7a; 9a). In his work, a comparative study has also been made between rectangular SRC profile fins when they are operated in wet conditions. Hong Web (996) calculated the fin efficiency for wet dry circular fins with a constant thickness. Kundu Barman () have studied a design analysis of annular fins under dehumidifying conditions with a polynomial relationship between humidity ratio saturation temperature by using differential transform method. In case of longitudinal fins of rectangular geometry, approximate analytic solution for performances has been demonstrated by Kundu (9b). Kundu Miyara (9) have established an analytical model for determination of the performance of a fin assembly under dehumidifying conditions. Kundu et al. (8) have described analytically to predict the fin performance of longitudinal triangular fins subject to simultaneous heat mass transfer. The heat mass transfer analysis for dehumidification of air on fin--tube heat exchangers was done experimentally by the few authors. The different techniques, namely, new reduction method, tinny circular fin method, finite circular fin method review of data reduction method used for analyzing the heat mass transfer characteristics of wavy fin--tube exchangers under dehumidifying conditions had been investigated by Pirompugd et al. (7a; 7b; 8; 9). The above investigations had been focused on determination of the optimum profile subjected to convective environment. However a thorough research works have already been devoted for analyzing the performance optimization of wet fins. To carryout these analyses, suitable fin geometry has been chosen a priory. However, the optimum profile fin may be employed in air conditioning apparatus, especially, in aircrafts where reduction of weight is always given an extra design attention. Kundu (8) determined an optimum fin profile of thin fins under dehumidifying condition of practical interest formulated with the treatment by a calculus of variation. Recently, Kundu () focused to determine the optimum fin profile for both fully partially wet longitudinal fins with a nonlinear saturation curve. In this book chapter, a mathematical theory has been developed for obtaining the optimum fin shape of three common types of fins, namely, longitudinal, spine annular fins by satisfying the maximizing heat transfer duty for a given either fin volume or both fin volume length. The analysis was formulated for the dry, partially fully wet surface conditions. For the analytical solution of a wet fin equation, a relationship between humidity ratio temperature of the saturation air is necessary it is taken a linear variation. The influence of

6 Heat Transfer - Mathematical Modelling, Numerical Methods Information Technology wet fin surface conditions on the optimum profile shape its dimensions has also been examined. From the analysis, it can be mentioned that whether a surface is dry, partially or fully wet at an optimum condition, the air relative humidity is a responsible factor. The optimum fin profile design variables have been determined as a function of thermopsychrometric parameters. The dry surface analysis can be possible from the present fully wet surface fin analysis with considering zero value of latent heat of condensation. From the analysis presented, it can be highlighted that unlike dry partially wet fins, tip temperature for fully wet fins is below the ambient temperature for the minimum profile envelop fin.. Variational formulations for the optimum fin shape For determination of the optimum fin shape, it can be assumed that the condensate thermal resistance to heat flow is negligibly small as the condensate film is much thinner than the boundary layer in the dehumidification process. Under such circumstances, it may follow that the heat transfer coefficient is not influenced significantly with the presence of condensation. The condensation takes place when fin surface temperature is below the dew point of the surrounding air for its calculation, specific humidity of the saturated air on the wet surface is assumed to be a linear function with the local fin temperature. This assumption can be considered due to the smaller temperature range involved in the practical application between fin base dew point temperatures within this small range, saturation curve on the psychometric chart is possible to be an approximated by a straight line (Wu Bong, 994; Kundu, ; Kundu, 7a; Kundu, 7b; Kundu, 8; Kundu, 9). Owing to small temperature variation in the fin between fin-base fin-tip, it can be assumed that the thermal conductivity of the fin material is a constant. The different types of fins, namely, longitudinal, spine anuular fin are commonly used according to the shape of the primary surfaces. Depending upon the fin base, fin tip dew point temperatures, fin-surface can be dry, partially fully wet. The analysis for determination of an optimum profile of fully partially wet fins for longitudinal, spine annular fin geometries are described separately in the followings:. Fully wet longitudinal fins The schematic diagram of an optimum shape of longitudinal fins is illustrated in Fig.. The governing energy equation for one-dimensional temperature distribution on fully wet surface fins can be written under steady state condition as follows: d dt h y = ( T Ta) + hm( ω ωa) hfg h dx dx k h m is the average mass transfer coefficient based on the humidity ratio difference, ω is the humidity ratio of saturated air at temperature T, ω a is the humidity ratio of the atmospheric air, h jg is the latent heat of condensation. For the mathematical simplicity, the following dimensionless variables parameters can be introduced: X = hx k ; Y = hy k; hl k θ = ; 3 = ; ( T T) ( T T ) a a b m p () e = h h C () where, e is the ewis number. The relationship between heat mass transfer coefficients can be obtained from the Chilton-Colburn analogy (Chilton Colburn, 934). The relationship between the saturated water film temperature T the corresponding

Optimum Fin Profile under Dry Wet Surface Conditions 7 saturated humidity ratio ω is approximated by a linear function (Wu Bong, 994; Kundu, ; Kundu, 7a; Kundu, 7b; Kundu, 8; Kundu, 9) in this study: ω = a+ bt (3) where, a b are constants determined from the conditions of air at the fin base fin tip. Eq. () is written in dimensionless form by using Eqs. () (3) as follows: wet y b y b dry x l A B l x l Fig.. Schematic diagram of an optimum longitudinal fin under dehumidifying conditions: A. Fully wet; B. Partially wet. where ( ) ( ) ddxyd dx = + bξ (4) = + ; θ = ( ω a bt ) ( T T ) ( + bξ) θ θ p p a a a b Eq. (4) is subjected to the following boundary conditions: ; 3 ξ = hfg Cpe (5) at X =, = + θp = (6a) at X =, Yd dx = (6b) For determination of the heat transfer duty through fins, Eq. (4) is multiplied by, then integrated, the following expression are obtained with the help of the corresponding boundary conditions: (7) [ ] ( ) Yd dx = Y d dx + + bξ dx The heat transfer rate through the fins can be calculated by applying the Fourier s law of heat conduction at the fin base: q Q= = [ Yd dx] = Y( d dx) + ( + bξ) dx kt (8) ( T) a b

8 Heat Transfer - Mathematical Modelling, Numerical Methods Information Technology The fin volume is obtained from the following expression: V h k U = = YdX (9) The profile shape of a fin has been determined from the variational principle after satisfying the maximization of heat transfer rate Q for a design condition. In the present study, either the fin volume or both the fin volume length are considered as a constraint condition. A functional F may be constructed from Eqs. (8) (9) by employing agrange multiplier λ: ( ) () F = Q λu = Y d dx + + bξ λ Y dx The relation between the variation of F that of Y is obtained from the above equation for maximum value of F, δf is zero for any admissible variation of δy. Thus () δf = Y Y d dx λ Y δy dx = From the above equation, the following optimality criteria are obtained: ( ) λ Y d dx Y = () From Eq. (), it is obvious that the temperature gradient in the longitudinal fin for the optimum condition is a constant... Optimum longitudinal fin for the volume constraint Here the fin length is not a constant thus it can be taken as a variable. From Eq. (), the variation of function F with is as follows: δf = Y d dx + + bξ λ Y δx = (3) At X =, the above term vanishes as δ. At, δx is not zero; therefore, at the tip, the following optimality conditions can be obtained: Y d dx + + bξ λ Y = (4) Combining Eqs. (4), (6), () (4), yields the tip condition =. The tip thickness of a fin may be determined from the tip condition the optimality criterion boundary condition (Eqs. () (6b)). It can be seen that the tip thickness is zero. The tip temperature for fully wet surface θ t = -θ p, which is slightly less than the ambient value this temperature is obvious as a function of psychometric properties of the surrounding air. From Eqs. (4) (6), () tip condition, the temperature distribution fin profile are written as follows:

Optimum Fin Profile under Dry Wet Surface Conditions 9 ( θ ) θ = + p X (5) ( + bξ ) Y = X (6) The optimum fin length opt can be obtained from Eqs. (9) (6). The maximum heat transfer rate through the fin can be written by the design variables as follows: 3 { 6U ( + bξ )} { 3 ( ξ) 4} opt = 3 Q opt U + b (7a) (7b).. Optimum longitudinal fin for both length volume constraints In fin design, some times the length of the fin is required to specify due to restricted space ease of manufacturing. Under this design consideration, both length (fixed ) volume may be adopted as a constraint. For obtaining the temperature distribution fin profile, Eqs. (6), (9) () can be combined: θ = α X (8) ( + bξ ) α Y = X X α (9) where ( + bξ ) 3 ( ξ ) 3 α = 6U+ + b Here, it may be noted that the optimum fin shape for dry surface fins can be determined by using the above formula.. Partially wet longitudinal fins There are two regions dry wet in partially wet fins shown in Fig. B. For partially wet longitudinal fins, the energy equations are in the followings: () d dt y h ( T T ) for dry surface T T a d dx dx > k = d dt h y {( T Ta) hm( ω ωa) hfg h + } for wet surface T Td dx dx k (a) (b) By using Eqs. () (3), Eq. () is made in normalized form it can be expressed as follows:

Heat Transfer - Mathematical Modelling, Numerical Methods Information Technology d dθ Y for dry domain d dx dx θ θ > θ = d d Y ( + bξ ) for wet domain θ θd dx dx (a) (b) The heat transfer through the tip is negligibly small in comparision to that through the lateral surfaces fin base temperature is taken as a constant. In addition, continuity of temperature heat conduction satisfies at the section where dry wet separates. Thus, for solving Eq. () the following boundary conditions are taken: at X =, = (3a) at X =, θ = θd dθ dx = d dx (3b) (3c) at X =, Ydθ dx = (3d) Eq. () are multiplied by respective variables θ, the following relationships are obtained by integration using boundary conditions: [ ] [ ] ( ) (4a) Yd dx = Yd dx + Y d dx + + bξ dx [ ] Combining Eqs. (4a) (4b), one can get (4b) Yθdθ dx = Y dθ dx + θ dx ϕ = + + + + θ d [ Yϕdϕ dx] Y ( dθ dx) θ dx Y ( dϕ dx) ( bξ) ϕ dx (5) d The heat transfer rate through the fins is calculated by applying the Fourier s law of heat conduction at the fin base it can be written by using Eq. (5) as d d q Q = = Y = Y ( dθ dx) + θ dx kt T dx θ a b d + Y( d dx) + ( + bξ) dx The fin volume per unit width can be obtained from the following expressions: (6)

Optimum Fin Profile under Dry Wet Surface Conditions (7) V h k U = = YdX+ YdX The optimum profile shape of a fin can be determined from the variational principle by constructing a functional F from Eqs. (6) (7) using agrange multiplier λ. d F = Q λu = Y( d dx) + ( + bξ) λy dx d + Y( dθ dx) + θ λθdy d dx θ For maximizing value of F, the following condition is obtained from Eq. (8). δf = Y Y ( d dx) λy δy dx d + Y Y( dθ dx) λθdy d δydx = θ d From Eq. (9), the optimality criterion is derived as follows: (8) (9) ( ) Y dθ dx λ θ Y θ > θ d d d = Y d dx λ Y θ θd (3a) (3b).. Optimum longitudinal fin for volume constraint The variation of F with a function of yields the following expressions from Eq. (9): { } δ F Y d dx + + bξ λ Y δx θ θ d = = δ F ϕd { Y( dθ dx) + θ λy θd d} δx θ > θd θ d (3a) (3b) At X =, the above term should vanish as δ. At, δx is nonzero; thus, the location for both dry wet surfaces coexist the fin tip satisfies are the optimality conditions: ( ) Y dθ dx + θ λθ (3a) dy d = at X Y d dx bξ λ = Y + + (3b) Y dθ dx + θ λθ Y = at X = (33) d d

Heat Transfer - Mathematical Modelling, Numerical Methods Information Technology Combining Eqs. (4), (6b), (3), (3) (33), the tip temperature vanishes. Using optimality criteria, the temperature distribution fin profile can be expressed as ( θ ) X ( X) ( ) d X θ = θd X ( + bξ ) ( θ ) { ( θd ) } Y = ( ) + X X X for d (34a) (34b) (35a) Y = ( X) for X (35b) The length of the wet region can be determined by using an energy balance at that length where dry wet sections live together. = θ (36) d Here is not a constraint. can be obtained from Eqs. (7), (35) (36). The optimum length the maximum heat transfer rate through a fin can be written as Q opt opt ( 6U) 3 = θd θd ξ θd θd ( 3 ) + ( + b )( ) ( 3 + ) ( b )( )( ) ( 6U) θd + + ξ + θd θd = θd θd ξ θd θd ( 3 ) + ( + b )( ) ( 3 + ) 3 3 3 (37a) (37b).. Optimum longitudinal fin for both length volume constraints The temperature distribution fin profile can be determined by using Eqs. (4), (6) (3): θ α X X θ = θ α X X (38a) d (38b) ( θd α ) α ( ξ) { ( ) α( )} Y= + + + b X X X α (39a) where Y = d + X X X α ( θ α ) α (39b)

Optimum Fin Profile under Dry Wet Surface Conditions 3 3θd ( ) + 3( + b ξ) ( ξ ) α = 6U+ + b + 3 + 4 3 3 3 3 (4) = θ α (4) d.3 Fully wet annular fins Figure a is drawn for a schematic representation of an optimum annular fin under condensation of saturated vapor on its surfaces. The energy equation for one-dimensional temperature distribution in fully wet annular fins can be written under steady state condition as d dt y ri x h ( ri x) ( T Ta) hm( ω ωa) hfg h dx + = + + dx k Eq. (4) is made in dimensionless form by using Eqs. () (3) as (4) x Dry surface Wet surface l l r i y b O y A y b y Wet surface Dry surface O x B l l Fig.. Typical configuration of wet fins: A. Annular fin; B. Spine

4 Heat Transfer - Mathematical Modelling, Numerical Methods Information Technology d d Y( Ri X) ( Ri X)( bξ) ; Ri hri k dx + = + + = dx The boundary conditions for annular fins for temperature distribution are considered same as taken longitudinal fins described in Eq. (6). The actual nondimensional heat transfer rate is calculated from the following formula: ( π ) d i i ( ) ( ξ) a b dx q h Q= = Y R + X = R + X Y d dx + + b dx 4k T T (44) The fin volume in dimensionless form can be written as 3 π (43) V h k U = = ( Ri + X) YdX (45) 4 For the application of variational principle, a function F can be constructed by using agrange multiplier λin the followings : ( ) i (46) F = Q λu = R + X Y d dx + + bξ λ Y dx For maximizing value of F, Eq. (46) is differentiated with respect to Y finally it becomes to zero. δ F = Y R + X Y d dx λ Yδ Y dx = i (47) From Eq. (47), the following optimality condition for fully wet annular fins is obtained. ( ) λ Y d dx Y = (48) From the above optimality condition, it is to mention that the temperature gradient for fully wet anuular fins is a constant it does not depend upon the wetness condition of the fin..3. Optimum annular fin for volume constraint For volume constraint, fin length is variable. The variation of function F with length is determined from this expression. ( R + X) i δf = Y d dx + + bξ λ Y δx = At X =, the following optimality conditions are obtained: (49) Y d dx + + bξ λ Y = (5)

Optimum Fin Profile under Dry Wet Surface Conditions 5 Using Eqs. (6), (43), (48) (5), the condition for temperature at the tip is =. The thickness of fin profile at the tip is also zero. The temperature distribution fin profile for the optimum annular fin for fully wet surface conditions is obtained from Eqs. (6), (43), (48) optimum tip conditions as follows: ( + bξ ) ( θ ) θ = + p X (5) Y = X + X + R X 6 3 3 3 3 i R + X i The optimum length the maximum heat transfer rate through the fin can be determined as a function of design variables in the followings: ( ξ ) ( )( + ξ) 4 3 opt opt + Ri opt U + b = = Q 3R b 6 (5) (53) opt opt i opt (54) The optimum length can be possible to calculate by a numerical technique. The Newton- Raphson method can be applied for the solution of Eq. (53) to determine the optimum fin length..3. Optimum annular fin for both length volume constraints The temperature fin profile for an optimum annular fin under length volume constraints is determined from Eqs. (6), (43), (45) (48) they can be expressed as θ = α X (55) where ( + bξ ) ( R + X) 3 3 ( i ) Y = 6R i X + 3 α R X α X 6α i ( Ri + ) ( + b ) 3 ξ α = 3 U+ + b 4R + 3 ( ξ ) i (56) (57).4 Partially wet annular fins The energy equations for dry wet regions in partially wet fins can be written separately in the followings: d dt y( ri x) h + ( r ) dry domain i x T T T T a d dx dx > + k = d dt h y( ri x ) ( ri x) {( T Ta) hm( ω ωa) hfg h + + + } wet domain T Td dx dx k (58a) (58b)

6 Heat Transfer - Mathematical Modelling, Numerical Methods Information Technology Eq. (58) can be expressed in normalized form as d dθ Y( Ri + X) ( Ri X) d dx dx + θ θ > θ = d d Y( Ri + X) ( Ri + X)( + bξ ) θ θd dx dx (59a) (59b) For solving Eq. (59), four boundary conditions are required these are already written in Eq. (3). Heat transfer equation is obtained from Eq. (59) boundary conditions in the followings: ( π ) d q h d Q = = Y ( Ri + X) = ( Ri + X) Y ( dθ dx) + θ dx 4k T T dx θ a b d + ( Ri + X) Y( d dx) + ( + bξ) dx (6) The dimensionless fin volume for the partially wet fin can be formulated in the following: 3 V h k U = = R + X YdX+ R + X YdX 4 ( i ) ( i ) (6) π X = X = To apply variational principle, a functional F is constructed by using heat transfer rate fin volume expressions. d F = Q λu = ( Ri + X) Y( d dx) + ( + bξ) λy dx d + ( Ri + X) Y( dθ dx) + θ λθdy d dx θ The optimality criterion can be derived by differentiating functional F with respect to Y then equating to zero. (6) δ F = Y ( Ri + X) Y ( d dx) λy δy dx d + Y ( Ri + X) Y( dθ dx) λϕθdy d δydx = θ d (63) Therefore from Eq. (63), the following optimality criterion is found: ( d dx) ( d dx) θ λ θ θ > θ d d d = λ θ θ d (64a) (64b)

Optimum Fin Profile under Dry Wet Surface Conditions 7.4. Optimum annular fin for volume constraint The following expression is obtained for variation of F with length from Eq. (6) ( Ri + X) δf { Y ( d dx ) ( b ξ) λ Y } δ X + + d θ θ = = δf ( Ri + X) d { Y( dθ dx) + θ λy θd d} δx θ > θd θ d (65a) (65b) From Eq. (65), the following conditions are determined: ( ) Y dθ dx + θ λθ (66a) dy d = at X Y d dx bξ λ = Y + + (66b) Y dθ dx + θ λθdϕy d = at X = (67) From the above analysis, one can determine the dimensionless tip temperature tip thickness, both are zero for the optimum design condition. The fin temperature fin profile are obtained as follows: ( θ ) X ( X) ( ) d X θ = θd X (68a) (68b) Y = 6R + 3 R 3 3 ( i) i 6( Ri + X) ( + bξ ) 3 3 ( ) { 6Ri } ϕ X 3 ϕ θd Ri X θd X θd ( X ) + + for Y = 6 i 3 for 6 R X + R X X X ( R + X) i 3 3 ( i) (69a) (69b) The length of the optimum annular fin can be determined from the following equation given below: { } θ θ ξ 4ϕ 4θ ϕ 3 θ θ θ θ 3 3 + 6θd ( θd) ( bξ)( 3ϕ θd )( θd) θ + + + + d R i U = ( b ) 3 4 d d + + d d d + d d (7)

8 Heat Transfer - Mathematical Modelling, Numerical Methods Information Technology Eq. (7) can be solved numerically. The Newton-Raphson iterative method can be employed to determine the optimum length of the fin after satisfying the necessary convergence criterion. After estimating opt value, one can calculate the maximum or optimum actual heat transfer rate which can be determined from the expression give below: 3 { } { } { } opt Qopt = 6Riθd + 3( opt Ri ) ( θd ) opt θd 6 ( θ )( bξ) 6R 3( R )( θ ) ( θ ) + + + d opt i opt i d opt d (7).4. Optimum annular fin for both length volume constraints The temperature distribution fin profile for the annular fin can be determined by using Eqs. (6), (59), (6) (64) as θ α X X θ = θ α X X (7a) d (7b) Y = 6R i d + + 3 d Ri + 6α ( R + X) i 3 3 ( θ α ) ( θ α α ) α 3 3 ( bξ) { 6 } Ri X 3 αri X α X X + + + ( X ) 3 3 6( i θd α i ) 3 ( α α i θd) α ( R + X) Y = 6 R + R X + α R + X X where i 3 3 6Riθd( ) + ( + bξ)( 3Ri + ) + 4θd( ) ( ξ )( i ) i ( i ) α = U+ + b 4R + 3 R 3 + 4R + 3 4 3 3 3 3 (73a) (73b) (74a) = θ α (74b) d Combining Eqs. (74a) (74b), the following transcendental equation is obtained: 4 3 d d i d d + bξ 3 3θ 4 + 3θ + + R θ 3+ bξ 3 + bξ + 3θ ( θd)( i ) 3 4 + U+ 4R + 3 = (75) In order to determine the wet length in the fin, Eq. (75) can be solved by using Newton- Raphson iterative technique.

Optimum Fin Profile under Dry Wet Surface Conditions 9.5 Fully wet pin fins A schematic diagram of a pin fin is shown in Fig. B. The energy equation for pin fins subject to condensation of vapor under fully wet condition is written below: d dt h dx dx k y = y ( T Ta) + hm( ω ωa) hfg h (76) Eq. (76) is expressed in nondimensional form as d d Y = Y( + bξ ) dx dx (77) For the solution of Eq. (77), the boundary conditions expressed in Eq. (7) are taken. The actual heat transfer rate is calculated from the following expression: ( π ) d q h Q= = Y = Y Y ( d dx) + ( + bξ) dx (78) k T T dx a b The fin volume of a pin fin is written in dimensionless form as 3 V h k U = = Y dx (79) π For determination of the optimum shape, a functional F is defined as (8) F = Q λu = Y Y d dx + + bξ λ Y dx The above equation gives a relationship between F Y. The optimum condition can be obtained by differentiating functional F with respect to Y. ( ) (8) δf = Y d dx + + bξ λ Y δydx = The following expression is obtained from Eq. (8). Y d dx + + bξ λ Y = (8) Multiplying on both sides in Eq. (77) by d/dx yields the following expression after some manipulations: ( ) ( ) ( ξ ) ( ) d Y d dx dx + d dx d Y dx = + b Y d dx (83) Eliminating from Eq. (8) by using Eq. (83) then integrating, the following expression is obtained.

Heat Transfer - Mathematical Modelling, Numerical Methods Information Technology ( ) λ 3Y d dx Y = C (84) where, C is an integration constant determined by using boundary conditions..5. Optimum fully wet pin fin for volume constraint The variation of F with a function of yields from Eq. (8) as Y δf = Y d dx + + bξ λ Y δx = (85) As δx is nonzero at X =, the following tip condition can be achieved : Y d dx + + bξ λ Y = (86) From Eqs. (6), (77), (8) (86) the dimensionless temperature at the tip vanishes it can be indicated that the same condition is obtained in the case of longitudinal annular fins with fully wet condition. Using Eqs. (6), (77), (8) (86), temperature distribution fin profile for the optimum fin can be written as ( θ ) θ = + p X (87a) ( + bξ )( X) Y = (87b) The optimum fin length heat transfer rate are determined as follows: Q opt 5 { U ( + bξ ) } 3 4 { 5U ( bξ) 8} = 5 max + (88a) (88b).5. Optimum fully wet pin fin for both length volume constraints The temperature distribution fin profile can be found from Eqs. (6), (77) (84) as = Yb Yt Y + Y Yt + Yb, Yt Y Yb (89) Y ( Y Yt ) dy ( )( + ) ( bξ ) ( X) = + (9) Y Y Y Y Y Yt t t In order to determine the temperature profile fin profile from the above two equations, the unknown variables Y t Y b are required to calculate apriori. These can be determined by constructing two constraint equations as follows:

Optimum Fin Profile under Dry Wet Surface Conditions U = Yb ( Y Yt ) Y dy ( + ξ ) ( )( + ) (9a) = b Y Y Y Y Y Y Yt t t Yb ( Y Yt ) dy ( )( + ) ( bξ ) = + (9b) = Y Y Y Y Y Y Yt t t A simultaneous solution of Eqs. (9a) (9b) is provided to get the fin thickness at the base tip. The above all integration can be performed by Simson s /3 rule..6 Partially wet pin fins The energy equations for partially wet pin fins are written as d dt y h yt ( T) for dry domain T T a dx dx > k = d dt h { ( ω ω ) } dx dx k y y T Ta hm a hfg h + for wet domain T Td d (9a) (9b) Eq. (9) can be expressed in dimensionless form as ( ϕ Y ) ( b ) d Y dθ dx dx dry domain (93a) Yθ θ > θ = d d Y d dx dx + ξ wet domain θ θd (93b) Boundary conditions are reguired to solve Eq. (93) which can be taken as longitudinal fins expressed in Eq. (3). The actual heat transfer rate is calculated from the following equation: ( π ) q h d d ( θ ) θ a b dx X θ = d Q= = Y = Y Y d dx + dx k T T + Y Y( d dx) + ( + bξ) dx The fin volume for a partially wet pin fin is written in an integral form as (94) 3 (95) V h k U = = Y dx+ Y dx π For the application of variational method, a functional F is constructed from heat transfer rate fin volume expressions:

Heat Transfer - Mathematical Modelling, Numerical Methods Information Technology d F = Q λu = Y Y ( d dx) + ( + bξ) λy dx d + Y Y( dθ dx) + θ λθdy d dx θ (96) For maximizing value of F condition, the following expression is obtained: δf = Y( d dx) + ( + bξ) λyδydx d + Y ( dθ dx) + θ λθdy d δy dx = θ d (97) Therefore, from Eq. (97), the following optimality criterion is derived: ( ) Y dθ dx + θ λ θ Y θ > θ d d d = Y d dx b + + ξ λ Y θ θd (98a) (98b) The following expressions is obtained from Eq. (93) after multiplying on both sides of Eqs. (93a) (93b) by dθ dx d dx, respectively. { ( ) } + ( ) { ( ) } + ( ) ( ) ( + ξ ) ( ) d Y dθ dx dx dθ dx d Y dx Yd θ dx θ > θd = d Y d dx dx d dx d Y dx b Yd dx θ θ d (99a) (99b) Eleminating obtained: θ from Eq. (99) with the help of Eq. (98), the following equations are ( ) 3Y dθ dx λ θdy d C (a) θ > θd = 3Y d dx λ C Y θ θd (b) The integration constants C C are determined from the boundary conditions..6. Optimum partially wet pin fin for volume constraint only The variation of functional F with length parameter for partially wet pin fins is given below: { } Y δ F Y d dx + + b ξ λ Y δ X θ θ d = = δ F Yd { Y( dθ dx) + θ λyϕ θd d} δx θ > θd θ d (a) (b)

Optimum Fin Profile under Dry Wet Surface Conditions 3 For, the above terms vanish as δ. For, δx is not zero; therefore, the location where dry wet surfaces present, the tip, the following optimality conditions are achieved: Y dθ dx + θ λθ (a) dy d = at X Y d dx bξ λ = Y + + (b) Y dθ dx + θ λθ Y = at X = (c) d The temperature distribution fin profile for an optimum pin fin for volume constraint only are same with that of the longitudinal fin. The fin length heat transfer rate for the optimum profile fin is obtained from the following equations : d where Q opt opt ( 6U) = 4 3θd ( 5 θd) + Z{ Z + θd + }( 3 + θd ) 5 35 ( b )( )( ) ( 6U) θd + + ξ ϕ + θd θd = 4 3θd ( 5 θd) + Z{ Z + θd + }( 3 + θd ) 5 35 (3a) (3b) Z Z = + bξ θd (4a) { } ( θ )( bξ) 5( θ )( 5 6 θ ) 5 ( 4 3θ 3) 3( θ ) = + + + + + (4b) d d d d d.6. Optimum partially wet pin fin for both length volume constraints The temperature distribution fin profile can be found from Eqs. (6), (93), (95) () for the partially wet conditions as follows: Y b dy + d θd Y + θdy t =, Y Y Y Y dyb + d θd Y + θdyt b (5a) θ = θ Y Y Y + Y Y + Y, Y Y Y (5b) d t t t ( Y Y ) dy ( )( + ) Y = ( + bξ ) ( X) (6a) Y Y Y Y Y Y

4 Heat Transfer - Mathematical Modelling, Numerical Methods Information Technology Y ( Y Yt ) dy ( )( + ) ( X) = (6b) Y Y Y Y Y Yt t t From the above equations for obtaining temperature distribution fin profile, it is worthy to mention that the design variables Y, Y t, Y b are required to determine first. The Y t Y b can be determined from the constraint equations taken in the design as Yb ( ) ( ξ ) ( )( ) ( t ) Y Yb Y Y Y dy Y Y Y dy U = + b Y Y Y Y Y Y Y Y Y Y (7a) Y= Y + + Y= Yt t t + ( ) ( )( ) ( t ) Y Y Y dy Y Y dy + = ( + bξ ) + ( ) (7b) Y= Y + Y= Yt t t + Y Y Y Y Y Y Y Y Y Y The parameters Y are determined from the dewpoint temperature section where dry wet part coexist one can take the continuity of temperature heat conduction at this section. 3. Results discussion The humid air is a mixture of dry air water vapor. Three properties of air are used to calculate a thermodynamic state point to show the influence of dehumidification of air on the optimizing shape of three common type of fins, namely, ongitudinal, annular pin fins. The air properties such as pressure, temperature, relative humidity of air are used to determine a pychrometric state. The effect of the psychometric properties of air on the optimization study is investigated. Two constraints, namely, fin volume, both fin volume length have been adopted to furnish the result for any designed application. Fig. 3 is depicted the profile shape temperature distribution in an optimum fin under different surface conditions as a function of fin length for a constraint fin volume. For the dry surface fin, the optimum fin shape can be determined from the fully wet fin analysis presented above by taking zero value of latent heat parameter ξ. To make a comparison between dry wet surface fins, the results for the dry surface fin has been plotted in the same figure. From the figure, it can be noticeable that the temperature distribution in the fin at an optimum condition under dry, partially wet, fully wet surfaces varies linearly (Fig. 3A). For the same fin base temperature for all the surface conditions adopted here, temperature on the fin surface of wet fins differs from that of the dry surface fin the difference increases gradually from the fin base to fin tip. The discrepancy in temperature occurs due to evolving of latent heat of condensation of moisture on the fin surface in the case of wet fins. This difference in temperature becomes maximum for the % relative humidity of air as the maximum value of latent heat released for this relative humidity. From the analysis, it can be highlighted that for the dry partially wet fin at an optimum condition, the tip temperature equalizes to the ambient temperature. However, for fully wet fin, the tip temperature is sightly less than the ambient value. Figure 3B is illustrated the optimum fin profile for dry, partially wet fully wet conditions. For the dry fin, the optimum length

Optimum Fin Profile under Dry Wet Surface Conditions 5 optimum base thickness are larger smaller, respectively than that for the any wet fin for the same fin volume the difference is maximum for the fully wet fin of % relative humidity of air. Nevertheless, the effect of relative humidity on the optimum fin profile for a wet fin is marginal. The same trend for temperature distribution the optimum fin profile of the annular pin fins, are found they are displayed in Figs. 4 5, respectively. In the case of annular fins, it can be mentioned that the optimum temperature distribution fin profile are also function of the tube outer radius parameter R i. θ..8.6.4. Fully wet (γ = %) Dry surface Wet surface U=. T a =3 C T b =3 C Fully dry surface Partially wet (γ = 7%)....4.6.8 A X Y.5.4.3.. Fully wet (γ=%)....4.6.8 B Dry surface Wet surface U=. T a =3 C T b =3 C Partially wet (γ=7%) X Fully dry Fig. 3. Variation of temperature distribution fin profile of an optimum longitudinal fin for the different surface conditions under a volume constraint: A. fin temperature; B. fin profile With satisfying either maximizing heat transfer rate for a given fin volume or minimizing fin volume for a given heat transfer duty, the optimization of any fin is studied. Depending upon the requirement of design, any one of these two constraints is used but they give the same result. The optimum profile shape is determined from the solution of the optimality criteria of the fin design the constraint condition taken. The result from the optimization study of wet fins is depicted in Figs. 6, 7 8 as a function of fin volume for longitudinal, annular spine, respectively. The optimum result for the dry surface condition of each fin, in comparison, is plotted in the corresponding figure. From these illustrations, it is understable that the optimum parameters, namely, heat transfer rate, fin length, fin thickness at the base, enhance continually with the fin volume. The maximum rate of heat transfer is not only as a function of the fin material but also as a function of the condition of the surface. A fully wet surface with % relative humidity predicts a maximum optimum heat transfer rate per unit volume in comparison with that transferred by any other surfaces. A partially wet surface transfers a less amount of heat per unit volume in comparison to that by fully wet surface fin heat transfer rate decreases gradually with the decremented relative humidity. A dry surface fin at an optimum condition transfers least amount of heat per unit fin volume in comparison with the wet surface fin. This is happen due to latent heat evolved in wet fin heat transfer mechanism. For a lower value of fin volume, difference in heat transfer among fully wet, partially wet dry surfaces may not so much important

6 Heat Transfer - Mathematical Modelling, Numerical Methods Information Technology θ..8.6.4 Dry Surface Wet surface Fully dry fin Partially wet fin (γ=7%) Y.4.3. Fully dry fin Wet surface Dry surface U=. R i =.5 T a =3 C Tb=3 C A. Fully wet fin (γ=%)...55..65 X. Partially Fully wet wet fin (γ=7%) fin (γ=%)...55..65 Fig. 4. Variation of temperature distribution fin profile of an optimum annular fin for the different surface conditions for a volume constraint: A. fin temperature; B. fin profile. θ..8.6.4 Dry surface Wet surface Fully dry Fully wet. (γ=%) Partially wet (γ=7%).....3 A X B Y B.5.4.3. Fully wet (γ=%). Fully T b =3 C dry.....3 X Wet surface Dry surface Partially wet (γ=7%) X U=. T a =3 C Fig. 5. Variation of temperature distribution fin profile of an optimum spine for the different surface conditions under a volume constraint: A. fin temperature; B. fin profile than that for a higher value of fin volume as shown in Fig. 6A. Irrespective of any surface condition, the length fin thickness at the base of an optimum fin increases with the fin volume. The optimum fin length for fully wet fins is always shorter than that for the partially wet as well as dry fins. The optimum length is a maximum for the dry surface condition under the same fin volume (Fig. 6B). Nevertheless, an opposite trend is noticed for

Optimum Fin Profile under Dry Wet Surface Conditions 7 the variation of fin thickness at the base with the fin volume in comparison with the variation of fin length with volume for different surface condition as shown in Fig. 6C. A similar exercise has been made for the annular fin spine by plotting Figs. 7 8. In the case of annular fin, the above parameters also function of the thermogeometric parameter R i. With the increase in R i, the optimum heat transfer rate increases as well as the optimum length optimum base thickness decreases, separately with the same fin volume. From the above optimum results, it can be emphasized to highlight that the optimum fin shape obtained from an optimization technique with the consideration of only one constraint either fin volume or heat transfer rate, is a complex in nature fragile shape at the tip as already shown in Fig. 3 hence, it is difficult surely in manufacturing process. To overcome this problem to restrict the length of the fin, fin length can be taken as an additional constraint with the fin volume. In this case, the shape of the fin profile fragile Q opt.4.3. γ=7% γ=% opt.4.3. γ=7% γ=% ( Y b ) opt..8.6.4 γ=7% Fully wet Partially wet Fully dry γ=%. Fully wet Partially wet Fully dry. Fully dry Partially wet Fully wet. A...3.6.9 U...3.6.9 B U...3.6.9 C U Fig. 6. Design parameters of an optimum longitudinal fin as a function of fin volume: A. maximum heat transfer rate; B. optimum length; C. optimum semi-fin thickness at the base Q opt.8.7.6.5.4.3.. R i =. R i =.5 Fully wet (γ=%) Partially wet (γ=7%) Fully dry opt.3.5..5..5 Fully Dry R i =.5 R i =. Fully wet Partially (γ=%) wet (γ=%) T a =3 C T b =3 C ( Y b ) opt..5..5 Fully wet ( γ=%) R i =.5 R i =. Partially wet (γ=7%) Fully dry...3.6.9 A U B...3.6.9 U...3.6.9 C U Fig. 7. Design parameters of an optimum annular fin as a function of fin volume: A. maximum heat transfer rate; B. optimum length; C. optimum semi-fin thickness at the base

8 Heat Transfer - Mathematical Modelling, Numerical Methods Information Technology.6.45..5 Fully wet Partially wet Fully dry.4.35. Q opt.4.3.. γ=% γ=7% T a =3 C T b =3 C opt.3.5..5. γ=7% γ=% Fully dry Partially wet Fully wet ( Y b ) opt.8.6.4. γ=% γ=7% Fully wet Partially wet Fully dry...3.6.9 A U.5..3.6.9 B U...3.6.9 Fig. 8. Design parameters of an optimum pin fin as a function of fin volume: A. maximum heat transfer rate; B. optimum length; C. optimum fin thickness at the base geometry at the tip can be improved significantly. To avoid the same nature of the result, the optimum result under both volume length constraint is illustrated only for the longitudinal fin. The variation of temperature fin profile is determined under the aforementioned constraint, which is displayed in Fig. 9. From the temperature distribution, it can be mentioned that temperature at the tip of an optimum fin does not vanish depends upon the magnitude of constraints chosen. For a fully wet surface, temperature at the tip may be closer to the ambient value in comparison with that for the partially dry surface conditions. With the increase in relative humidity of air, condensation of moisture increases as a result fin surface temperature increases. This observation can be found in Fig. 9A. The profile shape under both volume length constraints for various surface conditions is illustrated as a function of dimensionless fin length shown in Fig. 9B. From this figure, it is clear that the profile shape is improved significantly with respect to a profile obtained from only volume constraint chosen, with the consideration of a suitable compatibility in the manufacturing technique. However, there is slightly different in shape between dry wet surface optimum fin profiles under the same design constants. C U 4. Conclusions The fin surface may be dry, fully wet or partially wet depending upon the psychometric conditions of the surrounding air participated as well as the constraints taken in the design. Owing to mass transfer occured with the heat transfer mechanism, the wet surface fin differs from that of dry surface fin. The optimum-envelop shape of wet fins is different with respect to that of dry surface fins. This deviation may be increased with the increase in relative humidity of air. In this chapter, the optimum profile shape of different fins, namely, longitudinal, spine annular are evaluated for dry, fully wet, partially wet surface conditions using variational principle. The analysis has also included for different constraint conditions, namely, fin volume both fin length fin volume presented for possible requirements of an optimum design. From the results, the optimum design variable for wet fins is not only function of the design constraints but also is a function of psychometric properties of air. Unlike dry partially wet surface fins, tip temperature for the fully wet

Optimum Fin Profile under Dry Wet Surface Conditions 9 optimum fins under the volume constraint is less than the surrounding temperature. A significant change in optimum design variables has been noticed with the design constants such as fin volume surface conditions. In order to reduce the complexcity of the optimum profile fins under different surface conditions, the constraint fin length can be selected suitably with the constraint fin volume. θ..9.8.7.6.5.4 γ =% U=. =.5 γ =7% Y.5.4.3. Fully wet (γ = %) Fully wet (γ = 7%) Fully dry U=. =.5.3. Fully dry surface Fully wet surface. A.....3.4.5 X B.....3.4.5 Fig. 9. Variation of temperature fin profile in a longitudinal fin as a function of length for both volume length constraints: A. Temperature distribution; B. Fin profile 5. Acknowledgement The authors would like to thank King Mongkut s University of Technology Thonburi (KMUTT), the Thail Research Fund, the Office of Higher Education Commission the National Research University Project for the financial support. 6. Nomenclatures a constant determined from the conditions of humid air at the fin base fin tip b slop of a saturation line in the psychometric chart, K C non-dimensional integration constant used in Eq. (84) - - specific heat of humid air, J kg K C p F functional defined in Eqs. (), (8), (46), (6), (8) (96) h convective heat transfer coefficient, W m - K h m h fg k l mass transfer coefficient, - - kg m S latent heat of condensation, J kg - thermal conductivity of the fin material, fin length, m W m - K X

3 Heat Transfer - Mathematical Modelling, Numerical Methods Information Technology l wet length in partially wet fins, m dimensionless fin length, hl k dimensionless wet length in partially wet fins, hl k e ewis number q heat transfer rate through a fin, W Q dimensionless heat transfer rate r i base radius for annular fins, m R i dimensionless base radius, hri k T temperature, K U dimensionless fin volume, see Eqs. (9), (7), (45), (6), (79), (9a) (95) V fin volume (volume per unit width for longitudinal fins), m 3 x, y coordinates, see Figs., m X, Y dimensionless coordinates, hx k hy k, respectively y semi-thickness of a fin at which dry wet parts separated, m Y dimensionless thickness, hy k Z, Z dimensionless parameters defined in Eqs. (4a) (4b), respectively Greek etters α parameter defined in Eqs. (), (4), (57) (74a) λ agrange multiplier ω specific humidity of air, kg w. v. / kg. d. a. ξ atent heat parameter dimensionless temperature, θ + θp dimensionless temperature at the fin base, + θp θ dimensionless fin temperature, ( T T) ( T T ) a a b θ p dimensionless temperature parameter, see Eq. (5) γ Relative humidity Subscripts a ambient b base d dewpoint opt optimum t tip 7. References Chilton, T.H. & Colburn, A.P. (934). Mass transfer (absorption) coefficients prediction from data on heat transfer fluid friction. Ind. Eng. Chem., Vol. 6, 83. Duffin, R. J. (959). A variational problem relating to cooling fins with heat generation. Q. Appl. Math., Vol., 9-9. Guceri, S. & Maday, C. J. (975). A least weight circular cooling fin. ASME J. Eng. Ind., Vol. 97, 9-93. Hanin,. & Campo, A. (3). A new minimum volume straight cooling fin taking into account the length of arc. Int. J. Heat Mass Transfer, Vol. 46, 545-55.

Optimum Fin Profile under Dry Wet Surface Conditions 3 Hong, K. T. & Webb, R.. (996). Calculation of fin efficiency for wet dry fins. HVAC&R Research, Vol., 7-4. Kundu, B. & Das, P.K. (998). Profiles for optimum thin fins of different geometry - A unified approach. J. Institution Engineers (India): Mechanical Engineering Division, Vol. 78, No. 4, 5-8. Kundu, B. (). Analytical study of the effect of dehumidification of air on the performance optimization of straight tapered fins. Int. Comm. Heat Mass Transfer, Vol. 9, 69-78. Kundu, B. & Das, P.K. (4). Performance optimization analysis of straight taper fins with simultaneous heat mass transfer. ASME J. Heat Transfer, Vol. 6, 86-868. Kundu, B. & Das, P. K. (5). Optimum profile of thin fins with volumetric heat generation: a unified approach. J. Heat Transfer, Vol. 7, 945-948. Kundu, B. (7a). Performance optimization analysis of SRC profile fins subject to simultaneous heat mass transfer. Int. J. Heat Mass Transfer, Vol. 5, 645-655. Kundu, B. (7b). Performance optimum design analysis of longitudinal pin fins with simultaneous heat mass transfer: Unified comparative investigations. Applied Thermal Engg., Vol. 7, Nos. 5-6, 976-987. Kundu, B. (8). Optimization of fins under wet conditions using variational principle. J. Thermophysics Heat Transfer, Vol., No. 4, 64-66. Kundu, B., Barman, D. & Debnath, S. (8). An analytical approach for predicting fin performance of triangular fins subject to simultaneous heat mass transfer, Int. J. Refrigeration, Vol. 3, No. 6, 3-. Kundu, B. (9a). Analysis of thermal performance optimization of concentric circular fins under dehumidifying conditions, Int. J. Heat Mass Transfer, Vol. 5, 646-659. Kundu, B. (9b). Approximate analytic solution for performances of wet fins with a polynomial relationship between humidity ratio temperature, Int. J. Thermal Sciences, Vol. 48, No., 8-8. Kundu, B. & Miyara, A. (9). An analytical method for determination of the performance of a fin assembly under dehumidifying conditions: A comparative study, Int. J. Refrigeration, Vol. 3, No., 369-38. Kundu, B. (). A new methodology for determination of an optimum fin shape under dehumidifying conditions. Int. J. Refrigeration, Vol. 33, No. 6, 5-7. Kundu, B. & Barman, D. (). Analytical study on design analysis of annular fins under dehumidifying conditions with a polynomial relationship between humidity ratio saturation temperature, Int. J. Heat Fluid Flow, Vol. 3, No. 4, 7-733. iu, C. Y. (96). A variational problem relating to cooling fins with heat generation. Q. Appl. Math., Vol. 9, 45-5. iu, C. Y. (96). A variational problem with application to cooling fins. J. Soc. Indust. Appl. Math., Vol., 9-9. Maday, C. J. (974). The minimum weight one-dimensional straight fin. ASME J. Eng. Ind., Vol. 96, 6-65. McQuiston, F. C. (975). Fin efficiency with combined heat mass transfer. ASHRAE Transaction, Vol. 7, 35-355. Mokheimer, E. M. A. (). Performance of annular fins with different profiles subject to variable heat transfer coefficient. Int. J. Heat Mass Transfer, Vol. 45, 363-364.

3 Heat Transfer - Mathematical Modelling, Numerical Methods Information Technology Pirompugd, W., Wang, C. C. & Wongwises, S. (7a). Heat mass transfer characteristics of finned tube heat exchangers with dehumidification. J. Thermophysics Heat transfer, Vol., No., 36-37. Pirompugd, W., Wang, C. C. & Wongwises, S. (7b). A fully wet fully dry tiny circular fin method for heat mass transfer characteristics for plain fin--tube heat exchangers under dehumidifying conditions. J. Heat Transfer, Vol. 9, No. 9, 56-67. Pirompugd, W., Wang, C. C. & Wongwises, S. (8). Finite circular fin method for wavy fin--tube heat exchangers under fully partially wet surface conditions. Int. J. Heat Mass Transfer, Vol. 5, 4-47. Pirompugd, W., Wang, C. C. & Wongwises, S. (9). A review on reduction method for heat mass transfer characteristics of fin--tube heat exchangers under dehumidifying conditions. Int. J. Heat Mass Transfer, Vol. 5, 37-378. Razelos, P. & Imre, K. (983). Minimum mass convective fins with variable heat transfer coefficient. J. Franklin Institute, Vol. 35, 69-8. Schmidt, E. (96). Warmeubertragung durch Rippen. Z. Deustsh Ing., Vol. 7, 885-95. Solov ev, B.A. (968). An optimum radiator-fin profile. Inzhenerno Fizicheskii Jhurnal, Vol. 4, No. 3, 488-49. Threlkeld, J.. (97). Thermal environment engineering. Prentice-Hall, New York. Wilkins, J. E. Jr. (96). Minimum mass thin fins with specified minimum thickness. J. Soc. Ind. Appl. Math., Vol. 9, 94-6. Wu, G. & Bong, T. Y. (994). Overall efficiency of a straight fin with combined heat mass transfer. ASRAE Transation, Vol., No., 365-374. Zubair, S. M.; Al-Garni, A. Z. & Nizami, J. S. (996). The optimum dimensions of circular fins with variable profile temperature-dependent thermal conductivity. Vol. 39, No. 6, 343-3439.

Heat Transfer - Mathematical Modelling, Numerical Methods Information Technology Edited by Prof. Aziz Belmiloudi ISBN 978-953-37-55- Hard cover, 64 pages Publisher InTech Published online 4, February, Published in print edition February, Over the past few decades there has been a prolific increase in research development in area of heat transfer, heat exchangers their associated technologies. This book is a collection of current research in the above mentioned areas describes modelling, numerical methods, simulation information technology with modern ideas methods to analyse enhance heat transfer for single multiphase systems. The topics considered include various basic concepts of heat transfer, the fundamental modes of heat transfer (namely conduction, convection radiation), thermophysical properties, computational methodologies, control, stabilization optimization problems, condensation, boiling freezing, with many real-world problems important modern applications. The book is divided in four sections : "Inverse, Stabilization Optimization Problems", "Numerical Methods Calculations", "Heat Transfer in Mini/Micro Systems", "Energy Transfer Solid Materials", each section discusses various issues, methods applications in accordance with the subjects. The combination of fundamental approach with many important practical applications of current interest will make this book of interest to researchers, scientists, engineers graduate students in many disciplines, who make use of mathematical modelling, inverse problems, implementation of recently developed numerical methods in this multidisciplinary field as well as to experimental theoretical researchers in the field of heat mass transfer. How to reference In order to correctly reference this scholarly work, feel free to copy paste the following: Balaram Kundu Somchai Wongwises (). Optimum Fin Profile under Dry Wet Surface Conditions, Heat Transfer - Mathematical Modelling, Numerical Methods Information Technology, Prof. Aziz Belmiloudi (Ed.), ISBN: 978-953-37-55-, InTech, Available from: http:///books/heattransfer-mathematical-modelling-numerical-methods--information-technology/optimum-fin-profile-underdry--wet-surface-conditions InTech Europe University Campus STeP Ri Slavka Krautzeka 83/A 5 Rijeka, Croatia Phone: +385 (5) 77 447 InTech China Unit 45, Office Block, Hotel Equatorial Shanghai No.65, Yan An Road (West), Shanghai, 4, China Phone: +86--64898 Fax: +86--64898