Convection and conduction and lumped models

Similar documents
a. (1) Assume T = 20 ºC = 293 K. Apply Equation 2.22 to find the resistivity of the brass in the disk with

Energy Storage Devices

Dishonest casino as an HMM

R th is the Thevenin equivalent at the capacitor terminals.

Problem 1 Known: Dimensions and materials of the composition wall, 10 studs each with 2.5m high

The Components of Vector B. The Components of Vector B. Vector Components. Component Method of Vector Addition. Vector Components

Cooling of a hot metal forging. , dt dt

Lecture 4 ( ) Some points of vertical motion: Here we assumed t 0 =0 and the y axis to be vertical.

Kinematics Review Outline

Different kind of oscillation

2015 Sectional Physics Exam Solution Set

10.7 Temperature-dependent Viscoelastic Materials

D Solute mass diffusivity ( m ) D Mass diffusivity ( m s ) Gc Pr. Sc Sr Re g. q Radiative heat flux. U Reference velocity ( m s )

Physic 231 Lecture 33

Solution in semi infinite diffusion couples (error function analysis)

Notes on Inductance and Circuit Transients Joe Wolfe, Physics UNSW. Circuits with R and C. τ = RC = time constant

(V 1. (T i. )- FrC p. ))= 0 = FrC p (T 1. (T 1s. )+ UA(T os. (T is

Module 4. Analysis of Statically Indeterminate Structures by the Direct Stiffness Method. Version 2 CE IIT, Kharagpur

A Demand System for Input Factors when there are Technological Changes in Production

Online Supplement for Dynamic Multi-Technology. Production-Inventory Problem with Emissions Trading

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!

Energy & Work

SPH3U1 Lesson 06 Kinematics

2010 Sectional Physics Solution Set

Square law expression is non linear between I D and V GS. Need to operate in appropriate region for linear behaviour. W L

Subject: Turbojet engines (continued); Design parameters; Effect of mass flow on thrust.

Graduate Macroeconomics 2 Problem set 5. - Solutions

( )a = "t = 1 E =" B E = 5016 V. E = BHv # 3. 2 %r. c.) direction of induced current in the loop for : i.) "t < 1

Physics 20 Lesson 9H Rotational Kinematics

PRINCE SULTAN UNIVERSITY Department of Mathematical Sciences Final Examination Second Semester (072) STAT 271.

5.1 Angles and Their Measure

by Lauren DeDieu Advisor: George Chen

Wp/Lmin. Wn/Lmin 2.5V

ELEG 205 Fall Lecture #10. Mark Mirotznik, Ph.D. Professor The University of Delaware Tel: (302)

(,,, ) (,,, ). In addition, there are three other consumers, -2, -1, and 0. Consumer -2 has the utility function

2/20/2013. EE 101 Midterm 2 Review

Differentiation Applications 1: Related Rates

Chapter 6 : Gibbs Free Energy

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

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

, where. This is a highpass filter. The frequency response is the same as that for P.P.14.1 RC. Thus, the sketches of H and φ are shown below.

Digital Integrated CircuitDesign

IrrItrol Products 2016 catalog

Lecture 12. Heat Exchangers. Heat Exchangers Chee 318 1

convection coefficient. The different property values of water at 20 C are given by: u W/m K, h=8062 W/m K

( ) [ ] MAP Decision Rule

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4

Physics Courseware Physics I Constant Acceleration

Unit-I (Feedback amplifiers) Features of feedback amplifiers. Presentation by: S.Karthie, Lecturer/ECE SSN College of Engineering

Chapter 7 AC Power and Three-Phase Circuits

Motion of Wavepackets in Non-Hermitian. Quantum Mechanics

Conduction Heat Transfer

DESIGN OF TENSION MEMBERS

AP Physics Kinematic Wrap Up

RAMIFICATIONS of POSITION SERVO LOOP COMPENSATION

" 1 = # $H vap. Chapter 3 Problems

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

β A Constant-G m Biasing

H = d d q 1 d d q N d d p 1 d d p N exp

1.2.1 Vectors. 1 P age. Examples What is the reference vector angle for a vector that points 50 degrees east of south?

Chapter 9 Compressible Flow 667

Physics 107 HOMEWORK ASSIGNMENT #20

Copyright 2012 by the American Society for Pharmacology and Experimental Therapeutics.

Graphene nanoplatelets induced heterogeneous bimodal structural magnesium matrix composites with enhanced mechanical properties

ε (5) 4. Radiative Heat Transfer 4.1 Fundamentals of thermal radiation (1) W/m / K Stefan-Boltzmann-constant

Introduction ( Week 1-2) Course introduction A brief introduction to molecular biology A brief introduction to sequence comparison Part I: Algorithms

Chapter 3, Solution 1C.

ECE-320: Linear Control Systems Homework 1. 1) For the following transfer functions, determine both the impulse response and the unit step response.

AE/ME 339. K. M. Isaac. 8/31/2004 topic4: Implicit method, Stability, ADI method. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.

Ordinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s

Chapter 5: 8, 16, 25, 27, 28, 35, 36, 49, 52, 57, 72, 82 and 87.

PEGN 513 Reservoir Simulation I Fall 2009

Module 7: Solved Problems

Steady State Conduction

DEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL

Chapter 6. Dielectrics and Capacitance

Physics 201 Lecture 15

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD

7-84. Chapter 7 External Forced Convection

Problem 1. Refracting Surface (Modified from Pedrotti 2-2)

( ) () we define the interaction representation by the unitary transformation () = ()

Chapter 5: Force and Motion I-a

Circuits Op-Amp. Interaction of Circuit Elements. Quick Check How does closing the switch affect V o and I o?

4F-5 : Performance of an Ideal Gas Cycle 10 pts

Let s start from a first-order low pass filter we already discussed.

1/2 and e0 e s ' 1+ imm w 4 M s 3 πρ0 r 3 m. n 0 ktr. .Also,since n 0 ktr 1,wehave. 4 3 M sπρ 0 r 3. ktr. 3 M sπρ 0

Optimal Control. Lecture. Prof. Daniela Iacoviello

Name Student ID. A student uses a voltmeter to measure the electric potential difference across the three boxes.

Scattering at an Interface: Oblique Incidence

AP Statistics Notes Unit Two: The Normal Distributions

SUMMER REV: Half-Life DUE DATE: JULY 2 nd

Lecture 3: Resistive forces, and Energy

10. A.C CIRCUITS. Theoretically current grows to maximum value after infinite time. But practically it grows to maximum after 5τ. Decay of current :

Appendix H: Rarefaction and extrapolation of Hill numbers for incidence data

Observer Design for Nonlinear Systems using Linear Approximations

Lecture 17: Free Energy of Multi-phase Solutions at Equilibrium

CONVEX COMBINATIONS OF ANALYTIC FUNCTIONS

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

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

Control Systems. Mathematical Modeling of Control Systems.

Transcription:

MIT Hea ranfer Dynamc mdel 4.3./SG nvecn and cndcn and lmped mdel. Hea cnvecn If we have a rface wh he emperare and a rrndng fld wh he emperare a where a hgher han we have a hea flw a Φ h [W] () where he area f he rface and h he ceffcen f hea ranfer. a h Fgre : Hea cnvecn If he rrndng fld crclaed by a fan r pmp we are alkng f frced cnvecn. If he fld mn a rel f byancy frce, we are alkng f free (r naral) cnvecn. The hery fr he calclan f he ceffcen f hea ranfer he man bjec f he dcplne f Hea Tranfer. gd nrdcn can be een n Hea and Ma Tranfer by Incrpera and DeW, Wley,. Fr many prpe ne can e emaed vale - ypcal vale are gven n able : Prce h [W/(m K)] Free cnvecn Gae -5 Lqd 5- Frced cnvecn Gae 5-5 Lqd - Table : Typcal vale f he cnvecn hea ranfer ceffcen. Hea cndcn nder a ld cylnder a hwn n fgre. he ne end he emperare kep n and n he her end a and >. The rd nlaed n de,.e. hea nly ranmed n he hrnal drecn, frm he h end he clder end. k L Fgre : Hea cndcn n a ld cylnder (r fla plae wh hckne L)

MIT Hea ranfer Dynamc mdel 4.3./SG The hea flw can be calclaed a k Φ [W] () L where k he ceffcen f cndcvy fr he maeral ed, L he lengh and he cr ecnal area f he cylnder. Typcal vale fr k gven n able Maeral k [W/(m K)] Mneral fber.3-.5 Wd.-.6 Sanle Seel 3- Seel 4-6 lmnm bber 33 Table : Typcal vale f he cnvecn hea ranfer ceffcen 3. Thermal reance we wll ee laer n h paper cnvenen make an analgy elecrc crc. Eqan () can be rewren Φ a [W] (3) cnv where cnv he cnvecve reance hea ranfer. Y mgh remember Ohm Law ayng ha Vlage = eance me rren r U =. Or = U/ whch he ame a (3) f we be he hea flw wh he elecrc crren, he vlage drp by he emperare drp and he elecrc reance by he hermal reance. h cnv [K/W] (4) Fr cnvecn we have cnv = /(h ).e. he hgher h and/r he lwer reance, whch gve meanng. Fr he hea cndcn prblem we can make he ame analgy gvng he reance hea ranfer by cndcn L k cnd [K/W] (5) We wll e h analgy laer draw eqvalen hermal crc wh lmped capacance and reance elemen. 4. Mlple lmped parameer mdel Nw le ry cmbne he w prncple, ncldng he law f energy balance:

MIT Hea ranfer Dynamc mdel 3 4.3./SG Fgre 3 hw a mdel f a dble wall nlaed n he pace beween he w hermal heavy brck wall. h h Fgre 3: Dble wall (belw: n elecrc mdel f he wall) The elecrc analgy hw ha he prblem can be lked lke an elecrc crc wh 3 reance and w capacr We nw ame, ha he emperare n he w brck wall are alm even drbed, e. hee w wall can be lked lke lmped. Energy balance fr he nner wall, n. gve: Φ Φ (6) where Φ h (7) k Φ (8) ρ (9) c Energy balance fr he er wall, n. gve: Φ Φ () where Φ h () () ρ c Thee eqan can be lved n Smlnk, ee fgre 4

MIT Hea ranfer Dynamc mdel 4 4.3./SG Fgre 4: Mdel f he dble wall [ ] PHI [W] 3 Pl f emperare 5 5 5 5 5 Pl hea flw 5 5 5 me, a [h] PHI PHI PHI Fgre 5: el, baed n nal emperare = Daa fr he mdel: h_ = 8; % [W/(mK)] ef. f hea ranfer, nner h_ = ; % [W/(mK)] ef. f hea ranfer, er = ; % [m] rea k =.5; % [W/(m*K)] ef. f hea cndcn =.; % [m] Thckne, nlan _ = ; % [ ] Inner emperare _ = ; % [ ] Oer emperare c_ = 88; % [J/(kgK)] Spec. hea cap., nner wall c_ = 88; % [J/(kgK)] Spec. hea cap., er wall rh_ = 45; % [kg/m3] Deny, nner wall rh_ = 45; % [kg/m3] Deny, er wall _ =.; % [m] Thckne, nner wall _ =.; % [m] Thckne, er wall m_ m_ =rh_*_*; % [J/K] Hea capacy, nner wall =rh_*_*; % [J/K] Hea capacy, er wall n=; % Inal emperare, nner wall n=; % Inal emperare, er wall The ln f a gven nal vale prblem hwn n fgre 5. 5. Sae-Space frmlan pecal way f preenng he eqan called ae pace. On vecr frm lk lke h x x B y x D The bld ymbl ndcae array r marce

MIT Hea ranfer Dynamc mdel 5 4.3./SG Her x =, a clmn vecr repreenng he emperaren n he fve ne. a w dmennal vecr repreenng he amben emperare g. The y vecr an p vecr here he w emperare. Ung he example abve, wh he fve ne, we ge: If we defne h and k and h Then we have frm (6) (9) and frm () (): Then and B l n hen and D In h frmlan he Smlnk mdel mple: Fgre 6: Sae-Space frmlan

MIT Hea ranfer Dynamc mdel 6 4.3./SG The cnen n he Sae-Space blck a fllw _ = -(/(_*_)+/(_*_)); _ = /(_*_); _ = /(_*_); _ = -(/(_*_)+/(_*_)); B_ = /(_*_); B_ = ; B_ = ; B_ = /(_*_); =[ ; ] B =[B_ B_; B_ B_] =[ ; ] D =[ ; ] % Spdaa: mme =5*36; m('m wall', mme) fgre() h=pl(a/36, (:,), 'r-',a/36, (:,), 'b-'); e(h,'lnewdh',) grd % æer grd ax([-nf nf 3]) xlabel('a [h]') ylabel(' [ ]') legend('_','_') Pl % Sae-Space frmlan % 3.3./SG cle all clear clc 3 5 h_ = 8; % [W/(mK)] ef. f hea ranfer, nner h_ = ; % [W/(mK)] ef. f hea ranfer, er = ; % [m] rea k =.5; % [W/(m*K)] ef. f hea cndcn =.; % [m] Thckne, nlan _ = ; % [ ] Inner emperare _ = ; % [ ] Oer emperare c_ = 88; % [J/(kgK)] Spec. hea cap., nner wall c_ = 88; % [J/(kgK)] Spec. hea cap., er wall rh_ = 45; % [kg/m3] Deny, nner wall rh_ = 45; % [kg/m3] Deny, er wall _ =.; % [m] Thckne, nner wall _ =.; % [m] Thckne, er wall [ ] 5 5 5 5 5 a [h] m_ m _ =rh_*_*; % [J/K] Hea capacy, nner wall =rh_*_*; % [J/K] Hea capacy, er wall = c_*m_; % [J/K] Hea capacy, nner wall = c_*m_; % [J/K] Hea capacy, er wall _ = /(h_*); % [K/W] Thermal reance, nner cnvecn _ = /(k*); % [K/W] Thermal reance, nlan _ = /(h_*); % [K/W] Thermal reance, er cnvecn n=; % Inal emperare, nner wall n=; % Inal emperare, er wall