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.
Temos uma parede permeável a passagem de um dado gás, e queremos determinar a posição do gás ao longo da espessura da parede em função do tempo.

Passo 1

O enunciado: temos uma parede (uma lâmina) permeável à passagem de um gás e queremos determinar a posição/concentração do gás ao longo da espessura da parede em função do tempo.
Informações explícitas dadas pelo enunciado: a parede é permeável e a variável de interesse é a posição/concentração do gás ao longo da espessura e sua evolução temporal.
O que o enunciado pede: escolher a abordagem mais apropriada — parâmetros distribuídos (difusão) ou parâmetros agrupados (coeficiente de transferência de massa) — e, implicitamente, indicar como modelar o problema para obter a concentração ao longo da espessura em função do tempo.

Passo 2

O que significam as duas abordagens (em linguagem simples):
  • Parâmetros distribuídos (modelo de difusão): a concentração varia com a posição dentro da parede; tratamos a variável concentração como função C(x,t) e usamos uma equação diferencial parcial (PDE) que descreve como a concentração muda no espaço e no tempo. É preciso resolver para toda a distribuição espacial, isto é, obter a curva C para cada x e t.
  • Parâmetros agrupados (modelo lumped, coeficiente de transferência de massa k): assumimos que a concentração dentro da parede é praticamente uniforme (não há gradiente espacial significativo). Então descrevemos apenas a evolução média da concentração no tempo com uma equação diferencial ordinária (ODE) usando um coeficiente k que representa a resistência à troca entre a superfície e o ambiente. Resultado é C(t) médio — não temos perfil espacial.

Passo 3

Quando cada abordagem é adequada — explicação física e critério prático:
  • Use modelo distribuído quando existirem gradientes importantes de concentração dentro da parede, ou quando o problema pede explicitamente a informação espacial C(x,t) (como no enunciado: “posição do gás ao longo da espessura”). Aqui é natural esperar variação espacial, logo o modelo distribuído é o apropriado.
  • Use modelo agrupado (lumped) somente quando as diferenças de concentração dentro da parede forem desprezíveis em comparação com as diferenças entre a superfície e o fluido externo. Em termos práticos existe um critério adimensional, o número de Biot para massa, que compara a resistência externa (convecção) com a resistência interna (difusão).

Passo 4

Definição do número de Biot para transferência de massa e seu significado:
Definimos o número de Biot de massa como
\mathrm{Bi}=\frac{k\,L_c}{D},
onde:
  • k é o coeficiente de transferência de massa (unidade: \mathrm{m/s}); representa a resistência convectiva externa entre a superfície da parede e o ambiente;
  • D é a difusividade do gás na parede (unidade: \mathrm{m^2/s}); representa a capacidade do gás se difundir internamente;
  • L_c é um comprimento característico do corpo. Para uma lâmina de espessura e é comum usar L_c=V/A. Para uma placa plana uniforme com área A e volume V=A\,e, temos L_c=V/A=e se considerarmos toda a espessura ou L_c=e/2 se a formulação for simétrica (dependendo da convenção). A escolha habitual para critério é L_c = e/2 quando o problema é simétrico de ambos os lados; usar L_c=e também é aceitável se isso for explicitado.
Interpretação:
  • Se \mathrm{Bi}\ll 1 (regra prática: \mathrm{Bi}\lesssim 0{,}1), a resistência interna por difusão é pequena em comparação com a resistência externa; o interior se mistura rápido e a concentração dentro pode ser tomada como praticamente uniforme — então o modelo lumped é válido.
  • Se \mathrm{Bi}\gtrsim 0{,}1 (especialmente \gg 0{,}1), a resistência interna não é desprezível e surgirão gradientes de concentração dentro da parede — então é obrigatório usar o modelo distribuído (difusão).

Passo 5

Equações básicas do modelo distribuído (a “aula” sobre difusão que você pediu):
A difusão unidimensional e transiente é descrita pela equação de difusão (lei de Fick combinada com balanço de massa):
\frac{\partial C}{\partial t}=D\frac{\partial^2 C}{\partial x^2},
onde:
  • C(x,t) é a concentração do gás na parede em função da posição x (perpendicular às faces) e do tempo t;
  • D é a difusividade dentro do material.
A lei de Fick para o fluxo é
J(x,t)=-D\frac{\partial C}{\partial x},
onde J é o fluxo de massa (unidade: massa por área por tempo) direcionado segundo o gradiente de concentração (sinal negativo porque fluxo vai do maior para o menor).
Condições de contorno típicas (exemplos e significado físico):
  • Superfície com concentração imposta (Dirichlet): \;C(0,t)=C_s(t) e/ou \;C(e,t)=C_{s2}(t). Isso modela a superfície mantendo uma concentração conhecida (caso limite).
  • Superfície com resistência convectiva externa (Neumann/Robin mista): \;-D\left.\frac{\partial C}{\partial x}\right|_{x=0}=k\left(C(0,t)-C_\infty\right). Isso significa que o fluxo que sai da superfície por difusão é igual ao fluxo convectivo para o ambiente, com C_\infty a concentração no fluido longe da parede. É aqui que o coeficiente k aparece acoplado com a equação interna.
Condição inicial:
C(x,0)=C_0(x),
que tipicamente é C_0 constante (parede inicialmente com concentração uniforme).

Passo 6

Como resolver, em linhas gerais (dicas práticas para a prova):
  • Primeiro passo: calcule \mathrm{Bi}=\dfrac{k\,L_c}{D} com valores dados (ou estimados). Escolha L_c consistente (por exemplo, L_c=e/2 para uma placa).
  • Se \mathrm{Bi}\le 0{,}1 e o problema pede apenas a concentração média na parede, use o modelo lumped: derive a ODE de balanço e resolva (exponencial). Para uma placa com espessura e e concentração uniforme C(t), balanceando taxa de variação de massa com fluxo convectivo pela superfície obtém-se
\frac{\mathrm{d}C}{\mathrm{d}t}=\frac{k}{L_c}\bigl(C_\infty-C(t)\bigr),
com solução
C(t)=C_\infty+\bigl(C(0)-C_\infty\bigr)\exp\!\left(-\frac{k}{L_c}t\right).
  • Se \mathrm{Bi}\gtrsim 0{,}1 ou se você precisa do perfil C(x,t) (como no enunciado), monte a PDE de difusão com as condições de contorno apropriadas e resolva:
    • Casos simples permitem solução analítica: por exemplo, meio semi-infinito com superfície mantida em C_s tem solução com função erro; placa finita com contorno convectivo tem solução via séries (Crank) ou via transformadas.
    • Nos casos práticos com condições complicadas ou D dependendo de C, use método numérico: diferenças finitas explícitas/implícitas, método dos volumes finitos, ou elementos finitos. Para prova, é comum usar separação de variáveis (série de senos/cossenos) ou mencionar a solução tabelada de Crank para placa com convecção nas faces.

Passo 7

Relação com escalas temporais (ajuda a entender o que acontece):
  • Tempo característico de difusão interna:
t_d\sim\frac{L_c^{\,2}}{D}.
  • Tempo característico imposto pela convecção externa (lumped):
t_k\sim\frac{L_c}{k}.
Comparando os tempos,
\frac{t_k}{t_d}\sim\frac{D}{k\,L_c}=\frac{1}{\mathrm{Bi}}.
Se t_k\ll t_d ( \mathrm{Bi}\gg 1) a convecção é rápida comparada à difusão interna — surgem gradientes internos e é necessário o modelo distribuído. Se t_k\gg t_d ( \mathrm{Bi}\ll 1) o interior se equilibra rapidamente e o comportamento pode ser aproximado por uma ODE lumped.

Passo 8

Aplicando tudo ao enunciado específico (resposta prática e intuitiva):
  • O enunciado pede explicitamente a posição/concentração do gás ao longo da espessura em função do tempo — isto é exatamente informação espacial C(x,t). Mesmo sem números, o objetivo exige o perfil espacial: portanto a abordagem correta é a de parâmetros distribuídos (modelo de difusão transiente).
  • Na prática, você deve montar a PDE
\frac{\partial C}{\partial t}=D\frac{\partial^2 C}{\partial x^2},
com a condição inicial C(x,0)=C_0(x) e as condições de contorno adequadas (por exemplo, superfícies com convecção: -D\partial_x C|_{x=0}=k(C(0,t)-C_\infty), análogo na outra face). Depois solve analiticamente (se possível) ou numericamente.
  • Observação útil para prova: mesmo escolhendo o modelo distribuído, sempre calcule o número de Biot. Se for encontrado \mathrm{Bi}\ll 0{,}1 (o que poderia ocorrer para paredes muito finas e/ou difusividade muito alta), você pode justificar rapidamente uma simplificação lumped para obter só a média temporal — mas isso só é aceitável se a pergunta pedisse a concentração média; se a pergunta pede o perfil, mesmo com \mathrm{Bi}\ll 0{,}1 você pode apresentar que o perfil é praticamente plano e derivar a solução lumped como limite.

Resposta

Parâmetros distribuídos (difusão) — montar e resolver a equação de difusão transiente \dfrac{\partial C}{\partial t}=D\dfrac{\partial^2 C}{\partial x^2} com as condições de contorno apropriadas; usar o critério \mathrm{Bi}=\dfrac{kL_c}{D} para verificar se uma simplificação lumped (válida quando \mathrm{Bi}\lesssim 0{,}1) seria aceitável apenas quando for suficiente obter a concentração média e não o perfil espacial.

Exercícios relacionados