A MATHEMATICAL MODEL FOR NATURAL COOLING OF A CUP OF TEA

Similar documents
Institute of Actuaries of India

UNIT #5 EXPONENTIAL AND LOGARITHMIC FUNCTIONS

A THREE COMPARTMENT MATHEMATICAL MODEL OF LIVER

Boyce/DiPrima 9 th ed, Ch 2.1: Linear Equations; Method of Integrating Factors

Lecture 1: Numerical Integration The Trapezoidal and Simpson s Rule

Microscopic Flow Characteristics Time Headway - Distribution

An Indian Journal FULL PAPER. Trade Science Inc. A stage-structured model of a single-species with density-dependent and birth pulses ABSTRACT

Lecture 2: Current in RC circuit D.K.Pandey

Elementary Differential Equations and Boundary Value Problems

2.1. Differential Equations and Solutions #3, 4, 17, 20, 24, 35

Voltage v(z) ~ E(z)D. We can actually get to this wave behavior by using circuit theory, w/o going into details of the EM fields!

Spring 2006 Process Dynamics, Operations, and Control Lesson 2: Mathematics Review

General Article Application of differential equation in L-R and C-R circuit analysis by classical method. Abstract

Charging of capacitor through inductor and resistor

On the Derivatives of Bessel and Modified Bessel Functions with Respect to the Order and the Argument

CSE 245: Computer Aided Circuit Simulation and Verification

ERROR ANALYSIS A.J. Pintar and D. Caspary Department of Chemical Engineering Michigan Technological University Houghton, MI September, 2012

7.4 QUANTUM MECHANICAL TREATMENT OF FLUCTUATIONS *

Midterm exam 2, April 7, 2009 (solutions)

EE 434 Lecture 22. Bipolar Device Models

Logistic equation of Human population growth (generalization to the case of reactive environment).

The transition:transversion rate ratio vs. the T-ratio.

AR(1) Process. The first-order autoregressive process, AR(1) is. where e t is WN(0, σ 2 )

1. Inverse Matrix 4[(3 7) (02)] 1[(0 7) (3 2)] Recall that the inverse of A is equal to:

Applied Statistics and Probability for Engineers, 6 th edition October 17, 2016

Lecture 1: Growth and decay of current in RL circuit. Growth of current in LR Circuit. D.K.Pandey

On the Speed of Heat Wave. Mihály Makai

Boyce/DiPrima 9 th ed, Ch 7.8: Repeated Eigenvalues

4.1 The Uniform Distribution Def n: A c.r.v. X has a continuous uniform distribution on [a, b] when its pdf is = 1 a x b

Double Slits in Space and Time

whereby we can express the phase by any one of the formulas cos ( 3 whereby we can express the phase by any one of the formulas

a dt a dt a dt dt If 1, then the poles in the transfer function are complex conjugates. Let s look at f t H t f s / s. So, for a 2 nd order system:

Modelling of three dimensional liquid steel flow in continuous casting process

CHAPTER CHAPTER14. Expectations: The Basic Tools. Prepared by: Fernando Quijano and Yvonn Quijano

Transfer function and the Laplace transformation

UNSTEADY FLOW OF A FLUID PARTICLE SUSPENSION BETWEEN TWO PARALLEL PLATES SUDDENLY SET IN MOTION WITH SAME SPEED

5. An object moving along an x-coordinate axis with its scale measured in meters has a velocity of 6t

Lagrangian for RLC circuits using analogy with the classical mechanics concepts

10. If p and q are the lengths of the perpendiculars from the origin on the tangent and the normal to the curve

H is equal to the surface current J S

Economics 302 (Sec. 001) Intermediate Macroeconomic Theory and Policy (Spring 2011) 3/28/2012. UW Madison

Chapter 12 Introduction To The Laplace Transform

Wave Equation (2 Week)

Modeling and Experimental Investigation on the Internal Leakage in a CO2 Rotary Vane Expander

C From Faraday's Law, the induced voltage is, C The effect of electromagnetic induction in the coil itself is called selfinduction.

Control System Engineering (EE301T) Assignment: 2

Consider a system of 2 simultaneous first order linear equations

EXERCISE - 01 CHECK YOUR GRASP

The Relativistic Stern-Gerlach Force C. Tschalär 1. Introduction

Ch. 24 Molecular Reaction Dynamics 1. Collision Theory

Chapter 3: Fourier Representation of Signals and LTI Systems. Chih-Wei Liu

Poisson process Markov process

Math 3301 Homework Set 6 Solutions 10 Points. = +. The guess for the particular P ( ) ( ) ( ) ( ) ( ) ( ) ( ) cos 2 t : 4D= 2

Physics 160 Lecture 3. R. Johnson April 6, 2015

Impulsive Differential Equations. by using the Euler Method

DEPARTMENT OF ELECTRICAL &ELECTRONICS ENGINEERING SIGNALS AND SYSTEMS. Assoc. Prof. Dr. Burak Kelleci. Spring 2018

Let s look again at the first order linear differential equation we are attempting to solve, in its standard form:

Phys463.nb Conductivity. Another equivalent definition of the Fermi velocity is

CPSC 211 Data Structures & Implementations (c) Texas A&M University [ 259] B-Trees

14.02 Principles of Macroeconomics Problem Set 5 Fall 2005

Elements of Statistical Thermodynamics

Voltage v(z) ~ E(z)D. We can actually get to this wave behavior by using circuit theory, w/o going into details of the EM fields!

Text: WMM, Chapter 5. Sections , ,

4.3 Design of Sections for Flexure (Part II)

British Journal of Economics, Finance and Management Sciences 64 October 2011, Vol. 2 (1)

46. Let y = ln r. Then dy = dr, and so. = [ sin (ln r) cos (ln r)

Software Reliability using SPRT: Inflection S- shaped Model

Single Electron Devices for Logic Applications

Electrical Conductivity Measurement of Oxides Melts

Statistical Thermodynamics: Sublimation of Solid Iodine

Review Lecture 5. The source-free R-C/R-L circuit Step response of an RC/RL circuit. The time constant = RC The final capacitor voltage v( )

Midterm Examination (100 pts)

3(8 ) (8 x x ) 3x x (8 )

EE 350 Signals and Systems Spring 2005 Sample Exam #2 - Solutions

DE Dr. M. Sakalli

Coherence and interactions in diffusive systems. Cours 4. Diffusion + e-e interations

Quasi-Classical States of the Simple Harmonic Oscillator

CHAOS MODELS IN ECONOMICS

t=0 t>0: + vr - i dvc Continuation

The Variance-Covariance Matrix

are given in the table below. t (hours)

fiziks Institute for NET/JRF, GATE, IIT JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES MATEMATICAL PHYSICS SOLUTIONS are

Reliability Analysis of a Bridge and Parallel Series Networks with Critical and Non- Critical Human Errors: A Block Diagram Approach.

MEM 355 Performance Enhancement of Dynamical Systems A First Control Problem - Cruise Control

S.Y. B.Sc. (IT) : Sem. III. Applied Mathematics. Q.1 Attempt the following (any THREE) [15]

Summary: Solving a Homogeneous System of Two Linear First Order Equations in Two Unknowns

3.9 Carbon Contamination & Fractionation

0 +1e Radionuclides - can spontaneously emit particles and radiation which can be expressed by a nuclear equation.

k (but not necessarily much larger).

Arturo R. Samana* in collaboration with Carlos Bertulani*, & FranjoKrmpotic(UNLP-Argentina) *Department of Physics Texas A&M University -Commerce 07/

Introduction to Condensed Matter Physics

ANSWERS TO EVEN NUMBERED EXERCISES IN CHAPTER 11

Final Exam : Solutions

Higher order derivatives

Coherence and interactions in diffusive systems. Lecture 4. Diffusion + e-e interations

4. Which of the following organs develops first?

Instructors Solution for Assignment 3 Chapter 3: Time Domain Analysis of LTIC Systems

Almost power law : Tempered power-law models (T-FADE)

Ma/CS 6a Class 15: Flows and Bipartite Graphs

Themes. Flexible exchange rates with inflation targeting. Expectations formation under flexible exchange rates

Transcription:

MTHEMTICL MODEL FOR NTURL COOLING OF CUP OF TE 1 Mrs.D.Kalpana, 2 Mr.S.Dhvarajan 1 Snior Lcurr, Dparmn of Chmisry, PSB Polychnic Collg, Chnnai, India. 2 ssisan Profssor, Dparmn of Mahmaics, Dr.M.G.R Educaional and Rsarch Insiu, Chnnai, India. bsrac Hr i is proposd h rsuls of h simpl xprimn o masur h mpraur a any im in a cup of ho a. This aciviy offrs svral opporuniis o xplor mahmaics. In his invsigaion, h following mahmaical aciviis and concps ar applid: (i) horical principls modld by mahmaical xprssions (ii) xprssing mpirical rsuls in graphical form (iii) fiing horical prdicions o masurd daa (iv) comparing h qualiy of h horical prdicion wih masur daa. Hr w s our H o, h null hypohsis as; hr is no diffrnc bwn h rcordd valus and horically calculad valus. Kywords: Thrmodynamics, a, mpraur, ambin mpraur. Inroducion: n xprssion for h cooling of a cup of a is dvlopd from basic hrmodynamic principls. This hory is sd agains a simpl sysm, a cooling cup of a. In his xploraion an Excl spradsh was usd o graph h mpirical rsuls and h horical prdicion. Th acual daa and h mahmaical rprsnaion of h mpraur of a cooling of a 250ml cup of ho a is graphd vrsus im. = 0, h ha was rmovd from h a and h mpraur of h a was masurd a 1 minu inrvals. Room ambin mpraur was moniord a ½ minu inrvals o valida h assumpion ha hr was no significan im dpndnc in h room ambin mpraur. Th mach bwn h prdicd and acually valus suggss hr is somhing missing in h modl. Hnc a conjcur was mad as o h caus of his dispariy. I is proposd o us chi squar s o s h diffrnc bwn h obsrvd and horical prdicions. Thorical Modling: Th firs law can b xprssd as an quaion known as ra quaion. de W Qs I dscribs ha h ra of chang of nrgy in a sysm de is qual o h work, W, don on h sysm or powr nring h sysm minus ha ransfr ra Qs ou of h sysm. In our xprimn, 1

a =0, w assum ha all of h nrgy has bn impard o h sysm. Thrfor w s W=0 and his quaion bcoms: de Q s This quaion sas ha h ra of chang of nrgy in h sysm is qual o h ha ransfr ra. So a his poin, i is ncssary o raliz ha h nrgy of h sysm is dircly rlad o mpraur lvaion or chang and ha h ha ransfr ra can also b masur in rms of a mpraur diffrnc ovr im. Thrfor, his quaion can b wrin in rms of a mpraur diffrnc, T, vrsus im: dt () kt (), whr K is a consan Finding soluion: By solving h abov quaion, w g, Ingraing on boh sids, w g, dt () T () dt () T () k k log T( ) k C T () k C k T () T ( ) k, whr his im, i is convnin o impar physical maning o h consans usd in his soluion. lso no ha T() was dfind as a mpraur diffrnc, so ha h acual a mpraur is his mpraur diffrnc abov an ssnially im indpndn ambin mpraur: Thn h a mpraur T w can b xprssd as: C k T () T W Thr ar wo characrisics of h sysm and K. Th consan is rlad o h iniial mpraur ris ovr h ambin mpraur. Thrfor, T0 T, T0 is h iniial mpraur and T is h ambin mpraur (i..) room mpraur. Th consan K is rlad o h ha capaciy of h sysm and o h insulaing characrisics. So call h insulaing characrisics a hrmal rsisanc, q, C 2

and h ha capaciy of h a C. Sinc h grar h insulaion, or hrmal rsisanc, and h grar h ha capaciy, h slowr w would xpc h mpraur o chang, hs consans ar invrsly proporional o K, hus: c T () T T T W 0 Now ha a gnral xprssion for h a mpraur has bn dvlopd, do a quick saniy chck on how his quaion bhavs. = 0, h xponnial is uniy, and: T (0) T T T T W 0 0 This maks sns, sinc T O b h iniial mpraur of h ho a. s approachs zro and T ( ) w T, h room ambin mpraur as w would xpc. h xponnial Fiing Thorical Prdicions o Masurd Daa his poin w hav nough informaion from h mpirical daa o fill in h consans for our quaion for T W (). =T O -T = 69 C C = 38 C; his is h diffrnc bwn h saring mpraur and h final mpraur, which is h avrag valu a ambin h a mpraur is: w( ) 38 T o C =. So our quaion o prdic Th abl of daa of obsrvd and prdicions ar shown blow: SL NO Tim Room Tmpraur Obsrvd (O) Thorical (E) w( ) 38 T 1 0 29 69 69 2 5 64 62.98 3 10 32 59 57.92 4 15 55.5 53.65 5 20 50 50.07 6 25 29 48 47.05 7 30 45.2 44.51 8 35 42 42.37 9 40 32 41 40.57 3

10 45 38.4 39.05 11 50 32 37 37.78 12 55 35.4 36.7 13 60 32 33 35.8 Toal 617.5 617.45 vrags 47.5 47.496154 Th graph bwn im and obsrvd mpraur, im and xpcd mpraur and im and room mpraur is shown blow. T m p r a u r 80 70 60 50 40 30 20 10 0 Tim Vs Tmpraurs 0 10 20 30 40 50 60 70 Tim Room Tmpraur Obsrvd Tmprur Expcd Tmprur Comparing h Qualiy of h Thorical Prdicion wih Masurd: Hr, in ordr o compar h qualiy of h prdicd and obsrvd, i is usd paird s. calculaing valu for h abov daa, on can g calc 0.991126. Dgrs of frdom = (n-1) = (13-1) =12. By abl valu for 12 dgrs of frdom is 2.179. 0.991126 2.179 calc abl Hnc our null hypohsis is accpd. Thr is no diffrnc bwn h avrags of wo ss of valus. 4

Conclusion: In his aricl, w prdicd ha, hr is no significan diffrnc bwn h hory and obsrvd valus of a cooling sysm; w ook ha as Ho a for a cas sudy. Th mos obvious of his is h procss of vaporaion. This xprimn has bn conducd in Chnnai, India. Th acual rsuls may vary wih diffrn gographical rgions. Hr, i is akn undr crain room condiion mpraurs, in which hr was no allowanc of cooling sysms in h room. If so, h rsuls may vary. par from hs poins hr is apparnly an xcssivly rapid cooling, compard o h prdicion, followd by a slowing in h cooling. Hr w assum ha hr is no such rapidnss. Hr w omid fw variabls, lik h conainr. W ook h conainr as insulaor bu in som ohr siuaions, i may no hold. Bcaus of his, his rial is dubbd an opn sysm. Th sam concp can b xprimnd wih diffrn liquids o s h hory. Rfrncs: 1. Van Wyln, G. J. & Sonnag, R. E. (1973). Fundamnals of Classical Thrmodynamics (pp. 94). Nw York: John Wily and Sons, Inc. 2. Exrnal link: Chmical Thrmodynamics - Univrsiy of Norh Carolina 3. Hrbr B. Calln (1960). Thrmodynamics. Wily & Sons. Th clars accoun of h logical foundaions of h subjc. ISBN 0-471-13035-4. Library of Congrss Caalog No. 60-5597 Wb basd rfrncs: 4. hp://n.wikipdia.org/wiki/thrmodynamics 5