Natural Convection Experiment Measurements from a Vertical Surface

Similar documents
Thermal interaction between free convection and forced convection along a vertical conducting wall

Model Prediction of Heat Losses from Sirosmelt Pilot Plant

3B SCIENTIFIC PHYSICS

THERMODYNAMICS Lecture 15: Heat exchangers

Review for Exam #2. Specific Heat, Thermal Conductivity, and Thermal Diffusivity. Conduction

Maximizing Energy in Terahertz Pulse Radiation from a Switched Oscillator

Effects of Baffle on Entropy Generation in Separated Convection Flow Adjacent to Inclined Backward-Facing Step

Homework Set 4. gas B open end

Research on Static Tension Ratio Characteristic of Double-Vessel Friction Hoist System Components

STRESS ANALYSIS OF RUBBER BLOCKS UNDER VERTICAL LOADING AND SHEAR LOADING. A Dissertation. Presented to

Heat exchangers: Heat exchanger types:

2. Mass transfer takes place in the two contacting phases as in extraction and absorption.

FEM ANALYSES OF CUTTING OF ANISOTROPIC DENSELY COMPACTED AND SATURATED SAND

Role of Thermal Conductivity for Thermoelectrics with Finite Contacts

Experimental Investigation on the Effect of Fluid Flow Rate on the Performance of a Parallel Flow Heat Exchanger

Part G-4: Sample Exams

Chapter 3. Problem Solutions

A Study of Dufour and Soret Effect on MHD Mixed Convection Stagnation Point Flow towards a Vertical Plate in a Porous Medium

Vibration Control Using Heat Actuators

The Compton effect according to Schrödinger s theory

Millennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion

Determination of heat transfer intensity between free streaming water film and rigid surface using thermography

23.1 Tuning controllers, in the large view Quoting from Section 16.7:

International Journal of Advance Engineering and Research Development PERFORMANCE EVALUATION OF COMPOUND MULTILAYER INSULATION (77K-300K)

EFFECT OF PITCH NUMBER IN OVERALL HEAT TRANSFER RATE IN DOUBLE PIPE HELICAL HEAT EXCHANGER

The Second Law of Thermodynamics

Chapter 3 Lecture 7. Drag polar 2. Topics. Chapter-3

Process engineers are often faced with the task of

Four-dimensional equation of motion for viscous compressible substance with regard to the acceleration field, pressure field and dissipation field

Analytical Solution for Bending Stress Intensity Factor from Reissner s Plate Theory

Wave Propagation through Random Media

Lecture 27: Entropy and Information Prof. WAN, Xin

IMPACT MODELLING OF THE COEFFICIENT OF RESTITUTION OF POTATOES BASED ON THE KELVIN- VOIGHT PAIR

Flow over a hill covered with a plant canopy

NUMERICAL STUDY AND THERMAL PERFORMANCE OF THE FLAT PLATE SOLAR AIR HEATERS WITH AND WITHOUT THERMAL STORAGE

NUMERICAL SIMULATION OF ATOMIZATION WITH ADAPTIVE JET REFINEMENT

Lecture 27: Entropy and Information Prof. WAN, Xin

Measuring & Inducing Neural Activity Using Extracellular Fields I: Inverse systems approach

EF 152 Exam #3, Spring 2016 Page 1 of 6

Cavity flow with surface tension past a flat plate

COMPUTER METHODS FOR THE DETERMINATION OF THE CRITICAL PARAMETERS OF POLLUTED INSULATORS

Maximum work for Carnot-like heat engines with infinite heat source

Prediction of Coating Thickness

Why do moisture convergence deep convection schemes work for more scales than those they were in principle designed for?

Earlier Lecture. This gas tube is called as Pulse Tube and this phenomenon is called as Pulse Tube action.

Development of the Numerical Schemes and Iteration Procedures Nielsen, Peter Vilhelm

Observations on harmonic Progressions *

Optimal control of solar energy systems

Bottom Shear Stress Formulations to Compute Sediment Fluxes in Accelerated Skewed Waves

Design of a Wind Tunnel Apparatus to assist Flow and Aeroelastic Control via Zero Net Mass Flow Actuators

Where as discussed previously we interpret solutions to this partial differential equation in the weak sense: b

COMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION

Physical Laws, Absolutes, Relative Absolutes and Relativistic Time Phenomena

Main Menu. SEG Houston 2009 International Exposition and Annual Meeting

Physics 41 Chapter 22 HW

LATTICE BOLTZMANN METHOD FOR MICRO CHANNEL AND MICRO ORIFICE FLOWS TAIHO YEOM. Bachelor of Science in Mechanical Engineering.

Physics 231 Lecture 35

DEVELOPMENT OF A MULTI-FEED P-T WELLBORE MODEL FOR GEOTHERMAL WELLS

Study the Effect of Variable Viscosity and Thermal Conductivity of Micropolar Fluid in a Porous Channel

EFFECT OF SURFACE FINISH ON HEAT TRANSFER PERFORMANCE OF PLATE HEAT EXCHANGER

UNIT 1 OPEN CHANNEL FLOW 2 MARK QUESTIONS AND ANSWERS

arxiv: v1 [hep-ph] 5 Sep 2016

THE ESSENCE OF QUANTUM MECHANICS

f 2 f n where m is the total mass of the object. Expression (6a) is plotted in Figure 8 for several values of damping ( ).

arxiv:gr-qc/ v2 24 Jul 2002

3 Tidal systems modelling: ASMITA model

UNIVERSAL RELATIONSHIP BETWEEN COLLECTION EFFICIENCY AND THE CORONA POWER OF THE ELECTROSTATIC PRECIPITATOR

+Ze. n = N/V = 6.02 x x (Z Z c ) m /A, (1.1) Avogadro s number

Improvements in the Modeling of the Self-ignition of Tetrafluoroethylene

Wave-Particle Duality: de Broglie Waves and Uncertainty

EXPERIMENTAL STUDY ON BOTTOM BOUNDARY LAYER BENEATH SOLITARY WAVE

Supporting information

A = h w (1) Error Analysis Physics 141

NUMERICAL SIMULATION OF SEA ICE GROWTH AND DESALINATION

A NORMALIZED EQUATION OF AXIALLY LOADED PILES IN ELASTO-PLASTIC SOIL

Lecture 3 Heat Exchangers

An Experimental Study of Heat Transfer Enhancement in the Perforated Rectangular Fin

Stability of a Cold Core Eddy in the Presence of Convection: Hydrostatic versus Nonhydrostatic Modeling

CALCULATION OF THE HEAT TRANSFER AND TEMPERATURE ON THE AIRCRAFT ANTI-ICING SURFACE

EF 152 Exam #3, Fall, 2012 Page 1 of 6

ESCI 341 Atmospheric Thermodynamics Lesson 11 The Second Law of Thermodynamics

10.2 The Occurrence of Critical Flow; Controls

DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS

Effects of the swirl ratio on the turbulent flow fields of tornadolike vortices by using LES turbulent model

Developing Excel Macros for Solving Heat Diffusion Problems

ACE Engineering Academy. Hyderabad Delhi Bhopal Pune Bhubaneswar Lucknow Patna Bengaluru Chennai Vijayawada Vizag Tirupati Kukatpally Kolkata

Analysis of Leakage Paths Induced by Longitudinal Differential Settlement of the Shield-driven Tunneling

Analysis of discretization in the direct simulation Monte Carlo

The Verlet Algorithm for Molecular Dynamics Simulations

Simulation of hybrid Photovoltaic-Thermal Collector (PV-TC) Systems for domestic Heating and Cooling Case Study: Island of Rhodes

CONVECTION AT A.MODEL ICE EDGE

THERMAL MODELING OF PACKAGES FOR NORMAL CONDITIONS OF TRANSPORT WITH INSOLATION t

NUMERICAL SIMULATION OF SEA ICE GROWTH AND DESALINATION

Lightning electromagnetic environment in the presence of a tall grounded strike object

Chapter 2 Linear Elastic Fracture Mechanics

Green s function for the wave equation

Adaptive neuro-fuzzy inference system-based controllers for smart material actuator modelling

Duct Acoustics. Chap.4 Duct Acoustics. Plane wave

Physics 107 Problem 2.5 O. A. Pringle h Physics 107 Problem 2.6 O. A. Pringle

CRITICAL EXPONENTS TAKING INTO ACCOUNT DYNAMIC SCALING FOR ADSORPTION ON SMALL-SIZE ONE-DIMENSIONAL CLUSTERS

Transcription:

OBJECTIVE Natural Convetion Experiment Measurements from a Vertial Surfae 1. To demonstrate te basi priniples of natural onvetion eat transfer inluding determination of te onvetive eat transfer oeffiient. 2. To demonstrate te boundary layer arater of external natural onvetion. BACKGROUND Sine te density of most fluids vary wit temperature, temperature gradients witin a fluid medium will give rise to density gradients. If tese density gradients are su tat te fluid is in an unstable situation, eavy fluid on top of ligt fluid, te fluid will begin to move. Tis motion is termed natural onvetion. Newton's law of ooling governs tis pysial proess wi states tat te eat transfer from te surfae is diretly proportional to te temperature differene between te surfae and te fluid far away from te surfae or ( T - T ) q& (1) surfae fluid Introduing te onvetive eat transfer oeffiient as te onstant of proportionality, we ave ( T - T ) q& = (2) surfae fluid It is lear tat te major obstale in utilizing Eq. (2) for onvetive eat transfer alulations is te evaluation of te onvetive eat transfer oeffiient. Tere are tree standard metods used to evaluate. Te first involves a matematial solution to te onservation equations in differential form. For problems were tese equations are too ompliated to be solved analytially, we an employ te seond metod, a omputational solution. Finally for problems tat are so ompliated tat we annot even write te appropriate onservation equations, we must go into te laboratory and make measurements in employing an experimental solution. Before we ontemplate employing one of tese solution metods, it is useful to use our intuition to figure out upon wat te onvetive eat transfer oeffiient depends. Our intuition tells us tat te tree major fators assoiated wit te alulation of sould be (i) fluid meanis 1

(ii) (iii) fluid transport properties geometry Starting wit te fluid meanis, we reognize tere are a variety of ways to araterize te flow. We an onsider wat is driving te flow and lassify it as fored or natural onvetion. Next, we onsider ow many boundaries te fluid flow interats wit and lassify it as external or internal flow. Reall tat te differene between internal flows and external flows is often one of perspetive. A tird way to araterize te flow is by te presene or absene of turbulene. An important onsideration in te andling of onvetive eat transfer oeffiients is te notion of dynami similarity. It is found tat ertain systems in fluid meanis or eat transfer are found to ave similar beaviors even toug te pysial situations may be quite different. Reall te fluid meanis of flow in a pipe. Wat we are able to do is to take data as sown in Fig. 1 for different fluids and pipe diameters and by appropriately saling ollapse tese urves into one urve. Figure 1. Dynami Similarity for Pipe Flow Pressure Drop Frition Fator Veloity Reynolds Number In onvetive eat transfer we may apply dynamis saling to make a parallel transformation. 2

Figure 2. Dynami Similarity for Convetive Heat Transfer Convetive Heat Transfer Coeffiient Nusslet Number Temperature Differene Rayleig Number We ave defined a dimensionless onvetive eat transfer oeffiient alled te Nusselt number as Nu = L k (3) were : onvetive eat transfer oeffiient L: arateristi lengt k: termal ondutivity of te fluid. Te arateristi lengt is osen as te system lengt tat most affets te fluid flow. For flow along a flat plate our arateristi lengt is te plate lengt, L, and we write L Nu = D (4) k Te Rayleig number is indiative of te buoyany fore tat is driving te flow and is given by gβts T Ra = να f L 3 (5) 3

were g: aeleration due to gravity β: fluid termal expansion oeffiient T s : surfae temperature T f : fluid temperature L: arateristi lengt ν: fluid kinemati visosity α: fluid termal diffusivity. Anoter important feature introdued by te fluid meanis is te loal nature of te onvetive eat transfer oeffiient. If as te fluid flows over different regions of te surfae, te fluid meanis ange, ten te onvetive eat transfer oeffiient will ange. For natural onvetion over a vertial flat plate we will ave te boundary layer flow sown in Fig. 3. Figure 3. Flat Plate Boundary Layer Development δ x r g As te boundary layer tikness grows, te fluid meanis ange signifiantly, so tat te onvetive eat transfer oeffiient will vary along te lengt of plate and we will ave a loal onvetive eat transfer oeffiient, (x). Ten we may also define loal Nusselt and Rayleig numbers as Nu x = (x) x k (6) 4

gβts (x) T Ra = να f x 3 (7) Toug loal eat transfer onditions an be extremely important, an average eat transfer oeffiient over te entire surfae lengt is often desirable. By definition we ave 1 = L l, avg 0 (x) dx (8) wit an average Nusselt Number given as Nu L,avg =,avg k L (9) Te dimensionless parameter wi is used to represent te affet of fluid properties is te Prandtl number ν Pr = (10) α Te influene of geometry may be seen in a ouple of ways. First, for tose onfigurations tat ave two lengt dimensions, su as a ylinder, we introdue a dimensionless geometri parameter D Χ = (11) L Te seond way in wi we see geometrial influenes is troug te funtional form of te Nusselt number orrelation. In general we may write Nu = fn(ra,pr,x) (12) For simple situations tese may often be written as power law relationsips n Pr m Nu = a Ra (13) were te onstants a, m, and n will ange for different geometries. 5

In tis experiment te student will develop te relationsip among te Nusselt number and oter dimensionless parameters for natural onvetion from a vertial flat surfae. Te students will determine loal eat transfer oeffiients wi are indiative of a boundary layer penomena. Tese loal measurements will ten be averaged and ompared to publised orrelations. Natural onvetion eat transfer from a vertial surfae will be investigated using an eletrially eated flat plate. A semati of te apparatus is sown in Figure 4. Figure 4. Semati of Experimental Apparatus Rotary Seletion Swit Power Supply Flat Plate Voltage DMM Current DMM Test Setion T/C DMM Assuming te eating is uniform, te eat flux at any loation on te plate is given by, V I q& = (14) 2 w L were V and I are te eletri voltage and urrent supplied to te plate and w and L are te widt and lengt of te eated surfae. Temperature measurements are made along te lengt of te plate wit termoouple embedded beneat te eated surfae. 6

Te loal eat transfer oeffiient (te eat transfer oeffiient at a ertain distane from te leading edge) is ten alulated as q& (x) = T (x) - T s f (15) Sine in most engineering appliations te interest would be in an average eat transfer oeffiient for te entire surfae, we use te matematial definition 1 L = (x) dx (16) L 0 Now substituting and reognizing tat te termoouple is designed to read te temperature differene diretly, we ave q& L dx = L T(x) 0 (17) We now apply a trapezoidal rule integration to obtain q& N-1 1 1 = 0.5 + L i=1 Ti Ti-1 ( x - x ) i i-1 (18) were N is te number of temperature data points along te plate and i=0 orresponds to te point at te leading edge. Te average Nusslet number an ten alulated. PROCEDURE 1. Cek tat te power supply is unplugged and tat te transformer is set at zero. 2. Plug in te power supply and turn on te transformer. 5. Set te transformer at a power level of about 5% and allow te system to stabilize. You will want to ave measurements for at least five power settings. Sine natural onvetion is a relatively poor eat transfer proess you will want to monitor te maximum surfae temperature to prevent burnout of te system. It is reommended tat you do not exeed 25% power level on te transformer. 6. Read and reord te voltage measurements assoiated wit te eletrial voltage and urrent supplied for eating and te termoouples. 7

7. Inrease te power level, allow te system to stabilize, and reord appropriate measurements. Repeat te experiment as needed 8. Turn off te transformer and unplug te power supply. DATA ANALYSIS 1. Ea ase sould be reorded on te Exel spread seets provided in te lab. 2. From te experimentally determined values of te loal eat transfer oeffiient alulate te average eat transfer oeffiient and surfae temperature for ea power setting. 3. Plot te loal eat transfer oeffiient versus distane from te leading edge. 4. Calulate te average Nusselt number and Rayleig number for ea power setting. 5. Plot te Nusselt number versus Rayleig number for ea ase along wit te publised orrelations. 6. Estimate te unertainty error in your experimentally determined Nusselt number and Rayleig number. 7. Provide one sample and alulation of your data proessing. SUGGESTIONS FOR DISCUSSION 1. How does te experimental data ompare to te publised orrelations? 2. Wat are some possible errors in te experiment? 3. Wat an you tell about te struture of te boundary layer from te loal eat transfer oeffiient? 4. Can you identify a transition from laminar to turbulent flow? 8