Download mel242-20

Industrial Calculations For Transient Heat
Conduction
P M V Subbarao
Associate Professor
Mechanical Engineering Department
IIT Delhi
pplications where rate/duration of heating/cooling
is a Design Parameter……
General Form of Fourier’s Equation
Cartesian coordinates:
Cylindrical coordinates:
Spherical coordinates:
One-dimensional Transient Conduction
• One-dimensional transient conduction refers to a case where
the temperature varies temporally and in one spatial direction.
• For example, temperature varies with x and time.
• Three cases of 1-D conduction are commonly studied:
conduction through a plate, in a cylinder, and in a sphere.
• In all three cases, the surface of the solid is exposed to
convection.
• The exact analytical solutions to the three cases are very
complicated.
• Industry uses an approximate solution, obtained by using
graphical tools.
• The graphs allow you to find the centerline temperature at
any given time, and the temperature at any location based on
the centerline temperature.
One Dimensional Fourier’s Equations
Constant thermal conductivity & No heat generation:
Cartesian coordinates:
T
T
 2
t
x
2



T

T
1

T
Cylindrical coordinates:
 2 

t
r r 
 r
2
2

T
 T 2 T 
Spherical coordinates:
 2 

t
r r 
 r
Primitive Variables & derived Variables
Characteristic Space Dimension : L or ro
Characterisitc Time Dimension : ?
Characteristic Temperature Variables:
Initial Temperature : T0
Far Field Fluid Temperature : T
Characterisitc Medium Property : 
Extent of Solution Domain
Highest Measure of Excess Temperature : q0 =|T0 -T|
Instantaneous local Excess temperature: q  |T -T|
Cartesian coordinates:
q
q
 2
t
x
2



q

q
1

q
Cylindrical coordinates:
 2 

t
r r 
 r
2
2

q
 q 2 q 
Spherical coordinates:
 2 

t
r r 
 r
Nature of Solution
Cartesian coordinates:
q  q0  f ( x, t :  , L)
q  q0  X ( x :  , L)  (t :  , L)
Cylindrical coordinates:
q  q0  g (r , t :  , ro )
q  q0  R(r :  , ro )  (t :  , ro )
Spherical coordinates:
q  q 0  h(r , t :  , ro )
q  q0  R(r :  , ro )  (t :  , ro )
Dimensional Analysis & Similarity Variables
Buckingham Method
The Buckingham theorem states that the function relationship
among n quantities or variables whose units may be given in terms
of m fundamental units or dimensions may be written as (n – m)
dimensionless groups.
The variables in local instantaneous temperature excess (q)
calculation are: q0, x, L, t, h,k,C and r.
Total number of variables : 8
Fundamental variables : 4
Number of non-dimensional variables : 4
Final form of Equations:
The π groups must be independent of each other and no one
group should be formed by multiplying together powers of other
groups.
The recurring set must contain three variables that cannot
themselves be formed into a dimensionless group.
In this case there are two restriction :
Both q and q0 cannot be chosen as they can be formed into a
dimensionless group, (q0/q).
Similarly, x and L cannot be chosen.
Heisler Parameters
• Heisler divided the problem into two parts.
• Part 1 : Instantaneous center line temperature. Variables are q0,,L, t,
and .
• Part 2 : Spatial temperature distribution for a given center
line temperature. Variables are : qcenter,,x,L, and .
•
•
•
•
•
•
•
Two different charts were developed.
Three parameters are needed to use each of these charts:
First Chart :
Normalized centerline temperature,
the Fourier Number,
and the Biot Number.
The definition for each parameter are listed below:
Mid Plane Temperature of Slab
Reading the Graphs
• To find the centerline temperature, first calculate Fo and Bi-1.
• This chart has Fo on the horizontal axis, θ0* on the vertical,
and lines representing different values of Bi-1.
• Locate the line corresponding to the Bi-1 you have calculated.
• Next, locate the value of Fo on the x-axis, and draw a vertical
line from the axis so that it intersects the Bi-1 line.
• From the intersection point, draw a horizontal line.
• The value at which the horizontal line crosses the y-axis is
the normalized centerline temperature.
• To find the temperature at any location within the solid, we
now will use the second chart.
Mid Plane Temperature of Slab
Temperature Distribution in A Slab
• The second chart has Bi-1 on the horizontal axis, normalized
temperature on the vertical axis, and lines corresponding to normalized
distance, x/L or r/r0.
• Calculate the normalized distance for the point which you are
interested in, and locate the corresponding line on the graph.
• Using the calculated value of Bi-1, draw a vertical line through the
graph.
• Draw a horizontal line from the point of intersection between
normalized distance curve and the vertical line.
• The point at which the horizontal line intersects the y-axis is the
normalized temperature at that location.
Temperature Distribution in A Slab
Change in Internal Energy of A Slab
Centre Line Temperature of An Infinite Cylinder
Temperature Distribution in An Infinite Cylinder
Change in Internal Energy of An Infinite Cylinder
Centre Temperature of A Sphere
Temperature Distribution in A Sphere
Change in Internal Energy of A Sphere
Multi-dimensional Transient Conduction
Finite Cartesian Bodies:
q ( x, y, z, t )  q0  X ( x, t )  Y ( y, t )  Z ( z, t )
Finite Cylindrical Bodies:
q ( x, y, z, t )  q0  R(r , t )  Z ( z, t )
Multi-dimensional Conduction
• The analysis of multidimensional conduction is simplified by
approximating the shapes as a combination of two or more
semi-infinite or 1-D geometries.
• For example, a short cylinder can be constructed by
intersecting a 1-D plate with a 1-D cylinder.
• Similarly, a rectangular box can be constructed by intersecting
three 1-D plates, perpendicular to each other.
• In such cases, the temperature at any location and time within
the solid is simply the product of the solutions corresponding
to the geometries used to construct the shape.
• For example, in a rectangular box, T(x*,y*,z*,t) - the
temperature at time t and location x*, y*, z* - is equal to the
product of three 1-D solutions: T1(x*,t), T2(y*,t), and T3(z*,t).
Transient Conduction in A Finite Cylinder