Qual é a interpretação física do número de Fourier?
Que exigência deve ser satisfeita para se usar a aproximação de primeiro termo para determinar a resposta térmica transiente de uma parede piana, de um cilindro longo e de uma esfera nos quais há condução unidimensional devido a uma mudança nas condições na superfície? Em que estágio do processo transiente a exigência não é satisfeita?

Passo 1

O enunciado faz duas perguntas principais:
  1. Qual é a interpretação física do número de Fourier?
  2. Que exigência deve ser satisfeita para que a aproximação de primeiro termo possa ser usada para obter a resposta térmica transiente de uma parede plana, de um cilindro longo e de uma esfera (condução unidimensional causada por mudança nas condições na superfície)? Em que estágio do processo transiente essa exigência não é satisfeita?
Informações relevantes dadas pelo enunciado (resumidas, para não perder nada):
  • Trata-se de condução térmica transiente unidimensional em três geometrias: parede plana (slab), cilindro longo e esfera.
  • A situação física é uma mudança nas condições na superfície (por exemplo, mudança na temperatura ou convecção na superfície).
  • Queremos saber quando é válido usar a aproximação de primeiro termo (isto é, manter só o primeiro modo/autovalor da solução em série).

Passo 2

Definição e fórmula do número de Fourier. Vou apresentar a definição, explicar cada símbolo e dar a forma comum usada em problemas de condução transiente.
A definição usual do número de Fourier é:
\mathrm{Fo}=\frac{\alpha \, t}{L_c^{2}}
onde:
  • \alpha é a difusividade térmica do material, definida por \alpha=\dfrac{k}{\rho c_p}, com k condutividade térmica, \rho densidade e c_p calor específico (tudo em unidades do SI);
  • t é o tempo transcorrido desde a mudança nas condições de contorno;
  • L_c é um comprimento característico da geometria do corpo (explico exemplos abaixo).
Exemplos de L_c usados frequentemente:
  • Para uma parede plana de espessura total 2L que tem simetria no centro (convecção nas faces externas), usa-se L_c=L (a meia-espessura).
  • Para um cilindro longo de raio r_0, usa-se L_c=r_0.
  • Para uma esfera de raio r_0, usa-se L_c=r_0.
    (Alternativamente, para formas arbitrárias pode-se usar L_c=V/A — volume sobre área superficial.)
Interpretação física direta do número de Fourier:
  • \mathrm{Fo} é um tempo adimensional que compara o tempo real t com o tempo característico de difusão térmica t_{\text{diff}}=L_c^{2}/\alpha. Ou seja,
\mathrm{Fo}=\frac{t}{t_{\text{diff}}}.
  • Física: \mathrm{Fo} mede quantas “tempos de difusão” já se passaram. Para \mathrm{Fo}\ll 1 a difusão térmica ainda não “atingiu” grande parte do corpo (efeito superficial dominando); para \mathrm{Fo}\gtrsim 1 a difusão já percorreu uma distância da ordem de L_c e o corpo inteiro sente a mudança.
  • Outra intuição: a profundidade de penetração térmica cresce como \delta \sim \sqrt{\alpha t}. Escrevendo \delta em relação a L_c \delta/L_c\sim\sqrt{\mathrm{Fo}}. Assim \mathrm{Fo} também indica a fração da dimensão do corpo que já foi “afetada”.

Passo 3

Como surgem as séries e o significado da “aproximação de primeiro termo”. Vou explicar de onde vem a solução em série (sem desenvolver toda a dedução técnica), o papel dos autovalores e por que os termos mais altos decaem rapidamente.
Quando se resolve o problema transiente linear por separação de variáveis (ou por métodos equivalentes) para uma parede/cilindro/esfera com condições convectivas na superfície, a solução para a temperatura adimensionalizada \theta (por exemplo \theta=(T-T_\infty)/(T_i-T_\infty)) geralmente aparece como uma soma infinita:
\theta(\text{posição},t)=\sum_{n=1}^{\infty} A_n\,\Phi_n(\text{posição})\,\exp\!\big(-\lambda_n^{2}\,\mathrm{Fo}\big),çãçã
onde:
  • A_n são coeficientes que dependem da condição inicial e dos autovetores; esses coeficientes dependem do número de Biot ( \mathrm{Bi}) e da forma geométrica;
  • \Phi_n(\text{posição})çã são as funções espaciais (autofunções) associadas ao modo n;
  • \lambda_n são os autovalores positivos, ordenados tal que 0<\lambda_1<\lambda_2<\lambda_3<\dots;
  • o fator \exp(-\lambda_n^{2}\mathrm{Fo}) descreve o decaimento temporal de cada modo; note o quadrado de \lambda_n.
Observações importantes (didáticas):
  • Os autovalores \lambda_n são raízes de equações transcendentais que dependem do \mathrm{Bi}. Por exemplo, para a parede plana com convecção em ambas as faces os \lambda_n satisfazem a equação
\lambda_n \tan(\lambda_n)=\mathrm{Bi}.
  • Para cilindro e esfera os autovalores são raízes de equações envolvendo funções de Bessel ou funções especiais — o ponto crucial é que eles são discretos e crescentes.
  • Como \lambda_n cresce com n, o argumento do exponencial \lambda_n^{2}\mathrm{Fo} cresce muito com n, o que faz com que os termos de ordem alto decaiam muito mais rápido com o tempo.

Passo 4

Condição matemática para que a aproximação de primeiro termo seja válida. Vou escrever a condição geral e explicar cada parte.
A aproximação de primeiro termo consiste em truncar a soma e manter só o termo n=1:
\theta(\text{posição},t)\approx A_1\,\Phi_1(\text{posição})\,\exp\!\big(-\lambda_1^{2}\,\mathrm{Fo}\big).çãçã
Para que isso seja aceitável com uma precisão desejada \varepsilon (por exemplo \varepsilon=0{,}01 para 1%), os termos negligenciados (a soma dos termos n\ge2) devem ser muito menores que o primeiro termo. Tomando apenas o segundo termo como estimativa do erro, pedimos que
\left|\frac{A_2\,\Phi_2\,\exp(-\lambda_2^{2}\mathrm{Fo})}{A_1\,\Phi_1\,\exp(-\lambda_1^{2}\mathrm{Fo})}\right|\ll \varepsilon.
Isolando \mathrm{Fo} obtemos a forma útil:
\exp\!\big[-(\lambda_2^{2}-\lambda_1^{2})\,\mathrm{Fo}\big]\ll \varepsilon\,\frac{|A_1\Phi_1|}{|A_2\Phi_2|}.
Tomando logaritmo (e invertendo o sinal porque é exponencial decrescente),
\mathrm{Fo}\gg \frac{\ln\!\big(\dfrac{|A_1\Phi_1|}{\varepsilon\,|A_2\Phi_2|}\big)}{\lambda_2^{2}-\lambda_1^{2}}.
Interpretação e comentários:
  • Os termos A_n e \Phi_n dependem de \mathrm{Bi} e da geometria; portanto o valor mínimo de \mathrm{Fo} necessário para uma tolerância \varepsilon depende do \mathrm{Bi} e da geometria.
  • A forma acima dá um critério explícito: para uma dada tolerância \varepsilon calcule os dois primeiros autovalores \lambda_1,\lambda_2 e os coeficientes A_1,A_2 (há tabelas e fórmulas nos livros e gráficos de Heisler); então aplique a desigualdade para obter o \mathrm{Fo} mínimo necessário.
  • Se só queremos um critério qualitativo: como \lambda_2^{2}-\lambda_1^{2} é positivo e frequentemente de ordem 1 ou maior, basta que \mathrm{Fo} seja suficientemente grande para que o exponencial do segundo termo seja pequeno.

Passo 5

Regras práticas, regras de bolso e onde a exigência falha no processo transiente — explico como usar isso numa prova ou num problema prático.
Regras práticas (nível “como faço rápido na prova”):
  • Primeiro calcule o número de Biot:
\mathrm{Bi}=\frac{h\,L_c}{k},
onde h é o coeficiente convectivo e k é a condutividade do material. Isso indica se há grandes gradientes espaciais ( \mathrm{Bi}\gg 0{,}1) ou não ( \mathrm{Bi}\ll 0{,}1).
  • Se \mathrm{Bi}\le 0{,}1 frequentemente usa-se o modelo de corpo perfeitamente bem-misturado (lumped capacitance) e nem é preciso a série. Nesse caso a aproximação de primeiro termo por uma série espacial não é o caminho natural.
  • Se \mathrm{Bi} é moderado/alto (há gradientes), use a solução em série. A aproximação de primeiro termo será válida quando \mathrm{Fo} não for muito pequeno — ou seja, quando já passou o estágio inicial no qual apenas uma camada superficial foi afetada.
Regras de bolso numéricas (valores típicos usados em engenharia):
  • Tipicamente, para muitas geometrias e valores de \mathrm{Bi}, a aproximação de primeiro termo dá erro pequeno (por exemplo < a alguns por cento) quando
\mathrm{Fo}\gtrsim 0{,}2\\ \text{a}\\\ 0{,}5.
Isto é um valor prático: ao redor de \mathrm{Fo}\sim 0{,}1 ou menos o primeiro termo geralmente não é suficiente — é a fase inicial.
  • Estes números são orientativos; o critério rigoroso é a fórmula do Passo 4 usando \lambda_1,\lambda_2,A_1,A_2 (tabelados ou obtidos numericamente).
Em que estágio a exigência NÃO é satisfeita?
  • No estágio inicial do processo transiente, quando \mathrm{Fo} é pequeno (por exemplo \mathrm{Fo}\ll 0{,}1), a mudança térmica só penetrou uma pequena camada superficial de profundidade \delta\sim\sqrt{\alpha t}. A solução correta nessa fase tende ao comportamento de meio semi-infinito (solução com função erro / error function) e os modos mais altos da série são importantes. Portanto a aproximação de primeiro termo NÃO é satisfeita nos tempos iniciais logo após a mudança de contorno.
  • Em palavras mais físicas: quando a camada térmica ainda é muito pequena em comparação com o tamanho do corpo ( \delta\ll L_c ou \mathrm{Fo}\ll 1), a resposta local da superfície não representa o modo fundamental global; modos de ordem superior e soluções de curta-escala dominam. Só depois, quando a difusão já “vazou” pelo corpo (Fo não tão pequeno), o primeiro modo domina.

Passo 6

Resumo prático final (passos que você pode anotar na prova):
  • Calcule \mathrm{Bi} e escolha se o modelo de corpo bem misturado (lumped) é aplicável ( \mathrm{Bi}\le 0{,}1).
  • Se não for lumped, identifique L_c (meia-espessura para parede, raio para cilindro/esfera) e calcule \mathrm{Fo}=\dfrac{\alpha t}{L_c^{2}}.
  • Para precisão controlada: use a desigualdade do Passo 4 com \lambda_1,\lambda_2,A_1,A_2 apropriados (consulte tabelas/Heisler). Para regra rápida: se \mathrm{Fo}\gtrsim 0{,}2 0{,}5 a aproximação de primeiro termo costuma ser aceitável; se \mathrm{Fo}\ll 0{,}1 ela falha e então use aproximação de semi-infinito (função erro) ou mantenha mais termos da série.

Resposta

A interpretação física do número de Fourier é: \mathrm{Fo}=\dfrac{\alpha t}{L_c^{2}}, que representa o tempo adimensional (o quociente entre o tempo transcorrido t e o tempo característico de difusão térmica L_c^{2}/\alpha); fisicamente indica quantos “tempos de difusão” já ocorreram e está relacionado à profundidade de penetração térmica \delta\sim\sqrt{\alpha t}.
A exigência para usar a aproximação de primeiro termo é que os termos de ordem superior da solução em série sejam desprezíveis em relação ao primeiro termo; matematicamente, para tolerância \varepsilon,
\mathrm{Fo}\gg \frac{\ln\!\big(\dfrac{|A_1\Phi_1|}{\varepsilon\,|A_2\Phi_2|}\big)}{\lambda_2^{2}-\lambda_1^{2}},
onde \lambda_1,\lambda_2 são os dois primeiros autovalores (dependendo de \mathrm{Bi} e da geometria) e A_n,\Phi_n são os coeficientes e autofunções correspondentes. Em termos práticos, para parede plana (meia-espessura L), cilindro de raio r_0 e esfera de raio r_0 usa-se L_c=L,\;L_c=r_0,\;L_c=r_0 respectivamente, e uma regra prática é que a aproximação de primeiro termo costuma ser aceitável para \mathrm{Fo}\gtrsim 0{,}2 0{,}5 (verifica-se rigorosamente com a fórmula anterior ou tabelas).
A exigência NÃO é satisfeita no estágio inicial do processo transiente (tempos pequenos, \mathrm{Fo} pequeno, tipicamente \mathrm{Fo}\ll 0{,}1), quando a perturbação térmica está confinada a uma camada superficial (comportamento de meio semi-infinito) e os modos de ordem superior têm contribuição significativa.

Exercícios relacionados