Transient Conduction: Spatial Effects and the Role of Analytical Solutions

Similar documents
Conduction Heat Transfer

Chapter 3, Solution 1C.

3-42. Chapter 15 Steady Heat Conduction. Heat Conduction in Cylinders and Spheres

Department of Civil Engineering & Applied Mechanics McGill University, Montreal, Quebec Canada

CIRCLE YOUR DIVISION: Div. 1 (9:30 am) Div. 2 (11:30 am) Div. 3 (2:30 pm) Prof. Ruan Prof. Naik Mr. Singh

Physic 231 Lecture 33

Analysis The characteristic length of the junction and the Biot number are

Lecture 12. Heat Exchangers. Heat Exchangers Chee 318 1

Analytical Modeling of Natural Convection in Horizontal Annuli

LECTURER: PM DR MAZLAN ABDUL WAHID PM Dr Mazlan Abdul Wahid

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

QUESTION ANSWER. . e. Fourier number:

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

Exploiting vector space properties for the global optimization of process networks

A/2 l,k. Problem 1 STRATEGY. KNOWN Resistance of a complete spherical shell: r rk. Inner and outer radii

_J _J J J J J J J J _. 7 particles in the blue state; 3 particles in the red state: 720 configurations _J J J _J J J J J J J J _

Chapter 6 : Gibbs Free Energy

SIMULATION OF THREE PHASE THREE LEG TRANSFORMER BEHAVIOR UNDER DIFFERENT VOLTAGE SAG TYPES

Wp/Lmin. Wn/Lmin 2.5V

IGEE 401 Power Electronic Systems. Solution to Midterm Examination Fall 2004

Chapter 7. Systems 7.1 INTRODUCTION 7.2 MATHEMATICAL MODELING OF LIQUID LEVEL SYSTEMS. Steady State Flow. A. Bazoune

Approach: (Equilibrium) TD analysis, i.e., conservation eqns., state equations Issues: how to deal with

Time-Dependent Conduction :

Comparison of Building Codes and Insulation in China and Iceland

A New Method for Solving Integer Linear. Programming Problems with Fuzzy Variables

EE 221 Practice Problems for the Final Exam

Module 7: Solved Problems

CHAPTER 3 ANALYSIS OF KY BOOST CONVERTER

Design of Analog Integrated Circuits

Faculty of Engineering

Natural Convection in a Horizontal Annulus with Oscillating Inner Cylinder Using Lagrangian-Eulerian Kinematics

THEORY OF HYPERBOLIC TWO-TEMPERATURE GENERALIZED THERMOELASTICITY

Section 3: Detailed Solutions of Word Problems Unit 1: Solving Word Problems by Modeling with Formulas

Nanoscale Heat Transfer using Phonon Boltzmann Transport Equation

Problem Set 5 Solutions - McQuarrie Problems 3.20 MIT Dr. Anton Van Der Ven

Mode-Frequency Analysis of Laminated Spherical Shell

Water vapour balance in a building moisture exposure for timber structures

Physics 107 HOMEWORK ASSIGNMENT #20

V. Electrostatics Lecture 27a: Diffuse charge at electrodes

A BESTEST VALIDATION STUDY OF THE DYNAMIC GROUND-COUPLED HEAT TRANSFER MODEL USED IN ACCURATE. Dong Chen 1. PO Box 56, Highett. Vic.

ME2142/ME2142E Feedback Control Systems. Modelling of Physical Systems The Transfer Function

SUPPLEMENTARY INFORMATION

4DVAR, according to the name, is a four-dimensional variational method.

Tubular Flow with Laminar Flow (CHE 512) M.P. Dudukovic Chemical Reaction Engineering Laboratory (CREL), Washington University, St.

SELECTION OF MODEL PARAMETERS OF BIOGAS IC ENGINE. Karol Cupiał, Grzegorz Katolik

element k Using FEM to Solve Truss Problems

Chapter 4. Unsteady State Conduction

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

ASSUMPTIONS: (1) One-dimensional, radial conduction, (2) Constant properties.

55:041 Electronic Circuits

Chapter 4: Transient Heat Conduction. Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University

6. ELUTRIATION OF PARTICLES FROM FLUIDIZED BEDS

Numerical Transient Heat Conduction Experiment

XXIX CILAMCE November 4 th to 7 th, 2008 Maceió - Brazil

( ) = ( ) + ( 0) ) ( )

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

Lucas Imperfect Information Model

Fall 2010 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede. (n.b. for now, we do not require that k. vectors as a k 1 matrix: ( )

Pull-Out Strength of a Cast-In-Place Anchor Bolt in Concrete Exposed to High Temperature

PT326 PROCESS TRAINER

Transient Effects on High Voltage Diode Stack under Reverse Bias

Finite Element Solution of MHD Transient Flow past an Impulsively Started Infinite Horizontal Porous Plate in a Rotating Fluid with Hall Current

Department of Applied Mathematics, Tsinghua University Beijing , People's Republic of China Received 17 August 1998; accepted 10 December 1998

A HYDRAULIC OPEN LOOP SYSTEM FOR CONTROLLED EXCAVATION ALONG PRESCRIBED PATH. E. Bundy, W. Gutkowski

Spring 2002 Lecture #17

Chem 204A, Fall 2004, Mid-term (II)

2 Analysis of the non-linear aerodynamic loads of hypersonic flow. 1 General Introduction

Process Engineering Thermodynamics E (4 sp) Exam

Thermodynamics of Materials

15-69C Under the conditions of complete combustion with stoichiometric amount of air.

2 Finite difference basics

Monin Obukhov Similarity and Local-Free-Convection Scaling in the Atmospheric Boundary Layer Using Matched Asymptotic Expansions

Chapter 5 Time-Dependent Conduction

Lecture 3: Shannon s Theorem

Int. J. of Applied Mechanics and Engineering, 2014, vol.19, No.3, pp DOI: /ijame

Introduction to Electronic circuits.

Radiation Chapter 12 L8 (MMV031) Martin Andersson

External Forced Convection :

CTN 2/23/16. EE 247B/ME 218: Introduction to MEMS Design Lecture 11m2: Mechanics of Materials. Copyright 2016 Regents of the University of California

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient

MODULE 7 HEAT EXCHANGERS

Appendix B. The Finite Difference Scheme

Feedback Principle :-

CHAPTER 2 Algebraic Expressions and Fundamental Operations

Solution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method

10/9/2003 PHY Lecture 11 1

Bi-level Optimization Method of Air-conditioning System Based on Office Building Energy Storage Characteristics

Reprint (R34) Accurate Transmission Measurements Of Translucent Materials. January 2008

The Effect Of Type-III Antifreeze Proteins (AFPs) On CO2 Hydrate Slurry Formation

Phys 344 Ch 5 Lect 4 Feb 28 th,

BME 5742 Biosystems Modeling and Control

6. Stochastic processes (2)

6. Stochastic processes (2)

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models

Short notes for Heat transfer

Outline. Unit Eight Calculations with Entropy. The Second Law. Second Law Notes. Uses of Entropy. Entropy is a Property.

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

Specific yield for a two-dimensional flow

FIND: (a) Sketch temperature distribution, T(x,t), (b) Sketch the heat flux at the outer surface, q L,t as a function of time.

MECH 375, Heat Transfer Handout #5: Unsteady Conduction

Transcription:

Transent Cnductn: Spatal Effects and the Rle f Analytcal Slutns

Slutn t the Heat Equatn fr a Plane Wall wth Symmetrcal Cnvectn Cndtns If the lumped capactance apprxmatn can nt be made, cnsderatn must be gven t spatal, as well as tempral, varatns n temperature durng the transent prcess. Fr a plane wall wth symmetrcal cnvectn cndtns and cnstant prpertes, the heat equatn and ntal/bundary cndtns are: 2T 1 T = x2 α t T x,0 T T x x= 0 (5.26) = (5.27) = 0 T k = h T( L, t) T x x= L Exstence f seven ndependent varables: T T x,, t T, T, k, α, h (5.28) (5.29) = (5.30) Hw may the functnal dependence be smplfed?

Nn-dmensnalzatn f Heat Equatn and Intal/Bundary Cndtns: Dmensnless temperature dfference: θ * x Dmensnless crdnate: x * L αt Dmensnless tme: t * F L 2 F the Furer Number θ = θ T T T T hl The Bt Number: B k Exact Slutn: θ * = f x *, F, B sld cs * = C nexp n 2 F nx * n= 1 θ ζ ζ C = 4snζ n ζ tanζ 2ζ + sn 2 = ( ζ ) n n n n n B (5.39a) (5.39b,c) See Appendx B.3 fr frst fur rts (egenvalues ζ,..., ζ 1 4 ) f Eq. (5.39c)

The One-Term Apprxmatn ( F > 0.2) : Varatn f mdplane temperature (x * = 0) wth tme F : θ ( T T ) ( T T ) C exp ( ζ F) * 1 1 2 Table 5.1 C and ζ as a functn f B 1 1 Varatn f temperature wth lcatn (x * ) and tme F : θ = θ cs ζ ( x ) * * 1 * Change n thermal energy strage wth tme: Δ E = Q st (5.41) (5.40b) (5.43a) snζ 1 Q = Q 1 * θ (5.46) ζ 1 Q = ρc T T (5.44) Can the fregng results be used fr a plane wall that s well nsulated n ne sde and cnvectvely heated r cled n the ther? Can the fregng results be used f an sthermal cndtn Ts T s nstantaneusly mpsed n bth surfaces f a plane wall r n ne surface f a wall whse ther surface s well nsulated?

Graphcal Representatn f the One-Term Apprxmatn The Hesler Charts Mdplane Temperature:

Temperature Dstrbutn: Change n Thermal Energy Strage:

Radal Systems Lng Rds r Spheres Heated r Cled by Cnvectn. B = hr / k F= αt/ r 2 One-Term Apprxmatns: Lng Rd: Eqs. (5.49) and (5.51) Sphere: Eqs. (5.50) and (5.52) C, ζ Table 5.1 1 1 Graphcal Representatns: Lng Rd: Fgs. D.4 D.6 Sphere: Fgs. D.7 D.9

The Sem-Infnte Sld A sld that s ntally f unfrm temperature T and s assumed t extend t nfnty frm a surface at whch thermal cndtns are altered. Specal Cases: Case 1: Change n Surface Temperature (T s ) ( 0, ) (,0) T t = T T x = T s T x, t Ts x = erf T T 2 αt (5.57) s q = s ( T ) k T s παt (5.58)

Case 2: Unfrm Heat Flux ( q = q ) 1 2 α π 2 2 q t/ x T( x, t) T = exp k 4α t q x x erfc k 2 αt (5.59) s Case 3: Cnvectn Heat Transfer, ht T k = h T T t x x= 0 (, ) ( 0, ) T x t T x = erfc T T 2 αt 2 hx h αt x h αt exp + erfc k k 2 + 2 αt k (5.60)

Multdmensnal Effects Slutns fr multdmensnal transent cnductn can ften be expressed as a prduct f related ne-dmensnal slutns fr a plane wall, P(x,t), an nfnte cylnder, C(r,t), and/r a sem-nfnte sld, S(x,t). See Equatns (5.64) t (5.66) and Fg. 5.11. Cnsder superpstn f slutns fr tw-dmensnal cnductn n a shrt cylnder: (,, ) T r x t T T T = (, ) (, ) P x t x C r t (, ) T x t T T r,t T = x T T T T Plane Wall Infnte Cylnder

Prblem 5.66: Chargng a thermal energy strage system cnsstng f a packed bed f Pyrex spheres. KNOWN: Dameter, densty, specfc heat and thermal cnductvty f Pyrex spheres n packed bed thermal energy strage system. Cnvectn ceffcent and nlet gas temperature. FIND: Tme requred fr sphere t acqure 90% f maxmum pssble thermal energy and the crrespndng center and surface temperatures. SCHEMATIC:

ASSUMPTIONS: (1) One-dmensnal radal cnductn n sphere, (2) Neglgble heat transfer t r frm a sphere by radatn r cnductn due t cntact wth adjnng spheres, (3) Cnstant prpertes. ANALYSIS: Wth B h(r /3)/k = 75 W/m 2 K (0.0125m)/1.4 W/m K = 0.67, the lumped capactance methd s napprprate and the apprxmate (ne-term) slutn fr ne-dmensnal transent cnductn n a sphere s used t btan the desred results. T btan the requred tme, the specfed chargng requrement ( Q/ Q = 0.9) must frst be used t btan the dmensnless center temperature, θ *. Frm Eq. (5.52), 3 ζ1 Q 1 ( ζ1) ζ1 cs ( ζ1) Q θ = 3sn Wth B hr /k = 2.01, ζ 1 2.03 and C 1 1.48 frm Table 5.1. Hence, 3 ( ) θ 0.1 2.03 0.837 = 0.155 3 0.896 2.03 0.443 = 5.386 =

Frm Eq. (5.50c), the crrespndng tme s 2 r θ t = ln 2 αζ 1 C 1 3 7 2 α = k / ρc = 1.4 W / m K / 2225 kg /m 835J/kg K = 7.54 10 m /s, 2 ( 0.0375m) ln( 0.155/1.48) 7 2 2 7.54 10 m /s ( 2.03) t = = 1,020s Frm the defntn f θ *, the center temperature s T = T + 0.155 T T = 300 C 42.7 C = 257.3 C g, g, The surface temperature at the tme f nterest may be btaned frm Eq. (5.50b) wth r = 1, ( ζ ) 1 θ sn 0.155 0.896 Ts = Tg, + ( T Tg,) = 300 C 275 C = 280.9 C ζ 2.03 1 Is use f the ne-term apprxmatn apprprate?

Prblem: 5.82: Use f radatn heat transfer frm hgh ntensty lamps ( q 4 2 s = 10 W/m ) fr a prescrbed duratn (t=30 mn) t assess ablty f frewall t meet safety standards crrespndng t maxmum allwable temperatures at the heated (frnt) and unheated (back) surfaces. KNOWN: Thckness, ntal temperature and thermphyscal prpertes f cncrete frewall. Incdent radant flux and duratn f radant heatng. Maxmum allwable surface temperatures at the end f heatng. FIND: If maxmum allwable temperatures are exceeded. SCHEMATIC:

ASSUMPTIONS: (1) One-dmensnal cnductn n wall, (2) Valdty f semnfnte medum apprxmatn, (3) Neglgble cnvectn and radatve exchange wth the surrundngs at the rradated surface, (4) Neglgble heat transfer frm the back surface, (5) Cnstant prpertes. ANALYSIS: The thermal respnse f the wall s descrbed by Eq. (5.59) ( α π) 1/2 2 x 2q t/ q x x T( x, t) = T + exp erfc k 4α t k 2 α t 7 2 where, α = k/ ρc = 6.92 10 m /s and fr p ( α π) 1/2 t = 30 mn = 1800s, 2q t / / k = 284.5 K. Hence, at x = 0, T ( 0,30 mn) = 25 C + 284.5 C = 309.5 C < 325 C 2 At 1/2 x = 0.25m, x / 4αt = 12.54, q x / k = 1, 786K, and x / 2 αt = 3.54. Hence, 6 T 0.25m, 30 mn = 25 C + 284.5 C 3.58 10 1786 C ~ 0 25 C

Bth requrements are met. Is the assumptn f a sem-nfnte sld fr a plane wall f fnte thckness apprprate under the fregng cndtns? COMMENTS: The fregng analyss may r may nt be cnservatve, snce heat transfer at the rradated surface due t cnvectn and net radatn exchange wth the envrnment has been neglected. If the emssvty f the surface and the temperature f the surrundngs are assumed t be ε = 1 and T sur = 298K, radatn exchange at T s = 309.5 C wuld be 4 4 2 q = εσ T T = 6,080 W/m K, rad s sur whch s sgnfcant (~ 60% f the prescrbed radatn). Hwever, under actual cndtns, the wall wuld lkely be expsed t cmbustn gases and adjnng walls at elevated temperatures.