Questão 11
Derive uma expressão para a distribuição de temperatura em uma parede plana na qual as fontes de calor distribuídas variam de acordo com a relação linear.
\dot{q}=\dot{q}_{w}\left[1+\beta\left(T-T_{w}\right)\right]
onde \dot{q}_{w} é uma constante e igual ao calor gerado por unidade de volume na parede à temperatura da parede de T_{w}. Ambos os lados da parede são mantidos a T_{w}, e a espessura da placa é 2 L. Considere \beta um parâmetro constante.

Passo 1

O enunciado pede para achar a distribuição de temperatura numa placa plana de espessura 2L (ou seja, x\in[-L,L]), onde o calor gerado por unidade de volume varia linearmente com a temperatura segundo
\dot{q}=\dot{q}_{w}\left[1+\beta\left(T-T_{w}\right)\right],
sendo \dot{q}_{w} a geração volumétrica quando a temperatura local é T_{w}. Ambas as faces da placa (em x=\pm L) são mantidas à temperatura \,T_{w}.
Informações relevantes dadas:
\dot{q}=\dot{q}(T)=\dot{q}_{w}\bigl[1+\beta(T-T_{w})\bigr],
espessura total 2L, temperaturas de contorno T(\pm L)=T_{w}.
Vou assumir conduta usual: regime estacionário unidimensional e condutividade térmica constante k (se o problema não dá k, é padrão introduzi-la para relacionar fluxo e gradiente). A tarefa é encontrar T(x) para -L\le x\le L.

Passo 2

Escrevendo a equação de balanço de energia para condução unidimensional no regime estacionário com geração volumétrica:
\frac{d^{2}T}{dx^{2}}+\frac{\dot{q}(T)}{k}=0,
onde k é a condutividade térmica (unidade: W/m·K). Substituindo a expressão dada para \dot{q}(T):
\frac{d^{2}T}{dx^{2}}+\frac{\dot{q}_{w}}{k}\bigl[1+\beta(T-T_{w})\bigr]=0.
Agora vou manipular essa equação para colocá‑la em forma padrão de uma equação diferencial linear de segunda ordem com coeficientes constantes. Expandindo:
\frac{d^{2}T}{dx^{2}}+\frac{\dot{q}_{w}\beta}{k}\,T+\frac{\dot{q}_{w}}{k}\bigl(1-\beta T_{w}\bigr)=0.
Defino para simplificar:
A\equiv\frac{\dot{q}_{w}\beta}{k},\qquad B\equiv\frac{\dot{q}_{w}}{k}\bigl(1-\beta T_{w}\bigr).
A equação fica
\frac{d^{2}T}{dx^{2}}+A\,T+B=0.

Passo 3

Resolvendo a EDO. Esta é uma EDO linear não homogênea com coeficientes constantes. Método padrão: achar uma solução particular T_{p} (constantemente) e a solução geral da homogênea.
  1. Solução particular constante: suponho T_{p}= constante. Substituindo em \,A T_{p}+B=0, obtemos
A\,T_{p}+B=0\quad\Rightarrow\quad T_{p}=-\frac{B}{A}.
Voltando às definições de A e B:
T_{p}=-\frac{\dfrac{\dot{q}_{w}}{k}\bigl(1-\beta T_{w}\bigr)}{\dfrac{\dot{q}_{w}\beta}{k}} =-\frac{1-\beta T_{w}}{\beta}=T_{w}-\frac{1}{\beta}.
(Observação: isso é válido se \beta\neq 0. O caso \beta=0 será tratado depois.)
  1. Equação homogênea associada:
\frac{d^{2}T_{h}}{dx^{2}}+A\,T_{h}=0.
Se A>0 (isto é, \dot{q}_{w}\beta/k>0), defino
m\equiv\sqrt{A}=\sqrt{\frac{\dot{q}_{w}\beta}{k}},
e a solução homogênea é
T_{h}(x)=C_{1}\cos(mx)+C_{2}\sin(mx).
Se A<0 ( \beta<0 se \dot{q}_w/k>0), escrevo A=-\mu^{2} com \mu=\sqrt{-A}; então
T_{h}(x)=C_{1}\operatorname{cosh}(\mu x)+C_{2}\operatorname{sinh}(\mu x).
Esses dois casos são análogos (seno/cosseno versus seno hiperbólico/cosseno hiperbólico).
Portanto a solução geral é a soma:
T(x)=T_{p}+T_{h}(x).

Passo 4

Impondo a simetria e as condições de contorno. A geração e as condições são simétricas em x=0, então a solução física será par (simétrica); isto implica \displaystyle T(-x)=T(x) e, portanto, a componente ímpar da homogênea deve sumir. Consequentemente:
  • Para o caso A>0: C_{2}=0 e T(x)=T_{p}+C_{1}\cos(mx).
  • Para o caso A<0: C_{2}=0 e T(x)=T_{p}+C_{1}\operatorname{cosh}(\mu x).
Agora aplico as condições de contorno T(\pm L)=T_{w}. Como solução é par, basta aplicar em x=L.
Caso A>0 ( \beta>0 supondo \dot{q}_{w}/k>0):
T(L)=T_{p}+C_{1}\cos(mL)=T_{w}.
Isolando C_{1}:
C_{1}=T_{w}-T_{p}.
Mas T_{w}-T_{p}=T_{w}-\bigl(T_{w}-\tfrac{1}{\beta}\bigr)=\tfrac{1}{\beta}. Então
C_{1}=\frac{1}{\beta}.
Atenção: essa expressão para C_{1} vale somente se \cos(mL)=1. Porém no cálculo anterior eu precisei levar em conta a \emph{equação} correta:
corrigindo o passo: a equação correta é
T_{p}+C_{1}\cos(mL)=T_{w}\quad\Rightarrow\quad C_{1}=\frac{T_{w}-T_{p}}{\cos(mL)}=\frac{1/\beta}{\cos(mL)}.
Portanto a solução final para A>0 é
T(x)=T_{p}+\frac{1}{\beta}\frac{\cos(mx)}{\cos(mL)}.
Substituindo T_{p}=T_{w}-1/\beta,
T(x)=T_{w}-\frac{1}{\beta}+\frac{1}{\beta}\frac{\cos(mx)}{\cos(mL)} =T_{w}+\frac{1}{\beta}\left(\frac{\cos(mx)}{\cos(mL)}-1\right).
Caso A<0 ( \beta<0 supondo \dot{q}_{w}/k>0): defina \mu=\sqrt{-A}=\sqrt{\dfrac{-\dot{q}_{w}\beta}{k}}. Pela mesma lógica:
T(x)=T_{p}+C_{1}\operatorname{cosh}(\mu x),
com
C_{1}=\frac{T_{w}-T_{p}}{\operatorname{cosh}(\mu L)}=\frac{1/\beta}{\operatorname{cosh}(\mu L)}.
Logo
T(x)=T_{w}+\frac{1}{\beta}\left(\frac{\operatorname{cosh}(\mu x)}{\operatorname{cosh}(\mu L)}-1\right).

Passo 5

Verificações e o caso \beta=0:
  1. Caso \beta\to 0 (ou seja, geração independente de T): a expressão em série de \cos(mx)/\cos(mL) com m^{2}=\dfrac{\dot{q}_{w}\beta}{k} dá, no limite \beta\to 0, a clássica solução polinomial
T(x)=T_{w}+\frac{\dot{q}_{w}}{2k}\bigl(L^{2}-x^{2}\bigr),
que é o resultado esperado para geração volumétrica constante \dot{q}=\dot{q}_{w}. Aqui pode-se verificar tomando limite usando expansão de Taylor (mostrei a consistência).
  1. Observação sobre sinais e singularidades: a solução com termo 1/\beta exige cuidado se \beta for muito pequena numericamente; o limite \beta\to 0 é regular e dá a fórmula polinomial acima. Se \beta tem sinal negativo, usamos a forma com \operatorname{cosh}. As fórmulas acima são as soluções gerais e satisfazem T(\pm L)=T_{w} e a EDO original.

Resposta

\text{Para }\beta>0:\qquad T(x)=T_{w}+\frac{1}{\beta}\!\left(\frac{\cos\!\bigl(mx\bigr)}{\cos\!\bigl(mL\bigr)}-1\right), \quad m\equiv\sqrt{\frac{\dot{q}_{w}\beta}{k}}.
\text{Para }\beta<0:\qquad T(x)=T_{w}+\frac{1}{\beta}\!\left(\frac{\operatorname{cosh}\!\bigl(\mu x\bigr)}{\operatorname{cosh}\!\bigl(\mu L\bigr)}-1\right), \quad \mu\equiv\sqrt{\frac{-\dot{q}_{w}\beta}{k}}.
\text{Caso }\beta=0:\qquad T(x)=T_{w}+\frac{\dot{q}_{w}}{2k}\bigl(L^{2}-x^{2}\bigr).

Exercícios relacionados