Escoamento em meios porosos

Teoria Completa

Introdução

Antes de começar a ver as equações para meios porosos, queria revisar dois operadores com você, porque vamos precisar bastante deles. São eles: gradiente e divergente.

- Operador gradiente: aplicamos em um escalar e obtemos como resultado um vetor. Considerando eixos cartesianos e usando a temperatura como exemplo T , ele fica:

g r a d T = ∇ →   T = ∂ T ∂ x   e → x + ∂ T ∂ y   e → y + ∂ T ∂ z   e → z

Em que e → x ,   e → y ,   e → z são os vetores unitários de cada eixo.

Aplicando o gradiente em uma propriedade, obtemos a informação de como essa propriedade varia no espaço, em cada um dos eixos.

- Operador divergente: aplicamos em um vetor e obtemos um escalar. Usando a velocidade como exemplo e considerando eixos cartesianos:

d i v u → = ∇ → . u → = ∂ u → x ∂ x   + ∂ u → y ∂ y   + ∂ u → z ∂ z

Em que u → x ,   u → y ,   u → z são as componentes x ,   y ,   z da velocidade, respectivamente.

E ah! Repara que tem um pontinho entre ∇ → . u → , que simboliza o produto escalar.

Equações da continuidade e Movimento

Lá da mecânica dos fluidos, podemos escrever equação da continuidade (é um balanço de massa) e equação do movimento (é um balanço de forças) para o FLUIDO.

E q .   d a   c o n t i n u i d a d e :       ∂ ρ ∂ t + d i v ρ u → = 0

E q .   d o   M o v i m e n t o :     ρ ∂ u → ∂ t + g r a d   u → u → = - g r a d   P + d i v   τ ̿ + ρ g →

Em que:

u →   é   o   v e t o r   v e l o c i d a d e τ ̿   é   o   t e n s o r   t e n s ã o   ( r e l a c i o n a d o   a   f o r ç a s   v i s c o s a s )

Mas temos um problema: o meio poroso é todo louco, sem nenhuma simetria; o fluido vai indo por vários poros diferentes e aleatórios. Além disso, definir a velocidade ponto a ponto seria quase impossível, assim como definir as paredes. Precisamos achar outro jeito de continuar nosso desenvolvimento pra meios porosos.

Teoria das misturas da mecânica do contínuo

A solução é simples até: vamos reescrever essas equações, mas, agora, considerando não só o fluido mas sim FLUIDO E SÓLIDO, em uma proporção igual à porosidade ( ε ) . Tipo, se ε = 0,3 consideramos um balanço pegando 30 % de fluido e 70 % de sólido. As equações ficam:

E q .   d a   c o n t i n u i d a d e :       ∂ ( ε ρ ) ∂ t + d i v ε ρ u → = 0

E q .   d o   M o v i m e n t o :     ε ρ ∂ u → ∂ t + g r a d   u → u → = - g r a d   P + d i v   τ ̿ - M → + ρ g →

Em que:

M →   é   a   f o r ç a   q u e   o   l í q u i d o   e x e r c e   n a   m a t r i z   s ó l i d a

Daí, - M → é a reação que a matriz faz no fluido. E repara que escrevemos a equação do movimento considerando apenas as forças que atuam no fluido.

Mas ainda temos o problema de não conseguir definir u → ponto a ponto...e não tem muito o que fazer não; vamos adotar a estratégia de perder esse detalhamento pra conseguir resolver. Vamos definir uma velocidade superficial q → , considerando a área da do tubo (ou outro lugar) em que o meio poroso está inserido;

O módulo de q → vale:

q → = Q A

Em que:

Q   é   a   v a z ã o   v o l u m é t r i c a A   é   a   á r e a   t r a n v e r s a l   o c u p a d a   p e l o   m e i o   p o r o s o

q → é como se fosse a velocidade antes e depois do meio poroso. Por isso ele recebe esse nome de “velocidade superficial”. O fluido tem essa velocidade até a superfície do meio poroso. É uma condição de contorno. Podemos relacioná-la com u → da seguinte maneira:

q → = ε u →

Juntando as duas expressões:

Q A =   ε u →

Daí:

u → = Q ε A

Com base em q → , as equações ficam:

E q .   d a   c o n t i n u i d a d e :       ∂ ( ε ρ ) ∂ t + d i v ρ q → = 0

E q .   d o   M o v i m e n t o :     ρ ε ∂ ∂ t q → ε + g r a d   q → ε q → = - g r a d   P + d i v   τ ̿ - M → + ρ g →

E são essas equações que vamos, efetivamente, usar.

Equações para M →

Até agora sabemos da existência de uma forca entre líquido e sólido (força M → ), mas não conhecemos fórmulas pra ela. Bora ver duas correlações empíricas:

- Correlação de Forchheimer

Vou começar da correlação mais geral, ok? Para meio isotrópico (quando as propriedades físicas não variam com a posição), Forchheimer diz que:

M → = μ K 1 + C K ρ q → μ q →

Em que

K   é   a   p e r m e a b i l i d a d e   d o   m e i o   ( t e m   u n i d a d e   d e   á r e a ) C   é   o   f a t o r   s e m   n o m e ,   o u   f a t o r   d e   F o r c h h e i m e r ,   e   é   a d i m e n s i o n a l

Sendo que

R e M P = C K ρ q → μ

Em que R e M P é o Reynolds do meio poroso.

E falando em permeabilidade, deixa eu apresentar logo pra você uma unidade “especial”, que é o Darcy ( D ) .

1 D = 10 - 8   c m 2

- “Lei” de Darcy

A lei de Darcy é um caso especial da correlação que a gente acabou de ver, mas para situações onde (além do meio isotrópico e fluido Newtoniano) estamos falando de escoamento lento, com Reynolds baixo. Mais precisamente, a lei de Darcy é válida quando:

R e M P ≪ 1

E a “lei” fica com essa carinha:

M → = μ K q →

Determinação experimental de K e C

Para a determinação experimental de K e C , usamos um permeâmetro; nada mais é do que um tubo contendo o meio poroso; variamos a vazão Q de um fluido que passa pelo meio poroso e anotamos a queda de pressão ∆ P do fluido ao passar pelo meio poroso (MP) com essa tal vazão. Para isso, medimos a pressão antes e depois do meio poroso, né.

Como ao escoar pelo meio poroso o fluido perde carga:

P 0 > P L

Depois de fazer os experimentos e ter uma tabelinha relacionando Q e ∆ P , vamos usar a equação do movimento para determinar K e C . Acontece que o fluido pode ser incompressível ou compressível e o escoamento pode ser lento (darcyano) ou um escoamento qualquer (aí usa Forchheimer para M → ). Vou te mostrar agora como trabalhar com cada um desses casos; em todos eles, vamos assumir:

  • Fluido Newtoniano
  • Escoamento unidimensional, permanente (estado estacionário) e isotérmico
  • Meio poroso isotérmico (Temperatura constante)

  1. Escoamento de fluido incompressível
  2. Vamos pegar a equação do movimento e fazer as simplificações que der:

    1. Em estado permanente não há variação no tempo
    2. Gradiente nos diz como uma propriedade varia no espaço; assumimos meio poroso isotrópico daí NÃO há variação no espaço
    3. Não está listado nas hipóteses, mas vamos assumir que as forças viscosas são desprezíveis face a M →
    4. O escoamento é na horizontal e a gravidade está na vertical. É importante lembrar que isso é para ESSE exemplo que eu estou desenvolvendo com você, ok? Se o escoamento fosse na mesma direção que a gravidade, ia ter que considerar esse termo!

    Depois das simplificações, só sobra:

    - g r a d   P - M → = 0

    Mas assumimos escoamento unidimensional, lembra? Chamando o eixo da direção do escoamento de z , o gradiente se resume a d P d z . Daí:

    - d P d z = M →

    Mas esse escoamento incompressível pode ser lento (darcyano) ou um escoamento qualquer; bora ver cada caso:

    1. Escoamento Darcyano (lento)
    2. Se o escoamento é Darcyano, vale:

      M → = μ K q →

      Daí, voltando na equação do movimento com as simplificações:

      - d P d z = μ K q →

      Separando as variáveis e integrando do inicio ao fim do meio poroso, chamando o módulo de q → de q :

      - ∫ P 0 P L d P = μ K q ∫ 0 L d z

      Repara que μ ,   K e q saíram da integral porque são constantes; μ é constante porque assumimos isotérmico, K é constante porque assumimos isotrópico e q é constante porque a área transversal do meio poroso é constante.

      Resolvendo a integral:

      P 0 - P L = μ K q L

      Fazendo ∆ P = P 0 - P L , ou seja, maior menos menor, e rearranjando:

      ∆ P L = μ K q

      Isso é uma reta! Tendo os pontos relacionando ∆ P e q (conseguimos tirar q a partir de Q ), podemos fazer regressão linear:

    3. Escoamento qualquer

    Quando não pudermos assumir escoamento darcyano, temos que usar Forchheimer para M → :

    M → = μ K 1 + C K ρ q → μ q →

    Daí a equação do movimento com as simplificações e já integrada fica:

    ∆ P L = μ K q + μ K C K ρ q 2 μ