Compressible and immiscible fluids with arbitrary density ratio

Este artigo apresenta uma nova teoria baseada na termodinâmica para modelar fluidos compressíveis e imiscíveis com razões de densidade arbitrariamente altas, abordando as limitações das equações tradicionais de Navier-Stokes e Euler ao derivar equações de evolução de densidade que consideram adequadamente a verdadeira evolução do momento.

O essencial

Imagine que você esteja tentando descrever uma dança entre dois parceiros muito diferentes: um elefante pesado e de movimentos lentos (água) e uma pena leve e de movimentos rápidos (ar). No mundo da física de fluidos, eles são "imiscíveis", o que significa que não se misturam como óleo e água, mas possuem uma fronteira difusa onde se encontram.

Por muito tempo, cientistas lutaram para escrever as "regras da dança" (as equações) para sistemas onde a diferença de densidade é enorme, como ser 1.000 vezes mais pesado para o elefante do que para a pena. As regras antigas tinham uma falha importante: tratavam o peso do elefante como se fosse um número constante e imutável, mesmo justamente na fronteira difusa onde o elefante se transforma em uma pena. Isso é como tentar descrever uma pessoa se transformando em um fantasma dizendo: "Eles ainda são um humano sólido durante todo o tempo", o que não faz sentido.

Aqui está uma explicação simples do que este artigo propõe para corrigir esse problema.

1. O Problema com as Regras Antigas

A forma tradicional de calcular como os fluidos se movem (usando as equações de Navier-Stokes e Euler) baseia-se em um atalho chamado aproximação de Boussinesq.

• A Analogia: Imagine que você está empurrando um carrinho de compras. Se o carrinho estiver cheio de tijolos, ele é pesado. Se estiver vazio, é leve. As regras antigas assumiam que, se você estivesse empurrando um carrinho que está parcialmente cheio de tijolos e parcialmente cheio de ar, o peso do carrinho nunca mudaria enquanto você o empurra. Elas simplesmente assumem que o peso é uma média fixa.
• A Falha: Na realidade, conforme o carrinho se move e os tijolos se deslocam (difusão), o peso muda. As regras antigas ignoravam o fato de que o "momento" (massa $\times$ velocidade) muda porque a própria massa está mudando na fronteira. Elas também assumiam que o ar e a água ocupam uma quantidade perfeitamente previsível de espaço, ignorando que, quando eles se misturam na fronteira, o espaço que ocupam pode, na verdade, oscilar e mudar devido à pressão.

2. A Nova Abordagem: Minimização de Energia

Em vez de adivinhar como a densidade muda, o autor parte de um princípio fundamental: A Natureza sempre tenta usar a menor quantidade de energia possível.

• A Analogia: Pense em uma bola rolando montanha abaixo. Ela não se importa com os passos específicos que dá; ela apenas quer chegar ao fundo (menor energia). O autor usa esse conceito de "colina de energia" para derivar novas regras sobre como a água e o ar interagem.
• A Inovação Principal: O autor introduz um conceito chamado "Volume em Excesso".
  • Imagine que você tem um balde de água e um balde de ar. Se você os despejar juntos, pode esperar que o volume total seja exatamente a soma dos dois baldes. Mas, ao nível microscópico, quando eles se encontram, as moléculas podem se compactar mais ou menos, criando um pouco de espaço "extra" ou "faltante".
  • As regras antigas assumiam que esse espaço extra era zero em todos os lugares. Este artigo diz: "Não, esse espaço extra existe, ele muda de um lugar para outro e afeta a densidade".

3. As Novas "Regras da Dança" (Os Resultados)

Ao contabilizar essa mudança de "espaço extra" e usar a minimação de energia, o autor deriva um novo conjunto de equações que faz três coisas principais:

A. Uma Nova Forma de Ouvir o Som (Velocidade do Som)
O artigo mostra que a velocidade do som não é apenas um número aleatório; ela vem diretamente de como a energia do fluido muda conforme ele é espremido.

• A Metáfora: Pense no som como uma ondulação em uma multidão. A velocidade dessa ondulação depende de quão apertadas as pessoas (moléculas) estão compactadas e quanta energia elas têm. A nova fórmula calcula essa velocidade naturalmente, sem precisar que lhe digam qual é o valor previamente. Ela até sugere que, em um gás, a velocidade do som é aproximadamente a mesma velocidade média com que as moléculas do gás estão saltitando de um lado para o outro.

B. Uma Nova Regra para Pressão e Velocidade (Lei de Bernoulli)
Você provavelmente já ouviu o princípio de Bernoulli: "Quando um fluido se move mais rápido, sua pressão cai".

• A Reviravolta: A regra antiga funciona muito bem para a água fluindo em um tubo, mas ela falha quando há um salto enorme na densidade (como a água atingindo o ar). O autor cria uma Lei de Bernoulli Generalizada.
• A Metáfora: Imagine um rio fluindo para uma queda d'água. A regra antiga diz que a energia permanece a mesma. A nova regra diz: "Espere, conforme a água se transforma em névoa (ar), parte da energia é perdida ou transformada porque a 'densidade/estupidez' da água está mudando". A nova equação contabiliza essa mudança de energia, tornando-a precisa mesmo quando o fluido está mudando sua natureza de pesado para leve.

C. O "Calombo" na Densidade
Este é talvez o resultado mais visual.

• A Visão Antiga: Se você observar a fronteira entre a água e o ar, os modelos antigos diziam que a densidade apenas deslizaria suavemente de "água pesada" para "ar leve", como uma rampa descendente.
• A Nova Visão: A matemática do autor prevê um calombo. Ao cruzar a fronteira, a densidade na verdade sobe ligeiramente antes de cair para o nível do ar.
• A Metáfora: Imagine uma multidão de pessoas (água) tentando passar por uma porta para um corredor vazio (ar). Conforme elas se espremem pela porta, podem se compactar mais fortemente por um breve momento antes de se espalharem. A nova teoria prevê esse "calombo de compactação", o que coincide com o que simulações de computador avançadas (chamadas de Teoria do Funcional da Densidade) já observaram, mas os modelos simples antigos perderam.

Resumo

Este artigo propõe uma nova maneira de escrever as leis da física para fluidos que são muito diferentes em peso (como água e ar).

1. Ele deixa de fingir que o peso é constante na fronteira.
2. Ele admite que o espaço que as moléculas ocupam muda (volume em excesso).
3. Ele usa o princípio do "menor gasto de energia" para derivar novas regras que explicam como o som viaja, como a pressão muda e por que a densidade na fronteira entre água e ar possui um pequeno "calombo".

O autor afirma que este novo framework funciona para qualquer mistura de fluidos, não importa o quão diferentes sejam seus pesos, abrindo as portas para simulações de computador mais precisas de coisas como chuva caindo, ondas quebrando ou bolhas subindo.

Resumo técnico

Resumo Técnico: Fluidos Compressíveis e Imiscíveis com Razões de Densidade Arbitrárias: Teoria da Interface Difusa

Problema
Os modelos atuais de dinâmica de fluidos computacional (CFD) têm dificuldade em simular com precisão sistemas de fluidos imiscíveis com altas razões de densidade, como água-ar (razão $\approx 1000$). As abordagens existentes dependem de dois métodos principais, ambos os quais possuem limitações significativas:

1. A Aproximação de Boussinesq: Assume uma densidade ($\rho$) constante dentro da interface, utilizando a equação $\rho(du/dt) = -\nabla p + \dots$. Esta formulação falha em considerar a "verdadeira" evolução do momento $d(\rho u)/dt$, negligenciando mudanças de densidade e a resultante transferência de momento. Além disso, não consegue gerar naturalmente forças de gravidade e empuxo sem adições ad hoc, pois estas surgem de diferenças de densidade.
2. O Método de Interpolação: Define a densidade como uma função linear da concentração de volume ($\rho = \rho_w \phi + \rho_a(1-\phi)$). Esta abordagem cria um conflito fundamental entre a incompressibilidade ($d\rho/dt = 0$) e a difusão ($d\phi/dt \neq 0$). Além disso, assume um excesso de volume uniforme (frequentemente zero) através de todo o domínio, o que contradiz a realidade física onde o excesso de volume local depende da pressão local. Consequentemente, as equações de continuidade padrão e a interpolação linear falham em capturar perfis de densidade não monotônicos observados na interface água-ar.

Metodologia
Os autores propõem um novo arcabouço teórico fundamentado no princípio termodinâmico da minimização de energia, distinto da derivação tradicional das equações de Navier-Stokes e Euler. A metodologia envolve:

• Reformulação de Volume e Densidade: Em vez de assumir um excesso de volume zero, o modelo introduz um excesso de volume espacialmente variável ($v_e$) dentro da interface. Este volume é redistribuído entre os componentes fluídicos, levando a densidades parciais ($\rho'_w, \rho'_a$) que não são mais constantes, mas dependem da pressão local e da mistura. Isso exige duas equações de evolução acopladas: uma para a concentração de volume ($\phi$) e outra para a densidade ($\rho$).
• Princípio de Dissipação-Conservação: Os autores estabelecem um teorema ligando a primeira lei da termodinâmica (conservação) com a segunda lei (dissipação). Ao definir a energia total do sistema ($E$) como a soma da energia livre ($F$), energia de pressão ($P$), energia cinética ($K$) e energia livre de parede ($W$), eles derivam equações cinéticas minimizando a dissipação de energia.
• Relação de Gibbs-Duhem Modificada: Para resolver incompatibilidades matemáticas entre gradientes espaciais e derivadas funcionais em modelos existentes, uma relação de Gibbs-Duhem isotérmica modificada é derivada. Esta considera explicitamente a entropia dependente da velocidade e o trabalho inercial, levando a uma definição rigorosa de pressão que inclui contribuições da energia livre e da entropia de velocidade.
• Definição de Energia Cinética: A energia cinética é formulada como $K = \int \rho u \cdot u \, d\Omega$ (densidade de momento linear $\zeta = \rho u$) em vez do convencional $\frac{1}{2}\rho u^2$. Esta distinção permite a separação do balanço de momento (mudando tanto $u$ quanto $\rho$) da dissipação de velocidade (mudando apenas $u$).

Principais Contribuições
A derivação resulta em uma equação de evolução de densidade generalizada e três grandes avanços teóricos:

1. Lei de Bernoulli Generalizada: O princípio de Bernoulli padrão ($p + \frac{1}{2}\rho u^2 = \text{constante}$) mostra-se um caso especial. A forma generalizada contabiliza mudanças de entropia de velocidade e a evolução da densidade na interface: $p + \rho u^2 = \text{constante}$ (sob condições específicas de estado estacionário invariante), reconhecendo que saltos de energia ocorrem através das interfaces devido às diferenças de densidade.
2. Formulação Universal da Velocidade do Som: Uma expressão geral para a velocidade do som é derivada com base na derivada da energia livre em relação à densidade: $c = \sqrt{d\Psi/d\rho}$. Para gases ideais, isso aproxima a velocidade térmica média das moléculas.
3. Equação de Estado (EOS) Generalizada: A teoria demonstra que a equação de evolução da densidade gera naturalmente uma EOS para os fluidos (ex: $p \sim \rho \ln \rho$ para gases ideais) sem a necessidade de uma EOS prescrita, contrastando com os métodos tradicionais de interface difusa.

Resultos e Validação
O artigo valida a teoria através de três estudos de caso:

• Propagação do Som: A equação derivada da velocidade do som é consistente com a análise de perturbação e prevê corretamente a relação de dispersão das ondas sonoras no ar.
• EOS de Gás Ideal: No estado estacionário, o modelo reproduz a lei dos gases ideais, confirmando que a equação de evolução da densidade captura inerentemente as relações de estado termodinâmicas.
• Dinâmica da Interface Água-Ar: Uma solução analítica para uma interface plana água-ar em movimento num tubo demonstra que o modelo prevê um perfil de densidade não monotônico. Especificamente, a densidade aumenta do valor da água no bulk para um pico máximo no centro da interface antes de diminuir para o valor do gás no bulk. Este resultado alinha-se com as descobertas da Teoria do Funcional da Densidade (DFT) [11, 12] e contradiz a diminuição monotônica prevista pelos métodos de interpolação linear. O modelo também reproduz corretamente a equação de Young-Laplace para interfaces curvas.

Significância
O artigo afirma fornecer um arcabouço geral para fluidos imiscíveis com razões de densidade arbitrariamente altas, aplicável tanto aos regimes compressível quanto incompressível. Ao tratar a evolução da densidade como uma variável dinâmica acoplada ao excesso de volume e à dissipação viscosa, o modelo resolve as contradições inerentes à aproximação de Boussinesq e aos métodos de interpolação linear. O trabalho sugere que a dissipação viscosa altera inerentemente a densidade, desafiando a estrita suposição de incompressibilidade em fluxos dissipativos. Os autores postulam que esta teoria abre uma nova janela para o CFD, oferecendo uma descrição fisicamente consistente da transferência de momento e energia em interfaces fluido-fluido onde as razões de densidade são extremas.