(5.59 Incropera) Uma esfera de 80 mm de diâmetro ( k=50 \mathrm{~W} / \mathrm{m} . \mathrm{K} e \alpha=1,5.10^{-6} \mathrm{~m}^{2} / \mathrm{s} ) está inicialmente a uma temperatura uniforme elevada e é imersa em um banho de óleo mantido a 50^{\circ} \mathrm{C} . \mathrm{O} coeficiente de convecção para o processo de resfriamento é 1000 \mathrm{~W} / \mathrm{m}^{2} . \mathrm{K} . \mathrm{Em} um dado instante, a temperatura da superfície da esfera é de 150^{\circ} \mathrm{C}. Qual a temperatura correspondente no centro da esfera?
Passo 1
O que o enunciado dá e o que pede (resumindo e pegando todos os números):
Diâmetro da esfera: 80\\ \text{mm} → raio R=40\\ \text{mm}=0{,}04\\ \text{m}.
Regra prática: se \mathrm{Bi}<0{,}1 usa-se lumped capacitance. Aqui \mathrm{Bi}\approx0{,}267>0{,}1, logo a temperatura interna não é praticamente uniforme — não podemos usar o método lumped. Precisamos da solução de condução transiente no interior da esfera com convecção na superfície.
Passo 3
Apresentar a forma da solução (separação de variáveis) e obter a equação dos autovalores.
Para condução transiente em uma esfera com condição de convecção na superfície, através de separação de variáveis obtemos modos espaciais da forma (apenas a ideia necessária aqui):
\Phi_n(r)=\frac{\sin(\mu_n r/R)}{\mu_n r/R},
onde \mu_n são números adimensionais (autovalores) que satisfazem a condição de contorno na superfície. A condição de convecção -k\left.\dfrac{\partial T}{\partial r}\right|_{r=R}=h\left(T(R,t)-T_{\infty}\right) leva, para os autovalores, à equação transcendental
\mu_n\cot\mu_n=1-3\mathrm{Bi}.
Explicando os símbolos:
\mu_n é o n-ésimo autovalor adimensional (o primeiro autovalor chamamos \mu_1).
\mathrm{Bi} é o número de Biot já calculado.
A forma \Phi_n(r) dá o perfil espacial do modo n; em r=0 tem-se \Phi_n(0)=1 (limite bem definido) e em r=R tem-se \Phi_n(R)=\dfrac{\sin\mu_n}{\mu_n}.
Passo 4
Relação entre temperatura no centro e na superfície usando a aproximação de um único modo (primeiro termo).
A solução transiente é combinação de modos multiplicados por decaimentos exponenciais em função do tempo. Quando o primeiro modo domina (após tempos não muito pequenos) a solução pode ser aproximada por um único termo:
No centro (r=0) temos \Phi_1(0)=1. Na superfície (r=R) temos \Phi_1(R)=\dfrac{\sin\mu_1}{\mu_1}. Assim, dividindo as duas expressões (o fator temporal e A_1 se cancelam), obtemos a relação útil, independente do instante específico:
Observação didática: isso é uma consequência direta do fato de que o perfil espacial do primeiro modo tem valor 1 no centro e \sin\mu_1/\mu_1 na superfície; portanto, conhecendo T(R,t) e \mu_1 consegue-se T(0,t) sem precisar do tempo nem da temperatura inicial — desde que o primeiro modo domine razoavelmente (esta é a aproximação padrão em problemas como o do livro de Incropera).
Temos que achar \mu tal que \mu\cot\mu=0{,}2. Fazendo uma busca numérica (bisseção/interpolação) encontra-se o primeiro raiz positiva aproximadamente
\mu_1\approx1{,}432.
(Resumo rápido do procedimento numérico: avaliou-se a função f(\mu)=\mu\cot\mu-0{,}2 em vários pontos e interpolou-se entre os sinais opostos até chegar ao valor acima.)
Agora calculamos \sin\mu_1. Com \mu_1\approx1{,}432 temos numericamente \sin\mu_1\approx0{,}991 (valor usado abaixo com precisão compatível).
Arredondando de maneira razoável, a temperatura no centro é aproximadamente 195^\circ\mathrm{C}.
Observação final e encorajamento: usei a aproximação de primeiro modo — é a abordagem que o livro de Incropera e problemas de prova costumam pedir quando se dá a temperatura da superfície num instante (essa relação centro/superfície é padrão). Se quiser, eu posso fazer uma busca numérica mais refinada para \mu_1 e dar o valor com mais casas decimais, ou discutir a precisão dessa aproximação (quando os modos superiores são negligenciáveis).