Beyond binary scission: a generalized three-species cascade breakage model for wormlike micellar solutions

Sintesi Tecnica: Un modello di rottura a cascata a tre specie generalizzato per soluzioni di micelle similiformi (wormlike micelles)

Definizione del Problema
Le soluzioni di micelle similiformi (WLM) esibiscono comportamenti reologici complessi, tra cui viscoelasticità, shear thinning e shear banding, guidati dalla continua rottura e ricombinazione di reti micellari auto-assemblate. Gli attuali framework teorici, in particolare i modelli a due specie ampiamente utilizzati (ad esempio, il modello Vasquez-Cook-McKinley o VCM), si basano su una rappresentazione binaria grossolana della popolazione micellare (tipicamente distinguendo solo tra catene lunghe aggrovigliate e frammenti corti). Sebbene siano efficaci nel catturare determinati tratti macroscopici, questi modelli possiedono un limite fondamentale: non riescono a risolvere esplicitamente la rottura a cascata multiscala e i conseguenti regimi di rilassamento continui inerenti alle reali reti micellari. Questa approssimazione binaria filtra gli stati strutturali intermedi, portando all'incapacità di descrivere accuratamente gli spettri di rilassamento ampi, in particolare l'incremento ad alta frequenza dei moduli dinamici osservato negli esperimenti, e limita l'interpretazione fisica delle risposte transitorie e dei meccanismi di shear banding.

Metodologia
Per affrontare tali limitazioni, gli autori sviluppano un modello di rottura a cascata a tre specie che introduce uno stato micellare intermedio nel framework cinetico. Il modello assume che la rete micellare consista in tre stati strutturali discreti:

1. Specie A: Una rete di gel (primaria, altamente aggrovigliata).
2. Specie B: Catene lunghe (stato intermedio).
3. Specie C: Catene corte (frammenti disaggrovigliati).

La formulazione matematica si basa su un approccio di teoria cinetica utilizzando distribuzioni nello spazio delle fasi continuo. Gli autori derivano equazioni di Smoluchowski accoppiate per le funzioni di densità di probabilità delle tre specie, incorporando termini di convezione, diffusione e reazione. Le caratteristiche metodologiche chiave includono:

• Cinetica Gerarchica: Il modello definisce un percorso a cascata in cui la rete di gel (A) si rompe in catene lunghe (B), che successivamente si rompono in catene corte (C). Sono definiti anche percorsi di ricombinazione per garantire la conservazione della massa.
• Scissione Dipendente dal Flusso: I tassi di rottura per la rete di gel e le catene lunghe sono accoppiati al flusso macroscopico locale tramite un meccanismo di convezione di filamenti parzialmente estensibile (PEC). Ciò crea un meccanismo di rottura dipendente dallo stress, in cui la deformazione accelera la scissione.
• Equazioni Costitutive: Il modello deriva equazioni di evoluzione adimensionali per le densità numeriche e i tensori di conformazione, accoppiate all'equazione del momento di Cauchy macroscopica. Lo stress totale è la sovrapposizione degli stress elastici di tutte e tre le specie.
• Implementazione Numerica: Gli autori analizzano il modello sia sotto condizioni di flusso omogeneo (per studiare la viscoelasticità lineare e la risposta di taglio stazionaria) sia in un flusso di Couette cilindrico disomogeneo (per investigare lo shear banding, l'avvio transitorio e il rilassamento).

Contributi Chiave e Risultati

1. Risposta Viscoelastica a Tre Modi:
   Sotto flusso omogeneo, il modello produce analiticamente una struttura di Maxwell a tre modi discreti. A differenza dei modelli a due specie, questo framework riesce a catturare con successo l'incremento ad alta frequenza dei moduli di accumulo ($G'$) e di perdita ($G''$). Gli autori dimostrano che il modo di rilassamento aggiuntivo fornito dalla specie intermedia (B) e dalle catene corte (C) migliora significativamente l'adattamento alle curve master sperimentali per soluzioni CPyCl/NaSal, in particolare nel regime ad alta frequenza.

2. Curve di Flusso Non Monotoniche e Shear Banding:
   Il modello predice una curva costitutiva non monotona con una regione a pendenza strettamente negativa, un segno distintivo di instabilità costitutiva. Nelle simulazioni di flusso di Couette disomogeneo, questa instabilità porta alla formazione spontanea di bande di taglio coesistenti (regioni ad alto e basso tasso di taglio). Il modello a tre specie rivela che la specie intermedia (B) forma uno strato di transizione distinto all'interfaccia tra le bande, collegando la banda a basso shear rate dominata dalla rete di gel e la banda ad alto shear rate dominata dalle catene corte. Ciò fornisce una descrizione microstrutturale dell'interfaccia di banda più risolta rispetto a quanto consentito dai modelli binari.

3. Dinamica Transitoria ed Isteresi:
   Il modello cattura comportamenti transitori complessi:

• Overshoot di Stress: Durante l'avvio (startup), la competizione tra il carico elastico e la rottura accelerata a cascata porta a pronunciati overshoot di stress.
• Isteresi di Flusso: Il modello esibisce isteresi, in cui lo stato stazionario finale dipende dalla storia del flusso (avvio da fermo vs. riduzione a passo di alto shear).
• Rilassamento Multistep: Al termine del taglio, il modello predice un processo di rilassamento multistep. Invece di un singolo decadimento esponenziale, lo stress si rilassa attraverso stadi sequenziali: rilascio rapido dalle catene corte (C), aggiustamento intermedio dalle catene lunghe (B) e recupero lento della rete di gel (A). Questo risolve il "rilassamento prolungato della scala temporale" spesso osservato negli esperimenti ma difficile da catturare con modelli ridotti.

4. Consistenza Meccanica del Continuo:
   Il modello soddisfa la relazione di Lodge-Meissner ($N_1/\sigma_{r\theta} = \gamma_0$) durante grandi deformazioni a gradino di deformazione, confermando la sua consistenza tensoriale e validità fisica sotto deformazione finita.

Significatività e Rivendicazioni
L'articolo sostiene che il proposto framework a tre specie colma il divario tra i modelli di rete fenomenologici discreti e le teorie cinetiche continue (come la teoria originale di reptazione-reazione di Cates). Risolvendo esplicitamente la cinetica a cascata che governa gli stati strutturali intermedi, il modello fornisce un framework fisicamente interpretabile per lo shear banding delle micelle similiformi.

Gli autori affermano che questo approccio cattura la connessione intrinseca tra:

• Evoluzione microstrutturale a cascata (A $\to$ B $\to$ C).
• Dinamica di rilassamento ampia (tre distinti tempi di scala).
• Localizzazione del flusso macroscopico (shear banding).

La significatività risiede nella capacità del modello di mantenere la trattabilità computazionale macroscopica pur recuperando i modi di rilassamento intermedi mancanti. Ciò offre una base fisica più completa per l'interpretazione della ricca dinamica transitoria e stazionaria dei fluidi micellari complessi, potenzialmente aiutando l'interpretazione di futuri esperimenti microfluidici e applicazioni industriali che coinvolgono le WLM. Il documento non propone nuovi protocolli sperimentali, ma offre uno strumento teorico per interpretare meglio i fenomeni reologici esistenti.