Batterie al Sodio-Ione vs LFP nello Stoccaggio Energetico: Comportamento Termico e Curve di Scarica

Introduzione e Fondamenti Termodinamici ed Elettrochimici delle Tecnologie Na-ion e LiFePO₄

La transizione energetica globale richiede sistemi di accumulo dell'energia elettrochimica (BESS) che non solo massimizzino la densità energetica, ma garantiscano anche una resilienza termica superiore, una sicurezza intrinseca e una catena di approvvigionamento svincolata da materie prime critiche come il litio e il cobalto. In questo contesto, le batterie agli Ioni di Sodio (Na-ion) emergono come un'alternativa dirompente rispetto alle chimiche al Litio-Ferro-Fosfato (LiFePO₄ o LFP). Dal punto di vista dell'ingegneria elettrica di alta potenza e della progettazione di sistemi avanzati, è imperativo analizzare i fondamenti termodinamici e cinetici che governano entrambe le tecnologie per modellare correttamente il loro comportamento operativo, le loro curve di scarica e i loro profili di sicurezza in condizioni di stress severo.

Termodinamicamente, la differenza fondamentale risiede nel potenziale redox standard della coppia ionica e nel raggio atomico. Il sodio presenta un raggio ionico maggiore (1.02\, Å rispetto a 0.76\, Å del litio) e una massa molare superiore (22.99\, g/mol rispetto a 6.94\, g/mol ). Ciò si traduce in un potenziale standard di riduzione più elettronegativo per il litio (-3.04\, V vs. $SHE$ ) rispetto al sodio (-2.71\, V vs. $SHE$ ), il che stabilisce un limite superiore teorico per la tensione di cella nelle batterie Na-ion. Tuttavia, l'energia libera di Gibbs (\Delta G) delle reazioni di intercalazione definisce i profili di tensione operativa e la stabilità strutturale degli elettrodi anodici e catodici.

Nelle celle LiFePO₄, il catodo opera mediante un meccanismo di inserzione a due fasi (transizione di fase di primo ordine) tra $LiFePO _4$ e $FePO _4$ , generano un plateau di scarica straordinariamente piatto intorno a 3.2\, V - 3.3\, V . Al contrario, i sistemi Na-ion attuali utilizzano una gamma più ampia di materiali catodici, prevalentemente ossidi di metalli di transizione stratificati (Layered Transition Metal Oxides) e analoghi del blu di Prussia (Prussian Blue Analogues). Questi materiali mostrano frequentemente soluzioni solide monofasiche in ampi intervalli dello stato di carica (SoC), il che conferisce alle curve di scarica del Na-ion una pendenza caratteristica anziché un plateau piatto, complicando gli algoritmi di stima del SoC ma offrendo vantaggi termodinamici in termini di gestione del calore e prevenzione di sovratensioni localizzate.

L'energia interna del sistema e la sua variazione termica sono regolate dall'equazione fondamentale della calorimetria elettrochimica, che relaziona il calore generato con l'entropia della reazione:

\dot{Q} = I \left( U_{ oc } - V_{ term } \right) - I T \left( \frac{\partial U_{ oc }}{\nobreak\partial T} \right)

Dove \dot{Q} è il tasso di generazione di calore, $I$ è la corrente di funzionamento ( $I > 0$ in scarica), $U_{ oc }$ è la tensione a circuito aperto, $V_{ term }$ è la tensione ai terminali della cella, $T$ è la temperatura assoluta, e il termine \left(\frac{\partial U_{ oc }}{\partial T}\right) rappresenta il coefficiente entropico. Nelle celle LFP, il coefficiente entropico varia drasticamente vicino agli estremi di carica a causa delle transizioni di fase, provocando picchi locali di generazione di calore entropico. Nelle celle Na-ion, l'assenza di transizioni di fase repentine in determinati design di catodi stratificati attenua questo comportamento, riducendo la suscettibilità a hot spot interni indotti da gradienti entropici durante impulsi ad alta corrente.

Parametri Elettrici, Chimici e Normativi Comparativi

Al fine di stabilire una base di progettazione rigorosa per progetti di ingegneria elettrica su media e grande scala, è indispensabile confrontare i parametri intrinseci delle celle Na-ion e LFP ai sensi delle normative internazionali standardizzate (quali IEC 62619, IEEE 1679, UL 1973 e UL 9540A). La tabella seguente consolida le metriche critiche di prestazioni elettriche, i limiti normativi e le conseguenze operative.

Parametro Tecnico / Metrica	Tecnologia agli Ioni di Sodio (Na-ion)	Tecnologia al Litio-Ferro-Fosfato (LFP)	Normativa di Riferimento / Standard IEC/IEEE/UL	Conseguenza Operativa e Dielettrica
Tensione Nominale di Cella	3.0\, V - 3.3\, V	3.2\, V - 3.3\, V	IEC 62619 / IEEE 1679	Influenza la topologia degli inverter di potenza e il rapporto di conversione nei convertitori DC-DC.
Intervallo di Tensione Operativa	1.5\, V - 4.0\, V	2.5\, V - 3.65\, V	UL 1973 / IEC 62619	Richiede soglie di protezione BMS specifiche; la tecnologia Na-ion tollera la scarica a 0\, V senza danni strutturali critici.
Densità Energetica Gravimetrica	130 - 170\, Wh/kg	160 - 200\, Wh/kg	IEC 62485-5	La tecnologia Na-ion richiede un maggiore ingombro planimetrico ( $m ^2$ ) e un aumento del peso strutturale nei rack delle batterie.
Densità Energetica Volumetrica	250 - 320\, Wh/L	300 - 400\, Wh/L	IEC 62485-5	Impatta sul design della ventilazione forzata e sulla separazione spaziale nei container BESS di tipo ISO.
Resistenza Interna Relativa ( $R_{ int }$ )	Bassa (<0.5\, m \Omega per cella equivalente)	Media-Bassa (0.5 - 1.0\, m \Omega)	IEC 61951 / IEEE 1188	La tecnologia Na-ion offre minori perdite per effetto Joule ( $I^2 R$ ) interne durante scariche profonde ad alto tasso C.
Capacità di Scarica a Bassa Temperatura	Eccellente (>85\% a -20^\circ C )	Moderata-Bassa (50 - 65\% a -20^\circ C )	IEC 62619 / SAE J2464	La tecnologia Na-ion elimina la necessità di sistemi di preriscaldamento intensivo in ambienti artici o sottostazioni esposte.
Massimo C-Rate di Carica	$3C - 5C$ continuo	$1C - 2C$ continuo	IEC 62619	Riduzione significativa del tempo di ricarica nelle microreti e supporto per servizi ausiliari ad alta frequenza.
Cicli di Vita Utile (80% SOH)	$3000 - 5000$ cicli	$4000 - 6000$ cicli	IEC 62619 / IEEE 1679	Entrambe le tecnologie superano le soglie commerciali standard, richiedendo modelli di degradazione di Arrhenius modificati.

Analisi Termodinamica e Cinematica delle Curve di Scarica

La modellazione analitica delle curve di scarica è fondamentale per dimensionare correttamente i convertitori di potenza e prevedere la caduta di tensione dinamica in regimi di carico transitori. La curva di scarica di una batteria rappresenta la relazione funzionale tra la tensione ai terminali $V_{ term }$ e lo stato di scarica (Depth of Discharge, DoD) o la capacità estratta, sotto una determinata corrente $I$ e temperatura $T$ .

Per la tecnologia LFP, il modello semi-empirico di Shepherd o il modello avanzato di Nernst-Duhem vengono adattati incorporando termini che rappresentano la polarizzazione ohmica, di attivazione e di concentrazione. La tensione ai terminali si esprime come:

V_{ term }(t) = E_0 - R_{ int } I(t) - K \frac{Q}{Q - Q_{ ext }(t)} \left( I^*(t) + Q_{ ext }(t) \right) + A \exp\left( -B \cdot Q_{ ext }(t) \right)

Dove $E_0$ è la tensione a circuito aperto costante durante il plateau, $R_{ int }$ è la resistenza interna equivalente, $K$ è il coefficiente di polarizzazione, $Q$ è la capacità massima della cella, $Q_{ ext }$ è la carica estratta, I^*(t) è la corrente filtrata a bassa frequenza, e $A, B$ rappresentano la zona di caduta esponenziale iniziale e finale della curva LFP.

In netto contrasto, le celle Na-ion basate su catodi stratificati non presentano il plateau di tensione esteso tipico delle LFP. La loro curva di scarica è caratterizzata da una pendenza quasi lineare e continua dallo stato di carica al 100% fino allo 0%. Questa caratteristica viene modellata matematicamente incorporando un termine lineare dipendente dallo stato di carica ( $SoC$ ):

V_{ term, Na }(t) = U_{ 0,base } + m \cdot SoC (t) - I(t) R_{ int }(T, SoC ) - \eta_{ act }(I) - \eta_{ diff }(t)

Dove $U_{ 0,base }$ è la tensione di base di riferimento, $m$ esprime il coefficiente di pendenza caratteristico del materiale stratificato, \eta_{ act } è il sovrapotenziale di attivazione regolato dall'equazione di Butler-Volmer, e \eta_{ diff } è il sovrapotenziale di diffusione ionica dovuto ai gradienti di concentrazione nell'elettrolita e all'interno degli elettrodi.

Questa differenza nella pendenza della curva di scarica comporta profonde implicazioni per il Sistema di Gestione della Batteria (BMS). Mentre nella tecnologia LFP la pendenza quasi piatta rende estremamente complessa la stima precisa dello stato di carica mediante la tensione a circuito aperto (Open Circuit Voltage, OCV) — richiedendo algoritmi di conteggio dei coulomb (Coulomb Counting) altamente precisi con correzioni per la deriva termica e temporale —, nella tecnologia Na-ion la pendenza pronunciata consente una stima robusta e diretta del SoC basata sulla lettura della tensione ai terminali a riposo, riducendo drasticamente gli errori cumulativi degli algoritmi di stima come il Filtro di Kalman Esteso (EKF).

Comportamento Termico e Meccanismi di Degradazione

Il comportamento termico definisce i limiti operativi di sicurezza e la vita utile prevista dei sistemi BESS. La dissipazione termica e la generazione di calore interno sono regolate dall'equazione di conservazione dell'energia applicata a un volume di controllo che rappresenta una cella elettrochimica:

\rho C_p \frac{\partial T}{\partial t} = \nabla \cdot (k \nabla T) + \dot{q}_{ gen }

Dove \rho è la densità volumetrica del materiale della cella, $C_p$ è la capacità termica specifica, $k$ è il tensore di conducibilità termica (solitamente anisotropo, essendo molto più elevato nel piano degli elettrodi rispetto a quello trasversale agli strati), e \dot{q}_{ gen } è il tasso di generazione di calore volumetrico, composto dal calore Joule e dal calore entropico:

\dot{q}_{ gen } = \frac{I^2 R_{ int }}{V_{ cell }} + \frac{I T}{V_{ cell }} \left( \frac{\partial U_{ oc }}{\partial T} \right)

Le batterie Na-ion dimostrano un vantaggio termodinamico e cinetico cruciale a temperature estreme. A basse temperature (<0^\circ C ), l'energia di attivazione per la desolvatazione e la diffusione degli ioni di sodio nell'elettrolita e nelle interfasi di elettrolita solido (SEI) è spesso inferiore, oppure presenta un incremento minore della viscosità dell'elettrolita rispetto alle formulazioni basate su carbonati di litio. Ciò si traduce in un'impedenza d'interfaccia ( $R_{ SEI }$ ) che aumenta in modo molto più moderato nei sistemi Na-ion al diminuire della temperatura, consentendo correnti di carica e scarica elevate senza incorrere in cadute di tensione catastrofiche né nella temuta deposizione metallica di sodio dendritico, che corrisponderebbe al fenomeno del plating superficiale del litio nelle celle LFP.

La degradazione a lungo termine nelle celle LFP è dominata dalla dissoluzione del ferro nell'elettrolita, dalla crescita progressiva dello strato SEI e dal cracking meccanico delle particelle secondarie dovuto allo stress volumetrico durante le transizioni di fase (\sim 6.8\% di variazione volumetrica nella transizione LiFePO _4 \leftrightarrow FePO _4). D'altro canto, le celle Na-ion sperimentano sollecitazioni volumetriche differenti a seconda del materiale catodico utilizzato. Gli ossidi stratificati subiscono variazioni di volume durante le transizioni di fase ad alte tensioni, ma gli anodi in carbonio duro (hard carbon) utilizzati universalmente nei sistemi Na-ion operano mediante un meccanismo di accumulo misto di adsorbimento-intercalazione (meccanismo "adsorption-insertion" o "slope-plateau"), che assorbe le deformazioni meccaniche con una minore degradazione strutturale nel corso di migliaia di cicli.

Ingegneria della Sicurezza Intrinseca e Analisi Forense dei Guasti

Dal punto di vista della sicurezza industriale e della protezione antincendio negli impianti elettrici (in conformità con le norme NFPA 855 e IEC 62933), la valutazione della vulnerabilità alla fuga termica (Thermal Runaway) rappresenta il parametro di progettazione più critico. La fuga termica si innesca quando un evento esotermico iniziale (cortocircuito interno, sovraccarico di tensione, surriscaldamento esterno) innalza la temperatura della cella fino al punto di avviare reazioni secondarie a catena.

Le celle LFP sono ampiamente riconosciute per la loro stabilità termica superiore rispetto alle chimiche ricche di nichel (come NMC o NCA), grazie alla forte unione covalente del legame fosforo-ossigeno ( $P-O$ ) nel gruppo anionico $[ PO _4]^{3-}$ , che inibisce il rilascio di ossigeno gassoso ad alte temperature. Ciononostante, al raggiungimento di temperature critiche (>200^\circ C ), l'elettrolita organico si decompone in modo esotermico, generando gas infiammabili (carbonati alchilici) e pressioni interne elevate che portano all'apertura della valvola di sfiato (venting) e, in casi estremi, all'ignizione.

Le batterie Na-ion presentano caratteristiche di sicurezza intrinseca che superano sotto certi aspetti quelle delle LFP:

Termodinamica dell'Anodo in Carbonio Duro: Il sodio non forma leghe stabili con l'alluminio a temperatura ambiente (a differenza del litio, che si lega con il collettore di corrente in alluminio a tensioni prossime a 0\, V ). Ciò consente alle celle Na-ion di essere scaricate in sicurezza e immagazzinate completamente scariche a 0\, V senza compromettere l'integrità strutturale del collettore di corrente anodico né provocare micro-cortocircuiti dovuti alla dissoluzione dell'alluminio.
Reattività dell'Elettrolita e dei Sali: I sali di sodio utilizzati (come $NaPF _6$ o $NaClO _4$ in solventi specifici) combinati con il carbonio duro mostrano un'entalpia di reazione esotermica inferiore durante la decomposizione termica iniziale rispetto al sistema equivalente Li-ion.
Profilo Termico di Attivazione del Runaway: Prove calorimetriche a saggio accelerato (Accelerating Rate Calorimetry, ARC) dimostrano che l'inizio dell'autosurriscaldamento esotermico incontrollato ( $T_{ onset }$ ) nelle celle Na-ion commerciali si sposta verso temperature superiori (180^\circ C - 210^\circ C ) con tassi di generazione di calore di picco (\left(\frac{dT}{dt}\right)_{\max}) notevolmente inferiori rispetto alle celle LFP.

Nell'ambito dell'analisi forense di ingegneria elettrica a seguito di guasti catastrofici nei sistemi BESS, il verificarsi di un cortocircuito interno a bassa resistenza ( $R_{ sc }$ ) provoca una scarica locale massiva dell'energia immagazzinata nella cella interessata attraverso il punto di guasto, modellata da:

I_{ sc }(t) = \frac{U_{ oc }(t)}{R_{ sc } + R_{ int }(t)}

In un banco LFP, un cortocircuito interno può tradursi rapidamente in un aumento di temperatura localizzato di diverse migliaia di gradi al secondo nel nucleo della cella. In un sistema Na-ion, grazie alla minore energia immagazzinata per unità di volume e alla minore reattività chimica dei componenti a livello d'interfaccia, la velocità di propagazione del guasto termico alle celle adiacenti (Thermal Propagation) è significativamente inferiore, facilitando il contenimento mediante barriere in aerogel e sistemi di estinzione incendi basati su agenti puliti o nebulizzazione d'acqua con additivi.

Strategie Pratiche di Progettazione, Mitigazione e Analisi con Vexten Suite

Per illustrare l'applicazione rigorosa dei concetti sopra esposti nella progettazione ingegneristica reale, viene simulata l'analisi operativa e delle protezioni elettriche di un impianto BESS industriale da 5\, MW / 20\, MWh utilizzando la suite di calcolo Vexten Suite (basata sulle norme IEC 60909, IEC 60287 e IEEE 141).

Considerando un'architettura in media tensione a 33\, kV con molteplici sub-rack di batterie operanti a 1.5\, kV DC , il dimensionamento dei cavi di potenza tra i rack delle batterie e l'inverter centrale deve rispettare rigorosamente i fattori di correzione per temperatura e portata di corrente (ampacity) secondo la norma IEC 60287. La corrente continua massima di progetto per un inverter di P_{ inv } = 2.5\, MW con una tensione nominale del bus DC pari a V_{ dc } = 1200\, V si calcola come:

I_{ dc,max } = \frac{P_{ inv }}{\eta_{ inv } \cdot V_{ dc,min }} = \frac{2.5 \times 10^6\, W }{0.985 \cdot 1000\, V } \approx 2538\, A

Suddividendo questa potenza su $4$ stringhe parallele, la corrente per ciascun feeder principale del cablaggio DC è pari a I_{ feeder } = 634.5\, A . Applicando i fattori di correzione della norma IEC 60287 per cavi unipolari con isolamento in XLPE posati in passerelle perforate con raggruppamento di circuiti:

I_z = \frac{I_{ nominal }}{f_t \cdot f_g \cdot f_d}

Dove $f_t$ è il fattore di temperatura del suolo/aria, $f_g$ è il fattore di raggruppamento delle passerelle ( $0.75$ per $6$ circuiti adiacenti), e $f_d$ è il fattore di profondità/disposizione. Per la tecnologia Na-ion, il design termico del container BESS viene ottimizzato grazie alla minore resistenza termica interna delle celle, consentendo densità di compattazione leggermente superiori senza superare i gradienti termici massimi consentiti dalla norma IEC 62619 (\Delta T < 5^\circ C tra celle adiacenti).

Nell'analisi di cortocircuito secondo IEC 60909 eseguita in Vexten Suite per valutare la capacità di interruzione richiesta per gli interruttori automatici in corrente continua (DC Circuit Breakers, DCCB) e i fusibili ultrarapidi pirotecnici (Pyro-fuses), la corrente di cortocircuito simmetrica iniziale $I''_{k}$ fornita dal banco batterie viene determinata considerando la resistenza interna equivalente dell'intero banco:

I''_{k} = \frac{c \cdot U_n}{\sqrt{3} \cdot Z_k}

Donde $c$ è il fattore di tensione normativo ( $1.05$ per sistemi di generazione e accumulo), $U_n$ è la tensione nominale del sistema, e $Z_k$ è l'impedenza della maglia di cortocircuito che include la resistenza interna delle celle ( $R_{ int }$ ), la resistenza delle barre collettrici (busbars) e l'induttanza parassita dei cavi $L_{ cable }$ . Poiché le celle Na-ion mostrano una resistenza interna leggermente inferiore o paragonabile a quella delle celle LFP in determinate condizioni di temperatura, la corrente di cortocircuito transitoria iniziale può raggiungere picchi molto elevati, richiedendo l'uso di interruttori automatici con tecnologia di estinzione dell'arco tramite soffiaggio magnetico o interruttori a stato solido (Solid-State Circuit Breakers, SSCB) con tempi di estinzione inferiori a 2\, ms .

Inoltre, nell'analisi della distorsione armonica e della risonanza armonica nello stadio di conversione della potenza (inverter multilivello di tipo NPC o MMC), la modulazione di larghezza d'impulso (PWM) genera armoniche di commutazione che interagiscono con le impedenze parassite dei cavi e i filtri di uscita $LCL$ . La frequenza di risonanza del filtro $LCL$ si calcola mediante:

f_{ res } = \frac{1}{2\pi} \sqrt{\frac{L_1 + L_2}{L_1 L_2 C_f}}

Le caratteristiche di impedenza dinamica delle batterie Na-ion ad alta frequenza differiscono da quelle delle LFP a causa della struttura dei loro elettrodi e della conducibilità dell'elettrolita. I modelli a circuito equivalente di Randles utilizzati in Vexten Suite incorporano elementi a fase costante (Constant Phase Elements, CPE) per rappresentare la diffusione anomala e la capacità di doppio strato sulle interfasi. L'impedenza complessa della cella viene formulata come:

Z_{ cell }(\omega) = R_s + \frac{Rct}{1 + j\omega Cdl Rct} + Z_w(\omega)

Dove $R_s$ è la resistenza serie ohmica, $Rct$ è la resistenza di trasferimento di carica, $Cdl$ è la capacità di doppio strato, e $Z_w$ è l'impedenza di Warburg associata alla diffusione allo stato solido. La validazione di questi parametri all'interno di Vexten Suite consente agli ingegneri progettisti di garantire che le armoniche di corrente iniettate dall'inverter non coincidano con le frequenze di risonanza del sistema di accumulo, prevenendo così surriscaldamenti dielettrici negli isolamenti dei cavi e scatti intempestivi delle protezioni di sovracorrente.

Conclusioni e Criteri Finali di Selezione Tecnologica

La valutazione comparativa approfondita tra le tecnologie agli Ioni di Sodio e al Litio-Ferro-Fosfato dimostra che non esiste una soluzione universale ottimale, bensì un dominio applicativo ideale dettato dai requisiti tecnici specifici del progetto. Le celle LFP continuano a dominare in termini di densità energetica volumetrica e gravimetrica assoluta, oltre a beneficiare di uno storico di installazioni commerciali massicce e di una catena di fornitura industriale matura. Tuttavia, le batterie Na-ion rappresentano un progresso tecnologico insostituibile per applicazioni che richiedono:

Funzionamento in ambienti con temperature estreme senza complessi sistemi di condizionamento termico attivo.
Sicurezza intrinseca superiore contro gli eventi di fuga termica e la capacità operativa unica di essere scaricate in sicurezza a 0\, V per facilitare le operazioni di manutenzione e il trasporto sicuro (riduzione della classe Hazmat).
Indipendenza geopolitica da minerali critici e una potenziale riduzione dei costi su scala di gigafactory grazie all'abbondanza universale dei composti del sodio.
Tassi di carica e scarica elevati ( $>3C$ ) con minore degradazione dovuta a sollecitazioni termomeccaniche negli elettrodi in carbonio duro.

L'uso di strumenti analitici avanzati e suite di ingegneria come Vexten Suite è indispensabile per modellare con precisione le curve di scarica non lineari, dimensionare le protezioni contro i cortocircuiti secondo la norma IEC 60909 e garantire la stabilità termica ed elettrica dei sistemi BESS del futuro. L'adozione della tecnologia Na-ion segna una pietra miliare nell'evoluzione dell'ingegneria di potenza e dell'accumulo elettrochimico, consolidando un paradigma più sicuro, sostenibile e robusto per le reti elettriche intelligenti.