EJSS Cube Block Cooling Model

http://weelookang.blogspot.sg/2014/11/ejss-cube-block-cooling-model.html
increasing surface area of copper cube increases rate of heat loss to surrounding
model:https://dl.dropboxusercontent.com/u/44365627/lookangEJSworkspace/export/ejss_model_cooling/cooling_Simulation.xhtml
zip model: https://dl.dropboxusercontent.com/u/44365627/lookangEJSworkspace/export/ejss_model_cooling.zip
source: https://dl.dropboxusercontent.com/u/44365627/lookangEJSworkspace/export/ejss_src_cooling.zip
author: lookang, Christian wolfgang
http://weelookang.blogspot.sg/2014/11/ejss-cube-block-cooling-model.html
dull copper lose heat slower than shiny copper
model:https://dl.dropboxusercontent.com/u/44365627/lookangEJSworkspace/export/ejss_model_cooling/cooling_Simulation.xhtml
zip model: https://dl.dropboxusercontent.com/u/44365627/lookangEJSworkspace/export/ejss_model_cooling.zip
source: https://dl.dropboxusercontent.com/u/44365627/lookangEJSworkspace/export/ejss_src_cooling.zip
author: lookang, Christian wolfgang

Newton's Law of Cooling


The Newton's Law of Cooling model computes the temperature of an object of mass M as it is heated or cooled by the surrounding medium.

Assumption:

The model assumes that the temperature T within the object is uniform. 

Validity:

 This lumped system approximation is valid if the rate of thermal energy transfer within the object is faster than the rate of thermal energy transfer at the surface.

Convection-cooling "Newton's law of cooling" Model:

Newton assumed that the rate of thermal energy transfer at the object's surface is proportional to the surface area and to the temperature difference between the object and the surrounding medium.

δQδt=hA(T(t)Tbackground)

Q is the thermal energy in joules
h is the heat transfer coefficient (assumed independent of T here) (Wm2K)
A is the heat transfer surface area (m2)
T is the temperature of the object's surface and interior (since these are the same in this approximation)
Tbackground is the temperature of the surrounding background environment; i.e. the temperature suitably far from the surface is the time-dependent thermal gradient between environment and object.

Definition Specific Heat Capacity:

The specific heat capacity of a material on a per mass basis is

Q=mc(TfinalTinitial)
Q is heat energy
m  is the mass of the body
c specific heat capacity of a material
Tfinal is the Tbackground
Tinitial is the T(t)


combing the 2 equations

mc(TbackgroundT(t))δt=hA(T(t)Tbackground)

assuming mc is constant'

mcδ(TbackgroundT(t))δt=hA(T(t)Tbackground)

assuming Tbackground is a infinite reservoir

(Tbackground)δt=0
therefore
mc(δ(T(t)))δt=hA(T(t)Tbackground)

negative sign can be taken out of the differential equation.

mc(δT(t))δt=hA(T(t)Tbackground)

(T(t))δt=hAmc(T(t)Tbackground)

let κ=hAmc

(T(t))δt=κ(T(t)Tbackground)

If heating is added on,

heating=δQδt=mc(δTδt)

the final ODE equation looks like

(T(t))δt=κ(T(t)Tbackground)+heatingmc

Definition Equation Used:

V=mρ

V is volume of object
ρ is density of object


A=6(mρ)23

A surface area of object

assumption of increased surface are

Aincreasedsurfaceareaduetofins=(2)(6)(mρ)23

Materials added:

copper shiny cCu = 385  JkgK
ρCu = 8933  kgm3
heat transfer coefficient  hCu = 400 W(Km2)

copper dull cCu = 385  JkgK
ρCu = 8933  kgm3
heat transfer coefficient  hCu = 200 W(Km2)

aluminium shiny cAl = 903  JkgK
ρAl = 2702  kgm3
heat transfer coefficient  hAl = 400 W(Km2)

aluminium dull cAl =  903  JkgK
ρAl = 2702  kgm3
heat transfer coefficient  hAl = 200 W(Km2)

iron shiny cAl = 447  JkgK
ρAl = 7870  kgm3
heat transfer coefficient  hAl = 400 W(Km2)

iron dull cAl =  447  JkgK
ρAl = 7870  kgm3
heat transfer coefficient  hAl = 200 W(Km2)


Users can select the mass of the object and the material and the model computes the surface area assuming a cubic shape. The model plots the object's temperature as a function of time as the user heats and cools the object. A data-tool button on the temperature graph allows users fit the data to analytic functions.


Note: A typical (rough) heat transfer coefficient h for still air and iron is 6 W/(K m^2) and 400 W/(K m^2) . The Newton's Law of Cooling model assumes h=400 for all shiny and h=200 for dull materials. The actual value of h depends on many parameters including the material, the fluid velocity, the fluid viscosity and the condition of the object's surface.

References:


  1. "Measuring the Specific Heat of Metals by Cooling," William Dittrich, The Physics Teacher, (in press).

Credits:

  1. The Newton's Law of Cooling model was created by Wolfgang Christian using the Easy Java Simulations (EJS) version 4.2 authoring and modeling tool.
  2. EJSS Cube Block Cooling Model was created by Wolfgang Christian and recreated by lookang using the Easy Java Simulations (EJS) version 5.1 authoring and modeling tool


Posted by
Labels: ejss Java Script Thermodynamics