2. Fundamentação Teórica

a. Descreva as equações diferenciais que regem a condução de calor transiente unidimensional em sólidos semi-infinitos.
b. Apresente a solução analítica clássica do problema e a identidade matemática da modelagem (uso da função erro, condições de contorno típicas, hipóteses).
c. Explique os principais parâmetros físicos envolvidos: condutividade, difusividade térmica, temperatura inicial e de fronteira, etc.

Passo 1

O enunciado pede três coisas relacionadas ao problema de condução de calor transiente unidimensional em sólidos semi-infinitos:
  1. (a) Descrever as equações diferenciais que regem o problema (PDE, condições iniciais e de contorno).
  2. (b) Apresentar a solução analítica clássica e a identidade matemática usada (função erro / erfc), com as condições de contorno típicas e hipóteses do modelo.
  3. (c) Explicar os principais parâmetros físicos envolvidos ( k, \alpha, T_i, temperaturas de fronteira, etc.).
Informações relevantes implícitas no enunciado:
  • Geometria: semi-infinita no eixo x\ge 0 (ou seja, uma barra/material que se estende efetivamente até x\to\infty).
  • Problema transiente unidimensional: dependência T(x,t) com x e t.
  • Espera-se a solução clássica que envolve a função erro (erf) / função erro complementar (erfc).
  • Condições típicas: temperatura inicial uniforme e condições na superfície x=0 (temperatura prescrita, fluxo prescrito ou convecção).

Passo 2

Escrevendo formalmente a equação diferencial que rege o problema (a equação do calor em 1D, com propriedades constantes):
\frac{\partial T}{\partial t} = \alpha \frac{\partial^{2} T}{\partial x^{2}},\qquad x>0,\; t>0,
onde T(x,t) é a temperatura, \alpha é a difusividade térmica (ver definição depois). Essa é a equação básica: taxa de variação temporal de temperatura igual ao transporte por difusão térmica no espaço.
Condições iniciais e de contorno típicas (escrevo as três mais comuns que aparecem em provas):
  • Condição inicial (inicialização uniforme):
T(x,0)=T_i,\qquad x\ge 0,
onde T_i é a temperatura inicial uniforme do sólido.
  • Condição de contorno à superfície x=0:
    • Temperatura prescrita (Dirichlet):
    T(0,t)=T_s,\qquad t>0.
    • Flux de calor prescrito (Neumann), positivo entrando no sólido:
    -k\left.\frac{\partial T}{\partial x}\right|_{x=0}=q_0,\qquad t>0,onde k é a condutividade térmica e q_0 é o fluxo (W/m ^2).
    • Convecção (condição de Robin / Newton):
    -k\left.\frac{\partial T}{\partial x}\right|_{x=0}=h\bigl(T(0,t)-T_\infty\bigr),\qquad t>0,onde h é coeficiente de convecção e T_\infty temperatura do fluido externo.
  • Condição no infinito (compatível com semi-infinito):
\lim_{x\to\infty} T(x,t)=T_i,\qquad t\ge 0.

Passo 3

Hipóteses do modelo (é importante checar isso numa prova; se não valem, a solução clássica não é aplicável):
  • Propriedades térmicas homogêneas e independentes do tempo: k,\rho,c constantes.
  • Material isotrópico e sem fontes internas de calor.
  • Transferência unidimensional em x (variações em outras direções desprezíveis).
  • Semi-infinito: o campo térmico não alcança x\to\infty durante o tempo de interesse; numericamente, a profundidade de penetração \sqrt{\alpha t} é muito menor que a dimensão física do corpo que faria efeito.
  • Condição de fronteira aplicada repentinamente em t=0^+ (ex.: mudança instantânea de temperatura na superfície).
  • Regime linear (pequenas variações de temperatura que mantêm k,\alpha constantes).

Passo 4

Técnica padrão para resolver o caso clássico (superfície subitamente levada a temperatura fixa T_s): usar variável de similaridade. Vou fazer a derivação passo a passo porque é exatamente o que cai na prova.
Defina a temperatura adimensional
\theta(x,t)=\frac{T(x,t)-T_i}{T_s-T_i}.
Queremos resolver a PDE para \theta com \theta(x,0)=0, \theta(0,t)=1, \theta(\infty,t)=0.
Agora introduza a variável de similaridade
\eta=\frac{x}{2\sqrt{\alpha t}}.
Hipótese: a solução depende de x e t apenas através de \eta, ou seja \theta(x,t)=\Theta(\eta). Vou calcular as derivadas usando a regra da cadeia (mostro cada passo para você replicar na prova):
Cálculo de \partial \eta/\partial t (útil):
\eta=\frac{x}{2\sqrt{\alpha t}} \quad\Rightarrow\quad \frac{\partial\eta}{\partial t}=-\frac{x}{4\sqrt{\alpha}\,t^{3/2}}.
Mas x=2\eta\sqrt{\alpha t}, então
\frac{\partial\eta}{\partial t}=-\frac{\eta}{2t}.
Agora derivadas de \theta:
\frac{\partial \theta}{\partial t} = \Theta'(\eta)\frac{\partial\eta}{\partial t} = \Theta'(\eta)\left(-\frac{\eta}{2t}\right).
\frac{\partial \theta}{\partial x} = \Theta'(\eta)\frac{\partial\eta}{\partial x} = \Theta'(\eta)\frac{1}{2\sqrt{\alpha t}}.
\frac{\partial^2 \theta}{\partial x^2} = \Theta''(\eta)\frac{1}{4\alpha t}.
Substitua tudo na PDE \partial \theta/\partial t = \alpha\,\partial^2\theta/\partial x^2:
\Theta'(\eta)\left(-\frac{\eta}{2t}\right) = \alpha \,\Theta''(\eta)\frac{1}{4\alpha t}.
Multiplicando ambos os lados por 4t e simplificando:
-2\eta\,\Theta'(\eta)=\Theta''(\eta).
Reescrevendo:
\Theta''(\eta)+2\eta\,\Theta'(\eta)=0.

Passo 5

Integração da ODE reduzida (é um passo padrão e curto, faço devagar):
Defina u(\eta)=\Theta'(\eta). Então a ODE vira
u'(\eta)+2\eta u(\eta)=0.
Isto é uma EDO linear de primeira ordem para u. Separando variáveis:
\frac{u'}{u}=-2\eta \quad\Rightarrow\quad \ln |u| = -\eta^2 + C_1.
Então
u(\eta)=C_1 e^{-\eta^2}.
Agora integre para obter \Theta:
\Theta(\eta)=C_1\int e^{-\eta^2}d\eta + C_2.
A primitiva de e^{-\eta^2} é relacionada à função erro:
\int_0^{\eta} e^{-s^2} ds = \frac{\sqrt{\pi}}{2}\operatorname{erf}(\eta).
Portanto
\Theta(\eta)=C_1 \frac{\sqrt{\pi}}{2}\operatorname{erf}(\eta) + C_2.
Aplique as condições-limite:
  • Quando \eta\to\infty (isto é x\to\infty), \Theta\to 0. Como \operatorname{erf}(\infty)=1, obtemos C_1\frac{\sqrt{\pi}}{2}+C_2=0.
  • Em \eta=0 (isto é x=0), \Theta(0)=1. Como \operatorname{erf}(0)=0, isto dá C_2=1.
Resolvendo os dois: C_2=1 e então C_1 = -\dfrac{2}{\sqrt{\pi}}.
Substituindo de volta:
\Theta(\eta)=1-\operatorname{erf}(\eta)=\operatorname{erfc}(\eta).
Voltando para a variável física,
\boxed{\displaystyle \frac{T(x,t)-T_i}{T_s-T_i} = \operatorname{erfc}\!\left(\frac{x}{2\sqrt{\alpha t}}\right).}
Ou explicitamente,
\boxed{\displaystyle T(x,t)=T_i + (T_s-T_i)\,\operatorname{erfc}\!\left(\frac{x}{2\sqrt{\alpha t}}\right).}
Verificação rápida (coisa que o examinador gosta): para x=0 temos \operatorname{erfc}(0)=1 e T(0,t)=T_s. Para x\to\infty, \operatorname{erfc}\to 0 e T\to T_i. Tudo consistente.

Passo 6

Soluções para outras condições de contorno superficiais (resumo importante para prova):
  • Fluxo superficial constante q_0 (Neumann: -k\partial T/\partial x|_{0}=q_0, com T(x,0)=T_i):
\boxed{\displaystyle T(x,t)-T_i = \frac{2q_0\sqrt{\alpha t}}{k\sqrt{\pi}}\,e^{-\,\eta^{2}} \;-\; \frac{q_0 x}{k}\,\operatorname{erfc}(\eta),}
onde \eta=\dfrac{x}{2\sqrt{\alpha t}}. (Este é o resultado clássico; checando em x=0 obtemos crescimento em \sqrt{t} conforme esperado para um fluxo constante.)
  • Convecção na superfície (condição de Newton: -k\partial T/\partial x|_{0}=h(T(0,t)-T_\infty)): existe solução em termos de combinações de \operatorname{erfc} e exponenciais. Forma compacta (modelo de solução por similaridade):
T(x,t)-T_i = (T_\infty - T_i)\!\left[\operatorname{erfc}(\eta) - \exp\!\bigl(\eta^2 + \beta\eta\bigr)\operatorname{erfc}(\eta+\beta)\right],
onde \beta=\dfrac{h\sqrt{\alpha t}}{k}. (Essa forma aparece em livros como Carslaw & Jaeger; é mais técnica, mas útil quando h é finito.)
Observação prática: nas provas normalmente pedem ou a solução Dirichlet com \operatorname{erfc} ou a fórmula para fluxo constante; a forma para convecção é dada ou pedida com indicação do parâmetro \beta.

Passo 7

Definição e propriedades da função erro (parte do “identidade matemática”):
  • Definição da função erro:
\operatorname{erf}(z)=\frac{2}{\sqrt{\pi}}\int_{0}^{z} e^{-u^{2}}\,du.
  • Função erro complementar:
\operatorname{erfc}(z)=1-\operatorname{erf}(z)=\frac{2}{\sqrt{\pi}}\int_{z}^{\infty} e^{-u^{2}}\,du.
Propriedades úteis:
  • \operatorname{erf}(0)=0, \operatorname{erfc}(0)=1.
  • \operatorname{erf}(\infty)=1, \operatorname{erfc}(\infty)=0.
  • Para pequenos argumentos \operatorname{erfc}(z)\approx 1 - \dfrac{2}{\sqrt{\pi}}z (expansão inicial), e para grandes z decai rapidamente como e^{-z^{2}}.
Dica prática: na prova, reconhecer \eta=x/(2\sqrt{\alpha t}) e lembrar que integrais gaussianas conduzem a \operatorname{erf} e \operatorname{erfc} é a pista para a solução.

Passo 8

Explicação física e unidades dos principais parâmetros (faça sempre essa checagem mental numa prova):
  • Condutividade térmica k (unidade: \mathrm{W\,m^{-1}\,K^{-1}}): mede a capacidade do material de conduzir calor. Fluxo por gradiente é -k\partial T/\partial x.
  • Densidade \rho (unidade: \mathrm{kg\,m^{-3}}) e calor específico c (unidade: \mathrm{J\,kg^{-1}\,K^{-1}}).
  • Difusividade térmica \alpha (unidade: \mathrm{m^{2}\,s^{-1}}), definida por
\alpha=\frac{k}{\rho c}.
Interpretação: \alpha controla a velocidade com que a perturbação térmica se propaga por difusão. Quanto maior \alpha, mais rápido o calor se espalha.
  • Temperatura inicial T_i e temperatura de superfície imposta T_s (ou T_\infty para convecção): definem os contrastes térmicos que dirigem o fluxo.
  • Fluxo superficial q_0 (unidade: \mathrm{W\,m^{-2}}) quando a condição na superfície é de fluxo.
  • Coeficiente de convecção h (unidade: \mathrm{W\,m^{-2}\,K^{-1}}) para transferência entre superfície e fluido.
Noções físicas importantes:
  • Profundidade de penetração térmica característica:
\delta(t)\sim\sqrt{\alpha t}.
Isto dá a escala em x onde mudanças de temperatura são relevantes. Justifica a hipótese de semi-infinito se \delta(t) for muito menor que a dimensão física do corpo.
  • Para fluxo constante, a temperatura na superfície aumenta como \sqrt{t} (veja termo \propto \sqrt{t} na solução de Neumann).

Passo 9

Como usar isso na prova — checklist prático e dicas rápidas (pronto para colar na sua folha de raciocínio):
  1. Identifique geometria e se é razoável usar semi-infinito (compare L com \sqrt{\alpha t}).
  2. Escreva a PDE e as condições iniciais/contorno claramente.
  3. Verifique se as propriedades são constantes (se não forem, a solução clássica muda).
  4. Se a fronteira tem temperatura prescrita: use a solução com \operatorname{erfc} que dei no Passo 5.
  5. Se a fronteira tem fluxo constante: use a fórmula do Passo 6.
  6. Se for convecção, ou apresente a fórmula pertinente ou transforme em condição com parâmetro adimensional e aplique tabela de soluções (Carslaw & Jaeger).
  7. Faça checagens de limites: x=0, x\to\infty, t\to 0^+, e unidades.
  8. Se precisar de valor numérico de \operatorname{erfc}, em provas costuma-se permitir tabela/valores ou pedir comportamento assintótico; algumas provas esperam só escrita simbólica.
Motivação/ânimo: se bater desespero, lembra: reconhecer que é semi-infinito e tirar a variável \eta=x/(2\sqrt{\alpha t}) já resolve grande parte do problema. Respira e aplica a transformação com calma — é mecânico.

Passo 10

Exemplos rápidos (curtos) para fixar:
  • Exemplo 1 (temperatura fixa na superfície): se \alpha=10^{-6}\,\mathrm{m^2/s}, t=100\,\mathrm{s}, x=1\,\mathrm{mm}=10^{-3}\,\mathrm{m}, então \eta=\dfrac{10^{-3}}{2\sqrt{10^{-6}\cdot 100}}=\dfrac{10^{-3}}{2\cdot 10^{-2}}=0{,}025. Assim \operatorname{erfc}(0{,}025)\approx 1- \dfrac{2}{\sqrt{\pi}}(0{,}025)\approx 0{,}971 (aprox.). Assim T ainda está muito próximo de T_s.
  • Exemplo 2 (fluxo constante): observe que a presença do termo \dfrac{2q_0\sqrt{\alpha t}}{k\sqrt{\pi}} indica crescimento com \sqrt{t} na superfície.

Resposta

(a) \displaystyle \frac{\partial T}{\partial t}=\alpha\frac{\partial^{2}T}{\partial x^{2}},\quad T(x,0)=T_i,\quad \text{e na superfície }x=0\\ \text{condições típicas: }T(0,t)=T_s\\ \text{(Dirichlet)}\\ \text{ou}\\\ -k\frac{\partial T}{\partial x}\bigg|_{0}=q_0\\ \text{(Neumann)},\\ \lim_{x\to\infty}T(x,t)=T_i.íçõí
(b) Para superfície subitamente fixada em T_s:
\displaystyle T(x,t)=T_i+(T_s-T_i)\,\operatorname{erfc}\!\Bigl(\frac{x}{2\sqrt{\alpha t}}\Bigr),
onde \operatorname{erf}(z)=\dfrac{2}{\sqrt{\pi}}\int_{0}^{z}e^{-u^{2}}du e \operatorname{erfc}(z)=1-\operatorname{erf}(z). Para fluxo superficial constante q_0:
\displaystyle T(x,t)-T_i=\frac{2q_0\sqrt{\alpha t}}{k\sqrt{\pi}}e^{-\eta^{2}}-\frac{q_0 x}{k}\operatorname{erfc}(\eta),\quad \eta=\frac{x}{2\sqrt{\alpha t}}.
Hipóteses: propriedades constantes, unidimensional, sem gerações internas, semi-infinito ( \delta\sim\sqrt{\alpha t}\ll dimensão física), mudança de fronteira aplicada em t=0^+.
(c) Parâmetros físicos: k = condutividade ( \mathrm{W\,m^{-1}\,K^{-1}}), \rho = densidade ( \mathrm{kg\,m^{-3}}), c = calor específico ( \mathrm{J\,kg^{-1}\,K^{-1}}), difusividade \alpha=\dfrac{k}{\rho c} ( \mathrm{m^{2}\,s^{-1}}), T_i = temperatura inicial, T_s / T_\infty = temperaturas de fronteira, q_0 = fluxo superficial, h = coeficiente de convecção. A escala de penetração térmica é \delta(t)\sim\sqrt{\alpha t}, que justifica a aproximação semi-infinita.

Exercícios relacionados