Universidade Eduardo Mondlane Faculdade de Engenharia Departamento de Engenharia Química Curso de Engenharia Química
Universidade Eduardo Mondlane Faculdade de Engenharia Departamento de Engenharia Química Curso de Engenharia Química
Universidade Eduardo Mondlane Faculdade de Engenharia Departamento de Engenharia Química Curso de Engenharia Química
FACULDADE DE ENGENHARIA
Discente:
Sabela, Filipe Diamantino
Docentes:
Prof. Doutor Miguel M. Uamusse, Eng.
Prof. Doutor Alberto Tsamba, Eng.
𝑘 = 25𝑊/𝑚. ℃
Resolução
1 𝜕 𝜕𝑇 𝜕𝑇
. (𝐴. 𝑘. ) + 𝑔 = 𝑐𝑝. 𝜌. (𝑥, 𝑡)
𝐴 𝜕𝑥 𝜕𝑥 𝜕𝑡
𝜕𝑇
(𝑥, 𝑡) = 0 ; 𝑔 = 0
𝜕𝑡
1 𝜕 𝜕𝑇
. (𝐴. 𝑘. )=0
𝐴 𝜕𝑥 𝜕𝑥
𝜕 𝜕𝑇
( )=0
𝜕𝑥 𝜕𝑥
𝜕𝑇
𝐼1 : = 𝑐1
𝜕𝑥
𝐼2 : 𝑇(𝑥) = 𝑐1 𝑥 + 𝑐2
Condições de contorno:
1ª:
𝑑𝑇
𝑞 = −𝑘.
𝑑𝑥
𝑑𝑇
𝑞 ∗ 𝐴 = −𝑘 ∗ 𝐴 ∗
𝑑𝑥
𝑑𝑇
𝑄 = −𝑘 ∗ 𝐴 ∗
𝑑𝑥
𝑑𝑇
𝑠𝑒𝑛𝑑𝑜 = 𝑐1
𝑑𝑥
𝑄 = −𝑘 ∗ 𝐴 ∗ 𝑐1
1500 = −25 ∗ 12 ∗ 10−3 ∗ 𝑐1
𝑐1 = −5000
2ª
𝑞𝑢𝑎𝑛𝑑𝑜: 𝑥 = 0.8 ∗ 10−2 𝑚 ; 𝑇2 = 90℃
𝑇(𝑥) = −5000 ∗ 𝑥 + 𝑐2
𝑐2 = 130℃
𝑇(𝑥) = −5000 ∗ 𝑥 + 130
𝑘 = 2.4𝑊/𝑚. ℃ ; ℎ∞ = 24𝑊/𝑚2 . ℃ ; 𝐴 = 20 𝑚2
Resolução
Condições de contorno:
1ª
𝑑𝑇
𝑛𝑜 𝑝𝑜𝑛𝑡𝑜 𝑥 = 0 ; 𝑞 = −𝑘. 𝑜𝑛𝑑𝑒 𝑞 = ℎ∞ (𝑇𝑤 − 𝑇∞ )
𝑑𝑥
𝑇𝑤 − 𝑇0
ℎ∞ (𝑇𝑤 − 𝑇∞ ) = −𝑘.
𝐿
𝑇𝑤 − 80
24 ∗ (𝑇𝑤 − 15) = −2.4 ∗ → 𝑇𝑤 = 28℃
0.4
𝑑𝑇
ℎ∞ (𝑇𝑤 − 𝑇∞ ) = −𝑘 ∗
𝑑𝑥
−𝑘 ∗ 𝑐1 = ℎ∞ (𝑇𝑤 − 𝑇∞ )
ℎ∞ ∗ (𝑇𝑤 − 𝑇∞ ) 24 ∗ (28 − 15)
𝑐1 = − =
𝑘 2.4
𝑐1 = −130
2ª
𝑇(0) = 80℃
𝑇(𝑥) = 𝑐1 𝑥 + 𝑐2 → 𝑇(0) = 𝑐2
𝑐2 = 80℃
Equação da distribuição da temperatura
𝑇(𝑥) = −130. 𝑥 + 80
Resolução
𝑑𝑇 𝑇1 − 𝑇2
𝑞 = −𝑘 ∗ → 𝑞 − 𝑘. 𝑜𝑛𝑑𝑒: 𝐿
𝑑𝑥 ∆𝑥
(30 − 70)
𝑞 − 70. → 𝑞 = 140𝑊/𝑚2
0.02
𝑄 = 𝑞 ∗ 𝐴 → 𝑄 = 140 ∗ 0.6
𝑄 = 84𝑊
Resolução
1 𝜕 𝜕𝑇 𝜕𝑇
Partindo da eq. geral (para uma esfera) : 2
. (𝑟 2 . 𝑘. ) + 𝑔 = 𝑐𝑝. 𝜌. (𝑟, 𝑡)
𝑟 𝜕𝑟 𝜕𝑟 𝜕𝑡
𝜕𝑇
𝑔=0; (𝑟, 𝑡) ; 𝑘 − 𝑐𝑜𝑛𝑡
𝜕𝑡
𝜕 𝜕𝑇
(𝑟 2 . ) = 0
𝜕𝑥 𝜕𝑟
𝜕𝑇 𝜕𝑇 𝑐1
𝐼1 : 𝑟 2 . = 𝑐1 → = 2
𝜕𝑥 𝜕𝑥 𝑟
𝑐1
𝐼2 : = 𝑇(𝑟) = + 𝑐2
𝑟
Condições de contorno:
1ª
𝑄𝑢𝑎𝑛𝑑𝑜 𝑟 = 𝑟1 = 2𝑚 ; 𝑇 = 0℃
𝑐1 𝑐1
𝑇(0) = + 𝑐2 → + 𝑐2 = 0
2 2
2ª
𝑑𝑇
𝑛𝑜 𝑝𝑜𝑛𝑡𝑜 𝑟 = 2.1𝑚 ; 𝑞 = −𝑘. 𝑜𝑛𝑑𝑒 𝑞 = ℎ∞ (𝑇∞ − 𝑇𝑤 )
𝑑𝑟
𝑇𝑤𝑖 − 𝑇𝑤
ℎ∞ (𝑇𝑤 − 𝑇∞ ) = −𝑘 ∗
∆𝑟
𝑇𝑤 − 𝑡𝑒𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑑𝑎 𝑝𝑎𝑟𝑒𝑑𝑒 𝑒𝑥𝑡𝑒𝑟𝑖𝑜𝑟
𝑇𝑤𝑖 − 𝑡𝑒𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑑𝑎 𝑝𝑎𝑟𝑒𝑑𝑒 𝑖𝑛𝑡𝑒𝑟𝑖𝑜𝑟
0 − 𝑇𝑤
18 ∗ (𝑇𝑤 − 25) = −30 ∗
0.1
𝑇𝑤 = 1.42℃
𝑐1
𝑇(𝑟) = + 𝑐2 ; 𝑟 = 2.1𝑚; 𝑇(2.1) = 1.42℃
𝑟
𝑐1
− + 𝑐2 = 1.42
2.1
𝑐1
− + 𝑐2 = 0 𝑐1 = 59.22
{ 𝑐 2 →{ .
1
− + 𝑐2 = 1.49 𝑐2 = 29.61
2.1
59.22
𝑇(𝑟) = − + 29.61
𝑟
Resolução
Condições de contorno:
1ª
𝑑𝑇
𝑞0 = −𝑘.
𝐸𝑚 𝑥 = 0 { 𝑑𝑥
𝑇𝑤 |0 = 80℃
𝑇(𝐿) = −4℃
Resolução
𝑑𝑇
𝑞 = −𝑘.
𝑑𝑟
𝑑𝑇 𝑔0 𝑟0 2 2𝑟 𝑑𝑇 𝑔0 . 2𝑟
= . 2 → =
𝑑𝑟 𝑘 𝑟0 𝑑𝑟 𝑘
𝑔0 . 2𝑟
𝑞 = −𝑘. → 𝑞 = −𝑔0 . 2𝑟 → 𝑞 = −35.2.4
𝑘
𝑞 = 280𝑊/𝑐𝑚2
𝜕𝑇 𝑔0 . 𝑟 2 𝜕𝑇 𝑔0 . 𝑟 𝑐1
𝐼1 : 𝑟 =− + 𝑐1 → =− +
𝜕𝑥 2. 𝑘 𝜕𝑥 2. 𝑘 𝑟
𝑔0 . 𝑟 2
𝐼2 : 𝑇(𝑟) = − + 𝑐1 . ln 𝑟 + 𝑐2
4. 𝑘
𝑐𝑜𝑛𝑑𝑖çã𝑜 𝑒𝑠𝑝𝑒𝑐𝑖𝑎𝑙:
𝜕𝑇
= 0 ; 𝑐1 = 0
𝜕𝑥
𝑔0 . 𝑟 2
𝑇(𝑟) =− + 𝑐2
4. 𝑘
𝐶𝑜𝑛𝑑𝑖çõ𝑒𝑠 𝑑𝑒 𝑐𝑜𝑛𝑡𝑜𝑟𝑛𝑜:
𝑄𝑢𝑎𝑛𝑑𝑜: 𝑟 = 𝑏 ; 𝑇 = 𝑇2
𝑔0 . 𝑏2
𝑇2 = − + 𝑐2
4. 𝑘
𝑔0 . 𝑏2
𝑐2 = 𝑇2 +
4. 𝑘
𝑔0 . 𝑟 2 𝑔0 . 𝑏2
𝑇(𝑟) = − + 𝑇2 +
4. 𝑘 4. 𝑘
𝑔0
𝑇(𝑟) = . (𝑏2 − 𝑟 2 ) + 𝑇2
4. 𝑘
8. v
Dados e Esquema:
Pressupostos
∗ 𝑐𝑜𝑛𝑑𝑖çã𝑜 𝑢𝑛𝑖𝑑𝑖𝑚𝑒𝑛𝑠𝑖𝑜𝑛𝑎𝑙
∗ 𝑒𝑠𝑡𝑎𝑑𝑜 𝑒𝑠𝑡𝑎𝑐𝑖𝑜𝑛á𝑟𝑖𝑜
∗ 𝑘 − 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒
∗ 𝑔 − 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒
Resolução
1 𝜕 𝜕𝑇 𝜕𝑇
Partindo da eq. geral(para um cilindro): . (𝑟. 𝑘. ) + 𝑔 = 𝑐𝑝. 𝜌. (𝑟, 𝑡)
𝑟 𝜕𝑟 𝜕𝑟 𝜕𝑡
𝜕𝑇
𝑐𝑜𝑚𝑜 𝑒𝑠𝑡á 𝑛𝑜 𝑒𝑠𝑡𝑎𝑑𝑜 𝑒𝑠𝑡𝑎𝑐𝑖𝑜𝑛á𝑟𝑖𝑜: (𝑟, 𝑡) = 0
𝜕𝑡
𝑘 − 𝑐𝑜𝑛𝑠𝑡, 𝑑𝑖𝑣𝑖𝑑𝑖𝑛𝑑𝑜 𝑡𝑜𝑑𝑎 𝑎 𝑒𝑥𝑝𝑟𝑒𝑠𝑠ã𝑜 𝑝𝑜𝑟 𝑘
𝜕 𝜕𝑇 𝑔. 𝑟
(𝑟. )=−
𝜕𝑟 𝜕𝑟 𝑘
𝜕𝑇 𝑔0 . 𝑟 2 𝜕𝑇 𝑔0 . 𝑟 𝑐1
𝐼1 : 𝑟 =− + 𝑐1 → =− +
𝜕𝑥 2. 𝑘 𝜕𝑥 2. 𝑘 𝑟
𝑔0 . 𝑟 2
𝐼2 : 𝑇(𝑟) =− + 𝑐1 . ln 𝑟 + 𝑐2
4. 𝑘
𝑐𝑜𝑛𝑑𝑖çã𝑜 𝑒𝑠𝑝𝑒𝑐𝑖𝑎𝑙:
𝜕𝑇
= 0 ; 𝑐1 = 0
𝜕𝑥
𝑔0 . 𝑟 2
𝑇(𝑟) = − + 𝑐2
4. 𝑘
𝐶𝑜𝑛𝑑𝑖çõ𝑒𝑠 𝑑𝑒 𝑐𝑜𝑛𝑡𝑜𝑟𝑛𝑜:
𝑄𝑢𝑎𝑛𝑑𝑜: 𝑟 = 𝑏 ; 𝑇 = 0
𝑔0 . 𝑏2
0=− + 𝑐2
4. 𝑘
𝑔0 . 𝑏2
𝑐2 =
4. 𝑘
𝑔0 . 𝑟 2 𝑔0 . 𝑏2
𝑇(𝑟) = − +
4. 𝑘 4. 𝑘
𝑔0
𝑇(𝑟) = . (𝑏 2 − 𝑟 2 )
4. 𝑘
Temperatura no centroide
𝑔0 . 𝑏2
𝑇(0) =
4. 𝑘
9. aa
Dados e Esquema:
𝑟2
𝑔(𝑟) = 𝑔0 . (1 − )
𝑏2
Resolução
1 𝜕 𝜕𝑇 𝜕𝑇
. (𝑟. 𝑘. ) + 𝑔 = 𝑐𝑝. 𝜌. (2)
𝑟 𝜕𝑟 𝜕𝑟 𝜕𝑡
𝜕𝑇
𝑐𝑜𝑚𝑜 𝑒𝑠𝑡á 𝑛𝑜 𝑒𝑠𝑡𝑎𝑑𝑜 𝑒𝑠𝑡𝑎𝑐𝑖𝑜𝑛á𝑟𝑖𝑜: =0
𝜕𝑡
1 𝜕 𝜕𝑇
. (𝑟. 𝑘. ) + 𝑔 = 0
𝑟 𝜕𝑟 𝜕𝑟
𝜕 𝜕𝑇 𝑔. 𝑟
(𝑟. )=−
𝜕𝑟 𝜕𝑟 𝑘
𝜕 𝜕𝑇 𝑟 𝑟2
(𝑟. ) = − . 𝑔0 . (1 − 2 )
𝜕𝑟 𝜕𝑟 𝑘 𝑏
𝜕𝑇 𝑔0 . 𝑟 2 𝑔0 . 𝑟 3 𝜕𝑇 𝑔0 . 𝑟 2 𝑔0 . 𝑟 𝑐1
𝐼1 : 𝑟 = −( − ) + 𝑐1 → =( − )+
𝜕𝑥 2. 𝑘 3𝑏. 𝑘 𝜕𝑥 3𝑏. 𝑘 2. 𝑘 𝑟
𝑔0 . 𝑟 3 𝑔0 . 𝑟 2
𝐼2 : 𝑇(𝑟) = ( − ) + 𝑐1 . ln 𝑟 + 𝑐2
9𝑏. 𝑘 4. 𝑘
𝑔0 . 𝑏2 𝑔0 . 𝑏3
𝑐2 = −
4. 𝑘 9𝑏. 𝑘
𝑔0 . 𝑟 3 𝑔0 . 𝑟 2 𝑔0 . 𝑏2 𝑔0 . 𝑏3
𝑇(𝑟) = ( − )+ −
9𝑏. 𝑘 4. 𝑘 4. 𝑘 9𝑏. 𝑘
𝑔0 2 𝑔
𝑇(𝑟) = (𝑏 − 𝑟 2 ) + 0 (𝑟 3 − 𝑏3 )
4. 𝑘 9𝑏. 𝑘
Temperatura no centroide
𝑔0 . 𝑏2 𝑔0 . 𝑏3
𝑇(0) = −
4. 𝑘 9𝑏. 𝑘
10. Determinação da equação de queda de temperatura entre o eixo e a superfície de um
elemento combustível cilíndrico, de um reactor nuclear resfriado a gás, com geração
de calor dentro do elemento combustível.
Dados e Esquema:
𝑟2
𝑔(𝑟) = 𝑔0 . (1 − ) ; 𝑟 = 𝑏 → 𝑇 = 𝑇0
𝑏2
Resolução
1 𝜕 𝜕𝑇 𝜕𝑇
Partindo da eq. geral(para um cilindro): . (𝑟. 𝑘. ) + 𝑔 = 𝑐𝑝. 𝜌. (𝑟, 𝑡)
𝑟 𝜕𝑟 𝜕𝑟 𝜕𝑡
𝜕𝑇
(𝑟, 𝑡) = 0 ; 𝑘 − 𝑐𝑠𝑡
𝜕𝑡
𝜕 𝜕𝑇 𝑔. 𝑟
(𝑟. ) = −
𝜕𝑥 𝜕𝑥 𝑘
𝜕 𝜕𝑇 𝑔0 . 𝑟 𝑟2
(𝑟. ) = − . (1 − 2 )
𝜕𝑥 𝜕𝑥 𝑘 𝑏
𝜕 𝜕𝑇 𝑔0 𝑟3
(𝑟. ) = − . (𝑟 − 2 )
𝜕𝑥 𝜕𝑥 𝑘 𝑏
𝜕𝑇 𝑔0 𝑟 2 𝑟4
𝐼1 : 𝑟 = − .( − ) + 𝑐1
𝜕𝑥 𝑘 2 4. 𝑏2
𝜕𝑇 𝑔0 𝑟3 𝑐1
=− . (𝑟 − ) +
𝜕𝑥 2. 𝑘 2. 𝑏2 𝑟
𝑔0 𝑟 2 𝑟4
𝐼2 : 𝑇(𝑟) =− .( − ) + 𝑐1 . ln 𝑥 + 𝑐2
2. 𝑘 2 4.2. 𝑏2
𝑔0 𝑟4
𝑇(𝑟) = − . (𝑟 2 − ) + 𝑐1 . ln 𝑥 + 𝑐2
4. 𝑘 4. 𝑏2
𝐶𝑜𝑛𝑑𝑖çã𝑜 𝑒𝑠𝑝𝑒𝑐𝑖𝑎𝑙:
𝜕𝑇
= 0 ; 𝑐1 = 0
𝜕𝑥
𝑔0 𝑟4
𝑇(𝑟) = − . (𝑟 2 − ) + 𝑐2
4. 𝑘 4. 𝑏2
𝐶𝑜𝑛𝑑𝑖çõ𝑒𝑠 𝑑𝑒 𝑐𝑜𝑛𝑡𝑜𝑟𝑛𝑜:
𝑄𝑢𝑎𝑛𝑑𝑜: 𝑟 = 𝑏; 𝑇 = 𝑇0
𝑔0 𝑏4
𝑇0 = − . (𝑏2 − ) + 𝑐2
4. 𝑘 4. 𝑏2
𝑔0 𝑏2 3. 𝑔0 . 𝑏2
𝑐2 = 𝑇0 + . (𝑏2 − ) → 𝑐2 = 𝑇0 +
4. 𝑘 4 16. 𝑘
𝑔0 𝑟4 3. 𝑔0 . 𝑏2
𝑇(𝑟) = − . (𝑟 2 − ) + 𝑇0 +
4. 𝑘 4. 𝑏2 16. 𝑘
𝑔0 𝑟4 3. 𝑏2
𝑇(𝑟) = . (−𝑟 2 + + ) + 𝑇0
4. 𝑘 4. 𝑏2 4
∆𝑇 =?
∆𝑇 = 𝑇(0) − 𝑇(𝑏)
𝑔0 𝑟4 3. 𝑏2
𝑇(𝑟) = . (−𝑟 2 + + ) + 𝑇0
4. 𝑘 4. 𝑏2 4
𝑔0 04 3. 𝑏2
𝑇(0) = . (−02 + + ) + 𝑇0
4. 𝑘 4. 𝑏2 4
𝑔0 3. 𝑏2
𝑇(0) = .( ) + 𝑇0
4. 𝑘 4
𝑇(𝑏) = 𝑇0
𝑔0 3. 𝑏2
∆𝑇 = .( ) + 𝑇0 − 𝑇0
4. 𝑘 4
𝑔0 3. 𝑏2
∆𝑇 = .( )
4. 𝑘 4
b)
𝑏 = 2𝑐𝑚 = 0.02𝑚 ; 𝑔0 = 4. 107 𝑊/𝑚3 ; 𝑘 = 10 𝑊/𝑚. ℃
4. 107 3. 0.022
∆𝑇 = .( )
4.10 4
∆𝑇 = 300 ℃