A Volume of Fluid Immersed Boundary Method for Industrial Polymer Mixing

Este artigo apresenta um novo solucionador acoplado por blocos de Volume de Fluido com Fronteira Imersa (BC-VOF-IB) implementado no OpenFOAM que supera instabilidades numéricas causadas por altos contrastes de viscosidade para simular com precisão a mistura de polímeros com superfície livre em extrusoras industriais parcialmente preenchidas.

O essencial

Imagine que você está tentando misturar um tanque gigante de mel grosso e pegajoso com ar dentro de uma máquina giratória. Isso é essencialmente o que acontece na mistura industrial de polímeros, onde empresas como a Pirelli precisam combinar plástico derretido com aditivos para fabricar pneus, dispositivos médicos ou peças automotivas. O objetivo é misturar tudo perfeitamente para que o produto final seja forte e uniforme.

No entanto, simular esse processo em um computador é um pesadelo para matemáticos e engenheiros. Eis o porquê e como este artigo resolve isso, usando analogias simples:

O Problema: A Luta "Mel Grosso vs. Ar Fino"

Nessas máquinas, você tem dois fluidos muito diferentes:

1. Polímero Fundido: Extremamente grosso, pegajoso e de movimento lento (como mel frio).
2. Ar: Muito fino e de movimento rápido.

Quando você tenta simular como esses dois interagem dentro de uma máquina com parafusos giratórios, os programas de computador padrão ficam confusos. É como tentar calcular o movimento de um caracol e de um carro de corrida na mesma pista usando o mesmo conjunto de regras. O computador tenta dar passos minúsculos para manter o "caracol" (o plástico grosso) de se mover rápido demais, o que torna a simulação incrivelmente lenta — às vezes levando dias para finalizar alguns segundos de mistura em tempo real.

Além disso, as máquinas possuem partes complexas e giratórias (parafusos) que se movem dentro de um recipiente fixo. Tradicionalmente, para simular isso, você precisa construir uma malha digital (uma grade de caixas minúsculas) que se encaixe perfeitamente ao redor dos parafusos giratórios. À medida que os parafusos giram, essa grade precisa se remodelar constantemente, o que é como tentar tricotar um suéter enquanto a pessoa que o veste está correndo uma maratona. É bagunçado, difícil e propenso a erros.

A Solução: Uma Nova "Grade Inteligente" e uma "Abordagem em Equipe"

Os autores deste artigo desenvolveram uma nova maneira de executar essas simulações usando um software chamado OpenFOAM. Eles combinaram duas técnicas poderosas:

1. O Método da Fronteira Imersa (O Truque da "Parede Fantasma")
Em vez de remodelar a grade para se ajustar aos parafusos giratórios, eles mantiveram a grade fixa e rígida (como um bloco sólido de gelo). Em seguida, disseram ao computador: "Ei, há um parafuso giratório dentro deste bloco de gelo."

• A Analogia: Imagine uma piscina com uma grade fixa de azulejos no fundo. Em vez de mover os azulejos para se ajustar a um nadador, você apenas diz à água: "Não atravesse o nadador." O computador usa matemática para criar uma "parede fantasma" ao redor do parafuso, forçando o fluido a fluir ao seu redor sem nunca precisar reconstruir a grade. Isso torna o manuseio de formas complexas e em movimento muito mais fácil.

2. O Método do Volume de Fluido (VOF) (O Truque da "Tinta de Rastreamento")
Para ver onde o plástico grosso termina e o ar começa, eles usam uma "tinta" que preenche as células.

• A Analogia: Imagine que a grade do computador é um tabuleiro de xadrez 3D. Alguns quadrados são 100% plástico, alguns são 100% ar e alguns são uma mistura. O computador rastreia quanto "tinta de plástico" há em cada quadrado para ver a superfície do líquido.

3. O Esquema Acoplado por Blocos (O "Huddle" da Equipe)
Esta é a descoberta mais importante. Nas simulações padrão, o computador resolve a velocidade do fluido nas direções X, Y e Z uma por uma, como três pessoas falando por vez. Quando o fluido é supergrosso (como polímero), essa abordagem de "falar por vez" faz a simulação falhar ou ficar extremamente lenta, porque o fluido grosso acopla todas as direções juntas de forma rígida.

Os autores mudaram isso para uma abordagem Acoplada por Blocos.

• A Analogia: Em vez de três pessoas falando por vez, elas se reúnem em um círculo e resolvem o problema juntas exatamente ao mesmo tempo. Ao tratar o movimento em todas as direções como uma única equipe gigante e interconectada, o computador consegue lidar com a enorme diferença entre o plástico grosso e o ar fino sem ficar preso.

Os Resultados: De Horas para Minutos

A equipe testou seu novo método em dois cenários:

1. Um Canal em Forma de Osso de Cachorro: Um caso de teste onde plástico grosso é injetado em um canal estreito e sinuoso.

   • Antigo Método: O programa de computador padrão falhou ou levou 7 horas para simular alguns segundos porque foi forçado a dar passos minúsculos.
   • Novo Método: Seu novo método "Huddle da Equipe" concluiu o mesmo trabalho em apenas 16 minutos e não falhou, mesmo quando o plástico ficou extremamente grosso.

2. Máquinas Industriais Reais: Eles simularam extrusoras reais de parafuso único e parafuso duplo (as máquinas usadas para fazer grânulos de plástico).

   • Eles mostraram com sucesso como o plástico preenche a máquina, como a pressão se acumula e como o ar é empurrado para fora.
   • Eles provaram que seu método de "Parede Fantasma" funciona tão bem quanto o antigo e difícil método de remodelar a grade, mas muito mais rápido e fácil de configurar.

E Agora?

O artigo conclui que este é um grande passo à frente para a indústria. Ele preenche a lacuna entre a matemática acadêmica e as necessidades reais das fábricas. No entanto, os autores observam que seu modelo atual assume que a temperatura permanece constante (isotérmica). Na realidade, a mistura de plástico gera calor, o que altera a espessura do plástico. Adicionar efeitos de temperatura e comportamentos de plástico "elástico" mais complexos são os próximos passos para pesquisas futuras.

Em resumo: Eles criaram uma maneira mais rápida e estável de observar a mistura de plástico grosso com ar em máquinas giratórias em um computador, transformando um processo que antes levava horas em um que leva minutos, sem precisar reconstruir o mundo digital cada vez que um parafuso gira.

Resumo técnico

Resumo Técnico: Um Método de Fronteira Imersa com Volume de Fluido para Mistura Industrial de Polímeros

Enunciado do Problema
A simulação numérica de processos de mistura de polímeros apresenta desafios significativos devido à interação complexa entre reologia não newtoniana, dinâmica de superfície livre multifásica e geometrias complexas em movimento. Misturadores industriais, como extrusoras de rosca simples (SSE) e extrusoras de rosca dupla (TSE), frequentemente operam sob condições de preenchimento parcial, onde fusos de polímero altamente viscosos coexistem com ar. Simular esses cenários exige capturar a interface entre fases ao mesmo tempo em que se lidam com fronteiras sólidas rotativas.

Abordagens padrão de Dinâmica dos Fluidos Computacional (CFD) enfrentam dois obstáculos principais neste contexto:

1. Complexidade Geométrica: Gerar malhas ajustadas ao corpo (conformes) para parafusos rotativos com geometrias complexas é computacionalmente caro e frequentemente requer técnicas de malha dinâmica que podem levar a distorção da malha.
2. Instabilidade Numérica: O contraste extremo de viscosidade entre fusos de polímero e ar (frequentemente excedendo $10^5$ a $10^8$) causa instabilidades numéricas severas em solucionadores padrão segregados de Volume de Fluido (VOF). Especificamente, o tratamento explícito do gradiente de velocidade transposto no termo de difusão viscosa impõe restrições severas ao passo de tempo ($\Delta t \propto h^2$), tornando simulações de fluxos de alta viscosidade computacionalmente proibitivas ou instáveis.

Metodologia
Este trabalho propõe um novo framework numérico integrando uma abordagem de captura de interface por Volume de Fluido (VOF) com um método de Fronteira Imersa (IB) não conformante, implementado dentro da biblioteca de Volumes Finitos do OpenFOAM.

• Abordagem de Fronteira Imersa (IB): O método utiliza uma malha de fundo fixa onde as fronteiras sólidas (por exemplo, parafusos rotativos) são representadas por superfícies STL trianguladas. Um operador de interpolação por Mínimos Quadrados Ponderados (WLS) é empregado para impor condições de contorno (Dirichlet para velocidade, Neumann para fração volumétrica) em células intersectadas pela fronteira imersa. Isso elimina a necessidade de geração de malha conformante e permite o tratamento de linhas de contato móveis ao longo de geometrias complexas e rotativas.
• Formulação VOF: A interface é rastreada resolvendo uma equação de transporte para a fração volumétrica da fase ($\alpha$). O método emprega o algoritmo MULES (Limitador Universal Multidimensional com Solução Explícita) para manter a nitidez e limitação da interface, incluindo um termo de velocidade compressiva para contrapor a difusão numérica.
• Esquema Acoplado por Blocos (BC): Para abordar os problemas de estabilidade decorrentes de altos contrastes de viscosidade, os autores introduzem um esquema acoplado por blocos para a equação de momento. Diferentemente dos solucionadores segregados padrão que tratam componentes de velocidade sequencialmente, a abordagem BC trata o termo completo de difusão viscosa $\nabla \cdot (\mu(\nabla \mathbf{u} + \nabla^\top \mathbf{u}))$ implicitamente. Isso envolve o acoplamento das três componentes de velocidade dentro de um único sistema linear. A discretização do termo de gradiente tangencial utiliza um método de Área Finita (FA) combinado com interpolação WLS em estênceis estendidos (envolvendo vizinhos de ponto em vez de apenas vizinhos de face), adaptado de aplicações de elasticidade [11].
• Algoritmo de Solução: O solucionador emprega um algoritmo PIMPLE modificado (combinando SIMPLE e PISO) para lidar com o acoplamento pressão-velocidade. As restrições IB são integradas ao sistema algébrico, modificando os coeficientes da matriz e os termos fonte para satisfazer as condições de contorno nas superfícies imersas, mantendo a conservação de massa.

Principais Contribuições

1. Desenvolvimento de Solucionadores BC-VOF-IB: O artigo introduz o solucionador UCInterIbFoam, que combina o método IB não conformante com a abordagem VOF acoplada por blocos. Esta é a primeira implementação de tal método acoplado dentro do framework OpenFOAM-10 para fluxos bifásicos.
2. Melhoria de Estabilidade via Acoplamento por Blocos: O trabalho demonstra que o tratamento totalmente implícito da difusão viscosa relaxa significativamente as restrições de estabilidade do passo de tempo. Isso permite a simulação de fluxos de polímero de alta viscosidade com razões de viscosidade de até $10^8$ sem que a simulação falhe devido a violações do número de Courant.
3. Tratamento de Linhas de Contato Móveis: O framework aborda com sucesso as dificuldades geométricas e algorítmicas das linhas de contato móveis onde a interface bifásica intersecta fronteiras imersas, um requisito crítico para dispositivos de mistura parcialmente preenchidos.

Resultados Numéricos
Os métodos propostos foram validados e testados em duas categorias de problemas:

• Caso de Referência (Injeção de Molde em Forma de Ossos de Cão):

  • Uma comparação entre solucionadores segregados (interFoam), conformes acoplados por blocos (UCInterFoam) e não conformes acoplados por blocos (UCInterIbFoam) foi realizada em um caso de injeção de fluido de lei de potência de alta viscosidade.
  • O solucionador segregado padrão falhou em produzir resultados estáveis dentro de prazos razoáveis devido a restrições severas de passo de tempo ($\Delta t \approx 10^{-6}$ s).
  • Os solucionadores acoplados por blocos alcançaram soluções estáveis com passos de tempo de até $10^{-3}$ s, reduzindo o tempo computacional de 7 horas para 16 minutos (conforme) e 58 minutos (não conforme).
  • A solução IB não conformante mostrou concordância qualitativa com a solução conformante, com discrepâncias menores nos picos de viscosidade perto de restrições geométricas, atribuídas à natureza não conformante da malha.

• Aplicações Industriais (SSE e TSE):

  • Simulações foram conduzidas em Extrusoras de Rosca Simples (SSE) simplificadas e Extrusoras de Rosca Dupla (TSE) realistas sob condições tanto totalmente preenchidas quanto parcialmente preenchidas (alimentação por escassez).
  • Os solucionadores capturaram com sucesso a evolução da superfície livre, taxas de preenchimento e campos de pressão.
  • Os resultados alinharam-se com expectativas físicas: cenários totalmente preenchidos ($Q_{in} > Q_0$) exibiram gradientes de pressão negativos da entrada à saída, enquanto cenários parcialmente preenchidos ($Q_{in} < Q_0$) mostraram aumentos de pressão perto da saída.
  • O método IB provou ser eficaz no manuseio da geometria complexa de entrelaçamento da TSE sem exigir regeneração dinâmica da malha.

Significado e Alegações
O artigo alega que o framework BC-VOF-IB proposto representa um passo significativo em direção a fechar a lacuna entre a pesquisa acadêmica em CFD e as demandas industriais de processamento de polímeros. Ao superar as instabilidades numéricas associadas a altos contrastes de viscosidade e eliminar a necessidade de geração complexa de malhas ajustadas ao corpo para geometrias rotativas, o método permite previsões fisicamente consistentes dos campos de velocidade e pressão em dispositivos de mistura parcialmente preenchidos relevantes industrialmente.

Os autores reconhecem que o trabalho atual é limitado a processos isotérmicos e reologia newtoniana generalizada. Eles afirmam que desenvolvimentos futuros devem incluir efeitos térmicos (equação de energia) e modelos constitutivos viscoelásticos para capturar totalmente a física da mistura industrial de polímeros. No entanto, o framework atual estabelece uma base rigorosa para essas extensões e fornece uma ferramenta robusta para simular fluxos de superfície livre em geometrias complexas e em movimento.