Enunciado

Desenvolva a equação da difusão térmica, Equação2.27, em coordenadas esféricas, a partir do volume de controle diferencial mostrado na Figura 2.13.

Passo 1

E ai meus amores nessa questãozinha “Pegue o que puder, sem nada a devolver” okay?

Vamos aplicar o balanço de energia ao sistema

E ˙ i n + E ˙ g - E ˙ o u t = E ˙ s t

Onde, a entrada de energia é E ˙ i n

A saída de energia é E ˙ o u t

A taxa de geração de energia térmica é E ˙ g

E a taxa de armazenamento de energia é E ˙ s t

Considerando a expressão para o fluxo de calor, teremos

E ˙ i n = q r + q ϕ + q z

Onde, as taxas de calor por condução na direção r ,   ϕ   e   z da face de entrada do elemento são q r ,   q ϕ   e   q z respectivamente.

Para a saída de calor, teremos a seguinte expressão

E ˙ o u t = q r + d r + q ϕ + d ϕ + q z + d z

Onde, as taxas de calor por condução na direção r ,   ϕ   e   z da face de entrada do elemento são q r + d r ,   q ϕ + d ϕ   e   q z + d z respectivamente.

Passo 2

Vamos considerar a expressão para a taxa de geração de energia térmica

E ˙ g e n = q d V

Considerando a expressão para o armazenamento de energia

ρ c p ∂ T d t d V

Agora vamos substituir as expressões que encontramos no passo 2, show?

q r + q ϕ + q z + q ˙ d V - q r + d r + q ϕ + d ϕ + q z + d z = ρ c p ∂ T d t d V

Passo 3

Agora iremos escrever a taxa de calor por condução na direção de r em um material isotrópico pela lei de Fourier

q r = - k d A r ∂ T ∂ r

= - k r d θ r sin ⁡ θ d ϕ ∂ T ∂ r

= - k r 2 sin ⁡ θ d θ d ϕ ∂ T ∂ r

Onde a área da face na direção de r é d A r e o comprimento do volume elementar nas direções ϕ   e   θ , são r s i n   θ d ϕ   e   r d θ .

Vamos escrever a taxa de calor por condução na direção θ , em um material isotrópico pela lei de Fourier

q θ = - k d A θ ∂ T ∂ r θ

= - k d A r ∂ T ∂ θ

= - k r r sin ⁡ θ d θ d r ∂ T ∂ θ

Vamos escrever a taxa de calor por condução na direção ϕ , em um material isotrópico pela lei de Fourier

q ϕ = - k ( d A ) ∂ T d ( r sin ⁡ θ ϕ )

= - k d A ∂ T r sin ⁡ θ ∂ ϕ

= - k ( r d θ d r ) r s i n   θ ∂ T ∂ ϕ

Passo 4

Vamos escrever a expressão para o volume do elemento diferencial

d V = r d ϕ d r d z

  = ( d r ) ( r d θ   ) ( r   s i n θ d   ϕ )  

= r 2 sin ⁡ θ d r d θ d ϕ

Usaremos a expansão da série de Taylor na expressão a seguir, e desprezaremos os termos superiores.

E ˙ o u t = q r + d r + q ϕ + d ϕ + q z + d z

q r + d r + q ϕ + d ϕ + q z + d z = q r + ∂ q ∂ r d r + q ϕ + ∂ q ∂ ϕ d ϕ + q z + ∂ q ∂ z d z

Passo 5

Vamos diferenciar as seguintes equações, para obter as derivadas do fluxo de calor em diferentes direções.

= - k r 2 sin ⁡ θ d θ d ϕ ∂ T ∂ r

= - k r r sin ⁡ θ d θ d r ∂ T ∂ θ

= - k ( r d θ d r ) r s i n   θ ∂ T ∂ ϕ

Diferenciando

- ∂ q r d r = - ∂ ∂ r - k r 2 sin ⁡ θ d θ d ϕ ∂ T ∂ r

= sin ⁡ θ d ϕ d θ ∂ ∂ r k r 2 ∂ T ∂ r  

- ∂ q θ d θ = - ∂ ∂ θ - k r r sin ⁡ θ d ϕ d r ∂ T ∂ θ

= d ϕ d r ∂ ∂ θ k s i n   θ ∂ T ∂ θ  

- ∂ q ϕ ∂ ϕ = - ∂ ∂ ϕ - k d r d z sin ⁡ θ ∂ T ∂ ϕ

= d r d z sin ⁡ θ ∂ ∂ ϕ k ∂ T ∂ ϕ  

Passo 6

Manipulando essa outra expressão, teremos

q r + q ϕ + q z + q ˙ d V - q r + d r + q ϕ + d ϕ + q z + d z = ρ c p ∂ T d t d V

d r s i n   θ d θ d ϕ ∂ q ∂ r k r 2 ∂ T ∂ r + ( d ϕ ) ( d r d θ ) ∂ ∂ θ k s i n   θ ∂ T ∂ θ + d ϕ d r d z sin ⁡ θ ∂ ∂ ϕ k ∂ T ∂ ϕ   + q ˙ r 2 sin ⁡ θ d ϕ d r d θ

= ρ c p ∂ T ∂ t   r 2 sin ⁡ θ d ϕ d r d θ

Passo 7

Agora pra finalizar o exercício iremos substituir essas 3 equações diferenciadas que encontramos na seguinte fórmula

= ρ c p ∂ T ∂ t   r 2 sin ⁡ θ d ϕ d r d θ

Logo, vamos dividir essa expressão por r 2 sin ⁡ θ d ϕ d r d θ

1 r 2 ∂ ∂ r k r 2 ∂ T ∂ r + 1 r 2 sin 2 ⁡ θ ∂ ∂ ϕ k ∂ T ∂ ϕ + 1 r 2 sin ⁡ θ ∂ ∂ θ k s i n   θ ∂ T ∂ θ + q ˙ = ρ c p ∂ T ∂ t

Sendo assim, a equação de difusão de calor com geração de calor em coordenadas cilíndricas foi derivada.

Depois dessa questão, “Ninguém se mexe” Meu cérebro caiu.”

Haha bora pra proximaaaa

Resposta

Resposta nos passos.

Exercícios de Livros Relacionados

Desenvolva a equação da difusão térmica, Equação 2.24, em coordenadas cilíndricas, a partir do volume de controle diferencial mostrado na Figura 2.12.
Considere uma janela com 300m m × 300m m em um avião. Para uma diferença de temperaturas de 80 ° C entre as superfícies interna e externa da janela, calcule a perda térmica através de janelas com L =
Um aquecedor elétrico delgado, dissipando 4000W / m 2 , encontra-se imprensado entre duas placas, com espessura de 25 mm, cujas superfícies expostas trocam calor por convecção com um fluido a T ∞ = 20
Um aquecedor elétrico delgado, dissipando 4000W / m 2 , encontra-se imprensado entre duas placas, com espessura de 25 mm, cujas superfícies expostas trocam calor por convecção com um fluido a T ∞ = 20
Um anúncio de TV veiculado por um fabricante conhecido de isolamentos térmicos afuma que não é a espessura do material isolante que conta, mas sim o seu valor de R. O anúncio mostra que, para obter um