Um filme líquido de espessura constante igual a 1 cm desce por uma superfície plana (escoamento em regime laminar) em um experimento conhecido como “coluna de parede molhada”. A partir de uma posição situada em um plano
tal fluido é colocado em contato com um gás A, solúvel no líquido, de modo que na interface a condição de equilíbrio se estabelece de forma instantânea. Considerando que o tempo de exposição do fluido ao gás é de 55 s, calcular:
(a) Qual a concentração de A no líquido (em relação à de equilíbrio) em uma posição situada a 1 mm da superfície gás-líquido?
(b) Qual o comprimento do filme exposto ao gás A?
Dados:
Pa.s;
.
(a) Qual a concentração de A no líquido (em relação à de equilíbrio) em uma posição situada a 1 mm da superfície gás-líquido?
(b) Qual o comprimento do filme exposto ao gás A?
Passo 1
O enunciado diz que um filme líquido de espessura constante igual a
desce por uma superfície plana (escoamento laminar). A partir de uma posição em
o líquido é posto em contato com um gás
solúvel, e na interface a condição de equilíbrio se estabelece instantaneamente (ou seja, a concentração do dissolvido na superfície do líquido é igual à concentração de equilíbrio
). O tempo de exposição do filme ao gás é
.
As perguntas são:
(a) Qual a concentração de no líquido, em relação à de equilíbrio, numa posição situada a
da interface gás–líquido?
(b) Qual o comprimento do filme exposto ao gás ?
(a) Qual a concentração de
(b) Qual o comprimento do filme exposto ao gás
Dados fornecidos:
- Densidade do líquido:
. - Viscosidade dinâmica:
. - Difusividade molecular:
. - Espessura do filme:
. - Posição para (a):
. - Tempo de contato:
.
Observações iniciais (apenas para clareza do que vou assumir):
- Vou assumir que o filme está numa parede vertical (ângulo
, logo ), que é a configuração usual na “coluna de parede molhada”. - Vou assumir que, antes do contato com o gás, a concentração de
no líquido é (praticamente) zero e que a superfície do líquido é mantida instantaneamente na concentração de equilíbrio . Com isso o transporte no sentido perpendicular à interface é bem descrito por difusão transiente num meio semi-infinito (válido porque a espessura do filme é muito maior que a profundidade de penetração da difusão no tempo dado — isso será verificado mais adiante).
Passo 2
Modelagem do problema de difusão perpendicular (para a pergunta (a)) e solução conhecida.
Quando se impõe, de repente, uma concentração constante
na superfície
de um sólido (ou líquido) inicialmente sem soluto, e o transporte na direção
é por difusão molecular com difusividade
, a equação que descreve o problema unidimensional é a equação da difusão (equação de calor transiente):
Condições típicas usadas aqui:
-
para (inicialmente sem soluto), -
para (superfície mantida em equilíbrio), -
quando .
A solução clássica (método da função de similaridade) para um meio semi-infinito é
onde
é a função complemento do erro (comportamento:
e tende a
quando o argumento cresce). Portanto a concentração em relação à de equilíbrio é
Explicando de forma leve o que é
: a função erro
é uma integral ligada a distribuições gaussianas;
. Para valores do argumento maiores que 1 a
já é pequena (por exemplo,
é da ordem de
). Aqui ela mede quanto da concentração de equilíbrio já penetrou até a distância
em tempo
.
Antes de aplicar esta solução, preciso justificar a hipótese semi-infinita: a profundidade característica da penetração difusiva é do tipo
. Se
, o fundo do filme não afeta o perfil e a hipótese semi-infinita é razoável. Vou calcular
e comparar com
.
Passo 3
Cálculo da penetração difusiva e da concentração relativa em
.
Calcular
:
Ou seja,
Isso é bem menor que a espessura do filme
, então tratar o filme como semi-infinito para a penetração em
é razoável.
Agora calcule o argumento da
para
:
Usando uma aproximação assintótica ou tabelas/uma calculadora da função erro, acha-se
Portanto
ou seja, cerca de
da concentração de equilíbrio na posição a
da superfície, após
de exposição.
Vou explicar rapidamente como obtive esse número: para argumentos em torno de
–
,
é baixinho. Uma estimativa assintótica razoável é
para
; aplicando isso com
dá o valor aproximado acima, compatível com tabelas numéricas.
Comentário prático: quer dizer que a substância
só penetrou pouco no líquido em 55 s — a difusividade é baixa e a posição
está fora da “zona atingida” substancialmente.
Passo 4
Cálculo do comprimento do filme exposto ao gás
(pergunta (b)). Modelo do escoamento laminar de filme que desce por uma placa vertical (solução de Nusselt).
Para um filme em escoamento laminar sob ação da gravidade numa parede inclinada de ângulo
(vertical
), a equação de momento reduz a um balanço simples que leva à velocidade perfil parabólico:
onde
é a coordenada perpendicular à superfície,
na parede e
na interface. Para obter a velocidade média no filme (média sobre a espessura) calculamos o fluxo por unidade de largura
e dividimos por
:
A velocidade média
é
, portanto
É comum reescrever em termos da viscosidade cinemática
como
Assumindo parede vertical (
):
Comentário pedagógico curto: essa expressão vem de integrar a equação de Navier–Stokes em regime permanente, supondo que a variação ao longo da placa é suave e que as forças viscosas equilibram a componente tangencial da gravidade — a solução é clássica e chamada de perfil de Nusselt para filme fino.
Passo 5
Cálculos numéricos do comprimento exposto
.
Calcular a viscosidade cinemática:
Calcular
:
Numerador:
. Denominador:
. Portanto
Logo
.
Agora o comprimento do filme exposto é simplesmente o deslocamento médio do fluido no tempo
:
Ou seja,
, que são aproximadamente
.
Observações importantes sobre esse resultado:
- O valor numérico é extremamente grande e fisicamente surpreendente. Ele decorre diretamente dos dados fornecidos: a viscosidade dinâmica
combinada com a densidade relativamente baixa implica uma viscosidade cinemática muito pequena , o que fornece velocidades médias elevadas segundo a fórmula de Nusselt. Se na prática a viscosidade real for maior ou a placa não for perfeitamente vertical, o comprimento real será muito menor. - Se alguém obtiver esse exercício na prova, convém comentar (como fiz) que o número final deve ser verificado com sentido físico; mas matematicamente, com os dados fornecidos e as hipóteses usuais, este é o procedimento correto.
Resposta
(a)
(b)