An LES model with finite-rate phase change and subgrid spray based on a thermodynamically consistent four-equation multiphase model

Este trabalho propõe um modelo de LES com mudança de fase de taxa finita e spray de subgrade, baseado em um modelo multifásico de quatro equações termodinamicamente consistente, que utiliza uma nova formulação do modelo $\Sigma$ para prever a área superficial interfacial e garantir a precisão em simulações de combustão.

O essencial

O Mistério do Spray: Como simular gotas e vapor com precisão

Imagine que você está tentando filmar uma festa de aniversário com um canhão de confetes. Se você usar uma câmera lenta super potente, verá cada pedacinho de papel voando. Mas, se você estiver tentando filmar uma multidão correndo em um estádio, não dá para focar em cada pessoa individualmente; você vê apenas um "fluxo" de gente.

Na engenharia, temos o mesmo problema. Quando cientistas estudam como o combustível é injetado em um motor de carro ou em um foguete, eles lidam com um "spray" (um borrifo) de gotas minúsculas e vapor. Tentar calcular o movimento de cada molécula de cada gotinha seria tão difícil que nem o supercomputador mais rápido do mundo daria conta.

O que este artigo resolveu? Os pesquisadores criaram um "atalho inteligente" (um modelo matemático) que permite simular esse borrifo de forma rápida, mas sem perder a precisão científica.

---

1. O Problema do "Equilíbrio" (A Analogia da Festa)

Imagine que você tem uma sala com pessoas (o líquido) e um corredor vazio (o gás). Se a sala ficar muito cheia e quente, as pessoas começam a sair para o corredor. Isso é a evaporação.

Existem dois jeitos de simular isso:

1. O jeito "Instantâneo": Você assume que, assim que a sala esquenta, todo mundo sai correndo para o corredor ao mesmo tempo. Isso é rápido de calcular, mas é irreal, porque na vida real as pessoas levam um tempo para se mover.
2. O jeito "Lento": Você tenta calcular cada passo de cada pessoa. É perfeito, mas demora uma eternidade.

A solução do artigo: Eles criaram um modelo de "Taxa Finita". É como se eles dissessem ao computador: "Não assuma que todos saem de uma vez, mas também não perca tempo contando cada passo. Apenas use uma regra que diga que o movimento depende da temperatura e do espaço disponível". Eles garantiram que esse movimento nunca fosse "impossível" (como pessoas surgindo do nada), mantendo tudo dentro das leis da termodinâmica.

---

2. O Problema da "Área de Contato" (A Analogia do Açúcar)

Pense no seguinte: se você colocar um cubo de açúcar em um copo de café, ele demora para dissolver. Mas se você usar açúcar refinado (que tem muito mais "superfície" exposta), ele dissolve quase instantaneamente.

No spray de combustível, a velocidade com que o líquido vira gás depende de quanta "superfície" as gotas têm. O problema é que, em simulações de computador, as gotas são tão pequenas que o computador não consegue "vê-las" individualmente. Elas ficam invisíveis para a grade de cálculo.

A solução do artigo: Eles criaram um modelo chamado $\Sigma$ (Sigma). Imagine que, em vez de desenhar cada grão de açúcar, o computador apenas calcula uma "nuvem de área de contato". Ele estima: "Nesta região, embora eu não veja as gotas, eu sei que existe uma quantidade X de superfície de contato disponível para o combustível evaporar". Isso permite que o computador saiba o quão rápido o combustível vai queimar sem precisar desenhar cada gotinha.

---

3. Como eles testaram? (O Teste do Spray de Diesel)

Para saber se o "atalho" deles funcionava, eles testaram o modelo contra experimentos reais de um spray de n-dodecano (um tipo de combustível).

• Primeiro teste (Sem evaporação): Eles viram se o computador conseguia prever o formato do jato de líquido. O resultado? O computador "desenhou" o jato quase igual ao que aconteceu no laboratório.
• Segundo teste (Com evaporação): Eles esquentaram o ambiente para o combustível evaporar. O computador conseguiu prever exatamente o quanto o líquido avançava e o quanto o vapor se espalhava.

Resumo da Ópera

Os cientistas construíram uma ferramenta que é como um mapa inteligente: em vez de tentar desenhar cada árvore de uma floresta (o que levaria anos), eles criaram um mapa que mostra onde as florestas são densas e onde o vento sopra mais forte, permitindo prever o clima da floresta com uma precisão incrível e em uma fração do tempo.

Isso é fundamental para criar motores de carros mais eficientes, foguetes que gastam menos combustível e sistemas de combate a incêndios mais eficazes!

Resumo técnico

Resumo Técnico: Modelo LES com Mudança de Fase de Taxa Finita e Spray de Subgrade

Título Original: An LES model with finite-rate phase change and subgrid spray based on a thermodynamically consistent four-equation multiphase model

1. O Problema (Contexto e Motivação)

A modelagem de sistemas multifásicos multicomponentes em escala de engenharia (como motores de combustão interna, propulsão de foguetes e sistemas de supressão de incêndio) é extremamente complexa. Os desafios principais incluem:

• Custo Computacional: Modelos de não-equilíbrio total (que resolvem todas as variáveis de cada fase separadamente) são proibitivamente caros para simulações de larga escala.
• Limitações de Escala (LES): Em simulações de Grandes Escalas (LES), muitas estruturas de interface (gotículas e ligamentos) são menores que a malha computacional, exigindo modelos de subgrade para representar o spray e a área interfacial.
• Inconsistência Termodinâmica: Modelos de mudança de fase existentes (como o HEM - Homogeneous Equilibrium Model) frequentemente assumem equilíbrio instantâneo, o que superestima a taxa de transferência de massa em regimes onde a interface é finita e não resolvida. Além disso, modelos de taxa finita muitas vezes não respeitam os limites impostos pelo potencial químico, podendo gerar estados físicos impossíveis.

2. Metodologia

Os autores propõem um framework baseado em um modelo de quatro equações, que assume equilíbrio mecânico (pressão) e térmico (temperatura) entre as fases, reduzindo drasticamente o custo computacional sem perder a robustez. A metodologia divide-se em três pilares:

• Mudança de Fase de Taxa Finita Termodinamicamente Limitada: Em vez de assumir equilíbrio instantâneo, eles utilizam o modelo de Hertz-Knudsen para a taxa de transferência de massa. Para garantir a consistência, eles introduzem um novo mecanismo onde a taxa de mudança de fase é limitada pelo potencial químico (energia livre de Gibbs). Isso garante que o sistema nunca ultrapasse o estado de equilíbrio termodinâmico, evitando resultados não físicos.
• Modelo de Spray de Subgrade ($\Sigma$): Para representar a área interfacial não resolvida pela malha, utiliza-se o modelo de densidade de superfície ($\Sigma$). Os autores propõem uma forma "phase-confined" (confinada à fase) da equação de transporte de $\Sigma$, o que evita o vazamento não físico de massa/área através das interfaces de fase.
• Termos de Fonte de Spray: O modelo inclui termos de fonte para representar processos físicos de subgrade:
  • Regime Denso: Quebra turbulenta (turbulent breakup) e flashing (ebulição rápida por queda de pressão).
  • Regime Diluído: Quebra secundária de gotas e coalescência.
  • Indicadores de Regime: O modelo utiliza indicadores baseados na fração de volume e na escala de Hinze para alternar automaticamente entre os modelos de spray denso e diluído.

3. Principais Contribuições

1. Consistência Termodinâmica: Desenvolvimento de um modelo de mudança de fase de taxa finita que é matematicamente limitado pelo equilíbrio da energia livre de Gibbs.
2. Modelo $\Sigma$ Confinado: Uma nova formulação da equação de transporte de área interfacial que impede o vazamento de propriedades entre fases, essencial para a estabilidade em simulações LES.
3. Framework Integrado: A união de um modelo de quatro equações (eficiente) com modelos de subgrade de alta fidelidade (precisos), permitindo simular desde sprays densos até zonas de evaporação diluída.
4. Modelo de Flashing para LES: Uma abordagem para modelar a produção de área superficial durante eventos de flashing que também respeita os limites de transferência de massa termodinâmica.

4. Resultados e Validação

O modelo foi validado utilizando o caso de estudo ECN (Engine Combustion Network) Spray A (injeção de n-dodecano), comparando-o com dados experimentais em duas condições:

• Caso Não-Evaporante: O modelo de spray ($\Sigma$) previu com excelente precisão a densidade de massa projetada (PMD), a densidade de área superficial (PSA) e o diâmetro médio de Sauter (SMD), demonstrando que a dinâmica do spray e a quebra de gotas foram capturadas corretamente.
• Caso Evaporante: O modelo acoplado previu com sucesso a penetração do líquido e a extensão da nuvem de vapor. Embora tenha havido uma leve superestimação da penetração do gás (atribuída ao uso da lei dos gases ideais para fluidos supercríticos), os resultados mostraram concordância satisfatória com as medições experimentais de penetração de líquido e vapor.

5. Significância

Este trabalho representa um avanço significativo para a simulação de combustão e fluidos multifásicos. Ele oferece uma ferramenta que é, ao mesmo tempo, computacionalmente eficiente (devido ao modelo de quatro equações) e fisicamente rigorosa (devido ao tratamento termodinâmico da mudança de fase e ao confinamento da área de spray). Isso permite que engenheiros simulem processos de injeção de combustível em motores e sistemas de propulsão com uma fidelidade muito próxima da realidade, mas em escalas de tempo e tamanho viáveis para o design industrial.