A Volume of Fluid Immersed Boundary Method for Industrial Polymer Mixing

Questo articolo presenta un nuovo solver accoppiato a blocchi Volume of Fluid con Frontiera Immersa (BC-VOF-IB) implementato in OpenFOAM che supera le instabilità numeriche causate da forti contrasti di viscosità per simulare con precisione la miscelazione di polimeri con superficie libera in estrusori industriali parzialmente riempiti.

L'essenziale

Immagina di dover mescolare una gigantesca vasca di miele denso e appiccicoso con l'aria all'interno di una macchina rotante. Questo è essenzialmente ciò che accade nella miscelazione industriale dei polimeri, dove aziende come la Pirelli devono fondere la plastica con additivi per produrre pneumatici, dispositivi medici o componenti automobilistici. L'obiettivo è ottenere una miscelazione perfetta affinché il prodotto finale sia resistente e uniforme.

Tuttavia, simulare questo processo al computer è un incubo per matematici e ingegneri. Ecco il perché e come questo articolo risolve il problema, utilizzando semplici analogie:

Il Problema: La Lotta tra "Miele Denso e Aria Sottile"

In queste macchine, si hanno due fluidi molto diversi:

1. Fuso Polimerico: Estremamente denso, appiccicoso e lento (come il miele freddo).
2. Aria: Molto sottile e veloce.

Quando si tenta di simulare come questi due interagiscono all'interno di una macchina con viti rotanti, i programmi informatici standard si confondono. È come cercare di calcolare il movimento di una lumaca e di un'auto da corsa sulla stessa pista utilizzando lo stesso insieme di regole. Il computer cerca di compiere passi minuscoli, minuscoli per impedire che la "lumaca" (la plastica densa) si muova troppo velocemente, il che rende la simulazione incredibilmente lenta: a volte impiega giorni per completare pochi secondi di miscelazione in tempo reale.

Inoltre, le macchine hanno parti complesse e rotanti (viti) che si muovono all'interno di un contenitore fisso. Tradizionalmente, per simulare questo, è necessario costruire una mesh digitale (una griglia di piccoli cubi) che avvolga perfettamente le viti rotanti. Mentre le viti ruotano, questa griglia deve costantemente rimodellarsi, il che è come cercare di lavorare a maglia un maglione mentre la persona che lo indossa sta correndo una maratona. È disordinato, difficile e soggetto a errori.

La Soluzione: Una Nuova "Griglia Intelligente" e un "Approccio di Squadra"

Gli autori di questo articolo hanno sviluppato un nuovo modo per eseguire queste simulazioni utilizzando un software chiamato OpenFOAM. Hanno combinato due tecniche potenti:

1. Il Metodo della Frontiera Immersa (Il Trucco del "Muro Fantasma")
Invece di rimodellare la griglia per adattarla alle viti rotanti, hanno mantenuto la griglia fissa e rigida (come un blocco solido di ghiaccio). Hanno poi detto al computer: "Ehi, c'è una vite rotante dentro questo blocco di ghiaccio."

• L'Analogia: Immagina una piscina con una griglia fissa di piastrelle sul fondo. Invece di spostare le piastrelle per adattarle a un nuotatore, dici semplicemente all'acqua: "Non passare attraverso il nuotatore." Il computer utilizza la matematica per creare un "muro fantasma" attorno alla vite, costringendo il fluido a scorrere attorno ad essa senza dover mai ricostruire la griglia. Questo rende molto più semplice gestire forme complesse e in movimento.

2. Il Metodo del Volume di Fluido (VOF) (Il Trucco della "Pittura Tracciante")
Per vedere dove finisce la plastica densa e inizia l'aria, utilizzano una "pittura" che riempie le celle.

• L'Analogia: Immagina che la griglia del computer sia una scacchiera 3D. Alcuni quadrati sono al 100% plastica, altri al 100% aria e altri sono una miscela. Il computer traccia quanto "colore di plastica" c'è in ogni quadrato per vedere la superficie del liquido.

3. Lo Schema Accoppiato a Blocchi (L'"Huddle di Squadra")
Questa è la scoperta più importante. Nelle simulazioni standard, il computer risolve la velocità del fluido nelle direzioni X, Y e Z una alla volta, come tre persone che parlano a turno. Quando il fluido è super denso (come il polimero), questo approccio di "prendere il turno" fa crashare la simulazione o la rallenta fino a un'andatura di lumaca perché il fluido denso accoppia strettamente tutte le direzioni tra loro.

Gli autori hanno cambiato questo con un approccio Accoppiato a Blocchi.

• L'Analogia: Invece di tre persone che prendono il turno, si raggruppano tutte e risolvono il problema insieme esattamente nello stesso momento. Trattando il movimento in tutte le direzioni come un'unica grande squadra interconnessa, il computer può gestire la massiccia differenza tra la plastica densa e l'aria sottile senza bloccarsi.

I Risultati: Dalle Ore ai Minuti

Il team ha testato il loro nuovo metodo su due scenari:

1. Un Canale a Forma di Osso da Cane: Un caso di test in cui la plastica densa viene iniettata in un canale stretto e tortuoso.

   • Vecchio Metodo: Il programma informatico standard andava in crash o impiegava 7 ore per simulare pochi secondi perché era costretto a compiere passi minuscoli.
   • Nuovo Metodo: Il loro nuovo metodo "Huddle di Squadra" ha completato lo stesso lavoro in soli 16 minuti senza andare in crash, anche quando la plastica diventava estremamente densa.

2. Macchine Industriali Reali: Hanno simulato estrusori reali a vite singola e a doppia vite (le macchine utilizzate per produrre i pellet di plastica).

   • Hanno mostrato con successo come la plastica riempie la macchina, come la pressione si accumula e come l'aria viene spinta fuori.
   • Hanno dimostrato che il loro metodo "Muro Fantasma" funziona tanto bene quanto il vecchio e difficile metodo di rimodellare la griglia, ma è molto più veloce e facile da configurare.

Cosa Succede Ora?

L'articolo conclude che questo è un grande passo avanti per l'industria. Colma il divario tra la matematica accademica e le esigenze reali delle fabbriche. Tuttavia, gli autori notano che il loro modello attuale assume che la temperatura rimanga costante (isoterma). Nella realtà, la miscelazione della plastica genera calore, il che cambia la densità della plastica. Aggiungere gli effetti della temperatura e comportamenti della plastica più complessi e "elastici" sono i prossimi passi per la ricerca futura.

In breve: Hanno creato un modo più veloce e stabile per osservare al computer come la plastica densa si mescola con l'aria nelle macchine rotanti, trasformando un processo che richiedeva ore in uno che richiede minuti, senza dover ricostruire il mondo digitale ogni volta che una vite ruota.

Sommario tecnico

Sintesi Tecnica: Un Metodo Immersed Boundary Volume of Fluid per la Miscelazione Industriale dei Polimeri

Enunciato del Problema
La simulazione numerica dei processi di miscelazione dei polimeri presenta sfide significative a causa dell'interazione complessa tra reologia non newtoniana, dinamica multifase con superficie libera e geometrie complesse in movimento. I miscelatori industriali, come le estrusori a vite singola (SSE) e a vite doppia (TSE), operano spesso in condizioni di riempimento parziale, dove fusi polimerici ad alta viscosità coesistono con l'aria. La simulazione di tali scenari richiede la cattura dell'interfaccia tra le fasi gestendo al contempo confini solidi rotanti.

Gli approcci standard di Fluidodinamica Computazionale (CFD) affrontano due ostacoli principali in questo contesto:

1. Complessità Geometrica: La generazione di mesh adattate al corpo (conformi) per viti rotanti con geometrie complesse è computazionalmente costosa e richiede spesso tecniche di mesh dinamica che possono portare a distorsione della mesh.
2. Instabilità Numerica: L'estremo contrasto di viscosità tra fusi polimerici e aria (spesso superiore a $10^5$-$10^8$) causa gravi instabilità numeriche nei solutori VOF (Volume of Fluid) segregati standard. Nello specifico, il trattamento esplicito del gradiente di velocità trasposto nel termine di diffusione viscosa impone vincoli severi sul passo temporale ($\Delta t \propto h^2$), rendendo le simulazioni di flussi ad alta viscosità computazionalmente proibitive o instabili.

Metodologia
Questo lavoro propone un nuovo quadro numerico che integra un approccio di cattura dell'interfaccia Volume of Fluid (VOF) con un metodo Immersed Boundary (IB) non conforme, implementato all'interno della libreria a Volumi Finiti OpenFOAM.

• Approccio Immersed Boundary (IB): Il metodo utilizza una griglia di sfondo fissa in cui i confini solidi (ad esempio, le viti rotanti) sono rappresentati da superfici STL triangolate. Viene impiegato un operatore di interpolazione ai Minimi Quadrati Ponderati (WLS) per imporre le condizioni al contorno (Dirichlet per la velocità, Neumann per la frazione di volume) sulle celle intersecate dal confine immerso. Questo elimina la necessità di generazione di mesh conformi e permette la gestione di linee di contatto in movimento lungo geometrie complesse e rotanti.
• Formulazione VOF: L'interfaccia è tracciata risolvendo un'equazione di trasporto per la frazione di volume della fase ($\alpha$). Il metodo impiega l'algoritmo MULES (Multidimensional Universal Limiter with Explicit Solution) per mantenere la nitidezza e la limitatezza dell'interfaccia, includendo un termine di velocità compressiva per contrastare la diffusione numerica.
• Schema Accoppiato a Blocchi (BC): Per affrontare i problemi di stabilità derivanti da alti contrasti di viscosità, gli autori introducono uno schema accoppiato a blocchi per l'equazione della quantità di moto. A differenza dei solutori segregati standard che trattano le componenti di velocità sequenzialmente, l'approccio BC tratta implicitamente l'intero termine di diffusione viscosa $\nabla \cdot (\mu(\nabla \mathbf{u} + \nabla^\top \mathbf{u}))$. Ciò comporta l'accoppiamento delle tre componenti di velocità all'interno di un singolo sistema lineare. La discretizzazione del termine di gradiente tangenziale utilizza un metodo di Area Finita (FA) combinato con interpolazione WLS su stencil estesi (che coinvolgono punti vicini piuttosto che solo facce vicine), adattato da applicazioni di elasticità [11].
• Algoritmo di Soluzione: Il solutore impiega un algoritmo PIMPLE modificato (che combina SIMPLE e PISO) per gestire l'accoppiamento pressione-velocità. I vincoli IB sono integrati nel sistema algebrico, modificando i coefficienti della matrice e i termini sorgente per soddisfare le condizioni al contorno sulle superfici immerse mantenendo la conservazione della massa.

Contributi Chiave

1. Sviluppo di Solutori BC-VOF-IB: Il documento introduce il solutore UCInterIbFoam, che combina il metodo IB non conforme con l'approccio VOF accoppiato a blocchi. Questa è la prima implementazione di un tale metodo accoppiato all'interno del framework OpenFOAM-10 per flussi bifase.
2. Miglioramento della Stabilità tramite Accoppiamento a Blocchi: Il lavoro dimostra che il trattamento totalmente implicito della diffusione viscosa rilassa significativamente i vincoli di stabilità del passo temporale. Ciò permette la simulazione di flussi polimerici ad alta viscosità con rapporti di viscosità fino a $10^8$ senza che la simulazione si arresti a causa di violazioni del numero di Courant.
3. Gestione delle Linee di Contatto in Movimento: Il quadro affronta con successo le difficoltà geometriche e algoritmiche delle linee di contatto in movimento dove l'interfaccia bifase interseca i confini immersi, un requisito critico per dispositivi di miscelazione parzialmente riempiti.

Risultati Numerici
I metodi proposti sono stati validati e testati su due categorie di problemi:

• Caso di Riferimento (Stampaggio a Iniezione a Forma di Dog-bone):

  • È stata eseguita una comparazione tra solutori segregati (interFoam), accoppiati a blocchi conformi (UCInterFoam) e accoppiati a blocchi non conformi (UCInterIbFoam) su un caso di iniezione di fluido con legge di potenza ad alta viscosità.
  • Il solutore segregato standard non è riuscito a produrre risultati stabili entro tempi ragionevoli a causa di severe restrizioni sul passo temporale ($\Delta t \approx 10^{-6}$ s).
  • I solutori accoppiati a blocchi hanno raggiunto soluzioni stabili con passi temporali fino a $10^{-3}$ s, riducendo il tempo computazionale da 7 ore a 16 minuti (conforme) e 58 minuti (non conforme).
  • La soluzione IB non conforme ha mostrato accordo qualitativo con la soluzione conforme, con lievi discrepanze nei picchi di viscosità vicino a restrizioni geometriche, attribuite alla natura non conforme della griglia.

• Applicazioni Industriali (SSE e TSE):

  • Sono state condotte simulazioni su Estrusori a Vite Singola (SSE) semplificati ed Estrusori a Vite Doppia (TSE) realistici in condizioni sia completamente riempite che parzialmente riempite (alimentazione a fame).
  • I solutori hanno catturato con successo l'evoluzione della superficie libera, i rapporti di riempimento e i campi di pressione.
  • I risultati erano allineati con le aspettative fisiche: gli scenari completamente riempiti ($Q_{in} > Q_0$) hanno mostrato gradienti di pressione negativi dall'ingresso all'uscita, mentre gli scenari parzialmente riempiti ($Q_{in} < Q_0$) hanno mostrato aumenti di pressione vicino all'uscita.
  • Il metodo IB si è dimostrato efficace nella gestione della geometria complessa di ingranamento della TSE senza richiedere rigenerazione dinamica della mesh.

Significato e Affermazioni
Il documento afferma che il quadro BC-VOF-IB proposto rappresenta un passo significativo verso il colmare il divario tra la ricerca accademica in CFD e le esigenze di lavorazione industriale dei polimeri. Superando le instabilità numeriche associate agli alti contrasti di viscosità ed eliminando la necessità di una complessa generazione di mesh adattate al corpo per geometrie rotanti, il metodo permette previsioni fisicamente coerenti dei campi di velocità e pressione in dispositivi di miscelazione parzialmente riempiti di rilevanza industriale.

Gli autori riconoscono che il lavoro attuale è limitato a processi isotermi e reologia newtoniana generalizzata. Affermano che gli sviluppi futuri devono includere effetti termici (equazione dell'energia) e modelli costitutivi viscoelastici per catturare appieno la fisica della miscelazione industriale dei polimeri. Tuttavia, il quadro attuale stabilisce una base rigorosa per queste estensioni e fornisce uno strumento robusto per la simulazione di flussi con superficie libera in geometrie complesse e in movimento.