Daniel Nack Gisielly Schoeffel Linara Battisti Archer Maicon Luiz Hilgenstieler Sharlene Schmitt

DETERMINAO EXPERIMENTAL DO COEFICIENTE CONVECTIVO DE CALOR (h)

Trabalho apresentado para a disciplina de Laboratrio de Engenharia II do Curso de Engenharia Qumica do Centro de Cincias Tecnolgicas da Universidade Regional de Blumenau Professor: Larcio Ender

Blumenau, abril de 2003

INTRODUO
A cincia da Transferncia de Calor (T.C.) se refere anlise da taxa de transferncia de calor em um sistema. A energia transformada pelo fluxo de calor no pode ser diretamente, mas o conceito tem significado fsico porque relacionada grandeza mensurvel chamada Temperatura. O fluxo de calor ocorre sempre que h um gradiente de temperatura em um sistema, por isso, o conhecimento da distribuio de temperatura no sistema essencial nos estudos de transferncia de calor. O fluxo de calor a quantidade de calor transferido por unidade de rea por unidade de tempo. Existem trs modos de transferncia de calor: conduo, conveco e radiao. A distribuio de temperatura em um corpo controlada pelos efeitos dos trs modos de T.C., mas muitas vezes, para fins prticos, considera-se cada um separadamente.

FUNDAMENTAO TERICA
Para melhor entender o experimento deve-se deixar claro os conceitos de conduo e conveco. Conduo o modo de transferncia de calor em que a troca de energia tem lugar da regio da alta temperatura para a regio de baixa temperatura pelo movimento cintico ou pelo impacto direto de molculas, no caso de fluidos em repouso, e pelo movimento de eltrons no caso de metais. A lei emprica da conduo de calor Lei de Fourier estabelece qual o fluxo de calor por conduo, em uma dada direo do fluxo e ao gradiente de temperatura naquela direo. Se o fluxo de calor for na direo de x, a Lei de Fourier dada da seguinte forma: Qx = - K. dT dx ou qx = Qx = - K. dT A dx (2) (1)

Onde, Qx a taxa de fluxo de calor atravs da rea A no sentido positivo do eixo x; qx o fluxo de calor no sentido positivo do eixo x e K a condutividade trmica do material. Deste modo pode-se determinar a distribuio de temperatura nos slidos. O uso da Lei de Fourier tem como requisito o conhecimento da condutividade trmica (k). Esta propriedade, classificada como uma das propriedades de transporte da matria, fornece uma indicao da taxa segundo a qual a energia transferida pelo processo de difuso. Para um dado gradiente de temperatura, o fluxo trmico condutivo aumenta com o aumento da condutividade trmica. A condutividade trmica de um slido maior do que um lquido que por sua vez maior do que a de um gs. A conveco ocorre sempre que um fluido escoar sobre uma superfcie slida que est a uma temperatura diferente do fluido.

Se o movimento do fluido for induzido artificialmente, por alguma fora externa, diz-se que a transferncia de calor se processa por conveco forada. Se o movimento do fluido resultar dos efeitos da ascenso provocada pela diferena de densidade causada pela diferena de temperatura no fluido sobre a superfcie, a transferncia de calor chamada de conveco natural. A determinao da distribuio de temperatura e a transferncia de calor na conveco so complicadas. Nas aplicaes da engenharia, para simplificar os clculos da transferncia de calor entre uma superfcie quente, a Tw e um fluido que est escoando sobre ela a uma temperatura Tf, define-se um coeficiente de transferncia de calor como: q = h.(Tw Tf) fluido. Ou: q = h.(Tf Tw) quente para a parede fria. O coeficiente de transferncia de calor (h), varia com o tipo de fluxo (laminar ou turbulento), com a geometria do corpo e a rea de escoamento, com as propriedades fsicas do fluido, com a temperatura mdia e com a posio ao longo da superfcie do corpo, e ainda, depende do tipo de conveco (forada ou natural). A determinao do coeficiente de transferncia de calor a chave para a anlise da distribuio da temperatura em um corpo que est exposto a conveco. E justamente a determinao deste coeficiente o objetivo deste experimento. (4) A equao (4) utilizada quando a transferncia de calor se d do fluido (3) Onde q o fluxo de calor expresso em W/m 2 da parede quente para o

Sistemas com Conduo e Conveco O calor conduzido atravs de um corpo deve ser freqentemente removido (ou fornecido) por algum processo de conveco. O calor conduzido atravs do material e finalmente dissipado no ambiente por conveco. Para o estudo de sistemas que combinam conduo e conveco, deve-se fazer o uso do nmero de Biot.

Nmero de Biot O nmero de Biot (Bi) um parmetro admensional e representa a razo entre o coeficiente de transferncia convectiva de calor na superfcie do slido e a condutncia especfica do slido. A hiptese de temperatura uniforme no interior do slido vlida se a condutnicia especfica do slido for muito maior do que o coeficiente de transferncia convectiva de calor.

Bi = h.Ls / Ks

(5)

Onde: Ls = comprimento caracterstico; Ls = V/A (volume/rea) h = coeficiente convectivo de calor Ks = coeficiente condutivo de calor

O Biot usado para definir o mtodo a ser utilizado na soluo de problemas de Transferncia de calor transiente. Se Bi 0,1 : usa-se as cartas de temperatura transiente Se Bi < 0,1 : usa-se a anlise global, -dT = h.A. . ( To - T ) dt .Cp.V

Emprego das Cartas de Temperatura Transiente As cartas de temperatura transiente so empregadas quando os gradientes de temperatura no so desprezveis e no aplicvel a anlise global do sistema. Neste caso, a anlise dos problemas de conduo de calor envolve a determinao da distribuio de temperatura no interior do slido em funo do tempo e posio. Para a esfera usa-se as seguintes cartas: Anlise Global: Na anlise global, admite-se que a distribuio de temperatura dentro do slido, em qualquer instante seja suficientemente uniforme, de tal modo que a temperatura do slido possa a ser considerada funo exclusiva do tempo, isto , T(t). A equao da energia na transferncia de calor no slido pode ser escrita como: Taxa de T.C. por conveco no volume V do slido = Aumento da energia Interna do slido

h . A . ( T - To) (6)

.Cp.V

dT/dt

-dT(t) (7) dt

h.A.

To

.Cp.V

onde: To a temperatura inicial T a temperatura do fluido Cp a capacidade calorfica do slido

 a massa especfica do slido

METODOLOGIA
Para a realizao do experimento utilizaram-se duas esferas (uma de teflon e outra de alumino), um banho termostato, um multmetro, termopares e um cronmetro. As esferas foram primeiramente aquecidas no banho termostato (at uma temperatura prxima a 90o C) e quando toda a esfera estava com a mesma temperatura ela foi retira do banho e foi exposta ao ar livre. Foram medidas as temperaturas em alguns pontos da esfera em vrios intervalos de tempo. Com estes dados de temperatura versus tempo foram traadas as curvas de resfriamento de cada esfera.

Dados experimentais da esfera de Alumnio: t(s) 0 10 20 40 55 70 90 110 125 140 155 170 190 210 230 250 290 310 340 T (0) 88 86 85 84 84 82 82 82 82 81 81 81 80 80 80 80 79 79 79 T(r) 88 86 85 85 84 83 83 83 83 82 82 82 81 81 81 80 80 80 79 t(s) 370 400 430 460 490 520 590 650 710 1020 1320 1920 2520 3120 3720 5520 7320 9120 10920 T(0) 78 78 77 76 76 75 74 73 73 69 66 61 56 52 49 43 39 36 36 T(r) 79 79 78 77 76 76 75 74 73 70 67 61 57 53 49 43 39 36 36

Tabela 1 Dados experimentais da esfera de alumnio.

Curva de Resfriamento da Esfera de Alumnio

100 90 80 Temperatura (C) 70 60 50 40 30 20 10 0 0 5000 tempo (s) 10000 15000 T(0) C T(R) C

Grfico 1 Curva de resfriamento da esfera de alumnio. Atravs da anlise experimental, pode-se considerar que a transferncia de calor na esfera de alumnio ocorre uniformemente. Pode-se ento fazer uso da anlise global para determinar o valor do coeficiente convectivo de calor (h) para a esfera de alumnio. A temperatura final da esfera de alumnio a mesma que a do fluido, no caso o ar, utilizado para resfri-la. Dados: Cp = 919,28 J / kg C K = 239,4175 W / m C = 2702 kg / m

A = 0,0314 m R = 0,05 m Tempo de equilbrio = 10920 segundos T(t) = -1.10 18 t5 + 4.10-14 t4 6.10-10 t + 4.10-6 t - 0,0191 t + 85,233 T(10920) = C dT(t) / dt = -5.10 18 t4 + 16.10-14 t3 18.10-10 t2 + 8.10-6 t - 0,0191 dT (10920) / dt = -0,01018 dT(t) = h.A. . ( To - T ) dt .Cp.V -0,01018 = h.0, 0314 (88 36) 2702 . 919,28 . 4/3 . . (0,05)3 Isolando h temos: h = 8,1086 W/m2C (8)

Dados experimentais da esfera de Teflon: Tempo (s) 0 10 30 60 90 120 T ( R) 91 91 90 89 87 85 T (3/4 R) 91 91 91 91 91 90 T (R/2) 91 91 91 91 91 91 T(0) 91 91 91 91 91 91

150 83 90 91 180 82 89 91 210 81 89 91 240 80 89 91 270 79 88 90 300 78 88 90 330 77 88 90 390 75 87 90 510 72 86 90 630 70 85 89 750 69 83 88 1050 65 80 86 1350 62 77 83 1650 60 74 81 1950 58 71 77 2550 55 64 70 3150 52 59 64 3750 49 55 60 4350 47 51 55 4950 45 49 51 5550 43 46 49 6750 40 42 44 7950 39 40 41 9150 38 39 39 10350 37 38 38 11550 37 37 37 12750 36 37 37 14550 36 36 36 16350 36 36 36 18150 36 36 36 Tabela 2 Dados experimentais da esfera de teflon.

91 91 91 91 91 91 91 91 91 91 90 89 86 84 80 73 66 61 56 52 49 44 41 39 38 37 37 36 36 36

Curva de resfriamento da Esfera de teflon

100 Temperatura (C) 80 T ( R) 60 40 20 0 0 5000 10000 Tempo (s) 15000 20000 T (3/4 R) T (R/2) T(0)

Grfico 2 Curva de resfriamento da esfera de teflon. Como a esfera de teflon possui variao de temperatura com o tempo e com a posio, utilizam-se as cartas de temperatura transiente. Fazendo-se uso da primeira carta, acha-se o nmero de Biot no grfico e consequentemente o valor do coeficiente de conveco (h) da esfera de teflon. Para isso usa-se um tempo qualquer do resfriamento da esfera, no caso utilizamos o tempo de 210 segundos. Dados: Cp = 284076 J / kg C K = 0,41415 W / m C = 2200 kg/ m = 1.10-7 m2 / A = 0,0314 m R = 0,05 m Para achar o nmero de Biot no grfico necessita-se achar os valores de e de o.
s

(9)

o = To - T = 81 36 = 0,82
Ti - T 91 36

To a temperatura na superfcie da esfera Ti a temperatura no centro da esfera T a temperatura do fluido (ar)

= . t / r = 1.107 .210 / (0,005) = 0,0084

(10) Logo: 1/ Bi = 0,2 Conhecendo o valor de Biot poderamos usar as cartas de temperatura transiente, mas como esse valor foi maior que 0,1 no est dentro do compreendido nas cartas, o que nos fez fazer a anlise dos dados experimentais por outra forma. Para isto usamos a equao da conservao da energia trmica em escala macroscpica para uma esfera, onde obtivemos a seguinte equao: T = (r T/r) t r2 r (11)

A partir dessa equao, aplicando as condies de contorno, obtemos um sistema de Equaes Diferenciais Parciais de 2 Ordem: para t = 0 T(0,r) = To>T para r = 0, em qualquer tempo (T/r)r=0 = 0 para r = R, em qualquer tempo hconv . [ T(t,R) - T ] = - k (T/r )R

Conhecendo algumas equaes podemos ter a soluo exata:

 * = Cn. exp . (- n2 . Fo) . (1/ n . r*) . [sen ( n . r*)]

n=1

(12)

Onde: Cn = 4.[ sen ( n) - n . cos ( n)] (13)

2 . n sen (2 . n)

e, 1 - n . cotg ( n) = Bi Tendo as relaes: (14)

* = [T(r,t) - T] / (To - T)
r* = r/R Fo = . t / R2

(15) (16) (17)

Com essas equaes caimos numa equao transcedental, onde podemos dar um chute inicial para o valor de Biot, com isso acha os valores de n, depois de Cn, com isso achamos os valores de *. Mas para achar esses valores teramos que fazer um programa em alguma linguagem de programao, pois so muitos valores diferentes que devem ser chutados para Biot. Tentamos fazer esse programa em FORTRAN, mas ele no rodou. Ento tentamos fazer pelo mtodo grfico, encontrando uma curva para cada temperatura em cada tempo medido, com isso achendo a equao de cada curva e a derivada dessa equao, para ento tentar achar o valor de h pela equao (8), como mostrado a seguir:

Mas dessa forma tambm no conseguimos chegar a uma concluso. Acabamos no chegando a uma valor final do h convectivo no caso do resfriamento da esfera de teflon.

CONCLUSO
Para o caso da esfera de alumnio chegamos a um valor aceitvel, pois nas literaturas encontramos valores que variam um pouco, mas todos prximos de 10 W/m2C, e ns obtivemos um valor de h = 8,1086 W/m2C. Para o caso do resfriamento da esfera de teflon no chegamos a um valor final, mas com certeza o h de conveco ter que ser menor que o obtido no caso do alumnio, pois o tempo de resfriamento foi maior. Alguns erros podem ter ocorrido devido a correntes de ar presentes no laboratrio, variao de temperatura no decorrer do experimento, uso de silicone para fixar os termopares, o que pode ter influenciado na transferncia de calor, assim como alguns erros de leitura.

REFERNCIAS BIBLIOGRFICAS

BEJAN, Adrian. Heat Transfer. New York, John Wiley & Sons, Inc,1993 INCROPERA, Frank P. Fundamentops de Transferncia de Calor e de

massa. Rio de Janeiro: LTC, 1992 OZISIK, M. Necati. Transferncia de Calor: um texto bsico. Bogot: McGraw-Hill,1990. SCHMIDT, Frank W. Introduo s cincias trmicas. 2ed. So Paulo: Ed. Edgard Blcher,1996.