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:
(a) Descrever as equações diferenciais que regem o problema (PDE, condições iniciais e de contorno).
(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.
(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):
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):
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):
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):
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”):
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):
Identifique geometria e se é razoável usar semi-infinito (compare L com \sqrt{\alpha t}).
Escreva a PDE e as condições iniciais/contorno claramente.
Verifique se as propriedades são constantes (se não forem, a solução clássica muda).
Se a fronteira tem temperatura prescrita: use a solução com \operatorname{erfc} que dei no Passo 5.
Se a fronteira tem fluxo constante: use a fórmula do Passo 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).
Faça checagens de limites: x=0, x\to\infty, t\to 0^+, e unidades.
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:
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:
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.