Baterías de Sodio-Ion vs LFP en Almacenamiento Energético: Comportamiento Térmico y Curvas de Descarga en Campo

Introducción y Fundamentos Termodinámicos y Electroquímicos de las Tecnologías Na-ion y LiFePO₄

La transición energética global requiere sistemas de almacenamiento de energía electroquímica (BESS) que no solo maximicen la densidad energética, sino que garanticen una resiliencia térmica superior, seguridad intrínseca y una cadena de suministro desvinculada de materias primas críticas como el litio y el cobalto. Dentro de este contexto, las baterías de Sodio-Ion (Na-ion) emergen como una alternativa disruptiva frente a las químicas de Litio-Hierro-Fosfato (LiFePO₄ o LFP). Desde la perspectiva de la ingeniería eléctrica de alta potencia y el diseño de sistemas avanzados, es imperativo analizar los fundamentos termodinámicos y cinéticos que rigen ambas tecnologías para modelar correctamente su comportamiento operativo, sus curvas de descarga y sus perfiles de seguridad bajo condiciones de estrés severo.

Termodinámicamente, la diferencia fundamental radica en el potencial redox estándar del par iónico y el radio atómico. El sodio presenta un radio iónico mayor (1.02\, Å frente a los 0.76\, Å del litio) y una masa molar superior (22.99\, g/mol frente a 6.94\, g/mol ). Esto se traduce en un potencial estándar de reducción más electronegativo para el litio (-3.04\, V vs. $SHE$ ) en comparación con el sodio (-2.71\, V vs. $SHE$ ), lo que establece un límite superior teórico para la tensión de celda en las baterías Na-ion. Sin embargo, la energía libre de Gibbs (\Delta G) de las reacciones de intercalación define los perfiles de voltaje operativo y la estabilidad estructural de los electrodos anódicos y catódicos.

En las celdas LiFePO₄, el cátodo opera mediante un mecanismo de inserción de dos fases (transición de fase de primer orden) entre $LiFePO _4$ y $FePO _4$ , lo que genera una meseta de descarga extraordinariamente plana alrededor de 3.2\, V a 3.3\, V . Por el contrario, los sistemas Na-ion actuales utilizan una variedad más amplia de materiales catódicos, predominantemente óxidos de metales de transición层状 (Layered Transition Metal Oxides) y análogos de azul de Prusia (Prussian Blue Analogues). Estos materiales exhiben frecuentemente soluciones sólidas de fase única en amplios rangos de estado de carga (SoC), lo que confiere a las curvas de descarga del Na-ion una pendiente característica en lugar de una meseta plana, complicando los algoritmos de estimación de SoC pero ofreciendo ventajas termodinámicas en términos de gestión del calor y prevención de sobretensiones localizadas.

La energía interna del sistema y su variación térmica se rigen por la ecuación fundamental de la calorimetría electroquímica, la cual relaciona el calor generado con la entropía de la reacción:

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

Donde \dot{Q} es la tasa de generación de calor, $I$ es la corriente de operación ( $>0$ en descarga), $U_{ oc }$ es el voltaje de circuito abierto, $V_{ term }$ es el voltaje terminal de la celda, $T$ es la temperatura absoluta, y el término \left(\frac{\partial U_{ oc }}{\partial T}\right) representa el coeficiente entrópico. En las celdas LFP, el coeficiente entrópico varía drásticamente cerca de los extremos de carga debido a las transiciones de fase, provocando picos locales de generación de calor entrópico. En las celdas Na-ion, la ausencia de transiciones de fase abruptas en ciertos diseños de cátodos en capas suaviza este comportamiento, reduciendo la susceptibilidad a puntos calientes internos inducidos por gradientes entrópicos durante pulsos de alta corriente.

Parámetros Eléctricos, Químicos y Normativos Comparativos

Para establecer una base de diseño rigurosa en proyectos de ingeniería eléctrica de media y gran escala, es indispensable contrastar los parámetros intrínsecos de las celdas Na-ion y LFP bajo normativas internacionales estandarizadas (tales como IEC 62619, IEEE 1679, UL 1973 y UL 9540A). La siguiente tabla consolida las métricas críticas de desempeño eléctrico, límites normativos y consecuencias operativas.

Parámetro Técnico / Métrica	Tecnología Sodio-Ion (Na-ion)	Tecnología Litio-Hierro-Fosfato (LFP)	Normativa de Referencia / Estándar IEC/IEEE/UL	Consecuencia Operativa y Dieléctrica
Voltaje Nominal de Celda	3.0\, V - 3.3\, V	3.2\, V - 3.3\, V	IEC 62619 / IEEE 1679	Afecta la topología de los inversores de potencia y la relación de transformación en convertidores DC-DC.
Rango de Voltaje Operativo	1.5\, V - 4.0\, V	2.5\, V - 3.65\, V	UL 1973 / IEC 62619	Requiere umbrales de protección BMS específicos; Na-ion tolera descarga a 0\, V sin daño estructural crítico.
Densidad Energética Gravimétrica	130 - 170\, Wh/kg	160 - 200\, Wh/kg	IEC 62485-5	Na-ion exige mayor huella física ( $m ^2$ ) e incremento de peso estructural en racks de baterías.
Densidad Energética Volumétrica	250 - 320\, Wh/L	300 - 400\, Wh/L	IEC 62485-5	Impacta el diseño de ventilación forzada y separación espacial en contenedores BESS tipo ISO.
Resistencia Interna Relativa ( $R_{ int }$ )	Baja (<0.5\, m \Omega por celda equivalente)	Media-Baja (0.5 - 1.0\, m \Omega)	IEC 61951 / IEEE 1188	Na-ion ofrece menores pérdidas Joule ( $I^2 R$ ) internas bajo descargas profundas a alta tasa C.
Capacidad de Descarga a Baja Temperatura	Excelente (>85\% a -20^\circ C )	Moderada-Baja (50 - 65\% a -20^\circ C )	IEC 62619 / SAE J2464	Na-ion elimina la necesidad de sistemas de precalentamiento intensivo en entornos árticos o subestaciones expuestas.
Velocidad de Carga C-Rate Máxima	$3C - 5C$ continuo	$1C - 2C$ continuo	IEC 62619	Reducción significativa del tiempo de recarga en microrredes y soporte de servicios auxiliares de alta frecuencia.
Ciclos de Vida Útil (80% SOH)	$3000 - 5000$ ciclos	$4000 - 6000$ ciclos	IEC 62619 / IEEE 1679	Ambas tecnologías superan los umbrales comerciales estándar, requiriendo modelos de degradación de Arrhenius modificados.

Análisis Termodinámico y Cinemática de las Curvas de Descarga

El modelado analítico de las curvas de descarga es fundamental para dimensionar correctamente los convertidores de potencia y predecir la caída de tensión dinámica bajo regímenes de carga transitorios. La curva de descarga de una batería representa la relación funcional entre el voltaje terminal $V_{ term }$ y el estado de descarga (Depth of Discharge, DoD) o capacidad extraída, bajo una corriente $I$ y una temperatura $T$ dadas.

Para la tecnología LFP, el modelo semi-empírico de Shepherd o el modelo avanzado de Nernst-Duhem se adaptan mediante la incorporación de términos que representan la polarización óhmica, de activación y de concentración. El voltaje terminal se expresa como:

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)

Donde $E_0$ es el voltaje de circuito abierto constante durante la meseta, $R_{ int }$ es la resistencia interna equivalente, $K$ es el coeficiente de polarización, $Q$ es la capacidad máxima de la celda, $Q_{ ext }$ es la carga extraída, I^*(t) es la corriente filtrada de baja frecuencia, y $A, B$ representan la zona de caída exponencial inicial y final de la curva LFP.

En marcada contraposición, las celdas Na-ion basadas en cátodos en capas no presentan la meseta de voltaje extendida de LFP. Su curva de descarga se caracteriza por una pendiente casi lineal continua desde el estado de carga al 100% hasta el 0%. Esta característica se modela matemáticamente incorporando un término lineal dependiente del estado de carga ( $SoC$ ):

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

Donde $U_{ 0,base }$ es el voltaje base de referencia, $m$ es el coeficiente de pendiente característico del material en capas, \eta_{ act } es el sobrepotencial de activación regido por la ecuación de Butler-Volmer, y \eta_{ diff } es el sobrepotencial de difusión iónica debido a los gradientes de concentración en el electrolito y el seno de los electrodos.

Esta diferencia en la pendiente de la curva de descarga conlleva profundas implicaciones para el Sistema de Gestión de Baterías (BMS). Mientras que en LFP la pendiente casi plana dificulta la estimación precisa del estado de carga mediante tensión en circuito abierto (Open Circuit Voltage, OCV) —requiriendo contadores de culombios (Coulomb Counting) altamente precisos con correcciones por deriva—, en Na-ion la pendiente pronunciada permite una estimación robusta y directa del SoC basada en la lectura del voltaje terminal en reposo, reduciendo drásticamente los errores acumulativos de los algoritmos de estimación como el Filtro de Kalman Extendido (EKF).

Comportamiento Térmico y Mecanismos de Degradación

El comportamiento térmico define los límites operativos de seguridad y la vida útil proyectada de los sistemas BESS. La disipación térmica y la generación de calor interno se rigen por la ecuación de conservación de la energía aplicada a un volumen de control que representa una celda electroquímica:

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

Donde \rho es la densidad volumétrica del material de la celda, $C_p$ es la capacidad calorífica específica, $k$ es el tensor de conductividad térmica (usualmente anisotrópico, siendo mucho mayor en el plano de los electrodos que a través de las capas), y \dot{q}_{ gen } es la tasa de generación de calor volumétrica, compuesta por el calor joule y el calor entrópico:

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

Las baterías Na-ion demuestran una ventaja termodinámica y cinética crucial a temperaturas extremas. A bajas temperaturas (<0^\circ C ), la energía de activación para la desolvatación y difusión de los iones de sodio en el electrolito y las interfases de electrolito sólido (SEI) es frecuentemente inferior o presenta menor incremento en la viscosidad del electrolito en comparación con las formulaciones basadas en carbonatos de litio. Esto se traduce en una impedancia interfacial ( $R_{ SEI }$ ) que aumenta de manera mucho más moderada en Na-ion al descender la temperatura, permitiendo corrientes de carga y descarga elevadas sin incurrir en caídas de tensión catastróficas ni en la temida deposición metálica de sodio dendrítico que equivaldría a la litiación superficial (plating) en LFP.

La degradación a largo plazo en LFP está dominada por la disolución de hierro en el electrolito, el crecimiento progresivo de la capa SEI y el agrietamiento mecánico de las partículas secundarias debido a la tensión volumétrica durante las transiciones de fase (\sim 6.8\% de cambio volumétrico en la transición LiFePO _4 \leftrightarrow FePO _4). Por otro lado, las celdas Na-ion experimentan tensiones volumétricas diferentes según el material catódico utilizado. Los óxidos en capas sufren cambios de volumen durante las transiciones de fase a altos voltajes, pero las ánodas de carbono duro (hard carbon) utilizadas universalmente en Na-ion operan mediante un mecanismo de almacenamiento mixto de adsorción-intercalación (mecanismo "adsorption-insertion" o "slope-plateau"), el cual absorbe las deformaciones mecánicas con menor degradación estructural a lo largo de miles de ciclos.

Ingeniería de Seguridad Intrínseca y Análisis Forense de Fallas

Desde el punto de vista de la seguridad industrial y la protección contra incendios en instalaciones eléctricas (cumpliendo con NFPA 855 e IEC 62933), la evaluación de la vulnerabilidad ante la fuga térmica (Thermal Runaway) es el parámetro de diseño más crítico. La fuga térmica se desencadena cuando un evento exotérmico inicial (cortocircuito interno, sobretensión, sobrecalentamiento externo) eleva la temperatura de la celda hasta el punto de iniciar reacciones secundarias en cadena.

Las celdas LFP son ampliamente reconocidas por su estabilidad térmica superior en comparación con las químicas ricas en níquel (como NMC o NCA), debido a la fuerte unión covalente del enlace fósforo-oxígeno ( $P-O$ ) en el grupo aniónico $[ PO _4]^{3-}$ , el cual inhibe la liberación de oxígeno gaseoso a altas temperaturas. No obstante, al alcanzar temperaturas críticas (>200^\circ C ), el electrolito orgánico se descompone exotérmicamente, generando gases inflamables (carbonatos de alquilo) y presiones internas elevadas que conducen a la apertura de la válvula de alivio (venting) y, en casos extremos, a la ignición.

Las baterías Na-ion presentan características de seguridad intrínseca que superan en ciertos aspectos a las de LFP:

Termodinámica del Ánodo de Carbono Duro: El sodio no forma aleaciones estables con el aluminio a temperatura ambiente (a diferencia del litio, que forma aleaciones con el colector de corriente de aluminio a voltajes cercanos a 0\, V ). Esto permite que las celdas Na-ion se descarguen de forma segura y se almacenen completamente descargadas a 0\, V sin comprometer la integridad estructural del colector de corriente anódico ni provocar microcortocircuitos por disolución del aluminio.
Reactividad del Electrolito y Sales: Las sales de sodio utilizadas (como $NaPF _6$ o $NaClO _4$ en solventes específicos) combinadas con el carbono duro muestran una entalpía de reacción exotérmica menor durante la descomposición térmica inicial en comparación con el sistema Li-ion equivalente.
Perfil Térmico de Activación del Runaway: Pruebas calorimétricas de tasa acelerada (Accelerating Rate Calorimetry, ARC) demuestran que el inicio de la autocalentamiento exotérmico descontrolado ( $T_{ onset }$ ) en celdas Na-ion comerciales se desplaza hacia temperaturas superiores (>180^\circ C a 210^\circ C ) con tasas de generación de calor pico (\left(\frac{dT}{dt}\right)_{\max}) considerablemente más bajas que en LFP.

En el plano del análisis forense de ingeniería eléctrica ante fallas catastróficas en sistemas BESS, la ocurrencia de un cortocircuito interno de baja resistencia ( $R_{ sc }$ ) provoca una descarga local masiva de la energía almacenada en la celda afectada a través del punto de falla, modelada por:

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

En un banco LFP, un cortocircuito interno puede derivar rápidamente en un aumento de temperatura localizado de varios miles de grados por segundo en el núcleo. En un sistema Na-ion, debido a la menor energía almacenada por unidad de volumen y a la menor reactividad química de los componentes a nivel de interfaz, la velocidad de propagación de la falla térmica a las celdas adyacentes (Thermal Propagation) es significativamente menor, facilitando la contención mediante barreras aerogelíticas y sistemas de extinción de incendios basados en agentes limpios o inundación de agua con aditivos.

Estrategias Prácticas de Diseño, Mitigación y Análisis con Vexten Suite

Para ilustrar la aplicación rigurosa de los conceptos anteriores en el diseño de ingeniería del mundo real, se simula el análisis operativo y de protecciones eléctricas de una instalación BESS industrial de 5\, MW / 20\, MWh utilizando la suite de cálculo Vexten Suite (basada en normas IEC 60909, IEC 60287 y IEEE 141).

Considerando una arquitectura en media tensión a 33\, kV con múltiples sub-racks de baterías operando a 1.5\, kV DC , el dimensionamiento de los cables de potencia entre los racks de baterías y el inversor central debe cumplir estrictamente con los factores de corrección por temperatura y capacidad de conducción de corriente según IEC 60287. La corriente continua máxima de diseño para un inversor de P_{ inv } = 2.5\, MW con un voltaje de bus DC nominal V_{ dc } = 1200\, V se calcula como:

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

Dividiendo esta potencia en $4$ strings paralelos, la corriente por cada feeder principal de cableado DC es I_{ feeder } = 634.5\, A . Aplicando los factores de corrección de la norma IEC 60287 para cables unipolares con aislamiento de XLPE instalados en bandejas perforadas con agrupamiento de circuitos:

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

Donde $f_t$ es el factor de temperatura del suelo/aire, $f_g$ es el factor de agrupamiento de bandejas ( $0.75$ para $6$ circuitos adyacentes), y $f_d$ es el factor de profundidad/disposición. Para tecnología Na-ion, el diseño térmico del contenedor BESS se optimiza debido a la menor resistencia térmica interna de las celdas, permitiendo densidades de empaquetamiento ligeramente superiores sin exceder los gradientes térmicos máximos permitidos por la norma IEC 62619 (\Delta T < 5^\circ C entre celdas adyacentes).

En el análisis de cortocircuito según IEC 60909 realizado en Vexten Suite para evaluar la capacidad de interrupción requerida de los interruptores automáticos de corriente continua (DC Circuit Breakers, DCCB) y los fusibles ultrarrápidos pirotécnicos (Pyro-fuses), la corriente de cortocircuito simétrica inicial $I''_{k}$ aportada por el banco de baterías se determina considerando la resistencia interna equivalente del banco completo:

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

Donde $c$ es el factor de tensión normativo ( $1.05$ para sistemas de generación y almacenamiento), $U_n$ es el voltaje nominal del sistema, y $Z_k$ es la impedancia de la malla de cortocircuito incluyendo la resistencia interna de las celdas ( $R_{ int }$ ), la resistencia de las barras colectoras (busbars) y la inductancia parásita de los cables $L_{ cable }$ . Debido a que las celdas Na-ion exhiben una resistencia interna ligeramente inferior o comparable a LFP en ciertas condiciones de temperatura, la corriente de cortocircuito transitoria inicial puede alcanzar picos muy elevados, requiriendo el uso de interruptores automáticos con tecnología de extinción de arco por soplado magnético o semiconductores de estado sólido (Solid-State Circuit Breakers, SSCB) con tiempos de despeje inferiores a 2\, ms .

Asimismo, en el análisis de distorsión armónica y resonancia armónica en la etapa de conversión de potencia (inversores multinivel tipo NPC o MMC), la modulación de ancho de pulso (PWM) genera armónicos de conmutación que interactúan con las impedancias parásitas de los cables y los filtros de salida $LCL$ . La frecuencia de resonancia del filtro $LCL$ se calcula mediante:

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

Las características de impedancia dinámica de las baterías Na-ion a altas frecuencias difieren de las de LFP debido a la estructura de sus electrodos y la conductividad del electrolito. Los modelos de circuitos equivalentes de Randles utilizados en Vexten Suite incorporan elementos de fase constante (Constant Phase Elements, CPE) para representar la difusión anómala y la capacitancia de doble capa en las interfaces. La impedancia compleja de la celda se formula como:

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

Donde $R_s$ es la resistencia serie óhmica, $Rct$ es la resistencia de transferencia de carga, $Cdl$ es la capacitancia de doble capa, y $Z_w$ es la impedancia de Warburg asociada a la difusión en estado sólido. La validación de estos parámetros dentro de Vexten Suite permite a los ingenieros proyectistas asegurar que los armónicos de corriente inyectados por el inversor no coincidan con las frecuencias de resonancia del sistema de almacenamiento, previniendo así sobrecalentamientos dieléctricos en los aislamientos de los cables y disparos molestos de las protecciones de sobrecorriente.

Conclusiones y Criterios Finales de Selección Tecnológica

La evaluación comparativa exhaustiva entre las tecnologías de Sodio-Ion y Litio-Hierro-Fosfato demuestra que no existe una solución universal óptima, sino un dominio de aplicación óptimo dictado por los requerimientos técnicos del proyecto. Las celdas LFP continúan dominando en términos de densidad energética volumétrica y gravimétrica absoluta, así como en un historial de despliegue comercial masivo y una cadena de suministro industrial madura. Sin embargo, las baterías Na-ion representan un avance tecnológico insustituible para aplicaciones que exigen:

Operación en entornos con temperaturas extremas sin sistemas complejos de acondicionamiento térmico activo.
Seguridad intrínseca superior frente a eventos de fuga térmica y la capacidad operativa única de descargarse de forma segura a 0\, V para labores de mantenimiento y transporte seguro (Hazmat class reduction).
Independencia geopolítica de minerales críticos y una reducción potencial de costos a escala de gigafábrica debido a la abundancia universal de los compuestos sodíferos.
Tasas de carga y descarga elevadas ( $>3C$ ) con menor degradación por tensiones termomecánicas en los electrodos de carbono duro.

El uso de herramientas analíticas avanzadas y suites de ingeniería como Vexten Suite es indispensable para modelar con precisión las curvas de descarga no lineales, dimensionar las protecciones contra cortocircuitos según IEC 60909, y garantizar la estabilidad térmica y eléctrica de los sistemas BESS del futuro. La adopción de Na-ion marca un hito en la evolución de la ingeniería de potencia y el almacenamiento electroquímico, consolidando un paradigma más seguro, sostenible y robusto para las redes eléctricas inteligentes.