Condução térmica

Em transferência de calor, condução térmica ou difusão térmica (ou ainda condução ou difusão de calor) é a transferência de energia térmica entre átomos e/ou moléculas vizinhas em uma substância devido a um gradiente de temperatura. Noutras palavras, é um modo do fenômeno de transferência térmica causado por uma diferença de temperatura entre duas regiões em um mesmo meio, ou entre dois meios em contato, e sem perceber-se movimento global da matéria (na escala macroscópica) em oposição a convecção que é outra forma de transferência térmica.
Isso sempre ocorre a partir de uma região de maior temperatura para uma região de baixa temperatura, e atua equalizando as diferenças de temperatura. Genericamente, ocorre a propagação de calor sem transporte da substância formadora do sistema, ou seja, através de choques entre suas partículas integrantes ou intercâmbios energéticos dos átomos, moléculas, e elétrons. A condução térmica pode ser interpretada como a transmissão passo a passo de agitação térmica: um átomo (ou uma molécula) transfere parte de sua energia cinética ao átomo vizinho, sendo assim um fenômeno de transporte de energia interna devido à heterogeneidade da agitação molecular, sendo assim um fenômeno termodinamicamente irreversível.
Condução ocorre em todas as formas de matéria, sólidos, líquidos, gases e plasmas, mas não requer qualquer movimento de massa da matéria. Em sólidos, é devido à combinação das vibrações das moléculas em um retículo cristalino e o transporte de energia por elétrons livres. Os elétrons suscetíveis de mover denominam-se elétrons de condução no modelo do elétron livre.
Em fluidos (líquidos e gases) o transporte de energia é resultante da não-uniformidade do número de choques por unidade de volume, durante seu movimento aleatório, semelhante ao fenômeno da difusão. Em sólidos, a condução de calor é fornecida conjuntamente por condução de elétrons e vibração da rede cristalina (fônon)
Os metais e suas ligas, sejam sólidos ou líquidos, devido à elevada condutividade térmica, são excelentes meios de propagação de calor, normalmente associada à condutividade elétrica.Os gases e alguns sólidos, que possuem baixa condutividade térmica, são considerados péssimos meios de propagação de calor, sendo definidos como isolantes térmicos.
Calor pode ser transferido por radiação e/ou convecção, e frequentemente mais que um destes processos ocorrem em uma dada situação.

Lei de Fourier

A lei da condução térmica, também conhecida como lei de Fourier, estabelece que a taxa no tempo da transferência de calor através de um material é proporcional ao gradiente negativo na temperatura e na área em ângulos retos, para o gradiente, através do qual o calor está fluindo. Podemos afirmar, esta lei de duas formas equivalentes: a forma integral, em que olhamos para a quantidade de energia que flui para dentro ou para fora de um corpo como um todo, e na forma diferencial, em que olhamos para os fluxos de energia localmente.

Lei de Fourier.

Pode-se determinar o fluxo de calor transportado por condução pela Lei de Fourier:

$q^{\prime \prime} = \frac{q}{A} = -k \frac{\partial T}{\partial x}.$

A expressão acima aplica-se ao caso unidimensional, quando há gradiente de temperatura apenas na direção x.
Caso se conheça as temperaturas de duas superfícies específicas e queira-se calcular o fluxo de calor por condução entre elas, basta integrar a equação acima, que toma a forma:

$q^{\prime \prime} = \frac{k(T2 - T1)}{L}.$

A constante k, é a condutividade térmica do material. Entre duas substâncias, a que tiver condutividade maior conseguirá transferir uma quantidade maior de calor, para uma mesma diferença de temperatura.

Teorização completa

Modelo de um tubo de aquecimento, resfriado por hastes metálicas

Existem várias grandezas envolvidas, mas entre elas existem duas que são de muita importância de interesse prático no estudo de problemas de condução de calor. Estas grandezas são a razão de fluxo de calor e a distribuição da temperatura. As razões de fluxo de calor tratam da demanda de energia em um dado sistema, quando se requer uma distribuição de temperaturas conveniente para desenhar de maneira adequada no sistema, desde o ponto de vista dos materiais. Em um fenômeno qualquer, uma vez que seja conhecida a distribuição da temperatura é possível determinar as razões de fluxo de calor com ajuda da denominada Lei de Fourier (de 1822, estabelecida por Jean Baptiste Joseph Fourier).
A distribuição da temperatura é linear, e o fluxo de calor é constante de um extremo a outro de uma placa, para o caso da equação radial produzida.
E portanto a distribuição da temperatura apresenta-se em forma logarítmica:

$\ T = M.ln(r) + N$

O calor transferido $\dot{Q}$ é tratado pela lei de Fourier que descreve especificamente previsões (modelagens) de comportamento para o caso simples de um corpo sólido, com duas paredes paralelas^[5]^[6]:

$\dot{Q} = {\frac{\lambda}{d}} A (T_{W_{1}}-T_{W_{2}})$

A unidade de $\dot{Q}$ é o Watt (W), e sendo as grandezas:

$T_{W_{1}}$ a temperatura da superfície da parede mais quente
$T_{W_{2}}$ a temperatura da superfície da parede fria
$A$ da área através da qual o calor flui,
$λ$ a condutividade térmica, geralmente um parâmetro do material dependente da temperatura, e
$d$ a espessura do corpo, medido de parede a parede.

Atualmente a transferência de calor é descrita através do conceito mais rigoroso de fluxo de calor $\dot{\vec{q}}$ , em abordagens que visam reduzir-se aos tratamentos de Fourier e Newton. A notação $\dot{\vec{q}}$ é formulada apatir da derivada parcial no tempo do vetor fluxo de calor $\vec{q}$ . Aplica-se a seguinte definição:

$\dot{\vec{q}} = -\lambda \, \operatorname{grad} \, T$

Matematicamente, o fenômeno de "transferência de calor" é descrito por uma equação diferencial parcial, apresentando um padrão parabólico. Esta equação diferencial parcial, na forma especificada, apresenta a forma geral:

$\frac{\partial u(\vec r,t)}{\partial t} =a\, \Delta u(\vec{r},t)$

Sendo esta equação especial e chamada comumente equação do calor. Note-se que esta forma da equação do calor é válida somente para meios homogêneos e isotrópicos; noutras palavras, para meios que possuem a mesma composição em todos os lugares e nenhuma orientação preferencial (ocorrem orientações preferenciais, por exemplo, em fibras de materiais compostos, mas também por dilatações de grãos em chapas de aço laminadas, etc). Para estes casos - e apenas para isso , as propriedades materiais são adotadas com o objetivo de considerar apenas as grandezas dependentes da temperatura. Estritamente falando, a equação assim formulada não se aplica apenas quando o calor no corpo é introduzido ou removido por fenômenos estranhos à modelagem utilizada, sendo, neste caso, a fonte ou "fuga" um termo a ser adicionado ao equacionamento. Com estas restrições, segue-se a seguinte forma da equação do calor:

$\frac{\partial T(\vec r,t)}{\partial t} =a(T)\, . \Delta T(\vec{r},t)$

Sendo esta a equação diferencial que descreve os processos de transporte em geral (como o processo de difusão - que é um transporte material, devido a ser compreendido como uma diferença de concentração ou, no caso da equação do calor, apenas o "andar" da distribuição de temperatura em um corpo devido a um gradiente de temperatura). A solução analítica dessa equação não é em possível muitos casos práticos, sendo calculados, atualmente, com a ajuda de métodos tecnicamente relevantes de cálculos para fenômenos de transferência de calor como o método dos elementos finitos, obtendo-se como resultado a distribuição temporal da temperatura no espaço (um campo de temperaturas). Assim, pode-se obter influências, por exemplo, relacionadas ao comportamento de expansão espacial dos componentes (dilatação térmica), que por sua vez influenciará o estado de tensão local. Assim, o campo de temperaturas torna-se uma base importante para todas as tarefas de engenharia em que o componente de estresse térmico não pode ser negligenciado no projeto.

Forma diferencial

A forma diferencial da lei de Fourier da condução térmica mostra que o fluxo de calor local, $\overrightarrow{q}$ , é igual ao produto da condutividade térmica,

k

, e o gradiente de temperatura local negativo, $-\nabla T$ . O fluxo de calor é a quantidade de energia que flui através de uma superfície particular por unidade de área por unidade de tempo.

$\overrightarrow{q} = - k {\nabla} T$

onde (incluindo as unidades SI)

$\overrightarrow{q}$ é o fluxo de calor local, [W·m⁻²]

$\big.k\big.$ é a condutividade térmica do material, [W·m⁻¹·K⁻¹],

$\big.\nabla T\big.$ é o gradiente de temperatura, [K·m⁻¹].

A condutividade térmica,

k

, é frequentemente tratada como um constante, embora isto nem sempre seja verdade. Enquanto a condutividade térmica de um material geralmente varia com a temperatura, a variação pode ser pequena para uma significativa faixa de temperaturas para alguns materiais comuns. Em materiais anisotrópicos, a condutividade térmica varia tipicamente com orientação; neste caso

k

é representada por um tensor de segunda ordem. Em materiais não uniformes,

k

varia com a localização espacial.
Para muitas aplicações simples, a lei de Fourier é usada em sua forma unidimensional. Na direção x,

$q_x = - k \frac{d T}{d x}$

Forma integral

Pela integração da forma diferencial sobre a superfície total do material

S

, chega-se à forma integral da lei de Fourier:

$\frac{\partial Q}{\partial t} = -k \oint_S{\overrightarrow{\nabla} T \cdot \,\overrightarrow{dA}}$

onde (incluindo as unidades SI)

$\big. \frac{\partial Q}{\partial t}\big.$ é a quantidade de calor transferido por unidade de tempo [W] e

$\overrightarrow{dA}$ é um elemento de área de superfície orientado [m²]

A equação diferencial acima, quando integrada para um material homogêneo de geometria 1-D entre dois pontos a uma temperatura constante, dá a taxa de fluxo de calor como:

$\big. \frac{\Delta Q}{\Delta t} = -k A \frac{\Delta T}{\Delta x}$

onde

A é a área da superfície da seção transversal,

Δ T

é a diferença de temperatura entre os extremos,

Δ x

é a distância entre os extremos.

Esta lei forma a base para a derivação da equação do calor. A lei de Ohm é o análogo elétrico da lei de Fourier.

Condutância

Escrevendo-se

$\big. U = \frac{k}{\Delta x}, \quad$

onde U é a condutância, em W/(m² K).
Lei de Fourier pode também ser enunciada como:

$\big. \frac{\Delta Q}{\Delta t} = U A\, (-\Delta T).$

A recíproca da condutância é a resistência, R, dada por:

$\big. R = \frac{1}{U} = \frac{\Delta x}{k} = \frac{A\, (-\Delta T)}{\frac{\Delta Q}{\Delta t}}, \quad$

e é a resistência que é aditiva quando várias camadas codutivas situam-se entre as regiões quente e fria, porque A e Q são as mesmas para todas as camadas. Em uma partição de multicamadas, a condutância total está relacionada com a condutância de suas camadas por:

$\big. \frac{1}{U} = \frac{1}{U_1} + \frac{1}{U_2} + \frac{1}{U_3}+ \cdots$

Assim, quando trata-se de uma partição de multicamadas, a seguinte fórmula é geralmente usada:

$\big. \frac{\Delta Q}{\Delta t} = \frac{A\,(-\Delta T)}{\frac{\Delta x_1}{k_1} + \frac{\Delta x_2}{k_2} + \frac{\Delta x_3}{k_3}+ \cdots}.$

Quando o calor é conduzido a partir de um fluido para outro através de uma barreira, às vezes é importante considerar a condutância da fina película de fluido que permanece estacionária próximo a barreira. Esta fina película de líquido é difícil de quantificar, suas características, dependendo das complexas de turbulência e viscosidade, mas quando se lida com os obstáculos finos de alta condutância por vezes pode ser bastante significativa.

Representação de propriedade intensiva

As equações de condutância anteriores escritas em termos de propriedades extensivas, podem ser reformuladas em termos de propriedades intensivas. Idealmente, as fórmulas de condutância devem produzir uma grandeza com dimensões independentes da distância, como a lei de Ohm da resistência elétrica: $R = V/I\,\!$ , e condutância: $G = I/V \,\!$ .
A partir da fórmula elétrica: $R = \rho x / A \,\!$ , onde ρ é a resistividade, x = comprimento, A área de seção tranversal, temos $G = k A / x \,\!$ , onde G é condutância, k é condutividade, x = comprimento, A = área de seção tranversal.
Para o calor,

$\big. U = \frac{k A} {\Delta x}, \quad$

onde U é a condutância.
Lei de Fourier pode ser enunciada como:

$\big. Q = U \, \Delta T \quad$

análoga a lei de Ohm: $I = V/R \,\!$ ou $I = V G. \,\!$
A recíproca da condutância é resistência, R, dado por:

$\big. R = \frac{\, \Delta T}{Q}, \quad$

análoga a lei de Ohm: $R = V/I. \,\!$
A soma das condutâncias em série é ainda correta.

Cilindros

Condução através de cilindros pode ser calculada quando conhece-se variáveis tais como o raio interno

r 1

, o raio externo

r 2

, e o comprimento notado como $\ell$ .
A diferença de temperatura entre a parede interna e externa pode ser expressa como

T 2 - T 1

.
A área do fluxo de calor: $A_r= 2 \pi r \ell$
Quando a equação de Fourier é aplicada:

$Q = -k A_r \frac{\mathrm{d}T}{\mathrm{d}r} = -2 k \pi r \ell \frac{\mathrm{d}T}{\mathrm{d}r}$

Rearranjando-se:

$Q \int_{r_1}^{r_2} \frac{1}{r} \mathrm{d}r = -2 k \pi \ell \int_{r_1}^{r_2} \mathrm{d}T$

Portanto, a taxa de transferência de calor é

$Q = 2 k \pi \ell \frac{T_1 - T_2}{\ln r_2 - \ln r_1}$

A resistência térmica é

$R_c = \frac{\Delta T}{Q} = \frac{\ln r_2 - \ln r_1}{2 \pi k \ell}$

Tendo-se $Q = 2 \pi k \ell r_m \frac{T_1-T_2}{r_2-r_1}$ , onde $r_m = \frac{r_2-r_1}{\ln r_2 - \ln r_1}$ sendo importante notar que este é o raio log médio.

Condução em regime estacionário

Define-se um regime estacionário (ou permanente), quando as temperaturas não dependem do tempo. A temperatura depende apenas da disposição do ponto em que as medidas são tomadas ao longo do tempo. Para o restante desta divisão deste artigo, assume-se como estabelecido um estado de equilíbrio.

Superfície plana simples

O material é um meio condutor térmico delimitado por dois planos paralelos (o caso de uma parede). Cada plano tem uma temperatura T homogênea por toda a sua superfície. Consideramos que os planos têm duas dimensões infinitas para superar os efeitos de borda. Por conseguinte, o fluxo de entrada é igual ao de saída, não há perda de calor nas extremidades.

Perfil de temperatura em uma parede.

Note-se que

T 1

a temperatura do plano situado na abscissa

x 1

, e

T 2

a temperatura do plano situado na abscissa

x 2

. Note-se que

e = x 2 - x 1

é a espessura da parede. Em regime estacionário, T é uma função linear de

x

, assim:

$T = T_1 + \frac{x-x_1}{e}(T_2 - T_1)$

A densidade de fluxo de calor da superfície é escrita como:

$\varphi= -\lambda \frac{dT}{dx} = \frac{\lambda}{e} (T_1-T_2)$ .

O fluxo de calor através de uma superfície S é:

$\Phi = \frac{\lambda S}{e} (T_1-T_2)$

$\Phi = \frac{(T_1-T_2)}{\frac{e}{\lambda S}}$

Analogia com resistência elétrica.

Analogia elétrica
Por analogia com a eletricidade (lei de Ohm) no caso especial onde a área de contato entre cada material é constante (fluxo superficial $\varphi$ constante), podemos traçar um paralelo entre as duas expressões:

$U_1-U_2= RI\,$

$T_1-T_2= \frac{e}{\lambda S} \Phi\,$

Pode-se traçar um paralelo da primeira com a tensão e a temperatura, da outra parte a intensidade e o fluxo de calor :

$U_1-U_2 \leftrightarrow T_1-T_2,$

$I \leftrightarrow \Phi,$

Pode-se então definir uma resistência térmica, desempenhando a transferência de calor um papel semelhante ao da resistência elétrica.

$R \leftrightarrow R_{thc}= \frac{e}{\lambda S}\,$

onde

S

é a superfície do material e

e

a espessura. A resistência térmica

R t h c

é homogênea e expressa em K.W^-1

Surperfícies planas em série

Considere-se materiais A,B e C de espessura respectivamente

e A

e B

e C

e de condutividade térmica respetivamente

λ A

λ B

λ C

.
Os pressupostos são idênticos aos de uma única superfície plana. Consideramos que o contacto entre cada camada é perfeito, o que significa que a temperatura na interface entre dois materiais é idêntica em cada material (temperatura não salta na passagem por uma interface).
Finalmente, a área de contato entre cada material é constante, o que implica um fluxo de superfície $\varphi$ constante.
As resistências térmicas se adicionam:

$T_1- T_4= (\frac{e_A}{\lambda_A S}+ \frac{e_B}{\lambda_B S}+ \frac{e_C}{\lambda_C S}) \Phi\ = (R_{thA}+ R_{thB}+ R_{thC}) \Phi$

Demonstração

Globalmente, tem-se

$T_1- T_4= \frac{e}{\lambda S} \Phi\,$

Se decomposto

para a camada A : $T_1- T_2= \frac{e_A}{\lambda_A S} \Phi\,$

para a camada B : $T_2- T_3= \frac{e_B}{\lambda_B S} \Phi\,$

para a camada C : $T_3- T_4= \frac{e_C}{\lambda_C S} \Phi\,$

Nota: Dadas as hipóteses, o fluxo (ou a densidade de fluxo) permanece constante.
Com:

$T_1- T_4= (T_1-T_2)+(T_2-T_3)+(T_3-T_4)\,$

Portanto

$T_1- T_4= (\frac{e_A}{\lambda_A S}+ \frac{e_B}{\lambda_B S}+ \frac{e_C}{\lambda_C S}) \Phi\,$

$T_1- T_4= (R_{thA}+ R_{thB}+ R_{thC}) \Phi\,$

O perfil de temperaturas
Para cada material a variação da temperatura segue uma lei do tipo:

$T= T_1- \frac{e_X}{\lambda_X S}\Phi\,$

A variação de temperatura é linear na espessura do material considerado. A inclinação é inversamente proporcional a λ (condutividade térmica) característica de cada material. Quanto maior a condutividade térmica e menor o material (e, portanto, menor o isolamento), menor será a inclinação.

Analogia elétrica
Da mesma forma que as resistências elétricas em série são adicionadas, as resistências térmicas em série são somadas.

Superfícies planas em paralelo

Considere-se materiais planos justapostos lado a lado. Cada material é homogêneo e delimitado por dois planos paralelos. Este poderia ser o modelo de uma parede com uma janela. As hipóteses são idênticas as de uma superfície plana única. Além disso, considera-se que a temperatura da superfície é uniforme em cada elemento (T₁ et T₂). Considere-se S_A, S_B e S_C as respectivas superfícies dos elementos A, B e C.
Posteriormente, presume-se que o fluxo seja sempre perpendicular à parede composta; o que não é realista porque a temperatura da superfície de cada elemento da composição é diferente e há, portanto, um gradiente de temperatura lateral (a origem das pontes térmicas). Além disso, é necessário corrigir o fluxo de calor calculado na parede feita com os coeficientes de perda por unidade de comprimento, específicos para cada ramo de parede (e pode ser insignificante, cf. a regulamentação térmica RT 2000).
As condutâncias térmicas se adicionam:

$C_{th}= \frac{1}{R_{th}} = \frac{1}{\frac{e_A}{\lambda_A S_A}}+\frac{1}{\frac{e_B}{\lambda_B S_B}}+\frac{1}{\frac{e_C}{\lambda_C S_C}} = \frac{1}{R_{thA}} + \frac{1}{R_{thB}} + \frac{1}{R_{thC}}$

Demonstração

Para cada elemento , o fluxo é expresso em função da relação

$T_1- T_2= \frac{e_X}{\lambda_X S} \Phi\,$

Tendo-se em conta a analogia elétrica

$R_X= \frac{e_X}{\lambda_X S_X}\,$

como

X

é igual a

A

B

C

Tem-se, portanto

$\Phi_A= \frac{T_1- T_2}{R_A}\,$

$\Phi_B= \frac{T_1- T_2}{R_B}\,$

$\Phi_C= \frac{T_1- T_2}{R_C}\,$

O fluxo total é igual à soma dos fluxos em cada elemento

$\Phi= \Phi_A+ \Phi_B+ \Phi_C\,$

$\Phi= (T_1- T_2) (\frac{1}{R_A}+ \frac{1}{R_B}+ \frac{1}{R_C})\,$

Sendo S a surperfície total

$S= S_A+ S_B+ S_C\,$

O fluxo na superfície é então

$\varphi=\frac{\Phi}{S}\,$

Mais uma vez, por analogia com as leis elétricas, o inverso da resistência térmica é às vezes chamado condutividade térmica.

$C_{th}= \frac{1}{R_{th}} = \frac{1}{\frac{e_A}{\lambda_A S_A}}+\frac{1}{\frac{e_B}{\lambda_B S_B}}+\frac{1}{\frac{e_C}{\lambda_C S_C}} = \frac{1}{R_{thA}} + \frac{1}{R_{thB}} + \frac{1}{R_{thC}}$

Analogia elétrica
Também é possível fazer uma analogia entre um circuito elétrico de resistências em paralelo.


$I = (\frac{1}{R_1}+ \frac{1}{R_2}+ \frac{1}{R_3}) \Delta U\,$	$\Phi = (\frac{1}{R_{th1}}+ \frac{1}{R_{th2}}+ \frac{1}{R_{th3}}) \Delta T\,$

Surperfície cilíndrica simples

O tubo simples consiste de um material único e homogêneo. A temperatura é uniforme em cada superfície do tubo. Considera-se o tubo infinitamente longo para superar os efeitos da borda.
A variação de temperatura é escrito:

$\ T_1-T_2= \frac{\Phi}{2 \pi \lambda L } \ln \frac{R_2}{R_1}$

Demonstração

Considera-se uma variação dR dentro do material constituinte do tubo, a lei de Fourier é então expressa como:

$\Rightarrow \Phi= - \lambda S \frac{dT}{dR}\,$

Variação da temperatura na espessura do tubo

Evolução da temperatura na espessura de um tubo simples com T_interior > T_exterior.

Sendo S a superfície de um cilindro:

$S= 2 \pi R L\,$

Podemos escrever a lei de Fourier na forma:

$\Phi= - \lambda 2 \pi R L \frac{dT}{dR}\,$

$\frac{dR}{R}= - \frac{2 \pi \lambda L dT}{\Phi}\,$

$\int_{R_1}^{R_2} \frac{dR}{R}= - \frac{2 \pi \lambda L }{\Phi} \int_{T_1}^{T} dT\,$

$\ln \frac{R_2}{R_1}= \frac{2 \pi \lambda L }{\Phi} (T_1-T)\,$

A variação de temperatura no material é

$\ T= T_1 - \frac{\Phi}{2 \pi \lambda L } \ln \frac{R}{R_1}\,$

Em toda a espessura do tubo, a variação é

$\ T_1-T_2= \frac{\Phi}{2 \pi \lambda L } \ln \frac{R_2}{R_1}\,$

Superfícies cilíndricas concentricas

Esquema de um tubo concêntrico.

O tubo concêntrico é composto por tubos dispostos em camadas concêntricas . Considere-se que o contato é perfeito entre os tubos. A temperatura é uniforme em cada superfície do tubo. Considere-se que o tubo possui um comprimento interminável L para superar os efeitos de borda.
A resistência total do tubo é expressa por uma lei da tipo "série como a parede compondo uma série:

$\ R_{thT}= R_{thA} + R_{thB}$

Demonstração

Evolução da temperatura na primeira camada:

$\ T_1-T_2= \frac{\Phi}{2 \pi \lambda_A L } \ln \frac{R_2}{R_1}\,$

Evolução da temperatura na segunda camada:

$\ T_2-T_3= \frac{\Phi}{2 \pi \lambda_B L } \ln \frac{R_3}{R_2}\,$

Em toda a espessura do tubo:

$\ T_1-T_3= \frac{\Phi}{2 \pi L } \left ( \frac {\ln \frac {R_2}{R_1} }{\lambda_A} + \frac {\ln \frac {R_3}{R_2} }{\lambda_B}\right )\,$

A resistência térmica da camada A

$\ R_{thA}= \frac {\ln \frac {R_2}{R_1} }{\lambda_A}\,$

A resistência térmica da camada B

$\ R_{thB}= \frac {\ln \frac {R_3}{R_2} }{\lambda_B}\,$

Chegando-se a que resistência total do tubo é expressa por uma lei da tipo "série como a parede compondo uma série:

$\ R_{thT}= R_{thA} + R_{thB}\,$

Condução em regime dinâmico

A resolução da equação do calor em regime dinâmico é muito mais difícil. Ela usa os conceitos de transformada de Fourier de produto de convolução e distribuições. Apresentam-se alguns exemplos de resolução.

Caso de uma área ilimitado

Princípio geral

Escreve-se a equação do calor na forma:
$\frac{\partial T}{\partial t} - D \Delta T = P$ onde D é o coeficiente de difusividade térmica e P representa as fontes de calor. P pode ser uma função do tempo e da posição da fonte de calor, mas também uma distribuição. Por exemplo, a injeção instantânea e pontual de uma quantidade de calor pode ser representada pelo produto

δ(t)δ(x)

de uma distribuição de Dirac no instante t = 0 por uma distribuição de Dirac em x = 0, x é a abscissa em caso de problema unidimensional ou vetor de posição no caso geral.
Ela também dará o estado inicial do área

T 0 = T (0, x)

, que também pode ser uma função de x ou uma distribuição. Considera-se que T é nulo para

t < 0

.
O método de solução consiste em^[8] ^[9] :

Aplica-se uma transformada de Fourier sobre a variável x, a todos os termos da equação diferencial. Isto transforma a derivação com relação a x em um produto. Toma-se $F(T)(p, t) = \int T(x, t) \exp(-2i\pi px) dx$ , então a equação torna-se:

$\frac{\partial F(T)}{\partial t} + D 4\pi^2p^2F(T) = F(P)$ ou melhor, no sentido de distribuições:
$\frac{\partial F(T)}{\partial t} + D 4\pi^2p^2F(T) = F(P) + F(T_0)\delta(t)$ tendo-se em conta a condição inicial.

Reconhecendo-se nesta equação um produto de convolução :

$(\frac{\partial \delta(t)}{\partial t} + 4\pi^2Dp^2 \delta(t)) * F(T) = F(P) + F(T_0)\delta(t)$ O operador aplicado a F é um produto de convolução relativo à variável t.

Aplicando o operador inverso, prova-se que $H (t)exp( - 4π 2 D p 2 t)$ , onde H é a função de Heaviside, conduzindo a:

F (T) = F (P) * H (t)exp( - 4π 2 D p 2 t) + F (T 0) H (t)exp( - 4π 2 D p 2 t)

Se F(P) é uma função e não uma distribuição, esta relação se torna, para t > 0:
$F(T) = \int_0^t F(P)(\tau)\exp(- 4\pi^2Dp^2(t-\tau)) d\tau + F(T_0)\exp(- 4\pi^2Dp^2t)$

Toma-se a transformada de Fourier inversa para deduzir T.

Caso particular

Propagação por condução no plano a partir de um ponto quente. A altura de um dado ponto indica o valor da temperatura neste ponto.

Toma-se

T 0 = 0

P = δ(t)δ(x)

(injeção instantânea de calor em um determinado ponto), o método acima descrito conduz a:

F (P) = δ(t)

Assim, para t > 0 :

F (T) = exp( - 4π 2 D p 2 t)

cuja transformada de Fourier inversa é, para t > 0 :

$T = \frac{\exp(- \frac{x^2}{4Dt})}{2\sqrt{\pi tD}}$ no caso unidimensional.

$T = \frac{\exp(- \frac{r^2}{4Dt})}{8\sqrt{\pi tD}^3}$ no caso tridimensional.

Área ilimitado sem fonte de calor

Se toma-se somente a temperatura inicial

T 0

de uma fonte média independente de calor (P = 0), tem-se que:

$T = \frac{1}{2\sqrt{\pi tD}} \int_{-\infty}^{\infty} \exp(- \frac{(x - u)^2}{4tD}) T_0(u)\, du$ no caso unidimensional.

$T = \frac{1}{8\sqrt{\pi tD}^3} \int_{\mathbb R^3} \exp(- \frac{(r - s)^2}{4tD}) T_0(r)\, dx_s dy_s dz_s$ no caso tridimensional.

Caso de áreas limitadas, sem fonte de calor

Caso de um área limitada por um plano. O problema de Kelvin

Problema de Kelvin. O eixo x é orientado para a direita. A metade do semi-espaço x > 0, cuja temperatura inicial é T₀, tem um limite x = 0, onde a temperatura é sempre zero. A altura de um dado ponto indica o valor da temperatura neste ponto

Supondo-se o área limitada pelo plano x=0. Se coloca-se por condição os limites suplementares T(0,t) = 0 para todo t, então basta estender a distribuição inicial de temperatura

T 0

por uma função ímpar em x e aplicar-se o resultado anterior.
O caso mais célebre é o problema de Kelvin. Este considerou em 1860 que a Terra estava inicialmente a uma temperatura constante

T 0

da ordem de 3000° e que se resfriava por simples condução. Utilizando o valor atual do gradiente de temperatura em função da profundidade, ele deduziu uma estimativa para a idade da Terra. Pode-se aplicar o método de resolução acima para considerar a Terra como plana e infinitamente profunda, delimitado por um plano de sua superfície. O cálculo conduz a:
$T= \frac{T_0}{\sqrt{\pi t D}} \int_0^x \exp(- \frac{u^2}{4tD})\, du = T_0 \,{\rm erf}(\frac{x}{2\sqrt{Dt}})$ onde erf é chamada função erro de Gauss.
O gradiente de temperatura na superfície será:
$\frac{\partial T}{\partial x}= \frac{T_0}{\sqrt{\pi t D}}$ Conhecendo-se $\frac{\partial T}{\partial x}$ como sendo da ordem de 3 °C para cada 100 metros de profundidade e D estimado como $10^{-6} \,m^2s^{-1}$ , obtem-se que t vale 100 milhões de anos. Este resultado é largamente subestimado pois Kelvin ignora os fenômenos de convecção no interior do manto terrestre^[10] ^[11].

Caso de um área limitada por dois planos paralelos

Transferência de calor por condução em um área limitada 0 < x < L. Os dois limites da área são mantidos a uma temperatura constante . A altura de um dado ponto indica o valor da temperatura neste ponto

Considere-se uma área limitada por dois planos x=0 e x=L. Supondo-se que são dadas como condições de contorno os limites T(0,t) = T(L,t) = 0. Utiliza-se um métode de resolução baseado nas séries de Fourier, procurando-se T sob a forma:
$T = \sum_{n=1}^\infty b_n \sin(n\pi \frac{x}{L}) \exp(- \frac{Dn^2\pi^2t}{L^2})$ Esta expressão verifica tanto a equação do calor quanto as condições de contorno nos limites. Se é tomada a distribuição de temperatura inicial

T 0

, é suficiente expandir-se em uma série de Fourier para determinar-se os coeficientes

b n

.
Por exemplo, se tiver-se

T 0

constante, obtem-se:
$T = \frac{4T_0}{\pi} \sum_{n=0}^\infty \frac{1}{2n+1} \sin(\frac{(2n+1)\pi x}{L}) \exp(- \frac{D(2n+1)^2\pi^2t}{L^2})$ Ao deixar-se L ao infinito, encontra-se a solução de Kelvin do parágrafo anterior, a soma acima é considerado uma soma de Riemann convergindo para a integral.

Caso de uma área de geometria esférica

Transferência de calor por condução em uma área circular. O problema de Kelvin: a temperatura inicial é uniforme, a temperatura no limite do círculo é mantida como zero. A altura de um dado ponto indica o valor da temperatura neste ponto

Se o caso de propagação se sá em uma área esférica, e onde a temperatura não depede apenas das distância r ao centro, a equação de calor torna-se a expressão do laplaciano em uma esfera:
$\frac{\partial T}{\partial t} = D(\frac{2}{r} \frac{\partial T}{\partial r} + \frac{\partial^2 T}{\partial r^2})$ Se coloca-se

F = r T

, a equação torna-se:
$\frac{\partial F}{\partial t} = D \frac{\partial^2 F}{\partial r^2}$ Podemos, então, aplicar os métodos anteriores para determinar F, então deduzir T dividindo por r.
Assim, resolvendo-se o problema de Kelvin, no caso de uma esfera de raio R (temperatura inicial uniforme igual à

T 0

,a superfície é mantida a uma temperatura igual a zero) conduz à seguinte expressão para T:
$T(r,t) = 2T_0 \sum_{n=1}^\infty (-1)^{n+1} \, {\rm sinc}(n\pi \frac{r}{R}) \exp(- \frac{Dn^2\pi^2t}{R^2})$ onde sinc é a função sinus cardinal.

Caso de áreas limitadas, com fonte de calor

Considere-se a equação:
$\frac{\partial T}{\partial t} - D \Delta T = P$ com P não nulo. Em geral, procura-se uma solução particular para esta equação, de modo que, uma vez relacionada a T, pode ser reduzida a uma equação sem segundo membro. Apresentam-se alguns exemplos, onde P representa uma densidade de fonte de calor constante, independente da posição e do tempo.

Área limitada por dois planos paralelos

Transferência de calor por condução em um domínio limitado 0 < x < L. A área é aquecida uniformemente, enquanto que a temperatura inicial é zero, e que as duas bordas permanecem na temperatura zero. A altura de um dado ponto indica o valor da temperatura nesse ponto.

Considere-se uma área limitada pelos dois planos x=0 e x=L. Supondo-se que no momento inicial, a temperatura da área é igual a uma temperaturea de referência nula, e que as bordas da área se manterão permanentemente esta temperatura nula. T satisfaz, pois:

$\frac{\partial T}{\partial t} - D \frac{\partial^2 T}{\partial x^2} = P$

T(0,t) = T(L,t) = 0 para todo t positivo.

T(x,0) = 0 para todo x entre 0 e L.

A função $\frac{Px(L-x)}{2D}$ independente de t satisfaz as duas primeiras relações, de maneira que coloca-se que $G = T - \frac{Px(L-x)}{2D}$ , então G satisfaz :

$\frac{\partial G}{\partial t} - D \frac{\partial^2 G}{\partial x^2} = 0$

G (0, t) = G (L, t) = 0

$G(x, 0) = - \frac{Px(L-x)}{2D}$

Pode-se aplicar o método visto acima na busca de G sob a forma de uma série:
$G = \sum_{n=1}^\infty b_n \sin(\frac{n\pi x}{L}) \exp(- \frac{n^2\pi^2Dt}{L^2})$ que satisfaça as duas primeiras relações. Dado que, por razões de simetria, é esperado que

G (x) = G (L - x)

, pode-se supor que os coeficientes

b n

sejam nulos quando n é par, de maneira que:
$G = \sum_{n=0}^\infty b_{2n+1} \sin(\frac{(2n+1)\pi x}{L}) \exp(- \frac{(2n+1)^2\pi^2Dt}{L^2})$ Para t = 0, tem-se:
$- \frac{Px(L-x)}{2D} = \sum_{n=0}^\infty b_{2n+1} \sin(\frac{(2n+1)\pi x}{L})$ Obtem-se

b 2 n + 1

em desenvolver-se $- \frac{Px(L-x)}{2D}$ em série de Fourier. Encontrando-se:
$b_{2n+1} = - \frac{4PL^2}{(2n+1)^3\pi^3D}$ Assim G, finalmente:
$T = \frac{Px(L-x)}{2D} - \frac{4PL^2}{D\pi^3} \sum_{n=0}^\infty \frac{1}{(2n+1)^3} \sin(\frac{(2n+1)\pi x}{L}) \exp(- \frac{(2n+1)^2\pi^2Dt}{L^2})$ Quando t tende ao infinito, a temperatura da área tende a $\frac{Px(L-x)}{2D}$ , o aquecimento tende a se equilibrar com a perda de calor pelas duas bordas.

Área limitada por um plano

Propagação de calor por condução em uma área x > 0 limitada por uma borda. A área é aquecida uniformemente, enquanto que a temperatura inicial é nula, e que a borda continua a esta temperatura nula. A altura de um dado ponto indica o valor da temperatura nesse ponto.

A resolução do mesmo problema do caso em que x>0 consiste em determinar T tal que:

$\frac{\partial T}{\partial t} - D \frac{\partial^2 T}{\partial x^2} = P$

T(0,t) = 0 para todo t positivo.

T(x,0) = 0 para todo x > 0.

Pode ser obter a solução deixando L infinito na expressão dada no parágrafo anterior, igualando a série a uma soma de Riemann. Isto resulta na seguinte expressão:
$T = - \frac{Px^2}{2D} + \frac{Px^2}{2D} \,{\rm erf}(\frac{x}{2\sqrt{Dt}}) + \frac{Px\sqrt{t}}{\sqrt{D\pi}} \exp(- \frac{x^2}{4Dt}) + Pt \, {\rm erf}(\frac{x}{2\sqrt{Dt}})$ erf é a função chamada função erro de Gauss. Pode-se também encontrar esta expressão pela aplicação do método derivado do princípio geral relativo a uma área ilimitada, depois de estender a todo o espaço as funções T e P em duas funções ímpares em x, de fmodo que T se anula em x = 0.
Quando t tende ao infinito, T é de cerca de Pt, análogo aquela de uma área infinita. A borda única não é suficiente para dissipar o calor.

Área de uma geometria esférica

Propagação de calor por condução em uma área de borda circular. A área é aquecida uniformemente, enquanto que a temperatura inicial é zero, e que a borda continua a estar a esta temperatura zero. A altura de um dado ponto indica o valor da temperatura nesse ponto.

Dado o caso de uma área cuja borda seja uma esfera de raio R, utiliza-se a expressão do Laplaciano em esfera e tem-se de resolver:

$\frac{\partial T}{\partial t} = D\left(\frac{2}{r} \frac{\partial T}{\partial r} + \frac{\partial^2 T}{\partial r^2}\right) + P$

Para todo t, T(R,t) = 0

Para todo r, T(r,0) = 0

Pergunta-se se $G = rT + \frac{r^3P - rR^2P}{6D}$ , G satisfaz o sistema:

$\frac{\partial G}{\partial t} = D\frac{\partial^2 G}{\partial r^2}$

Para todo t, G(R,t) = 0

Para todo r, $G(r,0) = \frac{r^3P - rR^2P}{6D}$

O método das séries de Fourier sugere a busca de G sob a forma de uma série $\sum_{n=1}^\infty b_n \sin(\frac{n\pi r}{R}) \exp(- \frac{n^2\pi^2 Dt}{R^2})$ , onde os coeficientes

b n

encontram-se em desenvolvimento $\frac{r^3P - rR^2P}{6D}$ em série de Fourier. Obtendo-se:
$G = \frac{2P R^3}{D\pi^3} \sum_{n=1}^\infty \frac{(-1)^n}{n^3} \sin(\frac{n\pi r}{R}) \exp(- \frac{n^2\pi^2 Dt}{R^2})$ e, portanto:
$T = \frac{R^2P - r^2P}{6D} + \frac{2P R^2}{D\pi^2} \sum_{n=1}^\infty \frac{(-1)^n}{n^2} {\rm sinc}(\frac{n\pi r}{R}) \exp(- \frac{n^2\pi^2 Dt}{R^2})$ onde sinc é a função sinus cardinal.
Quando t tende ao infinito, a temperatura T tende à distribuição limite $\frac{R^2P - r^2P}{6D}$

estudar é preciso

quinta-feira, 23 de setembro de 2010

A condução térmica

Condução térmica

Lei de Fourier

Teorização completa

Forma diferencial

Forma integral

Condutância

Representação de propriedade intensiva

Cilindros

Condução em regime estacionário

Superfície plana simples

Surperfícies planas em série

Demonstração

Superfícies planas em paralelo

Demonstração

Surperfície cilíndrica simples

Demonstração

Superfícies cilíndricas concentricas

Demonstração

Condução em regime dinâmico

Caso de uma área ilimitado

Princípio geral

Caso particular

Área ilimitado sem fonte de calor

Caso de áreas limitadas, sem fonte de calor

Caso de um área limitada por um plano. O problema de Kelvin

Caso de um área limitada por dois planos paralelos

Caso de uma área de geometria esférica

Caso de áreas limitadas, com fonte de calor

Área limitada por dois planos paralelos

Área limitada por um plano

Área de uma geometria esférica

Nenhum comentário:

Postar um comentário