Heat Diffusion Equation

Steady State Conduction
=0
Heat Diffusion Equation
dE dE g dW
Energy balance equation:   qin  qout  E1  E 2 
dt dt dt
The above equation states:
Rate of change of total energy (dE/dt) = Rate of energy generated (dEg/dt) +
Rate of Heat Transfer in (qin) - Rate of Heat Transfer out (qout)
Apply this equation to a solid undergoing conduction heat
transfer as shown where:
dx E = mcpT= (rV)cpT = r(dxdydz)cpT
dy
T T
qin = q x   KA   K (dydz )
y x x x x
x q x
qout = q x  dx  q x  dq x  q x  dx
qx qx+dx x
Heat diffusion equation – mod 11/27/01

Heat Diffusion Equation (cont’d, page 2)
Substituting in the Energy Balance equation, we obtain:

EnergydE/dt Storage = EnergydEg/dtGeneration
qin -qout+ Net Heat Transfer

t
 r c pT  dxdydz  qdxdydz  qx  qx  dx
T qx
rcp dxdydz  qdxdydz  qx  ( qx  dx)
t x
 T
 qdxdydz  ( k ) dxdydz
x x
T  T  T  T
=> r c p  q  (k )  (k )  (k )
t x x y y z z
Note: partial differential operator

is used since T = T(x,y,z,t) Generalized to three-dimensions
Heat Diffusion Equation (cont’d, page 3)
T  T  T  T
rc p  q  (k )  (k )  (k )

t x x y y z z
Special case1: no generation q = 0
Special case 2: constant thermal conductivity k = constant
T 2T 2T 2T
rc p  k ( 2  2  2 )  q  k 2 T  q,
t x y z
 2
 2
 2
where  2  2  2  2 is the Laplacian operator
x y z

Special case 3:  0 and q = 0
t
 2
T  2
T  2
T
 T  2  2  2  0, The famous Laplace' s equation
2
x y z
1-D, Steady State Conduction
Assume steady and no generation, 1 - D Laplace' s equation
d 2T
2
 0, T ( x, y, x, t )  T ( x ), function of x coordinate alone
dx
Note: ordinary differential operator is used since T = T(x) only
The general solution of this equation can be determined by integting twice:
dT
First integration leads to  cons tan t  C1 .
Constant
dx
Integrate again T(x) = C1 x  C2
Second order differential equation: need two boundary conditions to
determine the two constants C1 and C 2 .
T Example:
T(x=0)=100°C=C2
100
T(x=1 m)=20°C=C1+C2, C1=-80°C
T(x)=100-80x (°C)
20
x
1-D Heat, Steady State Heat Transfer
Thermal Resistance (Electrical Circuit Analogy)
Recall Fourier' s Law:

dT
q = -kA , If the temperature gradient is a constant
dx
q = constent
Constant (heat transfer rate is a constant)
dT T2  T1
 , where T ( x  0)  T1 , T ( x  L)  T2
dx L
dT T1  T2 T1  T2
T1 T2 q   kA  kA 
dx L ( L / kA)
T1  T2 L
q , where R  :thermal resistance
R kA
x q (I)
This is analogous to an electrical circuit
L T2 (V2) I = (V1-V2)/R
T1 (V1)
R (R)
Thermal Resistance - Composite Wall Heat Transfer
Use of the thermal resistance concept make the analysis of complex geometries
relatively easy, as discussed in the example below.
However, note that the thermal heat resistance concept can only be applied for
steady state heat transfer with no heat generation.
Example: Consider a composite wall made of two different materials
R1=L1/(k1A) R2=L2/(k2A)
T T1
k2 T2
T1 k1 T2
T1  T2 T1  T2T T1  T2
q  
R R1  R2  L1   L2 
  
 1   2 
k A k A
T1  T  L1 
Also, q= , T  T1  qR1  T1  q  
L1 L2 R1  k1 A 
Composite Wall Heat Transfer (cont’d)
•Now consider the case where we have 2 different fluids on either sides of the wall at
temperatures, T,1 and T,2 , respectively.
• There is heat transfer by convection from the first fluid (on left) to material 1
and from material 2 to the second fluid (on right).
• Similar to conduction resistance, we can determine the convection resistance,
where Rconv = 1/hA
R1=L1/(k1A) R2=L2/(k2A)
Rconv,2= 1/(h2A)
Rconv,1= 1/(h1A)
T
T1 k1 k2 T,2
T2 T,1 T1 T T2
T , 2  T, 2 T , 2  T1 T1  T
where Q   
T,1 T,2 Rtot Rconv ,1 R1
T  T2 T 2  T , 2
 
R2 Rconv , 2
L1 L2 This approach can be extended to much more complex geometries
(see YAC, 8-4 and 8-5)
“R” value of insulations
In the US, insulation materials are often specified in terms of their thermal
resistance, denoted as R values, where Rvalue = L/k (thickness/thermal conductivity)
In the US, it has units of (hr ft2 °F)/Btu (Note: 1 Btu=1055 J)
R-11 for wall, R-19 to R-31 for ceiling.
Q: Why is the Rvalue given by L/k ?
A: From previous slide, we know that the thermal resistance to conduction is
defined as the temperature difference across the insulation by the heat flux going
through it, hence:
T T x
R  
q" k T k
x
Example:
The typical space inside the residential frame wall is 3.5 in. Find the R-value if
the wall cavity is filled with fiberglass batt. (k=0.046 W/m.K=0.027 Btu/h.ft.R)
x 0.292 ft
R   10.8( R. ft 2 . h / Btu)  R  11
k 0.027 Btu / h. ft. R

Heat Diffusion Equation

Uploaded by

Document Informationclick to expand document information

Copyright:

Available Formats

Heat Diffusion Equation

Uploaded by

Document Information

Original Description:

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Heat Diffusion Equation

Uploaded by

Copyright:

Available Formats

Steady State Conduction

Heat diffusion equation – mod 11/27/01

Substituting in the Energy Balance equation, we obtain:

Note: partial differential operator

Recall Fourier' s Law:

You might also like