Energy non-equipartition in vibrofluidized particles

Sintesi Tecnica: Non-equipartizione dell'energia in particelle vibrofluidizzate

Enunciato del Problema
Il teorema dell'equipartizione, che postula che l'energia cinetica sia distribuita equamente tra tutti i gradi di libertà in un fluido, devia frequentemente nei sistemi granulari. Sebbene studi precedenti abbiano documentato la non-equipartizione in letti vibrofluidizzati contenenti sia particelle lisce (portando all'anisotropia nell'energia cinetica fluttuante) sia particelle ruvide, l'influenza specifica dei parametri di interazione realistici delle particelle rimane oggetto di indagine. Nello specifico, questo lavoro affronta come il rapporto tra il coefficiente di rigidità della molla tangenziale e quello normale ($\kappa = k_t/k_n$) influenzi la ripartizione dell'energia tra i modi traslazionali e rotazionali. Lo studio mira a determinare se modelli realistici di contatto tra particelle, in particolare quelli che coinvolgono un alto attrito e rapporti di rigidità specifici, portino all'equipartizione dell'energia o a una non-equipartizione sostenuta.

Metodologia
Gli autori hanno impiegato il Metodo degli Elementi Discreti (DEM) utilizzando il software open-source LAMMPS per simulare un sistema granulare tridimensionale, verticalmente vibrato, composto da particelle sferiche.

• Modello di Interazione: Le forze di contatto sono state modellate utilizzando un modello molla-smorzatore lineare (LSD). La costante di rigidità della molla normale ($k_n = 10^8 mg/d$) è stata selezionata sufficientemente grande per garantire che le collisioni siano prevalentemente binarie ($t_c/t_f \ll 1$).
• Parametri: Lo studio ha variato il coefficiente di attrito ($\mu$) da quasi liscio ($\mu \approx 0$) a quasi perfettamente ruvido ($\mu = 10$). Sono stati investigati due valori principali per il rapporto di rigidità $\kappa$: $2/7$ e $3/4$. Ulteriori simulazioni sono state condotte per $0.67 \le \kappa < 1$.
• Protocollo di Simulazione: Le simulazioni sono state eseguite in due fasi: una fase iniziale di rilassamento ($7 \times 10^7 t_c$) per raggiungere uno stato stazionario, seguita dalla raccolta dati con un passo temporale più piccolo ($\Delta t = t_c/100$).
• Analisi: Le velocità lineari e angolari istantanee sono state analizzate per determinare i profili delle energie cinetiche medie fluttuanti traslazionali ($T$) e rotazionali ($T_R$). Le medie d'insieme sono state calcolate su più di $10^3$ configurazioni. Il rapporto $K = K_{ET}/K_{ER}$ (energia cinetica traslazionale su energia cinetica rotazionale) è stato utilizzato come metrica primaria per l'equipartizione.

Risultati Chiave

1. Dipendenza dal Rapporto di Rigidità ($\kappa$): Il comportamento del rapporto energetico $K$ dipende criticamente da $\kappa$.

   • Per $\kappa = 2/7$: Il rapporto $K$ diminuisce monotonicamente all'aumentare del coefficiente di attrito $\mu$. Ad alto attrito ($\mu = 10$), $K$ si avvicina all'unità, indicando che l'equipartizione dell'energia tra i modi traslazionali e rotazionali è raggiunta.
   • Per $\kappa = 3/4$: Il comportamento è non monotono. Dopo una riduzione iniziale con l'aumento di $\mu$, $K$ aumenta e si stabilizza a circa 5. Ciò indica una persistente mancanza di equipartizione dell'energia, anche a coefficienti di attrito molto elevati.
   • Le simulazioni per $0.67 \le \kappa < 1$ hanno confermato che la non-equipartizione persiste per particelle con alti coefficienti di attrito in questo intervallo.

2. Dinamica delle Collisioni e Scambio di Energia:

   • L'analisi delle collisioni individuali ha rivelato che per $\kappa = 2/7$, la frazione di contatti nel regime di "aderenza" (dove il coefficiente di restituzione rotazionale $\beta \approx 1$) aumenta con $\mu$. A $\mu = 1$, quasi il 75% dei contatti si trova in questo regime conservativo di energia, facilitando l'equipartizione.
   • Per $\kappa = 3/4$, la frazione di contatti nel regime stick-slip si stabilizza a circa il 55% per $\mu > 1$, con $\beta < 1$. Di conseguenza, il rapporto di scambio energetico $\Theta = E_{exch}/\Delta E$ rimane inferiore a 1 e diminuisce per particelle molto ruvide. La dissipazione supera lo scambio di energia tra i modi, impedendo l'equipartizione.

3. Validazione della Teoria: I risultati delle simulazioni per la temperatura granulare di equilibrio ($T_o$) hanno mostrato un eccellente accordo con le previsioni teoriche (basate su espansioni di momenti di ordine principale) sia per particelle perfettamente lisce sia per particelle ruvide in regimi diluiti, a condizione che venissero utilizzati i parametri di scala appropriati (coinvolgenti $\beta$).

Significato e Affermazioni
Il lavoro afferma che la selezione del rapporto di rigidità $\kappa$ è un fattore cruciale nella previsione del comportamento macroscopico del flusso di particelle granulari realistiche. Lo studio dimostra che l'equipartizione dell'energia non è un risultato garantito nei sistemi granulari ruvidi vibrofluidizzati; piuttosto, è condizionata dalla specifica interazione tra attrito e il rapporto di rigidità tangenziale-normale.

Nello specifico, gli autori concludono che:

• L'equipartizione è osservata per $\kappa = 2/7$ ad alto attrito a causa del dominio delle collisioni di aderenza conservatrici di energia.
• La non-equipartizione è osservata per $0.67 \le \kappa < 1$ (incluso $\kappa = 3/4$) perché le collisioni stick-slip in questi regimi non conservano l'energia, portando a una significativa disparità tra le energie cinetiche traslazionali e rotazionali.

Questi risultati suggeriscono che la modellazione realistica della reologia granulare richiede un'attenta considerazione del rapporto di rigidità, poiché altera fondamentalmente i meccanismi di distribuzione dell'energia all'interno del sistema. Il lavoro non propone nuovi setup sperimentali o applicazioni, ma fornisce un'analisi sistematica di come i parametri di interazione microscopica dettino la distribuzione macroscopica dell'energia.