Numerical Transient Heat Conduction Experiment

Size: px
Start display at page:

Download "Numerical Transient Heat Conduction Experiment"

Transcription

1 Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use of fnte dfference to solve transent heat conducton problems. BACKGROUND he prmary law that descrbes conducton heat transfer s Fourer's Law of Heat Conducton. Fourer's Law s based upon the observaton that the conductve heat flux s drectly proportonal to the negatve of the temperature gradent or q & x - x (1) Introducng the thermal conductvty as the constant of proportonalty, we then have the equaton q & x = - k x () In ths experment we wll utlze Fourer's Law to study the problem of transent, one dmensonal heat conducton n a cylnder and to use the law n determnng the thermal conductvty of a sold. A cylndrcal element whch s embedded wth several thermocouples s heated at one end by an electrc hot plate and cooled at the other end by flowng water. he sde of the cylnder s very well nsulated so that the heat conducton s assumed to be one dmensonal. A schematc of the apparatus s shown below.

2 Fgure 1. Expermental Apparatus for Cylndrcal Element Water Flow hermocuples Hot Plate We consder the transent problem n whch the cylnder begns at some constant, unform temperature and then suddenly the hot plate s turned on so that a heat flux s mposed at the lower boundary. Our descrbng equaton for conservaton of energy balances the nternal energy change wth the axal heat conducton. Hence, n dfferental form we wrte = α t z (3) where α s the thermal dffusvty of the cylnder materal. One way to solve ths equaton s to dscretze the space doman, wrte t n fnte dfference form, and solve the subsequent system of algebrac equatons for the dscrtzed temperatures. If we use a second order correct approxmaton, the fnte dfference form of Eq. (3) becomes - t (j-1) = α (?) +1 + (?) ( z) - (?) -1 (4) he subscrpt on the represents the space (or z) dscetzaton whle the superscrpt represents the tme dscetzaton. Note that n Eq. (4) we have not specfed the tme dscretzaton (j or j-1) for the spatal dervatve. here are two choces for the tme dscretzaton for the temperatures n the spatal dervatve. hey could be evaluated at the prevous tme step, j-1, or at the current tme step, j. he algorthm s called explct f the temperatures are evaluated at the prevous tme step, and clearly ths makes the algebrac system very easy to solve snce Eq. (4) can be solved drectly for

3 wth no couplng to the other spatal nodes (+1 and -1). However, the explct approach does not always lead to a stable soluton (can you say blow-up) and n fact stablty s only guaranteed when α t ( z) < 0.5 Snce for most materals α s of the order of 10-5, ths stablty crtera often leads to unacceptable tme steps or spatal grds. he mplct algorthm, when the spatal dervatve s evaluated at the current tme step, does not have ths stablty problem, but does requre smultaneous soluton of the spatal node equatons. Fortunately, Mcrosoft Excel s a powerful spreadsheet tool that can carry out these smultaneous calculatons. Hence, for a spatal doman wth N spatal nodes, we would have N smultaneous equatons of the form - t (j-1) = α +1 + ( z) - -1 (5) to be solved for the 's. In steady state the axal temperature profle should be lnear whch confrms Fourer's Law. he steady state heat transfer s determned by measurng the mass flow rate and temperature change of a coolant stream whch passes over one end of the element, or q& = m& c p ( - ) out n coolant (6) hen the thermal conductvty can be calculated by q/a & cross-secton,element k = slope of vs z graph (7) PROCEDURE #1 (For sectons dong heatng up experment) 1. Makng sure that all coolant sample valves are closed turn on the water supply to the apparatus table.

4 . Open the valves to provde coolng to the hot plate assembles. Valves should be turned so as to pont at ether coolant n or coolant out. 3. Check to make sure all panel swtches and the hot plate swtches are n the off poston. Plug n the cable and turn on the power to the apparatus table. 4. urn the panel swtch for the hot plates (3 UNI 4) to the on poston. 5. Record the temperatures for the ten thermocouples for Unt 4. hese wll serve as your ntal temperatures for the transent conducton process. 6. Immedately set the hot plate swtch for Unt 4 to approxmately 50 C and start the stop watch. 7. At ten mnute ntervals record the temperatures for the ten thermocouples for Unt 4. Also record the tme requred to record these temperatures. After one hour the system should have reached steady state, whch can be confrmed by the lnearty of the temperature data wth poston or lttle change n the slope calculaton. 8. At steady state the energy delvered to the coolng water wll be determned. Place the empty beaker on the scale and zero the scale. Usng the coolant sample valve allow water to flow nto the beaker, measurng the tme wth a stop watch. When the beaker s nearly full turn the coolant sample valve to off. Record the tme. Place the flled beaker on the scale and record the mass. 9. At steady state record thermocouple readngs from channels 4&5 on unt 5. PROCEDURE # (For sectons dong coolng down experment) 1. If your secton s runnng a coolng down transent then the system s already at the heated steady state. At steady state the energy delvered to the coolng water wll be determned. Place the empty beaker on the scale and zero the scale. Usng the coolant sample valve allow water to flow nto the beaker, measurng the tme wth a stop watch. When the beaker s nearly full turn the coolant sample valve to off. Record the tme. Place the flled beaker on the scale and record the mass.. At steady state record thermocouple readngs from channels 4&5 on unt Record the temperatures for the ten thermocouples for Unt 4. hese wll serve as your ntal temperatures for the transent conducton process as well as the steady state temperatures from whch to determne the thermal conductvty. 4. Immedately turn the hot plate swtch for Unt 4 off and start the stop watch.

5 5. At ten mnute ntervals record the temperatures for the ten thermocouples for Unt 4. Also record the tme requred to record these temperatures. After one hour the system should have completely cooled down, whch can be confrmed by a constant temperature along the element. 6. urn the panel swtch for the hot plates (3 UNI 4) to the off poston. 7. Check to make sure all panel swtches and the hot plate swtches are n the off poston. Unplug n the cable and turn off the power to the apparatus table. 8. Close the valves to provde coolng to the hot plate assembles. 9. urn off the water supply to the apparatus table. Your data wll be entered on an Excel spreadsheet smlar to that shown below. ME 41 One Dmensonal, ransent Heat Conducton Experment Materal hermal Dffusvty : 1.17E-04 m^/s /C Dstance: m emperature emperature Locaton Expermental Numercal Locaton Expermental Numercal me = 0 sec. me = 600 sec. d/dx= d/dx= Steady State Calculaton of hermal Conductvty Coolant Mass: gm Coolant n: C me: sec Coolant out: C Mass Flow: kg/sec Coolant avg: C Water Cp: J/kg K Heat Flow: W hermal Conductvty: W/m K Whte cells ndcate data nput by the student. Lghtly shaded cells have calculaton equatons suppled.

6 DAA ANALYSIS 1. Graph both the expermental temperature and the numercal temperature versus poston at three dfferent tmes (early tme, moderate tme, and steady state). Show a sample calculaton for the heat suppled at steady state by the hot plate to the element. 3. Show a sample calculaton for the thermal conductvty of the element. SUGGESIONS FOR DISCUSSION 1. How do the numercal and expermental temperatures compare? Explan any dfferences or trends n the comparson.. What sort of shape does the temperature dstrbuton have at early tmes? How does ths compare wth what s predcted by analytcal solutons? 3. Based on the thermal conductvty, what would you suppose the element to be made of? 4. How could you check the one-dmensonalty of the experment?

7 able 1 Dameter (n.) Length (n.) 11 Poston of frst thermocouple from hot plate (n.) 1 Dstance between thermocouples (n.) Cylndrcal Element Unt Unt 5 hermocouple Channels Coolant n: Channel 5 Coolant out: Channel 4

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018 MATH 5630: Dscrete Tme-Space Model Hung Phan, UMass Lowell March, 08 Newton s Law of Coolng Consder the coolng of a well strred coffee so that the temperature does not depend on space Newton s law of collng

More information

2 Finite difference basics

2 Finite difference basics Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T

More information

Gravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11)

Gravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11) Gravtatonal Acceleraton: A case of constant acceleraton (approx. hr.) (6/7/11) Introducton The gravtatonal force s one of the fundamental forces of nature. Under the nfluence of ths force all objects havng

More information

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

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850) hermal-fluds I Chapter 18 ransent heat conducton Dr. Prmal Fernando prmal@eng.fsu.edu Ph: (850) 410-6323 1 ransent heat conducton In general, he temperature of a body vares wth tme as well as poston. In

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

Appendix B. The Finite Difference Scheme

Appendix B. The Finite Difference Scheme 140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton

More information

Finite Element Modelling of truss/cable structures

Finite Element Modelling of truss/cable structures Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures

More information

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient 58:080 Expermental Engneerng 1 OBJECTIVE Lab 2e Thermal System Response and Effectve Heat Transfer Coeffcent Warnng: though the experment has educatonal objectves (to learn about bolng heat transfer, etc.),

More information

Chapter Newton s Method

Chapter Newton s Method Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve

More information

FTCS Solution to the Heat Equation

FTCS Solution to the Heat Equation FTCS Soluton to the Heat Equaton ME 448/548 Notes Gerald Recktenwald Portland State Unversty Department of Mechancal Engneerng gerry@pdx.edu ME 448/548: FTCS Soluton to the Heat Equaton Overvew 1. Use

More information

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc

More information

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

Principles of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer Prncples of Food and Boprocess Engneerng (FS 31) Solutons to Example Problems on Heat Transfer 1. We start wth Fourer s law of heat conducton: Q = k A ( T/ x) Rearrangng, we get: Q/A = k ( T/ x) Here,

More information

Additional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty

Additional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,

More information

DUE: WEDS FEB 21ST 2018

DUE: WEDS FEB 21ST 2018 HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant

More information

The Analysis of Convection Experiment

The Analysis of Convection Experiment Internatonal Conference on Appled Scence and Engneerng Innovaton (ASEI 5) The Analyss of Convecton Experment Zlong Zhang School of North Chna Electrc Power Unversty, Baodng 7, Chna 469567@qq.com Keywords:

More information

Experiment 1 Mass, volume and density

Experiment 1 Mass, volume and density Experment 1 Mass, volume and densty Purpose 1. Famlarze wth basc measurement tools such as verner calper, mcrometer, and laboratory balance. 2. Learn how to use the concepts of sgnfcant fgures, expermental

More information

The Finite Element Method

The Finite Element Method The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of

More information

THE CURRENT BALANCE Physics 258/259

THE CURRENT BALANCE Physics 258/259 DSH 1988, 005 THE CURRENT BALANCE Physcs 58/59 The tme average force between two parallel conductors carryng an alternatng current s measured by balancng ths force aganst the gravtatonal force on a set

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

Announcements. Lecture #2

Announcements. Lecture #2 Announcements Lectures wll be n 4 LeConte begnnng Frday 8/29 Addtonal dscusson TA Denns Chang (Sectons 101, 105) Offce hours: Mo 2-3 PM; Th 5-6 PM Lab sectons begn Tuesday 9/2 Read Experment #1 onlne Download

More information

6.3.4 Modified Euler s method of integration

6.3.4 Modified Euler s method of integration 6.3.4 Modfed Euler s method of ntegraton Before dscussng the applcaton of Euler s method for solvng the swng equatons, let us frst revew the basc Euler s method of numercal ntegraton. Let the general from

More information

A Numerical Study of Heat Transfer and Fluid Flow past Single Tube

A Numerical Study of Heat Transfer and Fluid Flow past Single Tube A Numercal Study of Heat ransfer and Flud Flow past Sngle ube ZEINAB SAYED ABDEL-REHIM Mechancal Engneerng Natonal Research Center El-Bohos Street, Dokk, Gza EGYP abdelrehmz@yahoo.com Abstract: - A numercal

More information

Introduction to circuit analysis. Classification of Materials

Introduction to circuit analysis. Classification of Materials Introducton to crcut analyss OUTLINE Electrcal quanttes Charge Current Voltage Power The deal basc crcut element Sgn conventons Current versus voltage (I-V) graph Readng: 1.2, 1.3,1.6 Lecture 2, Slde 1

More information

ONE-DIMENSIONAL COLLISIONS

ONE-DIMENSIONAL COLLISIONS Purpose Theory ONE-DIMENSIONAL COLLISIONS a. To very the law o conservaton o lnear momentum n one-dmensonal collsons. b. To study conservaton o energy and lnear momentum n both elastc and nelastc onedmensonal

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

ME 501A Seminar in Engineering Analysis Page 1

ME 501A Seminar in Engineering Analysis Page 1 umercal Solutons of oundary-value Problems n Os ovember 7, 7 umercal Solutons of oundary- Value Problems n Os Larry aretto Mechancal ngneerng 5 Semnar n ngneerng nalyss ovember 7, 7 Outlne Revew stff equaton

More information

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity 1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum

More information

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00 ONE IMENSIONAL TRIANGULAR FIN EXPERIMENT Techncal Advsor: r..c. Look, Jr. Verson: /3/ 7. GENERAL OJECTIVES a) To understand a one-dmensonal epermental appromaton. b) To understand the art of epermental

More information

Analytical Chemistry Calibration Curve Handout

Analytical Chemistry Calibration Curve Handout I. Quck-and Drty Excel Tutoral Analytcal Chemstry Calbraton Curve Handout For those of you wth lttle experence wth Excel, I ve provded some key technques that should help you use the program both for problem

More information

Boise State University Department of Electrical and Computer Engineering ECE 212L Circuit Analysis and Design Lab

Boise State University Department of Electrical and Computer Engineering ECE 212L Circuit Analysis and Design Lab Bose State Unersty Department of Electrcal and omputer Engneerng EE 1L rcut Analyss and Desgn Lab Experment #8: The Integratng and Dfferentatng Op-Amp rcuts 1 Objectes The objectes of ths laboratory experment

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

More information

LAB 4: Modulus of elasticity

LAB 4: Modulus of elasticity LAB 4: Modulus of elastcty 1. Preparaton: modulus of elastcty (chapter15, p.79) Hook s law graphcal determnaton of modulus of elastcty (p.8) determnaton of modulus of elastcty n tenson and flexural stress

More information

Week 9 Chapter 10 Section 1-5

Week 9 Chapter 10 Section 1-5 Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,

More information

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

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

Mathematics Intersection of Lines

Mathematics Intersection of Lines a place of mnd F A C U L T Y O F E D U C A T I O N Department of Currculum and Pedagog Mathematcs Intersecton of Lnes Scence and Mathematcs Educaton Research Group Supported b UBC Teachng and Learnng Enhancement

More information

Total solidification time of a liquid phase change material enclosed in cylindrical/spherical containers

Total solidification time of a liquid phase change material enclosed in cylindrical/spherical containers Appled Thermal Engneerng 25 (2005) 488 502 www.elsever.com/locate/apthermeng Total soldfcaton tme of a lqud phase change materal enclosed n cylndrcal/sphercal contaners Levent Blr *, Zafer _ Ilken Department

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

Chapter 4 The Wave Equation

Chapter 4 The Wave Equation Chapter 4 The Wave Equaton Another classcal example of a hyperbolc PDE s a wave equaton. The wave equaton s a second-order lnear hyperbolc PDE that descrbes the propagaton of a varety of waves, such as

More information

Lecture 21: Numerical methods for pricing American type derivatives

Lecture 21: Numerical methods for pricing American type derivatives Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)

More information

Boise State University Department of Electrical and Computer Engineering ECE 212L Circuit Analysis and Design Lab

Boise State University Department of Electrical and Computer Engineering ECE 212L Circuit Analysis and Design Lab Bose State Unersty Department of Electrcal and omputer Engneerng EE 1L rcut Analyss and Desgn Lab Experment #8: The Integratng and Dfferentatng Op-Amp rcuts 1 Objectes The objectes of ths laboratory experment

More information

SIMPLE REACTION TIME AS A FUNCTION OF TIME UNCERTAINTY 1

SIMPLE REACTION TIME AS A FUNCTION OF TIME UNCERTAINTY 1 Journal of Expermental Vol. 5, No. 3, 1957 Psychology SIMPLE REACTION TIME AS A FUNCTION OF TIME UNCERTAINTY 1 EDMUND T. KLEMMER Operatonal Applcatons Laboratory, Ar Force Cambrdge Research Center An earler

More information

Math1110 (Spring 2009) Prelim 3 - Solutions

Math1110 (Spring 2009) Prelim 3 - Solutions Math 1110 (Sprng 2009) Solutons to Prelm 3 (04/21/2009) 1 Queston 1. (16 ponts) Short answer. Math1110 (Sprng 2009) Prelm 3 - Solutons x a 1 (a) (4 ponts) Please evaluate lm, where a and b are postve numbers.

More information

1-Dimensional Advection-Diffusion Finite Difference Model Due to a Flow under Propagating Solitary Wave

1-Dimensional Advection-Diffusion Finite Difference Model Due to a Flow under Propagating Solitary Wave 014 4th Internatonal Conference on Future nvronment and nergy IPCB vol.61 (014) (014) IACSIT Press, Sngapore I: 10.776/IPCB. 014. V61. 6 1-mensonal Advecton-ffuson Fnte fference Model ue to a Flow under

More information

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur Module Random Processes Lesson 6 Functons of Random Varables After readng ths lesson, ou wll learn about cdf of functon of a random varable. Formula for determnng the pdf of a random varable. Let, X be

More information

PART 8. Partial Differential Equations PDEs

PART 8. Partial Differential Equations PDEs he Islamc Unverst of Gaza Facult of Engneerng Cvl Engneerng Department Numercal Analss ECIV 3306 PAR 8 Partal Dfferental Equatons PDEs Chapter 9; Fnte Dfference: Ellptc Equatons Assocate Prof. Mazen Abualtaef

More information

2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12

2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12 REVIEW Lecture 11: 2.29 Numercal Flud Mechancs Fall 2011 Lecture 12 End of (Lnear) Algebrac Systems Gradent Methods Krylov Subspace Methods Precondtonng of Ax=b FINITE DIFFERENCES Classfcaton of Partal

More information

Lab session: numerical simulations of sponateous polarization

Lab session: numerical simulations of sponateous polarization Lab sesson: numercal smulatons of sponateous polarzaton Emerc Boun & Vncent Calvez CNRS, ENS Lyon, France CIMPA, Hammamet, March 2012 Spontaneous cell polarzaton: the 1D case The Hawkns-Voturez model for

More information

Definition. Measures of Dispersion. Measures of Dispersion. Definition. The Range. Measures of Dispersion 3/24/2014

Definition. Measures of Dispersion. Measures of Dispersion. Definition. The Range. Measures of Dispersion 3/24/2014 Measures of Dsperson Defenton Range Interquartle Range Varance and Standard Devaton Defnton Measures of dsperson are descrptve statstcs that descrbe how smlar a set of scores are to each other The more

More information

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

( ) = ( ) + ( 0) ) ( ) EETOMAGNETI OMPATIBIITY HANDBOOK 1 hapter 9: Transent Behavor n the Tme Doman 9.1 Desgn a crcut usng reasonable values for the components that s capable of provdng a tme delay of 100 ms to a dgtal sgnal.

More information

Numerical Solution of Boussinesq Equations as a Model of Interfacial-wave Propagation

Numerical Solution of Boussinesq Equations as a Model of Interfacial-wave Propagation BULLETIN of the Malaysan Mathematcal Scences Socety http://math.usm.my/bulletn Bull. Malays. Math. Sc. Soc. (2) 28(2) (2005), 163 172 Numercal Soluton of Boussnesq Equatons as a Model of Interfacal-wave

More information

Lecture 10 Support Vector Machines II

Lecture 10 Support Vector Machines II Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed

More information

Using Spectrophotometric Methods to Determine an Equilibrium Constant Prelab

Using Spectrophotometric Methods to Determine an Equilibrium Constant Prelab Usng Spectrophotometrc Methods to Determne an Equlbrum Constant Prelab 1. What s the purpose of ths experment? 2. Wll the absorbance of the ulbrum mxture (at 447 nm) ncrease or decrease as Fe soluton s

More information

MODULE 2: Worked-out Problems

MODULE 2: Worked-out Problems MODUE : Worked-out Problems Problem : he steady-state temperature dstrbuton n a one dmensonal wall of thermal conductvty 5W/m and thckness 5 mm s observed to be ( C) abx, where a C, B- c/ m, and x n meters

More information

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY. Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR

More information

Numerical modelization by finite differences of a thermoelectric refrigeration device of double jump". Experimental validation.

Numerical modelization by finite differences of a thermoelectric refrigeration device of double jump. Experimental validation. Numercal modelzaton by fnte dfferences of a thermoelectrc refrgeraton devce of double jump". Expermental valdaton. A. Rodríguez, J.G. Ván, D. Astran, Dpto. Ingenería Mecánca, Energétca y de Materales.

More information

Thermal State Investigations of Two-stage Thruster with Anode Layer

Thermal State Investigations of Two-stage Thruster with Anode Layer Thermal State Investgatons of Two-stage Thruster wth Anode Layer IEPC-2007-127 Presented at the 30 th Internatonal Electrc Propulson Conference, Florence, Italy Mtroshn A.S. *, Garkusha V.I., Semenkn A.V.

More information

Solution Thermodynamics

Solution Thermodynamics Soluton hermodynamcs usng Wagner Notaton by Stanley. Howard Department of aterals and etallurgcal Engneerng South Dakota School of nes and echnology Rapd Cty, SD 57701 January 7, 001 Soluton hermodynamcs

More information

Chapter 8. Potential Energy and Conservation of Energy

Chapter 8. Potential Energy and Conservation of Energy Chapter 8 Potental Energy and Conservaton of Energy In ths chapter we wll ntroduce the followng concepts: Potental Energy Conservatve and non-conservatve forces Mechancal Energy Conservaton of Mechancal

More information

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam. ME 270 Sprng 2017 Exam 1 NAME (Last, Frst): Please revew the followng statement: I certfy that I have not gven unauthorzed ad nor have I receved ad n the completon of ths exam. Sgnature: Instructor s Name

More information

Canonical transformations

Canonical transformations Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,

More information

Physics 141. Lecture 14. Frank L. H. Wolfs Department of Physics and Astronomy, University of Rochester, Lecture 14, Page 1

Physics 141. Lecture 14. Frank L. H. Wolfs Department of Physics and Astronomy, University of Rochester, Lecture 14, Page 1 Physcs 141. Lecture 14. Frank L. H. Wolfs Department of Physcs and Astronomy, Unversty of Rochester, Lecture 14, Page 1 Physcs 141. Lecture 14. Course Informaton: Lab report # 3. Exam # 2. Mult-Partcle

More information

Electrostatic Potential from Transmembrane Currents

Electrostatic Potential from Transmembrane Currents Electrostatc Potental from Transmembrane Currents Let s assume that the current densty j(r, t) s ohmc;.e., lnearly proportonal to the electrc feld E(r, t): j = σ c (r)e (1) wth conductvty σ c = σ c (r).

More information

Constitutive Modelling of Superplastic AA-5083

Constitutive Modelling of Superplastic AA-5083 TECHNISCHE MECHANIK, 3, -5, (01, 1-6 submtted: September 19, 011 Consttutve Modellng of Superplastc AA-5083 G. Gulano In ths study a fast procedure for determnng the constants of superplastc 5083 Al alloy

More information

Second Order Analysis

Second Order Analysis Second Order Analyss In the prevous classes we looked at a method that determnes the load correspondng to a state of bfurcaton equlbrum of a perfect frame by egenvalye analyss The system was assumed to

More information

Advanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)

Advanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017) Advanced rcuts Topcs - Part by Dr. olton (Fall 07) Part : Some thngs you should already know from Physcs 0 and 45 These are all thngs that you should have learned n Physcs 0 and/or 45. Ths secton s organzed

More information

Feb 14: Spatial analysis of data fields

Feb 14: Spatial analysis of data fields Feb 4: Spatal analyss of data felds Mappng rregularly sampled data onto a regular grd Many analyss technques for geophyscal data requre the data be located at regular ntervals n space and/or tme. hs s

More information

Section 8.3 Polar Form of Complex Numbers

Section 8.3 Polar Form of Complex Numbers 80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the

More information

Four Parameter Heat Transfer Turbulence Models for Heavy Liquid Metals

Four Parameter Heat Transfer Turbulence Models for Heavy Liquid Metals J. Energy Power Sources Vol., No., 015, pp. 63-74 Receved: January 9, 015, Publshed: February 8, 015 Journal of Energy and Power Sources www.ethanpublshng.com Four Parameter Heat Transfer Turbulence Models

More information

MAGNUM - A Fortran Library for the Calculation of Magnetic Configurations

MAGNUM - A Fortran Library for the Calculation of Magnetic Configurations CRYO/6/34 September, 3, 6 MAGNUM - A Fortran Lbrary for the Calculaton of Magnetc Confguratons L. Bottura Dstrbuton: Keywords: P. Bruzzone, M. Calv, J. Lster, C. Marnucc (EPFL/CRPP), A. Portone (EFDA-

More information

Implicit Integration Henyey Method

Implicit Integration Henyey Method Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

FUZZY FINITE ELEMENT METHOD

FUZZY FINITE ELEMENT METHOD FUZZY FINITE ELEMENT METHOD RELIABILITY TRUCTURE ANALYI UING PROBABILITY 3.. Maxmum Normal tress Internal force s the shear force, V has a magntude equal to the load P and bendng moment, M. Bendng moments

More information

How Differential Equations Arise. Newton s Second Law of Motion

How Differential Equations Arise. Newton s Second Law of Motion page 1 CHAPTER 1 Frst-Order Dfferental Equatons Among all of the mathematcal dscplnes the theory of dfferental equatons s the most mportant. It furnshes the explanaton of all those elementary manfestatons

More information

CHAPTER 14 GENERAL PERTURBATION THEORY

CHAPTER 14 GENERAL PERTURBATION THEORY CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves

More information

), it produces a response (output function g (x)

), it produces a response (output function g (x) Lnear Systems Revew Notes adapted from notes by Mchael Braun Typcally n electrcal engneerng, one s concerned wth functons of tme, such as a voltage waveform System descrpton s therefore defned n the domans

More information

Section 3.6 Complex Zeros

Section 3.6 Complex Zeros 04 Chapter Secton 6 Comple Zeros When fndng the zeros of polynomals, at some pont you're faced wth the problem Whle there are clearly no real numbers that are solutons to ths equaton, leavng thngs there

More information

CFD VALIDATION OF STRATIFIED TWO-PHASE FLOWS IN A HORIZONTAL CHANNEL

CFD VALIDATION OF STRATIFIED TWO-PHASE FLOWS IN A HORIZONTAL CHANNEL CFD VALIDATION OF STRATIFIED TWO-PHASE FLOWS IN A HORIZONTAL CHANNEL 1. Introducton Chrstophe Vallée and Thomas Höhne In dfferent scenaros of small break Loss of Coolant Accdent (SB-LOCA), stratfed twophase

More information

MAE140 - Linear Circuits - Winter 16 Final, March 16, 2016

MAE140 - Linear Circuits - Winter 16 Final, March 16, 2016 ME140 - Lnear rcuts - Wnter 16 Fnal, March 16, 2016 Instructons () The exam s open book. You may use your class notes and textbook. You may use a hand calculator wth no communcaton capabltes. () You have

More information

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law: CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and

More information

Visco-Rubber Elastic Model for Pressure Sensitive Adhesive

Visco-Rubber Elastic Model for Pressure Sensitive Adhesive Vsco-Rubber Elastc Model for Pressure Senstve Adhesve Kazuhsa Maeda, Shgenobu Okazawa, Koj Nshgch and Takash Iwamoto Abstract A materal model to descrbe large deformaton of pressure senstve adhesve (PSA

More information

AP Physics 1 & 2 Summer Assignment

AP Physics 1 & 2 Summer Assignment AP Physcs 1 & 2 Summer Assgnment AP Physcs 1 requres an exceptonal profcency n algebra, trgonometry, and geometry. It was desgned by a select group of college professors and hgh school scence teachers

More information

Formal solvers of the RT equation

Formal solvers of the RT equation Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan

More information

χ 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

χ 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 Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown

More information

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

4DVAR, according to the name, is a four-dimensional variational method. 4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The

More information

13. One way of expressing the power dissipated by a resistor is P = ( V)

13. One way of expressing the power dissipated by a resistor is P = ( V) Current and esstance 9. One way of expressng the power dsspated by a resstor s ( ). Thus, f the potental dfference across the resstor s doubled, the power wll be ncreased by a factor of 4, to a value of

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

TEST 5 (phy 240) 2. Show that the volume coefficient of thermal expansion for an ideal gas at constant pressure is temperature dependent and given by

TEST 5 (phy 240) 2. Show that the volume coefficient of thermal expansion for an ideal gas at constant pressure is temperature dependent and given by ES 5 (phy 40). a) Wrte the zeroth law o thermodynamcs. b) What s thermal conductvty? c) Identyng all es, draw schematcally a P dagram o the arnot cycle. d) What s the ecency o an engne and what s the coecent

More information

Numerical Simulation of Wave Propagation Using the Shallow Water Equations

Numerical Simulation of Wave Propagation Using the Shallow Water Equations umercal Smulaton of Wave Propagaton Usng the Shallow Water Equatons Junbo Par Harve udd College 6th Aprl 007 Abstract The shallow water equatons SWE were used to model water wave propagaton n one dmenson

More information

Supplementary Materials for

Supplementary Materials for advances.scencemag.org/cg/content/full/2/7/e1600304/dc1 Supplementary Materals for Interface-drven topologcal Hall effect n SrRuO3-SrIrO3 blayer Jobu Matsuno, Naok Ogawa, Kenj Yasuda, Fumtaka Kagawa, Wataru

More information

Important Instructions to the Examiners:

Important Instructions to the Examiners: Summer 0 Examnaton Subject & Code: asc Maths (70) Model Answer Page No: / Important Instructons to the Examners: ) The Answers should be examned by key words and not as word-to-word as gven n the model

More information

Thermal Conductivity of Cryogenic Deuterium

Thermal Conductivity of Cryogenic Deuterium hermal onductvty of ryogenc euterum Alan he Pttsford Mendon Hgh chool Pttsford, Y Advsor: r. R.. raxton Laboratory for Laser nergetcs Unversty of Rochester Rochester, Y ovember 007 Alan he Abstract Accurate

More information

Circuit Variables. Unit: volt (V = J/C)

Circuit Variables. Unit: volt (V = J/C) Crcut Varables Scentfc nestgaton of statc electrcty was done n late 700 s and Coulomb s credted wth most of the dscoeres. He found that electrc charges hae two attrbutes: amount and polarty. There are

More information

Note 10. Modeling and Simulation of Dynamic Systems

Note 10. Modeling and Simulation of Dynamic Systems Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada

More information

Determining the Temperature Distributions of Fire Exposed Reinforced Concrete Cross-Sections with Different Methods

Determining the Temperature Distributions of Fire Exposed Reinforced Concrete Cross-Sections with Different Methods Research Journal of Envronmental and Earth Scences 4(8): 782-788, 212 ISSN: 241-492 Maxwell Scentfc Organzaton, 212 Submtted: May 21, 212 Accepted: June 23, 212 Publshed: August 2, 212 Determnng the Temperature

More information

Chapter 13: Multiple Regression

Chapter 13: Multiple Regression Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to

More information