heat equations of change ii
DESCRIPTION
Heat Equations of Change II. Outline. 7. Heat Equations of Change 7.1. Derivation of Basic Equations 7.1.1. Differential Equation for Heat Conduction 7.1.2. Energy Equation 7.1.3. Buckingham Pi Method 7.2. Unsteady-state Conduction - PowerPoint PPT PresentationTRANSCRIPT
![Page 1: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/1.jpg)
Heat Equations of Change II
![Page 2: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/2.jpg)
7. Heat Equations of Change
7.1. Derivation of Basic Equations
7.1.1. Differential Equation for Heat Conduction
7.1.2. Energy Equation
7.1.3. Buckingham Pi Method
7.2. Unsteady-state Conduction
7.2.1. Gurney-Lurie Charts
7.2.2. Lumped Systems Analysis
Outline
![Page 3: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/3.jpg)
Quiz β 2014.02.14
An electrically heated resistance wire has a diameter of 2 mm and a resistance of 0.10 ohm per foot of wire. The thermal conductivity of the wire is 20 W/mΒ·K. At a current of 100 A, calculate the steady state temperature difference between the center and surface of the wire.
TIME IS UP!!!
![Page 4: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/4.jpg)
From the previous lectureβ¦
A 10-cm diameter nickel-steel sphere has a thermal conductivity, k = 10 W/m-K. Within the sphere, 800 W/m3 of heat is being generated. The surrounding air is at 20Β°C and the heat transfer coefficient from the surroundings to the surface of the sphere is 10 W/m2-K. What is the temperature at the center of the sphere?
Exercise!
![Page 5: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/5.jpg)
A Comparison Between Methods
Energy Equation:Solution!
We are left with: 0=π[ 1π2 πππ (π2 ππππ )]+πBecause T depends on r only: 0=
πππ (π2 ππππ )+ππ
2
π
ππππ·ππ·π‘ =ππ»2π +π
ππΆππ·ππ·π‘ =β (π» βπ )β (π βπ» π )
In spherical coordinates:
From the Energy Equation
For solids with generation
+π
![Page 6: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/6.jpg)
1π2
πππ (π 2 ππππ )+ 1
π2 sinππππ (sinπ ππππ )+ 1
π 2sin 2ππ2πππ2
+ππ=
1πΌππππ‘
A Comparison Between Methods
Differential Equation of Heat Conduction:
Solution! Another way!!
We are left with: 0= 1π2
πππ (π2 ππππ )+ππ
Because T depends on r only:
π»2π+ππ=
1πΌππππ‘
In spherical coordinates:
From the Diff. Eqn of Heat Conduction0= π
ππ (π2 ππππ )+ππ2
π
![Page 7: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/7.jpg)
A Comparison Between Methods
Overall Shell Heat Balance:
Solution! Another way again!! The Shell:
Rate of Heat IN:
Rate of Heat OUT:
Generation:
(4π π2ππβ)|π(4π π2ππβ)|π+βπ(4ππ 2βπ )π
(π 2ππβ)|π+βπβ (π 2ππβ )|πβπ =ππ2
Adding the terms and dividing :
πππ (π 2ππ )=ππ2
![Page 8: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/8.jpg)
A Comparison Between Methods
Overall Shell Heat Balance:
Solution! Another way again!! The Shell:
Rate of Heat IN:
Rate of Heat OUT:
Generation:
(4π π2ππβ)|π(4π π2ππβ)|π+βπ(4ππ 2βπ )π
β πππ (π2π ππππ )=ππ2
Inputting Fourierβs Law:
πππ (π 2ππ )=ππ2
![Page 9: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/9.jpg)
A Comparison Between Methods
Overall Shell Heat Balance:
Solution! Another way again!! The Shell:
Rate of Heat IN:
Rate of Heat OUT:
Generation:
(4π π2ππβ)|π(4π π2ππβ)|π+βπ(4ππ 2βπ )π
Inputting Fourierβs Law:
From the Overall Shell Heat Balance0= π
ππ (π2 ππππ )+ππ2
πβ πππ (π2π ππππ )=ππ2
![Page 10: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/10.jpg)
A Comparison Between Methods
Solution! Apparentlyβ¦
0= πππ (π2 ππππ )+ππ
2
πFrom the Overall
Shell Heat Balance
From the Diff. Eqn of Heat Conduction0= π
ππ (π2 ππππ )+ππ2
π
0= πππ (π2 ππππ )+ππ
2
πFrom the
Energy Equation
Using any method below, we can obtain the same ODE to solve!!
Increasing range of
applicability of the method
![Page 11: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/11.jpg)
A Comparison Between Methods
Solution!
0= πππ (π2 ππππ )+ππ
2
πFrom the Overall
Shell Heat Balance
From the Diff. Eqn of Heat Conduction0= π
ππ (π2 ππππ )+ππ2
π
0= πππ (π2 ππππ )+ππ
2
πFrom the
Energy Equation
The next step is to solve the ODE and define boundary conditions.
Integrating twice:
π (π )=β ππ2
6π βπΆ1
π +πΆ2
Boundary conditions:
ππ‘ π=0 , ππ ππ πππππ‘π
ΒΏππ‘ π=π , π ππππ =h (π πππβπ (π ) )
![Page 12: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/12.jpg)
A Comparison Between Methods
Solution!
Applying the boundary conditions:
π (π )=π πππ+ππ 2
6 π [1β( ππ )2]+ππ 3h
π=10 πππΎh=10 π
π2πΎ
π =0.05ππ=800π
π3
π πππ=20 Β°πΆ
Final Answer:
21.37Β°C
![Page 13: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/13.jpg)
7. Heat Equations of Change
7.1. Derivation of Basic Equations
7.1.1. Differential Equation for Heat Conduction
7.1.2. Energy Equation
7.1.3. Buckingham Pi Method
7.2. Unsteady-state Conduction
7.2.1. Gurney-Lurie Charts
7.2.2. Lumped Systems Analysis
Outline
![Page 14: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/14.jpg)
Example!
1-D Transient Heat Conduction
Consider a flat slab of thickness L with an initial temperature of T0. It is later submerged in a large fluid with temperature T1 (> T0). Neglecting surface resistance (infinite h), determine the temperature at the center of the slab at any time t.
![Page 15: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/15.jpg)
Example!
1-D Transient Heat Conduction
What do we expect to happen?
β’ The surface temp. is held at T1. (because of infinite h)
β’ Temperature profile at any time is symmetrical
β’ Unidirectional (x only) heat flow only
β’ The temp. at the center, T0, slowly approaches T1.
T0
T1
![Page 16: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/16.jpg)
Example! Solution:
1-D Transient Heat Conduction
T0
T1
Differential Equation of Heat Conduction:
π»2π+ππ=
1πΌππππ‘
With the assumptions:ππππ‘ =πΌ π
2πππ₯2
Initial and Boundary conditions:
Initial:
Boundary:
![Page 17: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/17.jpg)
Example! Solution:
1-D Transient Heat Conduction
ππππ‘ =πΌ π
2πππ₯2
To facilitate solving, we will make this dimensionless!
Let: Y = dimensionless temp.
π=π1βππ 1βπ 0
*The choice of definition of Y is arbitrary.
ππππ‘ =
β1π 1βπ0
ππππ‘
Implications:
ππππ₯ =
β1π 1βπ0
πππ π₯
π2ππ π₯2
=β1
π 1βπ0π2πππ₯2
β (π1βπ 0 ) ππππ‘ =βπΌ (π 1βπ 0 ) π2ππ π₯2
![Page 18: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/18.jpg)
1-D Transient Heat Conduction
ππππ‘ =πΌ π
2ππ π₯2
To facilitate solving, we will make this dimensionless!
Let: X = dimensionless length
π=π₯πΏ
*The choice of definition of X is arbitrary.
π π=1πΏ ππ₯
Implications:
ππππ‘ =
πΌπΏ2π2ππ π 2
π π 2=( 1πΏππ₯)2
=1πΏ2ππ₯2
Example! Solution:
![Page 19: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/19.jpg)
1-D Transient Heat Conduction
ππππ‘ =
πΌπΏ2π2ππ π 2
To facilitate solving, we will make this dimensionless!
Let: Ο = dimensionless time
π=πΌπ‘πΏ2
*The choice of definition of Ο is arbitrary.
ππ= πΌπΏ2ππ‘
Implications:
ππππ =
π2ππ π 2
Example! Solution:
![Page 20: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/20.jpg)
1-D Transient Heat Conduction
ππππ =
π2ππ π 2
DIMENSIONLESS FORM:
Previous IC & BC: New IC & BC:
Easier to solve and no worries with units!
π=1 , π=0 , π=ππ=0 , π=π , π=0π=0 , π=π , π=1
π=πΌπ‘πΏ2
π=π₯πΏπ=
π1βππ 1βπ 0
Example! Solution:
![Page 21: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/21.jpg)
A new dimensionless numberβ¦
1-D Transient Heat Conduction
π=πΌπ‘πΏ2
π=π₯πΏπ=
π1βππ 1βπ 0
Dim. Group Ratio EquationFourier, Fo Rate of heat conduction/
Rate of heat storage
Rate of heat conduction:
ππΏ3
Rate of heat storage: ππππΏ3/π‘
![Page 22: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/22.jpg)
Going backβ¦
1-D Transient Heat Conduction
ππππ =
π2ππ π 2
DIMENSIONLESS FORM:
New IC & BC:
π=1 , π=0 , π=ππ=0 , π=π , π=0π=0 , π=π , π=1
π=πΌπ‘πΏ2
π=π₯πΏπ=
π1βππ 1βπ 0
T0
T1
![Page 23: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/23.jpg)
Going backβ¦
1-D Transient Heat Conduction
ππππ =
π2ππ π 2
DIMENSIONLESS FORM:
New IC & BC:π=1 , π=0 , π=ππ=0 , π=π , π=0π=0 , π=π , π=1
π=πΌπ‘πΏ2
π=π₯πΏπ=
π1βππ 1βπ 0
Solution: (A Fourier series)
π (π ,π )= 4πβπ=0
β πβ (2π+1)2 π2π
2π+1sin [ (2π+1 ) π π ]
![Page 24: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/24.jpg)
Gurney-Lurie Charts
Gurney-Lurie Charts
π=πΌπ‘πΏ2
π=π1βππ 1βπ 0
- plots of the dimensionless temperature Y against Fo with varying Bi and X for different geometries.
- Each point in the curves are solutions to the PDE involving heat conduction + convection.
π=π₯πΏ
![Page 25: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/25.jpg)
For flat plates with convection
Gurney-Lurie Charts
π=πΌπ‘πΏ2
π=π1βπ
π1βπ0
Bi
Geankoplis, Figure 5.3-6
![Page 26: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/26.jpg)
For long cylinders with convection
Gurney-Lurie Charts
π=πΌπ‘πΏ2
π=π1βπ
π1βπ0
Bi
Geankoplis, Figure 5.3-8
![Page 27: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/27.jpg)
For solid spheres with convection
Gurney-Lurie Charts
π=πΌπ‘πΏ2
π=π1βπ
π1βπ0
Bi
Geankoplis, Figure 5.3-10
![Page 28: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/28.jpg)
Gurney-Lurie Charts
How do we compute the Biot Numbers for different geometries?
Recall: Biot Number π΅π=h π₯1π
Characteristic length: Volume/Surface Area
![Page 29: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/29.jpg)
From a while agoβ¦
A 10-cm diameter nickel-steel sphere has a thermal conductivity, k = 10 W/m-K. Within the sphere, 800 W/m3 of heat is being generated. The surrounding air is at 20Β°C and the heat transfer coefficient from the surroundings to the surface of the sphere is 10 W/m2-K. What is the temperature at the center of the sphere?
Recall this exercise:
![Page 30: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/30.jpg)
From a while agoβ¦
A 10-cm diameter nickel-steel sphere has a thermal conductivity, k = 10 W/m-K. Within the sphere, 800 W/m3of heat is being generated.The sphere, initially at T0 = 30Β°C, is suddenly submerged into a fluidβ¦The surrounding air is at and the heat transfer coefficient from the surroundings to the surface of the sphere is 10 W/m2-K. What is the temperature at the center of the sphere across time?
Recall this exercise:
![Page 31: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/31.jpg)
From a while agoβ¦
To solve the new problem:Differential Equation of Heat Conduction:
π»2π+ππ=
1πΌππππ‘
1π2
πππ (π2 ππππ )+ 1
π2 sinππππ (sinπ ππππ )+ 1
π 2sin2ππ2πππ2
+ππ=
1πΌππππ‘
In spherical coordinates:
Steps:1. Turn the remaining PDE into
dimensionless form.2. Solve analytically for T(r, t).
But there is an easier way!
![Page 32: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/32.jpg)
Lumped Systems Analysis
Letβs assume that the sphere is too small for conduction to matter. The temperature distribution inside the sphere can, therefore, be assumed uniform!
Heat Balance: ( π ππ‘πππ hπππ‘ ππππ€πππ‘ππ ππππππ π£πππ’πππ h hπ‘ πππ’ππππ’πππππ¦ π π’ππππππ π΄ )=(π ππ‘πππ πππππππ ππππππ‘πππππππππππ¦ ππ
π ππππππ π£πππ’πππ )hπ΄ (π ββπ (π‘ ) )=ππππ
ππ (π‘)ππ‘
Initial Condition:
π (π‘ )=π 0 , π‘=0ππ (π‘)ππ‘ = hπ΄
ππππ(π ββπ (π‘ ) )
![Page 33: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/33.jpg)
Lumped Systems Analysis
ππ (π‘)ππ‘ = hπ΄
ππππ(π ββπ (π‘ ) )Initial Condition:
π (π‘ )=π 0 , π‘=0ππ (π‘ )
π ββπ (π‘)= hπ΄ππππ
ππ‘
β ln(π ββπ (π‘)πββπ0 )= hπ΄ π‘
ππππIntegrating and plugging the IC:
Rearranging: π (π‘ )βπβπ 0βπ β
=exp ( β hπ΄ π‘ππππ )
![Page 34: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/34.jpg)
Lumped Systems Analysis
Rearranging: π (π‘ )βπβπ 0βπ β
=exp ( β hπ΄ π‘ππππ )
(βh (π / π΄ )π )( π
ππππ‘
(π / π΄ )2 )
(βπ΅π ) (πΉπ )Temperature Profile at the center of the sphere across time:
π (π‘ )βπβπ 0βπ β
=πβπ΅ππΉπ
![Page 35: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/35.jpg)
Lumped Systems Analysis
Temperature Profile at the center of the sphere across time:
π (π‘ )βπβπ 0βπ β
=πβπ΅ππΉπ
This Lumped System Analysis is sufficiently accurate only when
Bi < 0.1.
π=10 πππΎh=10 π
π2πΎ
π =0.05ππ 0=30 Β°πΆπ β=20 Β°πΆ
π΅π=10 ππ2π
10 πππΎ( 0.05π3 )=0.0167<0.1
Given:
Valid!!Checking:
![Page 36: Heat Equations of Change II](https://reader035.vdocument.in/reader035/viewer/2022062501/56816034550346895dcf576c/html5/thumbnails/36.jpg)
Lumped Systems Analysis
For review:A person is found dead at 5 PM in a room where T = 20Β°C. The temperature of the body was measured at 25Β°C when found. The heat transfer coefficient is estimated to be 0.8 W/m2K. Estimate the time of death assuming:
(1) The body can be modeled as a 30-cm diameter, 1.7-m long cylinder.(2) The thermal properties of the body and the heat transfer are
constant.(3) The body temperature was 37Β°C at the time of death.(4) Since the human body has almost the same properties of water at
the average temperature of 31Β°C: k = 0.617 W/mK, density = 996 kg/m3, and cp = 4178 J/kgK.
(5) Radiation effects are negligible.