Cooling A Pipe Filled With Water
Cooling A Pipe Filled With Water
Cooling A Pipe Filled With Water
This application calculates the cooling down time of water in a carbon steel pipe.
Sheets 1 and 2, present the necessary theory, while sheets 3 and 4 present pipe
applications with and without insulation.
Sheet 5Calculation method of the National Mechanical Insulation Committee (NMIC)
Sheet 6 applies the developed routine to an example of the NMIC
Sheets 7 to 9 present routines for the calculation of convection in air
Sheet 10 presents a graphic of cool down time versus pipe diameter.
Sheet 11 presents references
VBA window
Committee (NMIC)
Resume of calculations and information
Sheet 1
1.- Coolling down of a pipe
This sheet contains:
- Derivation of the cool down equation
It is shown that the ASHRAE equation is a special case of the derived equation.
- Derivation of the equation for the water temperature as a function of time
Sheet 2
2.- t_average
- Derivation of the equation for the average water temperature
Sheet 3
3.- Cooling down of water til 0°C, in a pipe with insulation
Spreadsheet for the calculation of the cool down time to 0°C, for an insulated pipe
Application of this equation is checked with an application of the National Mechanical
Insulation Commitee (NMIC), Sheet 6.
Sheet 4
4.- Cooling down of water til 0°C, in a pipe without insulation
Spreadsheet for the calculation of the cool down time to 0°C, for an uninsulated pipe
Application cannot be checked with NMIC method, since this equation is not adequate
for this type of case.
Rev. cjc. 30.01.201
Sheet 5
5.- NMIC
NMIC Equations (National Mechanical Insulation Commitee)
Calculation method from the NMIC
ed equation.
Sheet 6
6.- NMIC example
Application of the equation and comparison with results of a NMIC table.
Sheet 7
7.- Outside tube combined natural and forced convection
Calculation routine for natural convection around an horizontal pipe in air
Sheet 9
9.- Cooling down time curve
uninsulated pipe Grphic of cooling down time (h) versus pipe diemeter (in)
on is not adequate for the given conditions
Sheet 10
10.- U-factor
Heat transfer coefficient for a pipe with insulation
Ref
References
Rev. cjc. 30.01.2014
1.- Cooling down of a pipe filled with water
d1
sI, kI
sp, kp
L
τf d 1 ¿ ρ⋅C P t
τ|0 =− ¿ ln ( t−t a )|t f
4⋅U 1 i
1.2 Global heat transfer coefficient [1]
1
q= A1⋅U 1⋅Δt
U 1=
Δt
( ) ( )
r2 r3
r 1⋅ln
r1
r 1⋅ln
r2 r1
(Eq. 1.4) q=
1
h1
+
kp
+
kI
+
h3⋅r 3
RT
1
= A 1⋅U 1
1 2⋅π⋅L⋅k
If the interior convection has a very
RT =
()
large value RT r3
1 ln
U 1= For the case of Equation 1.7 r2
() ()
r2 r3
r 1⋅ln r 1⋅ln (Eq. 1.5)
1
r1 r2 r1 U 1=
()
+ + r3
( )
kp kI h3⋅r 3 r3
r 1⋅ln ln
r2 r2
If in addition the pipe conductivity
k RT =
kp has a very large value 2⋅π⋅k⋅L
k
1 U 1=
()
U 1= r3 Valid for hi, kp y ki
()
r3 r 1⋅ln
r 1⋅ln (Eq. 1.6) with very large values
r2 r r2
+ 1
kI h3⋅r 3 A 1⋅k
1
=A 1⋅U 1 =
()
If in addition the exterior convection RT r3
has a very large value r 1⋅ln
r2
1
U 1=
1 2⋅π⋅r 1⋅L⋅k
()
r3
r 1⋅ln
(Eq. 1.7)
=
()
r2 RT r3
kI r 1⋅ln
r2
1.3 Heat loss per metre of pipe lenght (A3: Exterior surface of insulated pipe, per metre of lenght)
A3 ti t e W
q m
r3 1 r3 r2 r3 r3 1
ln ln
r1 hi k1 r1 k 2 r2 hrc
q= A3 * (ti-te) / ( (r3/r1)*(1/hi) + (r3/k1)*Ln(r2/r1) + (r3/k2)*Ln(r3/r2) + (1/hrc) )
( ) ( )
r2 r3
r 1⋅ln r 1⋅ln
1 r1 r2 r1
+ + +
h1 kp kI h3⋅r 3
dτ =−
d 1⋅ρ⋅C P
4⋅U 1
⋅
dt
( t −t a )
ln
[ ]
t −t a
t i −t a
=−
4⋅U 1
d 1⋅ρ⋅C P
⋅τ
4⋅U 1
re-arranging
t −t a −
d 1⋅ρ⋅C P
⋅τ
=e
dt 4⋅U 1 t i −t a
=− ⋅dτ
t−t a d 1⋅ρ⋅C P 4⋅U 1
− ⋅τ
d 1⋅ρ⋅C P
t −t a =( t i −t a ) e
the integral will be
t
dt 4⋅U 1 τ
∫ t−t =− d ¿ ρ⋅C ¿∫ dτ −
4⋅U 1
⋅τ
t i a 1 P 0 d 1⋅ρ⋅C P
t =t a + ( t i−t a ) e
solving
t 4⋅U 1 t
ln ( t−t a )|t i =− ¿ τ|0
d1 ¿ ρ⋅C P
t 4⋅U 1 t
ln ( t−t a )|t =− ¿ τ|0
i d1 ¿ ρ⋅C P
evaluating
4⋅U 1
ln ( t −t a ) −ln ( t i−t a ) =− ⋅τ
d 1⋅ρ⋅C P
Rev. cjc. 30.01.2014
d 1⋅ρ⋅C P
4⋅U 1 [
τ f =− ⋅ ln ( t f −t a )−ln ( t i−t a ) ]
dt d 1⋅ρ⋅C P
=U 1⋅π⋅d 1⋅L⋅( t −t a ) τf= ⋅[ ln ( t i −t a )− ln ( t f −t a ) ]
dτ 4⋅U 1
(Eq. 1.1)
τf=
d 1⋅ρ⋅C P
4⋅U 1
⋅ ln
[( t i−t a
( t f −t a ) )] (Eq. 1.2)
⋅C P
U1
¿∫
t
dt
tf
( t −t a )
i
4⋅U 1
d 1⋅ρ⋅c P
⋅τ f =ln
t i−t a
t f −t a ( ) (Eq. 1.3)
t
ln ( t−t a )|t f
i
The ASHRAEequation is
1 2⋅π⋅L⋅k
=
()
RT r3
ln
r2
ϑ=ρ⋅C P⋅π⋅
D1 2
2 ()
⋅L⋅RT⋅ln
t i−t a
t f −t a [ ] (Eq. 1.9)
RT =
ln
()
r3
r2 (Eq. 1.8)
with
RT =
1
A1⋅U 1
2⋅π⋅k⋅L
ϑ=
ρ⋅C P⋅π⋅( ) [ ]
D1 2
2
t −t
⋅L⋅ln i a
t f −t a
alid for hi, kp y ki A1⋅U 1
th very large values
ϑ=
ρ⋅C P⋅π⋅D1 ¿L⋅ln
2
[ ]
t i−t a
t f −ta
(Eq. 1.10)
4⋅π⋅D1 ¿ L⋅U 1
ϑ=
ρ⋅C P⋅D1⋅ln
[ ]
t i −t a
t f −t a
(Eq. 1.11)
(The factor "12" from Ashrae equation
4⋅U 1 is due to the use of imperial units)
1
U 3=
() ()
r2 r3
r 3⋅ln r 3⋅ln
r1 r2 1
+ +
kP kI hrc
4⋅U 1
⋅τ
d 1⋅ρ⋅C P
4⋅U 1
⋅τ
⋅ρ⋅C P
4⋅U 1
− ⋅τ
d 1⋅ρ⋅C P
)e
4⋅U 1
− ⋅τ (Eq. 1.12)
d 1⋅ρ⋅C P
e
Rev. cjc. 30.01.2014
1
U 1=
( ) ( )
r2 r3
r 1⋅ln r 1⋅ln
1 r1 r2 r1
+ + +
h1 kp kI h3⋅r 3
1
U 3=
( ) ( )
r2 r3
r 3⋅ln r 3⋅ln
r1 r2 1
+ +
kP kI hrc
q=U 3⋅A3⋅Δt
2.- Average inside water temperature thus, the area below the curve is
[ ]
4⋅U 1
− ⋅τ f
1 d 1⋅ρ⋅C P
The time function of the water temperature A =t a⋅τ f + ( t i−t a )⋅ ⋅ 1−e
4⋅U 1
is given by (Eq. 1.12) d 1⋅ρ⋅C P
4⋅U 1
− ⋅τ
d 1⋅ρ⋅C P The average water temperature tint_ave results from
t =t a + ( t i−t a ) e
A =t int_ ave⋅( τ f −0 )
The area below the temperature curve is A =t int_ ave⋅τ f
τf
A=∫0 t⋅dτ t int_ ave =
A
τf
[ ] [ ]
4⋅U 1 4⋅U 1
τf − ⋅τ 1 − ⋅τ
d ⋅ρ⋅C f
A =∫0 t a + ( t i−t a ) e
d 1⋅ρ⋅C P t int_ ave =t a + ( t i−t a )⋅ ⋅ 1−e 1 P
⋅dτ 4⋅U 1
⋅τ
d 1⋅ρ⋅C P f
[ ( )]
τf τf
A =∫0 t a⋅dτ +∫0 ( ti −t a ) e
d 1⋅ρ⋅C P
⋅dτ 4⋅U 1 t i−t a
⋅τ f = ln
4⋅U 1 d 1⋅ρ⋅C P ( t f −t a )
τf τ f − d ⋅ρ⋅C ⋅τ
A =t a⋅∫0 dτ +( t i −t a)⋅∫ e ⋅dτ
[ ]
1 P 4⋅U 1
0 ( t i−t a ) −
d 1⋅ρ⋅C P
⋅τ f
t int_ ave =t a + ⋅ 1−e
( )
4⋅U 1 (Eq. 2.1)
t i−t a
τ f − d ⋅ρ⋅C ⋅τ ln
A =t a⋅τ|0f + ( t i−ta ) ¿ ∫ e
τ
0
1 P
¿dτ t f −t a
[ ]
4⋅U 1 tf
− ⋅τ
τ −1 d 1⋅ρ⋅C P
A =t a⋅τ|0f + ( t i−t a ) ¿ e 0
4⋅U 1
d1 ¿ ρ⋅C P
[ ]
4⋅U1
− ⋅τ
−1 d 1⋅ρ⋅C P f
A =t a⋅τ f + ( t i−t a )⋅ ⋅e −1
4⋅U 1
d 1⋅ρ⋅C P
A =t a⋅τ f + ( t i−t a )⋅
−1
4⋅U 1
[ − ⋅τ
d ⋅ρ⋅C f
⋅ e 1 P −1 ]
d 1⋅ρ⋅C P
m
a
e
p
t
r
e
TC
a
u
)
°
(
r
tf = 73962 s 8
tint_ave = 7.43 °C
6
]
4⋅U1 ti = 20 ºC 2
− ⋅τ f
d 1⋅ρ⋅C P
tf = 0 ºC
0
tint_m = 10 °C
0 2 4 Tiemp6o(h) 8 10 12
(Eq. 2.1)
Rev. cjc. 30.01.2014
4⋅U 1
− ⋅τ
d 1⋅ρ⋅C P
−t a )⋅e
3.- Cooling down time "t" of water til 0°C, in a pipe with insulation
Cells in yellow are for input values. A cell in orange indicates a result. A cell in magenta indicates an assumed valu
The average exterior surface temperature (t3_ass) must be initially assumed and later its value replaced by the calc
button to perform this task.
The exterior film coefficient "h combined" is a function of the natural and forced convection contributions. This comb
calculated according reference [5]. A different value for this exterior convection (h e) may be input.
Since there is not a relationship available for natural convección inside a tube, a relation for laminar forced convec
the Reynolds number tends to zero and the Nussel for the case of constant surface temperature becomes constant. N
To check this calculation sheet, use sheet "6. NMIC example". To neglect film resistances, use very heigh value fo
For a more finer calculation, the materials and fluid properties should conform to the resulting temperatures (in gen
Material: Carbon steel Wind velocity U1 =
dn = 4 in v= 5 km/h U 1 : Overal heat transfer co
Sch STD Assumed exterior average surface h1 =
Pipe material thermal conductivity temperature d1 =
kp = 60.5 W/(m*K) t3_ass = -4.8 °C (Note 1) r1 =
Absolute rugosity Combined exterior convection coefficient kp =
Rabs = 0.1 mm For function used, see Note 2 d2 =
Insulation thickness and conductivity hcombined = #VALUE! W/(m²*K) r2 =
sI = 50 mm Average value of inside temp. kI =
kI = 0.043 W/(m*K) ti = 5.6 °C d3 =
Water properties tf = 0 °C r3 =
r= 1000 kg/m 3
tave= 2.8 °C h3 =
Cp = 4200 J/(kg*K) Water conductivity U1 =
Fluid initial temperature kw = #VALUE! W/(m*K) Cool down time
ti = 5.6 ºC Nusselt number [1], page 269 The equation for the cool d
Fluid final temperature Nu = 3.66 - is (from sheet 1)
tf =
ta =
0
Exterior ambient temperature
-5 ºC
ºC
h1 =
Inside film coefficient "h1"
h1 = Nu * kw / d
#VALUE! W/(m²*K)
τ=
d 1⋅ρ⋅C P
4⋅U 1
⋅ ln
[(
Pipe_Imp_CS_Dint_dn_sch(dn, SCH)
Pipe inside diameter Pipe outside diameter Insulation outside diamete
d1 = Pipe_Imp_CS_Dint_dn_sch d2 = Pipe_Imp_CS_Dext_dn d3 =
d1 = #VALUE! mm d2 = #VALUE! mm d2 =
A1 = #VALUE! m²/m A2= #VALUE! m²/m Note. For the deduction of used e
Note 1. This temperature will have a value near to the exterior ambient temperature.
Note 2. Combined_Convection_Coefficient_horizontal_tube_Material_dn_Vkmh_tsurf_tamb(Mat, dn, Vkm
Average water temperature tint_ave = ta + ( (ti-ta ) / ln( (ti-ta) / (tf-ta) ) ) * ( 1 - Exp( (-4*U1*tf) / (d1*r*cp)) )
ta = -5 °C Heat loss per meter of pipe Heat loss per square mete
ti = 5.6 °C qL = U1 * A1* (tint_m - te) surface
U1 = #VALUE! W/(m2*K) U1 = #VALUE! W/(m2*K) q3 =
d1 = #VALUE! m A1 = #VALUE! m²/m qL =
r= 1000 kg/m3 tint_m = #VALUE! °C A3=
c p= 4200 J/(kg*K) te = -5 °C q3 =
tf = #VALUE! s qL = #VALUE! W/m
tint_ave = #VALUE! °C
ta = -5 ºC
ti = 5.6 ºC
U1 = #VALUE! W/(m2*K)
Cool d
d1 = #VALUE! m
11
r= 1000 kg/m3 s min h
cp = 4200 J/(kg*K) 0 0 0.00
60 1 0.02
120
10 2 0.03
180 3 0.05
240 4 0.07
300 5 0.08
600
9 10 0.17
900 15 0.25
1200 20 0.33
1500 25 0.42
t (°C)1800 30 0.50
8
2300 38 0.64
2800 47 0.78
3300 55 0.92
3800
7
63 1.06
4300 72 1.19
4800 80 1.33
5300 88 1.47
5800 97 1.61
6
6300 105 1.75
6448 107 1.79
5
0 2 4
dn SCH Material
1/2 5 Carbon steel
3/4 10 Stainless steel
1 20 HDPE PE100
1 1/2 30
2 40
3 60
4 80
5 100
6 120
8 140
10 160
12 STD
14 XS
16 XXS
18
20
22
24
26
28
30
32
34
36
38
40
42
44
46
48
Rev. cjc. 30.01.2014
τ=
d 1⋅ρ⋅C P
4⋅U 1
⋅ ln
[( t i−t a
( t f −t a ) )]
Dt3 = #VALUE!
#VALUE!
°C
Nu =
Gz =
4 6 8 10 12
Time (h)
rced convection inside tubes [1], page 269
Re < 2,300
r laminar flow sufficiently far from the entrance of a tube, where the flow is
drodynamically fully developed and has the caracteristic parabolic velocity profile
Poiseuille, the friction factor is given by
Eq. (4.39)
verage Nusselt
verage value for a tube of length "L"
aminar convection Re≤2300
], page 272, Eq. (4.50)
0. 065⋅Gz
Nu=3 . 66+ 2
3
1+ 0 .04⋅Gz
Gz=Re⋅Pr⋅ () d
L
3.66+(0.065*Gz ) / (1+0.04*Gz^(2/3) )
Re * Pr * (d/L)
usselt for thermally fully developed flow, with constant surface temperature
the wall tremperature is uniform, for example, if steam is condensing on the outside
the tube wall, then sufficiently far downstream of whre heating starts, the flow
ecome thermally fully developed, the shape of the temperature profile is unchanging,
nd the Nusswelt number has a constat value
he thermal entrance length "Leh" is defined as the distance required for the Nusselt
umber to decrease to whithin 5% of its fully developed value "Nu ∞ " [1], page 271.
[5]
Heat and mass transfer. Fundamentasl and aplications
Yunus A. Cengel and Afshin J. Ghajar
Equation (9-41), page 488
Cells in yellow are for input values. A cell in orange indicates a result. A cell in magenta indicates an assumed value.
The average exterior surface temperature (t3_ass) must be initially assumed and later its value replaced by the calculat
button to perform this task.
The exterior film coefficient "h combined" is a function of the natural and forced convection contributions. This combine
calculated according reference [5]. A different value for this exterior convection (h e) may be input.
Since there is not a relationship available for natural convección inside a tube, a relation for laminar forced convection
the Reynolds number tends to zero and the Nussel for the case of constant surface temperature becomes constant. Nu =
ta = -5 ºC
ti = 5.6 ºC
U1 = #VALUE! W/(m2*K)
Cool dow
d1 = #VALUE! m
11
r= 1000 kg/m3 s min h
cp = 4200 J/(kg*K) 0 0 0.00
60 1 0.02
120
10 2 0.03
180 3 0.05
240 4 0.07
300 5 0.08
600
9 10 0.17
900 15 0.25
1200 20 0.33
1500 25 0.42
t (°C)1800
8 30 0.50
2300 38 0.64
2800 47 0.78
3300 55 0.92
3800
7
63 1.06
4300 72 1.19
4800 80 1.33
5300 88 1.47
5800 97 1.61
6
6300 105 1.75
6448 107 1.79
5
0 2 4
6
5
0 2 4
dn SCH Material
1/2 5 Carbon steel
3/4 10 Stainless steel
1 20 HDPE PE100
1 1/2 30
2 40
3 60
4 80
5 100
6 120
8 140
10 160
12 STD
14 XS
16 XXS
18
20
22
24
26
28
30
32
34
36
38
40
42
44
46
48
Rev. cjc. 30.01.2014
64
( 1/h1 + ( r1/kP) *ln(r2/r1) + ( r1/kI) *ln(r3/r2) + (1/h3) *(r1/r3) )^-1 f=
Re
: Overal heat transfer coefficient t = (d1 * r * Cp ) / (4 * U1) * ln((ti - ta) / (tf - ta) )
#VALUE! W/(m2*K) d1 = #VALUE! m
#VALUE! mm r= 1000 kg/m3 Constant wall surface tempe
#VALUE! m Cp = 4200 J/(kg*K)
60.5 W/(m*K) U1 = #VALUE! W/(m2*K) Average Nusselt
#VALUE! mm ti = 5.6 °C Average value for a tube of le
#VALUE! m ta = -5 °C Laminar convection
0.043 W/(m*K) tf = 0 °C [1], page 272, Eq. (4.50)
#VALUE! mm t= #VALUE! s
0. 065⋅Gz
#VALUE! m t= #VALUE! h Nu=3 . 66+ 2
3
#VALUE! W/(m2*K) The assumed temperature t3_ass shall 1+ 0 .04⋅Gz
#VALUE! W/(m2*K) be iterated until its value is close
ool down time
he equation for the cool down time
enough to the calculated value.
t3_ass = -1.0 °C
Gz=Re⋅Pr⋅ ()d
L
(from sheet 1) t3 = #VALUE! °C
τ=
d 1⋅ρ⋅C P
4⋅U 1 [(
⋅ ln
t i−t a
( t f −t a ) )] Dt3 = #VALUE!
#VALUE!
°C
Nu =
Gz =
4 6 8 10 12
Time (h)
4 6 8 10 12
Time (h)
[5]
rced convection inside tubes [1], page 269 Heat and mass transfer. Fu
Yunus A. Cengel and Afshin
Equation (9-41), page 488
Re < 2,300
r laminar flow sufficiently far from the entrance of a tube, where the flow is
drodynamically fully developed and has the caracteristic parabolic velocity profile
Poiseuille, the friction factor is given by [1], page 272, Eq. (4.50)
Nusselt at bulk temperature
Eq. (4.39) consideration of properties v
with temperature
Nub = 3.66+(0.065*Gz ) / (1+0.0
onstant wall surface temperature Gz =
Nub =
verage Nusselt
verage value for a tube of length "L"
aminar convection Re≤2300
], page 272, Eq. (4.50)
0. 065⋅Gz
Nu=3 . 66+ 2
3
1+ 0 .04⋅Gz
Gz=Re⋅Pr⋅ () d
L
3.66+(0.065*Gz ) / (1+0.04*Gz^(2/3) )
Re * Pr * (d/L)
usselt for thermally fully developed flow, with constant surface temperature
the wall tremperature is uniform, for example, if steam is condensing on the outside
the tube wall, then sufficiently far downstream of whre heating starts, the flow
ecome thermally fully developed, the shape of the temperature profile is unchanging,
nd the Nusswelt number has a constat value
he thermal entrance length "Leh" is defined as the distance required for the Nusselt
umber to decrease to whithin 5% of its fully developed value "Nu ∞ " [1], page 271.
eat and mass transfer. Fundamentasl and aplications
unus A. Cengel and Afshin J. Ghajar
quation (9-41), page 488
q= A1⋅U 1⋅Δt
Δt Δt=ti −t a
q=
RT
1
=A ⋅U
RT 1 1
Global heat transfer coefficient (Note 1)
1
U 1= L
() ()
r2 r3
r 1⋅ln r 1⋅ln (Ec. 2.1)
1 r1 r2 r1
+ + +
h1 kp kI h3⋅r 3
And if in addition, the exterior film c
If the interior convection has a very large very large valu (h3 → ∞)
value (h1 → ∞) 1
U 1=
1
( )
r3
U 1= r 1⋅ln
() ()
r2
r2 r3 (Ec. 2.2)
r 1⋅ln r 1⋅ln kI
r1 r2 r
+ + 1
kp kI h3⋅r 3
Global heat transfer coefficient refered to surfa
And if in addition the pipe thermal conductivity
has a large value ( 1/ kp → 0) This case is the one handled by the publication
(Ec. 2.3) Mechanical Insulation Design Guide - Desgin O
1 (Note 2) by the National Mechanical Insulation Commite
U 1= Last updated 07-26-2011
()
r3
r 1⋅ln http://www.wbdg.org/design/midg_design_tf
r2 r1
+
kI h3⋅r 3 This factor considers the case of a a pipe expo
"ta" and surrounded by surfaces at the same te
interchange is considered.
Pipe
dn = in 6
Sch STD
kp = W/(m*K) 60.5
Insulation thickness and conductivity
sI = mm 50
kI = Btu*in / (h * ft²*°F) 0.30
kI = W/(m*K) 0.043
Water properties
r= kg/m3 1000
Cp = J/(kg*K) 4200
Fluid initial temperature
ti = °F 42
ti = °C 5.6
Fluid final temperature
tf = °C 0
Exterior ambient temperature
ta = °F -18
ta = °C -27.8
Inside film coefficient
h1 = W/(m2*K) 10,000
Outside film coefficient
h3 = W/(m2*K) 10,000
Solution (from Sheet 6)
t= 6.9 h
ta
ti d1 d2
d3
sI, kI
sI, kI
sp, kp
/design/midg_design_tffc.php
http://www.wbdg.org/design/midg_design_tffc.php
0.084 m d1 = 0.15408 m
0.043 W/(m*K) r= 1000 kg/m3
268.3 mm Cp = 4200 J/(kg*K)
0.134 m U1 = 1.20 W/(m2*K)
100000 W/(m2*K) ti = 5.6 ºC
1.20 W/(m2*K) ta = -27.8 ºC
tf = 0 ºC
t= 24811 s
t= 6.9 h
lls, interior water film and exterior air films resistances are usually neglected.
onductivity is let with its real value, as it has a very low influence.
no checking without insulation has been done.
Stepwise solution
Combined natural and forced convection
3.- Combined natural and forced convection
mm NuCombined = Abs( (NuNat)^n ± (NuForc)^n )^(1/n)
mm with n ≈ 3÷ 4
(The Nusselt-value is not very sensible to the n-vale
Convection coefficient
he_Conv = NuCombined * k / d
ion coefficient and Re. NuCombined = #VALUE! W/(m2*K)
d, see Note 2 k= #VALUE! W / (m*K)
- d= #VALUE! m
W/(m²*k) he_Conv = #VALUE! W/(m2*K)
- Combined convection
Eq. (9-41), [5], page 487
1
n, Vkmh, ts, te) Nu combined =( Nuforced + Nunatural n
n n
)
dn, Vkmh, ts, te)
Mat, dn, Vkmh, ts, te) With the exponent "n" with a value in the
range of 3 and 4. Used value, n = 4
Microsoft Editor de
[5] ecuaciones 3.0
273.15 K
g= 9.80665 m/s²
flow in smooth pipes.
and Boelter.[2] Eq. 6-4, on page 274
Forced_Laminar_Interior_Array_Water_Material_dn_SchOrPN_Lm_twall_tin_Qm3h
Forced_Laminar_Interior_Nusselt_Water_Material_dn_SchOrPN_Lm_twall_tin_Qm3h
Forced_Laminar_Interior_Coefficient_Water_Material_dn_SchOrPN_Lm_twall_tin_Qm3h
Module: e1_For_Int_Laminar_Water
Re = #VALUE! -
−h⋅π⋅d⋅L Lh_Lam = #VALUE! m
A=
m⋅Cp
t out =t wall −( t wall −t in )⋅Exp( A ) - Thermally fully developed Eq. (8-12)
Lt_Lam = Lh_Lam * Pr
Temperature of fluid leving the tube Lh_Lam = #VALUE! m
tout = twall -(twall - tin ) * Exp(A) Pr = 3.0 -
A= -(h*p*d*L ) / ( m*Cp ) Lt_Lam = #VALUE! m
h= #VALUE! W/(m²*k)
d= #VALUE! m Tube length
L= 20 m L= 20.0 m
m= #VALUE! kg/s
Cp = #VALUE! J/(kg*K) Since #VALUE!
A= #VALUE! The laminar flow is
tout = twall -(twall - tin ) * Exp(A) #VALUE!
twall = 100 °C
tin = 15 °C Since #VALUE!
A= #VALUE! The laminar flow is
tout = #VALUE! °C #VALUE!
Note.
Nusselt numbers and thus convection heat transfer
coefficients are much higher in the entrance region.
Microsoft Editor de
ecuaciones 3.0
Datos
Cañería dn Tiempo
t = de enfriamiento ha
dn = Variable in in diámetro
h de la cañería, para
Sch STD 2 0.5
Rabs = 0.1 mm 3 12 1.1
kp = 60.5 W/(m*K) 4 1.8Cañería sin aislamiento
Espesor y conductividad de la aislación 6 4 (Pipe without insulation)
10
sI = 0.00001 mm (1) 8 6.8
kI = 0.035 W/(m*K) 10 10.5
Propiedades del agua 8
r= 1000 kg/m 3
Tiempo enfriamiento
Cp = 4200 J/(kg*K)
6 (h)
Temperatura inicial del agua
ti = 10 ºC
Temperatura final del agua 4
tf = 0 ºC
Temperatura del ambiente exterior
2
ta = -5 ºC
Velocidad del viento
v= 162 km/h 0
Flujo másico de agua 0 2 4 6
Diametro nominal
m= 1 kg/min (2)
Los puntos de la curva se calcularon con la hoja de cálculo "4.- Pipe without insulation"
Las ecuaciones utilizadas se muestran en la hoja de cálculo "1.- Cool down of a pipe"
Rev. cjc. 30.01.2014
sin aislamiento
ut insulation)
4 6 8 10 12
Diametro nominal cañería (in)
References
Introducing the Mechanical Insulation Design Guide From the National Insulation Association and the National Ins
http://www.insulation.org/articles/article.cfm?id=PR080101
[4] http://www.tak2000.com/data/handbookx.pdf
[6] http://www.youtube.com/watch?v=afv2p1gJKbA
[12] http://wwwme.nchu.edu.tw/Enter/html/lab/lab516/Heat%20Transfer/chapter_8.pdf
[17] Advanced Heat and Mass Transfer by Amir Faghri, Yuwen Zhang, and John R. Howell
5.3 Hydrodynamically and Thermally Fully Developed Laminar Flow
http://www.thermalfluidscentral.org/e-resources/download.php?id=??
5.4 Hydrodynamically fully developed and Thermally Fully Developing Lamina
http://www.thermalfluidscentral.org/e-resources/download.php?id=86
Tau
Microsoft Equation
3.0
Q
Microsoft Equation
3.0
U3 Microsoft Equation
3.0
e+tuberias+ante+la+congelacion.+Armaflex.+Inf.+Tec.+N%C2%B0+9&pf=p&output=search&sclient=psy-ab&oq=Proteccion+de+tuberias+a
te%20la%20congelacion.%20armaflex.%20inf.%20tec.%20n%C2%B0%209&cp=34&gs_id=i&xhr=t&q=Mechanical%20Insulation%20Design
U1
Microsoft Equation
3.0
q=Proteccion+de+tuberias+ante+la+congelacion.+Armaflex.+Inf.+Tec.+N%C2%B0+9&gs_l=&pbx=1&bav=on.2,or.r_gc.r_pw.r_qf.,cf.osb&fp
nical%20Insulation%20Design%20Guide&pf=p&sclient=psy-ab&oq=Mechanical+Insulation+Design+Guide&gs_l=&pbx=1&bav=on.2,or.r_gc
2,or.r_gc.r_pw.r_qf.,cf.osb&fp=2789c0abca5a439a&biw=1024&bih=674
s_l=&pbx=1&bav=on.2,or.r_gc.r_pw.r_qf.,cf.osb&fp=2789c0abca5a439a&biw=1024&bih=674&bs=1
Heat transfered trough a pipe and the insulation
( t 1 −t 2 )
Q= A 1⋅2⋅k P⋅ d3
ln 2
r
r1 [ ] (b)
kP
kI
Heat transfer by conduction trough the pipe material
( t 2 −t 3 )
Q= A 2⋅2⋅k I⋅
ln
r3
r2 [ ] (c)
Adding the four equations (e.) to (h)
[ ]
Heat transfer by convection and radiation from the r2
exterior surface to the ambient ln ln
1 r1
t i−t e=Q⋅ + Q⋅ +Q⋅
Q= A 3⋅he⋅(t 3 −t e ) (d) A 1⋅hi A 1⋅2⋅k P A
[ [] [
The temperature differential in each case are
r2 r3
1 ln ln
(ti −t 1 )=Q⋅
A1⋅hi
(.e)
1 r1 r2
t i−t e=Q⋅ + +
A 1⋅hi A 1⋅2⋅k P A 2⋅2⋅
( t 1 −t 2 )=Q⋅
ln
[ ]
r2
r1 (f)
A 1⋅2⋅k P t i −t e
Q=
( t2 −t 3 )=Q⋅
ln
[]
r3
r2
A 2⋅2⋅k I
(g)
1
+
ln
[ ] [ ]
r2
r1
+
ln
+
1
r3
r2
A 1⋅hi A 1⋅2⋅k P A 2⋅2⋅k I A 3⋅h
1
(t 3 −t e )=Q⋅ (h)
A 3⋅he
Rev. cjc. 30.01.2014
t i −t e
Q= A 3⋅
ti A3
+ A 3⋅
ln
r1 [ ]
r2
+ A 3⋅
r2
ln
A
+ 3
[]
r3
[ [] [] ]
r2 r3 1
ln ln U 3=
Q⋅
1
+ + +
r11
A 1⋅hi A 1⋅2⋅k P A 2⋅2⋅k I A 3⋅he
r2
d3 1 d 3
⋅ + ⋅
ln
d1 [ ]
d2
d3
+ ⋅
d2
+
1
ln
[ ]
d3
d1 hi d 1 2⋅k P d 2 2⋅k I he
t i −t e
+
ln
[ ] [ ]
r2
r1
+
ln
r3
r2
+
1
Q=π⋅d 3⋅L⋅U 3⋅( t i−t e )
Q
q= =π⋅d 3⋅U 3⋅( t i−t e ) [ W / m ]
hi A 1⋅2⋅k P A 2⋅2⋅k I A 3⋅he L
Multiplying numerator and denominator by "A1"
t −t
Q= A 1⋅ i e
3
ln
⋅
[]
r2
r3
A
+ 3
A1 1 A1
⋅ + ⋅
r1
ln
A
+ 1⋅
[ ]
r2
r2
ln
A 1
+ 1⋅
[ ]
r3
⋅
ln
[ ]
r3
r2
+
A3 1 d1 1 d1
⋅ + ⋅
ln
d2
d1 d1
+ ⋅
[ ]
ln
d3
d2 d1 1
+ ⋅
[ ]
2⋅k I A 3 ¿ he d 1 hi d 1 2⋅k P d 2 2⋅k I d3 he
t −t
Q= A 1⋅ i e
ln
[ ]
d3
d2
+
1 1
+
ln
d2
d1 [ ]
d1
+ ⋅
ln
d3
d2 d1 1
+ ⋅
[ ]
2⋅k I he hi 2⋅k P d 2 2⋅k I d3 he
Q= A 1⋅U 1⋅( t i −t e )
1
U 1=
1
ln
+
[ ]
d2
d1 d1
+ ⋅
ln
d2 [ ]
d3
d1 1
+ ⋅
hi 2⋅k P d 2 2⋅k I d 3 he
Microsoft Editor de
ecuaciones 3.0
Microsoft Equation
3.0
Microsoft Equation
3.0
Microsoft Equation
3.0
Microsoft Editor de
ecuaciones 3.0
Microsoft Editor de
ecuaciones 3.0