A Robust Truncated-Domain Approach for Cone--Jet Simulations in Electrospinning and Electrospraying

Questo studio presenta un nuovo approccio di dominio troncato per le simulazioni elettrofluidodinamiche del cono-getto nell'elettrofilatura e nell'elettrospruzzatura, che utilizza simulazioni elettrostatiche a dominio completo per definire condizioni al contorno precise, eliminando la necessità di parametri empirici e riducendo significativamente i costi computazionali mantenendo un'alta accuratezza predittiva.

L'essenziale

Immagina di dover creare fibre microscopiche o goccioline minuscole per creare farmaci, filtri o sensori. Per farlo, gli scienziati usano una tecnica chiamata elettrofilatura (o electrospinning). È un po' come se avessi un ago da siringa molto fine, lo riempi di un liquido appiccicoso e gli dai una scossa elettrica molto forte.

Ecco cosa succede: il liquido all'estremità dell'ago si deforma, formando una punta a forma di cono (chiamata "cono di Taylor") e da lì viene spruzzato un getto sottilissimo che si trasforma in fibre o gocce.

Il Problema: Il "Dilemma della Scala"
Il problema per i computer che simulano questo processo è enorme. Immagina di dover disegnare su un foglio:

1. La punta dell'ago è grande come un capello (pochi micron).
2. La distanza fino al "pavimento" dove le fibre atterrano è come la lunghezza di una stanza (20 centimetri).

C'è una differenza di scala di circa 100 volte. Per simulare tutto correttamente, il computer dovrebbe usare una griglia di punti così fitta da coprire l'intera stanza, ma con la precisione di un capello. Questo richiederebbe una potenza di calcolo mostruosa, come cercare di contare ogni singolo granello di sabbia in una spiaggia usando un microscopio. È troppo costoso e lento.

La Soluzione Vecchia: "Indovinare la Mappa"
Per risparmiare, gli scienziati hanno provato a simulare solo la parte vicina all'ago (il "cono") e a ignorare il resto della stanza. Ma per farlo, dovevano inventare delle regole matematiche (formule analitiche) per dire al computer come si comporta l'elettricità ai bordi di questo piccolo mondo tagliato.
Purtroppo, queste vecchie regole erano come una mappa disegnata a mano da un artista: non erano precise. Sottostimavano la forza dell'elettricità vicino alla punta e richiedevano di essere "aggiustate" (tarate) guardando prima i risultati di esperimenti reali o simulazioni complete. Se non avevi già fatto l'esperimento, non potevi usare la simulazione per prevedere cosa sarebbe successo. Era un circolo vizioso.

La Nuova Soluzione: "La Foto Precisa"
Gli autori di questo articolo (Ghanashyam K. C., Satyavrata Samavedi e Harish N. Dixit) hanno trovato un modo intelligente e robusto per risolvere il problema.

Ecco la loro idea, spiegata con una metafora:

Immagina che l'elettricità sia come il vento che soffia intorno a un albero (l'ago).

1. Il vecchio metodo: Cercava di calcolare il vento usando una formula generica ("il vento soffia sempre così"). Spesso sbagliava, specialmente vicino all'albero.
2. Il nuovo metodo: Prima di simulare il movimento del liquido, fanno una simulazione solo del vento (l'elettricità statica) su tutta la stanza. È un calcolo molto più semplice e veloce perché non devono preoccuparsi del liquido che si muove.
3. La mossa geniale: Una volta che hanno la mappa esatta del vento (il campo elettrico) vicino all'ago, la "fotocopiano" e la usano come regola per i bordi della loro piccola simulazione tagliata.

In pratica, dicono al computer: "Non preoccuparti di tutta la stanza. Simula solo il cono vicino all'ago, ma sappi che il vento ai bordi è esattamente come lo abbiamo misurato nella foto completa."

I Risultati: Perché è meglio?

• Nessun "aggiustamento": Non serve più indovinare parametri o guardare esperimenti precedenti. Il metodo funziona da solo.
• Precisione: Riproduce perfettamente la forma del cono, la velocità del liquido e la distribuzione della carica elettrica, esattamente come farebbe una simulazione completa (ma molto più veloce).
• Risparmio: Permette di usare domini di simulazione molto più piccoli, riducendo i tempi di calcolo da giorni a ore, senza perdere accuratezza.

In sintesi
Questo studio offre un nuovo modo per "tagliare" il problema in due: prima si calcola con precisione dove va l'elettricità (che è facile), e poi si usa quella informazione per simulare il movimento del liquido in una zona piccola. È come se, invece di dover calcolare il traffico di tutta una città per sapere come muoversi in un incrocio, guardaste prima una mappa satellitare aggiornata del traffico e poi simulaste solo il vostro incrocio sapendo esattamente cosa succede intorno.

È un approccio più intelligente, più veloce e più affidabile per progettare le tecnologie del futuro, dalle medicine ai nuovi materiali.

Sommario tecnico

Titolo: Un approccio robusto a dominio troncato per le simulazioni cono-getto nell'elettrofilatura e nell'elettrospruzzatura

1. Il Problema

Le simulazioni numeriche dirette (DNS) dei processi di elettrofilatura e elettrospruzzatura sono estremamente onerose dal punto di vista computazionale a causa della vasta separazione di scale spaziali coinvolta.

• Disparità di scale: La distanza tra l'ago e il collettore ($H'$) è tipicamente dell'ordine di 15-20 cm, mentre la regione del "cono-getto" (dove avviene la formazione del getto) si estende solo per 0,6-1 cm. Inoltre, il diametro interno dell'ago richiede risoluzioni di griglia dell'ordine di micron (1-5 µm).
• Costo computazionale: Simulare l'intero dominio fisico con la risoluzione necessaria per catturare le interfacce liquide e i dettagli del campo elettrico vicino all'ago rende le simulazioni tridimensionali proibitive per configurazioni sperimentali realistiche.
• Limiti degli approcci attuali: Le simulazioni a dominio troncato (TDS) sono un'alternativa valida, ma gli approcci esistenti basati su espressioni analitiche per il potenziale elettrico (come il modello di Jones e Thong) presentano due difetti critici:
  1. Sottostimano sistematicamente il campo elettrico vicino alla punta dell'ago, riducendo le forze elettrostatiche calcolate.
  2. Richiedono una regolazione empirica dei parametri basata su conoscenze a priori della configurazione cono-getto (ottenuta da esperimenti o simulazioni a dominio completo), limitando così la capacità predittiva del metodo per nuove configurazioni.

2. Metodologia

Gli autori propongono un nuovo quadro generale per le simulazioni elettro-idrodinamiche (EHD) a dominio troncato che supera le limitazioni dei metodi analitici. La metodologia si articola nei seguenti passaggi:

1. Simulazione Elettrostatica a Dominio Completo (a basso costo): Viene eseguita una simulazione puramente elettrostatica su un dominio completo (che include l'ago e il collettore). Poiché non coinvolge la dinamica dei fluidi, questo passaggio è computazionalmente economico.
2. Estrazione dei Dati: Da questa simulazione elettrostatica, vengono estratti i valori esatti del campo elettrico e del potenziale elettrico lungo i confini superiori e laterali del futuro dominio troncato. Sulla parte inferiore (vicino all'ago), viene registrata solo la componente normale del gradiente del campo elettrico.
3. Approssimazione e Adattamento: I dati estratti vengono lisciati e adattati a profili gaussiani a tre termini per garantire continuità e differenziabilità.
4. Simulazione EHD a Dominio Troncato: Questi profili adattati vengono impostati come condizioni al contorno per una simulazione EHD completa (che include la dinamica dei fluidi e l'interfaccia liquido-gas) eseguita esclusivamente sulla regione locale del cono-getto (dominio troncato).
5. Strumenti Numerici: Le simulazioni sono state condotte utilizzando OpenFOAM, estendendo il solver interFoam per incorporare il modello dielettrico perdente di Taylor-Melcher. L'interfaccia è catturata tramite il metodo VOF (Volume of Fluid) con il limite MULES e raffinamento adattivo della mesh.

3. Contributi Chiave

• Eliminazione dei parametri regolabili: Il metodo proposto non richiede la calibrazione empirica di coefficienti (come il coefficiente $K_V$ nel modello di Jones-Thong) basata su risultati sperimentali o simulazioni complete precedenti.
• Accuratezza del campo elettrico: Utilizzando i dati esatti della simulazione elettrostatica completa, il metodo garantisce la corretta distribuzione del campo elettrico vicino alla punta dell'ago, indipendentemente dalle dimensioni del dominio troncato.
• Indipendenza dalla configurazione: Il metodo è predittivo e non richiede la conoscenza a priori della forma del cono-getto, rendendolo applicabile a regimi di parametri inesplorati.
• Generalità: Sebbene dimostrato in configurazione assialsimmetrica, il metodo è intrinsecamente estendibile a geometrie tridimensionali complete e a fluidi non newtoniani o viscoelastici, poiché la procedura di ricostruzione elettrostatica è indipendente dal modello reologico del fluido.

4. Risultati

I risultati sono stati confrontati con simulazioni a dominio completo (FDS) e dati sperimentali (utilizzando 1-ottanolo come fluido di lavoro):

• Convergenza a domini più piccoli: Il metodo proposto raggiunge la convergenza della forma del cono-getto con dimensioni di dominio troncato significativamente più piccole (circa $6R_i \times 6R_i$) rispetto ai metodi analitici, che richiedono domini più grandi per stabilizzarsi.
• Accuratezza delle quantità fisiche: Oltre alla forma dell'interfaccia, il metodo riproduce con alta precisione quantità fisiche critiche difficili da misurare sperimentalmente, tra cui:
  • Distribuzioni di densità di carica.
  • Campi di velocità assiale.
  • Tensioni di Maxwell.
  • Correnti di conduzione e convezione superficiale.
• Confronto con il metodo Jones-Thong: Mentre il metodo di Jones-Thong può produrre forme di interfaccia simili se i parametri vengono "aggiustati" manualmente, fallisce nel riprodurre accuratamente i gradienti di campo e le distribuzioni di corrente senza tale calibrazione. Il metodo proposto, invece, ottiene un accordo quantitativo superiore con le simulazioni a dominio completo senza alcun aggiustamento.
• Riduzione dei costi: La strategia riduce drasticamente il costo computazionale mantenendo l'accuratezza fisica, eliminando la necessità di simulare l'intera distanza tra ago e collettore nella dinamica dei fluidi.

5. Significato e Implicazioni

Questo lavoro fornisce un quadro robusto e predittivo per lo studio dei flussi elettro-idrodinamici nel regime cono-getto.

• Ottimizzazione dei processi: Permette di esplorare parametri operativi (tensione, portata, proprietà del fluido) e ottimizzare i processi di elettrofilatura e elettrospruzzatura senza i costi proibitivi delle simulazioni a dominio completo.
• Affidabilità: Rimuove l'incertezza legata alla scelta empirica delle condizioni al contorno nei domini troncati, rendendo le simulazioni uno strumento affidabile per la progettazione di sistemi di produzione di fibre e particelle sub-micrometriche.
• Flessibilità: La metodologia è versatile e può essere applicata a configurazioni complesse (3D, flussi instabili, fluidi complessi), aprendo la strada a studi più approfonditi sulle instabilità del getto e sulla dinamica delle gocce.

In sintesi, l'approccio proposto sostituisce le approssimazioni analitiche con dati elettrostatici esatti e a basso costo, risolvendo il compromesso tra accuratezza fisica e efficienza computazionale nelle simulazioni di elettrospinning ed electrospraying.