Para cada situação, descreva qual abordagem mais apropriada a utilizar, por parâmentros distribuídos (difusão) ou por parâmetros agrupados (coeficiente de transferência de massa) nos processos de transferência de massa descritos.
Na secagem de uma placa de madeira, queremos determinar a concentração da água presente em um dado ponto da placa.
Passo 1
O enunciado: temos a secagem de uma placa de madeira e queremos determinar a concentração da água presente em um dado ponto da placa.
Informações relevantes que o enunciado dá (ou que normalmente são implicadas numa questão desse tipo):
A peça é uma placa (geometria do tipo “lâmina” com uma espessura bem definida).
Há água no interior da madeira (um campo de concentração espacial C(\mathbf{x},t)).
A placa está sendo exposta ao ambiente (ar) que remove umidade pela superfície — logo há transporte na superfície entre madeira e ar.
O objetivo explícito é saber a concentração da água em um ponto específico dentro da placa (não apenas a umidade média).
Pergunta que temos que responder: qual abordagem é mais apropriada — modelo por parâmetros distribuídos (difusão no interior, equações parciais) ou modelo por parâmetros agrupados (modelo lumped com coeficiente de transferência de massa) — para achar a concentração em um ponto da placa?
Passo 2
Explicação conceitual das duas abordagens (isso é essencial para escolher direito):
Modelo por parâmetros agrupados (lumped parameter, coeficiente de transferência de massa): trata-se a peça como se a concentração dentro dela fosse praticamente uniforme ao longo do volume. Então representamos a troca com o ambiente por um único coeficiente convectivo k_c na superfície e escrevemos um balanço de massa do tipo ODE para a concentração média \bar C(t). Uma forma típica da equação é
\frac{d\bar C}{dt} = -\frac{k_c A}{V}\bigl(\bar C - C_\infty\bigr),
onde A é a área externa, V o volume, e C_\infty a concentração de vapor/umidade no ar. Esse modelo dá apenas a evolução da concentração média \bar C(t), não a distribuição C(x,t).
Modelo por parâmetros distribuídos (difusão): representa explicitamente a variação espacial da concentração dentro da placa, usando uma equação de difusão (uma PDE). Para um meio homogêneo e isotrópico, a equação é
onde D_{\text{eff}} é o coeficiente efetivo de difusão da água na madeira, C_s é a concentração na superfície, e \partial/\partial n derivada normal. Este modelo fornece C(\mathbf{x},t), ou seja, a concentração em qualquer ponto.
Passo 3
Critério prático e objetivo para escolher o modelo: número de Biot de massa \mathrm{Bi}_m.
Definição:
\mathrm{Bi}_m=\frac{k_c\,L}{D_{\text{eff}}},
onde:
k_c é o coeficiente convectivo de massa (unidade: \mathrm{m/s}), que caracteriza a resistência externa entre a superfície e o ar;
D_{\text{eff}} é o coeficiente efetivo de difusão no sólido (unidade: \mathrm{m^2/s});
L é a dimensão característica do sólido para difusão. Para uma placa de espessura 2e costuma-se usar L=e (meia-espessura) ou L= espessura/2 conforme convenção.
Interpretação física:
Se \mathrm{Bi}_m \ll 0{,}1: a resistência interna (difusão) é pequena comparada à resistência externa; o conteúdo dentro da placa fica praticamente homogêneo durante a secagem → modelo lumped (parâmetros agrupados) é aceitável.
Se \mathrm{Bi}_m \gtrsim 0{,}1: há gradientes internos significativos; a difusão interna não é rápida o suficiente para manter uniformidade → é necessário modelo distribuído (difusão).
Passo 4
Aplicação do critério ao caso concreto (placa de madeira, concentração em um ponto):
Importante: a pergunta pede a concentração num ponto da placa. Mesmo que \mathrm{Bi}_m seja pequeno, o modelo lumped só fornece \bar C(t) (média). Se realmente precisamos do valor localizado C(x_0,t), o modelo lumped não entrega isso com fidelidade — salvo fazer a forte aproximação C(x_0,t)\approx\bar C(t) quando se justificar. Portanto a escolha depende de dois itens:
o valor de \mathrm{Bi}_m (para saber se a suposição de uniformidade interna é válida);
o grau de precisão exigido para o valor no ponto (se a prova/pedido aceita aproximar C(x_0,t) por \bar C(t) ou exige a distribuição exata).
Em madeira, geralmente:
D_{\text{eff}} costuma ser relativamente baixo (madeiras são porosas com transporte lento por difusão líquida/evaporação);
a espessura de placas não é desprezível;
além disso, D_{\text{eff}} pode depender muito do teor de umidade e da direção (anisotropia).
Esses fatores tendem a aumentar a importância da difusão interna, levando a \mathrm{Bi}_m tipicamente não muito pequeno — logo gradientes internos aparecem.
Passo 5
Procedimento prático que eu seguiria numa prova para decidir e resolver (passo a passo aplicável em exame):
Identificar/estimar parâmetros: pegue ou estime k_c, D_{\text{eff}} e a meia-espessura L.
Se valores não são dados, justifique estimativas ou explique que sem valores não se pode computar numericamente \mathrm{Bi}_m.
Calcule
\mathrm{Bi}_m=\frac{k_c L}{D_{\text{eff}}}.
Compare com o critério:
Se \mathrm{Bi}_m < 0{,}1 e aceitar aproximar C(x_0,t)\approx\bar C(t), use o modelo lumped e escreva a ODE para \bar C. Resolva a ODE e, se necessário, diga que a concentração no ponto será aproximada pela média.
Se \mathrm{Bi}_m \ge 0{,}1 ou se a pergunta exige o valor pontual com precisão, adote o modelo distribuído: escreva a equação de difusão com as condições de contorno convectivas e condições iniciais, e escolha solução analítica (caso de coeficientes constantes e geometria de lâmina há soluções em séries de Fourier — ver solução para “slab with convective boundary”) ou solução numérica (diferença finita, elementos finitos) se D_{\text{eff}} variar com C ou anisotropia estiver presente.
Passo 6
Como montar o modelo distribuído para a placa (exposição didática passo a passo para resolver na prova):
Geometria simples: tomo eixo x ao longo da espessura, placa simétrica com espessura 2e; uso x=0 no centro e x=\pm e nas superfícies.
normalmente C_0 constante se a placa começou homogênea.
Solução: para C_0 constante e D_{\text{eff}} constante existe solução por séries de senos/cossenos (série de Fourier) que fornece C(x,t). Na prova, você pode:
escrever a forma da série (autovalores \lambda_n obtidos da transcendental que envolve \mathrm{Bi}_m);
se pedirem só procedimento, descreva que se usa separação de variáveis, aplica BCs e obtém uma série convergente;
para cálculos práticos, frequentemente o primeiro termo da série já dá boa aproximação em regimes de tempo médio/alto.
Passo 7
Como montar o modelo lumped (caso a escolha caiba) — e por que ele NÃO é ideal aqui:
Balanço de massa médio:
V\frac{d\bar C}{dt} = -k_c A\bigl(\bar C - C_\infty\bigr)
ou reescrevendo
\frac{d\bar C}{dt} = -\frac{k_c A}{V}\bigl(\bar C - C_\infty\bigr).
Limitação importante: esta equação dá apenas \bar C(t). Se a prova pede a concentração num ponto específico e houver gradientes internos, aproximar C(x_0,t) por \bar C(t) pode ser grosseiro e errado.
Passo 8
Regra prática/intuicional para madeira (resumo motivacional e encorajador):
Se a placa for fina e o ar for muito pouco convectivo (pequeno k_c) — ou seja, se \mathrm{Bi}_m \ll 0{,}1 — você pode usar lumped, mas lembre-se: isso dá média.
Na maioria das situações reais de secagem de madeira (placas com espessura razoável, D_{\text{eff}} baixo, anisotropia), haverá gradientes de concentração ao longo da espessura. Como você quer a concentração em um ponto, o caminho seguro e recomendado é o modelo distribuído por difusão.
Então, não se desespere: na prova, primeiro calcule ou estime \mathrm{Bi}_m. Se ele indicar que há gradiente, escreva a PDE de difusão, as condições e resolva analiticamente (se possível) ou justifique a necessidade de solução numérica. Isso mostra entendimento e é a resposta correta.