Aula 12 - Sistemas Térmicos

Prévia do material em texto

Modelagem Matemática de Sistemas Térmicos 1
1 INTRODUÇÃO
Sistemas térmicos são sistemas nos quais estão envolvidos o armazenamento e o fluxo de calor
por condução, convecção ou radiação. A rigor, sempre estão envolvidas simultaneamente as três
formas de transferência de calor. Entretanto, na prática, tem-se em geral a preponderância de
uma forma sobre as demais ou então a preponderância de duas formas sobre a terceira, o que é
mais comum. Exemplos clássicos de sistemas térmicos são o sistema de arrefecimento do motor
de um automóvel, o refrigerador doméstico, o sistema de condicionamento de ar de um
escritório, etc.
Há três maneiras pelas quais o calor pode fluir de uma substância para outra: condução,
convecção e radiação. Na transferência de calor por condução ou convecção, o fluxo de calor q,
em kcal/s, é dado por
(1) θ∆= Kq
onde θ∆ = diferença de temperatura, em K
K = coeficiente de proporcionalidade, em kcal/s.K, dado por
(2) 
X
kAK ∆= na condução
(3) hAK = na convecção
onde k = condutividade térmica, em kcal/m.s.K
A = área normal ao fluxo de calor, m2
∆X = espessura do condutor, em m
h = coeficiente de transferência de calor por convecção, kcal/m2.s.K
Na transferência de calor por radiação, o fluxo de calor q, em kcal/s, é dado por
(4) )(Kq 42
4
1r θ−θ=
onde Kr = coef 4 de,
 tam
θ1 = temp
θ1 = temp
12 Modelagem Matemática de Sistemas
Térmicos
iciente de proporcionalidade, em kcal/s.K , que depende da emissivida
anho e configuração da superfície
eratura do emissor, em K
eratura do emissor, em K
1
Modelagem Matemática de Sistemas Térmicos
2
2
Na modelagem que será feita a seguir, serão consideradas apenas as transferências de calor por
condução e por convecção, desprezando-se o efeito da radiação.
2 VARIÁVEIS TÉRMICAS
As variáveis usadas para descrever o comportamento de um sistema térmico são:
θ = temperatura em Kelvins [K] o C = K - 273,15
q = fluxo de calor em Watts [W] 1 W = 1 J/s
As temperaturas em vários pontos de um corpo variam com a localização, o que significa que o
sistema térmico é inerentemente um sistema com parâmetros distribuídos. Em conseqüência, os
modelos matemáticos são constituídos por equações diferenciais parciais, pois as propriedades
são distribuídas e não concentradas. Na modelagem e na análise, entretanto, para simplificar o
problema, é conveniente admitir que um sistema térmico possa ser representado por um modelo
de parâmetros concentrados, no qual as substâncias que são caracterizadas pela resistência ao
fluxo de calor têm capacitância térmica desprezível e que as substâncias que são representadas
pela capacitância térmica têm resistência desprezível ao fluxo de calor. Isso nos conduzirá a
modelos regidos por equações diferenciais ordinárias, com as suas já conhecidas vantagens.
3 NÚMERO DE BIOT
Existe um parâmetro adimensional, denominado Número de Biot, que serve de critério para
definir se um sistema térmico pode ser admitido como de parâmetros concentrados. Ele é
definido como
(5) 
k
hLBi c=
onde h e k já foram definidos e onde Lc é o comprimento característico do sólido, definido por
(6) 
s
c A
VL =
onde V = volume do sólido, em m3
As = é a área da superfície de contato entre sólido e fluido, no caso de transferência de
 calor por convecção, em m2
 Evidentemente, Lc depende da forma do sólido.
Assim, para esferas de raio r, temos:
Modelagem Matemática de Sistemas Térmicos
3
3
r
3
1
r4
r
3
4
L 2
3
c =π
π
=
Para cilindros maciços de raio r e comprimento L:
 
)Lr(2
rL
r2rL2
LrL 2
2
c +=π+π
π=
Para cubos de aresta L:
 
6
L
L6
LL 2
3
c ==
Um critério aceitável para que a temmperatura no interior de um sólido não varie com a
localização é que
(7) 1,0
k
hLBi c <=
4 VARIÁVEIS INCREMENTAIS
Para a maioria dos sistemas térmicos existe uma condição de equilíbrio que define o ponto de
operação do sistema. Assim, podemos definir uma temperatura incremental e um fluxo de calor
incremental como
(8) 
−θ−θ=θ )t()t(^
(9) 
−−= q)t(q)t(q^
onde 
-
q e 
−θ são os valores das variáveis no ponto de operação.
5 CAPACITÂNCIA TÉRMICA
Existe uma relação entre a temperatura de um corpo físico e o calor nele armazenado. Não
havendo mudança de fase e desde que a faixa de temperaturas não seja excessiva, tal relação
pode ser considerada linear. Assim, sendo qi(t) o fluxo de calor que entra em um corpo e qo(t) o
fluxo de calor que sai do mesmo corpo, o calor líquido (no sentido contábil) armazenado no corpo
entre dois instantes de tempo t0 e t é dado por
 λλ−λ∫ d)](q)(q[ ott i0
onde λ é uma variável muda usada na integração.
Modelagem Matemática de Sistemas Térmicos
4
4
Vamos assumir que o calor armazenado durante esse intervalo de tempo seja igual a uma certa
constante C multiplicada pela variação de temperatura, ou seja
 )]t()t([Cd)](q)(q[ 0o
t
t i0
θ−θ=λλ−λ∫
onde θ(t0) é a temperatura do corpo no instante de referência t0. Podemos rescrever a equação
acima como
(10) λλ−λ+θ=θ ∫ d)](q)(q[C1)t()t( ott i0 0
onde a constante C é definida como a capacitância térmica do corpo, dada em [J/K]. Para um
corpo de massa M e calor específico c, a capacitância térmica é dada por C = Mc, para M em [kg]
e c em [J/kg.K].
Diferenciando a eq. (10), obtemos
(11) )]t(q)t(q[
C
1)t( oi
. −=θ
equação que é muito usada quando o sistema é modelado no espaço de estados.
6 RESISTÊNCIA TÉRMICA
No caso de transferência de calor por condução, a Lei de Fourier estabelece que o fluxo de calor
q(t) entre dois corpos com temperatura θ1(t) > θ2(t), separados por um meio condutor, é dado por
d
)t()t(
A)t(q 21
θ−θα=
onde α = condutividade térmica do material condutor [J/m.s.K] ou [W/m.K] (tabelada)
A = área normal ao fluxo de calor [m2]
d = espessura do condutor [m]
Podemos rescrever a equação acima como
(12) )]t()t([
R
1)t(q 21 θ−θ=
onde R é definida como a resistência térmica e é função do material e das dimensões do meio
condutor, sendo dada por
(13) α= A
dR
Modelagem Matemática de Sistemas Térmicos
5
5
Só podemos usar a eq. (12) quando não há armazenamento de energia térmica no meio condutor.
Caso isso não aconteça, devemos então incluir a capacitância térmica do meio condutor no
modelo.
Resistências térmicas em série
Consideremos dois corpos com temperaturas θ1(t) > θ2(t), separados por duas resistências
térmicas em série R1 e R2, conforme ilustra a fig. 1(a):
 Fig. 1
Sendo q(t) o fluxo de calor através das mesmas e estando as resistências perfeitamente isoladas
termicamente, queremos achar uma resistência térmica equivalente Req, conforme a fig. 1(b).
Chamando θB a temperatura na interface das duas resistências, podemos escrever a eq. (12) duas
vezes:
 )(
R
1q B1
1
θ−θ=
)(
R
1q 2B
2
θ−θ=
Eliminando θB nas equações acima, chegamos a
 )(
RR
1q 21
21
θ−θ+=
que, comparada com a eq. (12), permite que escrevamos
 Req = R1 + R2
donde podemos concluir que existe uma analogia com as resistências elétricas em série.
Podemos estender o resultado para n resistências térmicas em série:
(14) ∑
=
=
n
1i
ieq RR
Modelagem Matemática de Sistemas Térmicos
6
6
Resistências térmicas em paralelo
Aproveitando a analogia citada, podemos estabelecer uma expressão para n resistências
térmicasem paralelo:
(15) 
∑
=
= n
1i i
eq
R
1
1R
7 FONTE TÉRMICA
A fonte térmica ideal adiciona ou retira energia térmica do sistema. No primeiro caso, o fluxo de
calor qi(t) é positivo e, no segundo caso, qi(t) é negativo. A fonte térmica ideal é representada
pela fig. 2:
 Fig. 2
Vamos estudar, a seguir, a modelagem matemática de alguns sistemas térmicos através de
exemplos.
8 MODELOS MATEMÁTICOS DE SISTEMAS TÉRMICOS
Exemplo 1
A fig. 3 mostra uma capacitância térmica C isolada do ambiente por uma resistência térmica
equivalente R. A temperatura interna é θ, considerada uniforme, enquanto que a temperatura
ambiente é θa, também uniforme. Calor é adicionado ao interior do sistema com um fluxo qi(t). No
ponto de operação, os valores de qi e θ são 
-
i e q θ
−
, respectivamente. Desenvolver um modelo
matemático para o sistema em termos das variáveis incrementais.
 Fig. 3
Modelagem Matemática de Sistemas Térmicos
7
7
Solução
Aplicando a eq. (12): ])t([
R
1)t(q ao θ−θ=
Substituindo na eq. (11): ]})t([
R
1)t(q{
C
1)t( ai
. θ−θ−=θ
ou 
onde reconhecemos uma EDOL de 1a ordem com coeficientes constantes, não homogênea, com
duas entradas qi(t) e θa e saída θ(t). A constante de tempo do sistema é dada por τ = RC.
Em termos das variáveis incrementais 
−θ−θ=θ )t()t(^ e
 
−−= q)t(q)t(q^
podemos obter um modelo matemático substituindo θ(t), qi(t) e suas derivadas na EDOL acima,
chegando a
Vemos, agora, que temos um sistema com apenas uma entrada e uma saída.
Exemplo 2 - Termômetro de mercúrio
A fig. 4 ilustra um sistema térmico constando de um termômetro de mercúrio que está,
inicialmente, à temperatura ambiente 
−θ e é mergulhado em um reservatório cujo líquido está a
uma temperatura 
−θ + θb, isto é, θb acima da temperatura ambiente.
 Fig. 4
O reservatório tem capacitância térmica C e o termômetro tem resistência térmica R.
Desenvolver um modelo matemático para o sistema em termos das variáveis incrementais.
ai
.
)t(Rq)t()t(RC θ+=θ+θ
)t(qR)t()t(RC i
^^
.
^ =θ+θ
Modelagem Matemática de Sistemas Térmicos
8
8
 Solução
Aplicando a eq. (12) para o termômetro: )]t([
R
1)t(q bi θ−θ=
Substituindo na eq. (11): )]}t([
R
1{
C
1)t( b
. θ−θ=θ
ou 
(16)
onde reconhecemos uma EDOL de 1a ordem com coeficientes constantes, não homogênea, com
entrada θb e saída θ(t). A constante de tempo do sistema é dada por τ = RC.
Comparando a eq. (16) com a EDOL modelo matemático do circuito elétrico RC paralelo mostrado
na fig. 5, dada por
eedt
edRC ioo =+
 Fig. 5
vemos que existe uma analogia entre o sistema térmico e o sistema elétrico, denominada analogia
eletrotérmica, dada pela tabela seguinte:
Sistema elétrico Sistema térmico
voltagem e temperatura θ
corrente elétrica i fluxo de calor q
resistência elétrica R resistência térmica R
Capacitância C capacitância térmica C
Exemplo 3
A fig. 6 mostra um vaso indeformável de volume V, no qual um líquido de massa específica ρ e
calor específico c escoa através dele. Um "mixer" assegura que a temperatura do líquido
permaneça uniforme em todo o reservatório e igual a θ(t). O líquido entra no reservatório com
uma vazão volumétrica constante 
−
w à temperatura θi(t). Ele sai do reservatório com a mesma
vazão volumétrica à temperatura θo(t), considerada igual à temperatura do líquido θ(t), devido à
mistura perfeita feita pelo "mixer". A resistência térmica do vaso é R e a temperatura ambiente
é constante e igual a θa.
b
.
)t()t(RC θ=θ+θ
Modelagem Matemática de Sistemas Térmicos
9
9
 Fig. 6
É adicionado um fluxo de calor qh(t) ao líquido por meio de um aquecedor. Desenvolver um modelo
matemático para o sistema.
Solução
Calor que entra no vaso: )t(cw)t(q)t(q ihi θρ+=
−
Calor que sai do vaso: [ ] )t(cw)t(
R
1)t(q
_
a0 θρ+θ−θ=
Capacitância térmica: C = Mc = ρVc
Levando na eq. (11):
 )]}t(cw))t((
R
1[)]t(cw)t(q{[
cV
1)]t(q)t(q[
C
1)t( aihoi
. θρ+θ−θ−θρ+ρ=−=θ
−−
Rearrumando a equação acima, chegamos à EDOL de 1a ordem
onde a constante de tempo é dada por 
RC
1
V
w
1
+
=τ − .
Podemos observar que temos três entradas, θi(t), qh(t) e θa, e apenas uma saída, θ(t).
Exemplo 4
Uma esfera de cobre (ρ = 8954 kg/m3, c = 383,1 J/kg.0C e k = 385 W/m. 0C), de diâmetro 0,06
m, é subitamente colocada em um reservatório que contem um líquido quente a uma temperatura
θo. Em conseqüência, a temperatura da esfera, θ(t), cresce com o tempo. O coeficiente de
transferência de calor por convecção é h = 25 W/ m. 0C. Pedem-se:
ahi
.
RC
1)t(q
C
1)t(
V
w)t()
RC
1
V
w()t( θ++θ=θ++θ
−−
Modelagem Matemática de Sistemas Térmicos
10
10
(a) É válido usar um modelo com parâmetros concentrados? Justifique;
(b) Caso positivo, obter um modelo matemático para o sistema;
(c) Calcular a constante de tempo do sistema.
Solução
(a) m 01,0
2
06,0.
3
1r
3
1Lc ===
1,010x49,6
385
01,0x25
k
hLBi 4c <=== −
Logo, é possível.
(b) Aplicando a eq. (11):
 )]t(q)t(q[
C
1)t( oi
. −=θ
onde C = Mc = ρVc
 )]t([hAq osi θ−θ=
 0qo =
Logo: { }0)]t([hA
Vc
1)t( os
. −θ−θρ=θ
 o
.
s
)t()t(
hA
Vc θ=θ+θρ
(c) h 87,22s 1372
25
01,0x1,383x8954
h
cL
hA
Vc c
s
===ρ=ρ=τ
Modelagem Matemática de Sistemas Térmicos 11
EXERCÍCIOS
1 A figura mostra um vaso fechado, isolado, cheio de líquido e contendo um aquecedor
elétrico imerso no líquido. A resistência elétrica do aquecedor, por sua vez, está colocada
dentro de uma jaqueta metálica de resistência térmica RHL. A resistência térmica do vaso
e de seu isolamento é RLa. O aquecedor tem uma capacitância térmica CH e o líquido uma
capacitância térmica CL. A temperatura do aquecedor é θH e a do líquido é θL, a qual é
considerada uniforme devido ao "mixer".
O aquecedor elétrico e o líquido estão inicialmente à temperatura ambiente θa, estando o
aquecedor desligado. No instante t = 0, o aquecedor é ligado, de modo que o fluxo de calor
fornecido ao sistema é qi(t). Pedem-se:
(a) modelo matemático no espaço de estados, usand tado θH(t) e θL(t), as quais
podem ser obtidas diretamente a partir da eq. 
(b) usando o VisSim, graficar as temperaturas θH( ntradas θa = 300 K e qi(t)
sendo um degrau de amplitude 1,5 x 104 W;
(c) a partir do gráfico do item (b), achar o tempo ara atingir a temperatura
desejada θd = 365 K.
Obs.: para os itens (b) e (c) usar os dados numérico da figura.
2 
 
 
 
Da
Re
 Dados numéricos:
CH = 20 x 103 J/K
CL = 1 x 106 J/K
RHL = 1 x 10-3 K/W
RLA = 5 x 10-3 K/W
θa = 300 K
Uma esfera de alumínio de diâmetro 0,08 m
200 0C. Ela é retirada do forno e colocada
20 0C. Conhecendo as propriedades do alum
 (a) É válido usar um modelo com parâmetros
 (b) Caso positivo, obter um modelo matemá
 (c) Calcular a constante de tempo do sistem
dos do Alumínio: ρ = 2707 kg/m3; c = 8
h = 3,5 W/ m. 0C
sp.: (a) Sim; (b) ar
.
s
)t()t(
hA
Vc θ=θ+θρ
o as variáveis de es
(9);
t) e θL(t) para as e
 que leva o líquido p
s mostrados ao lado
11
 encontra-se em um forno à temperatura de
 ao ar livre que se encontra à temperatura de
ínio, dadas abaixo, pedem-se:
 concentrados? Justifique;
tico para o sistema;
a.
96 J/kg.0C; k = 204 W/m. 0C;
(c) 2,56 h
Modelagem Matemática de Sistemas Térmicos
12
12