Capturing the calendering U-shape in lithium-ion electrode thermal conductivity

Resumen Técnico: Captura de la forma en U de la calandrado en la conductividad térmica de electrodos de iones de litio

Planteamiento del Problema
La gestión térmica en las baterías de iones de litio depende de modelos precisos de la conductividad térmica efectiva a través del plano ($\lambda_{eff}$) para cada capa porosa. Existe una brecha crítica en el modelado de la dependencia de $\lambda_{eff}$ respecto al proceso de calandrado (laminado por rodillos). Los datos experimentales revelan una dependencia no monotónica, de "forma en U", de la conductividad respecto a la porosidad durante el calandrado: la conductividad inicialmente cae antes de volver a subir. Los cierres analíticos existentes, como los modelos basados únicamente en la porosidad o las teorías de medio efectivo de contacto estático (p. ej., Bruggeman, Maxwell-Eucken, Zehner–Bauer–Schlünder estándar), no logran reproducir esta forma en U, prediciendo únicamente un aumento monotónico. Además, los modelos estándar suelen ignorar el efecto Knudsen (conducción de gas rarefacto en poros micrométricos) y la evolución de las redes de contacto de partículas o la reorientación cristalográfica durante el procesamiento. Esta falta de un cierre analítico que explique la forma en U introduce errores dependientes del proceso que pueden aproximarse a un orden de magnitud.

Metodología
Los autores proponen un cierre sensible al calandrado basado en la teoría de medio efectivo de Zehner–Bauer–Schlünder (ZBS), mejorada con tres capacidades específicas:

1. Base Corregida por Knudsen: El modelo incorpora la corrección de Smoluchowski/Knudsen para la conductividad del gas en los poros, la cual es significativa en poros de escala micrométrica (especialmente bajo la purga de helio utilizada en el Análisis de Flash Láser, LFA).
2. Término de Contacto Sensible al Calandrado: Para abordar la forma en U, los autores introducen un término de contacto indexado por compresión, $\phi(\Pi)$, donde $\Pi$ es la tasa de compresión. Este término modela la evolución de la red de puentes sólido-sólido. Está formulado como una función cuadrática, $\phi(\Pi) = \max(0, \phi_0 + a\Pi + b\Pi^2)$, para capturar el daño inicial a la red de aglutinante (cizallamiento, $a < 0$) seguido de la recuperación mediante el entrelazamiento de partículas ($b \geq 0$) a mayor compresión.
3. Reorientación Anisotrópica (Ánodo de Grafito): Para los ánodos de grafito, el modelo incluye un mecanismo de reorientación acotado. Las escamas de grafito son altamente anisotrópicas ($\lambda_a \gg \lambda_c$). El modelo postula que el calandrado alinea las escamas en el plano, forzando la trayectoria de calor a través del plano a cruzar el débil eje $c$. Esto se modela mediante un factor de orientación de Hermans $S(\Pi)$ que evoluciona con la compresión, reduciendo la conductividad sólida efectiva $\lambda_s^\perp$ hacia el límite del eje $c$.

El modelo se implementa como un cierre diferenciable de baja dimensión utilizando JAX, lo que permite tanto la predicción directa como la inferencia inversa (p. ej., estimar la porosidad a partir de la conductividad). La calibración se centra en la conductividad aparente de la capa a través del plano derivada de mediciones LFA tras restar la contribución del colector de corriente.

Contribuciones Clave

• Contribución Primaria: Un cierre sensible al calandrado que reproduce con éxito la dependencia en forma de U de la conductividad y la porosidad en cuatro familias de electrodos (dos ánodos de grafito, dos cátodos de NMC). Este es el primer variante de baja dimensión encontrado que captura este comportamiento, reduciendo el Error Porcentual Absoluto Medio (MAPE) intra-muestra del 31.1% (solo porosidad) al 4.5%.
• Contribución Secundaria (Específica para Ánodos): Una cuenta de reorientación acotada para ánodos de grafito que explica la anomalamente baja conductividad sólida medida al alinear las escamas hacia el eje $c$. Aunque los datos de conductividad por sí solos no pueden distinguir de forma única esto de un ajuste de solo contacto, el modelo de reorientación es más parsimonioso y está físicamente acotado.
• Inverso Diferenciable: La implementación permite la inversión del cierre para inferir la porosidad del recubrimiento o la calidad de la interfaz a partir de mediciones térmicas, apoyando el control de calidad en línea.
• Análisis de Identificabilidad: El artículo aborda rigurosamente la identificabilidad de los parámetros dado el pequeño conjunto de datos (6–8 estados por familia). Reporta cantidades agrupadas (p. ej., la conductancia del puente $G_0 = \phi_0 \lambda_b$) y valles de verosimilitud de perfil en lugar de reclamar estimaciones puntuales únicas para parámetros confundidos como $\lambda_s$ y $\phi_0$.

Resultos

• Precisión: El cierre calibrado logra un MAPE intra-muestra del 4.5% en 27 estados de calandrado, coincidiendo con la incertidumbre estimada del objetivo LFA reconstruido. El error de "dejar un estado fuera" (leave-one-state-out) es del 7.8%.
• Discriminación de Mecanismos:
  • Cátodos: La forma en U es impulsada enteramente por la evolución de la red de contacto (daño y recuperación). La reorientación es inerte para partículas de NMC cuasi-isotrópicas.
  • Ánodos: Tanto un modelo de solo contacto como un modelo acoplado a la reorientación ajustan los datos igualmente bien (MAPE ~1.8–5.4%). El modelo de reorientación es favorecido por su parsimonia y su capacidad de explicar la baja $\lambda_s$ mediante la alineación hacia el eje $c$, pero sigue siendo degenerado con el modelo de contacto hasta que se mida la textura por XRD de la misma lámina.
• Consistencia Externa: El modelo supera varios controles de consistencia:
  • Predice el signo correcto del coeficiente de temperatura (negativo, indicando control por sólido) sin ajustar datos de temperatura.
  • Predice la conductividad en el plano ($\lambda_{\parallel}$) consistente con mediciones fototérmicas independientes, a pesar de haber sido ajustado solo a datos a través del plano.
  • Se extrapola a un stack de celdas 18650 comercial con un error de ~15% cuando se infiere la resistencia de la interfaz.
• Viabilidad de Inversa: En datos sintéticos, la inversa diferenciable recupera la porosidad con un RMSE de 0.008. En datos reales, el error es más grueso (0.05) debido a residuos de la forma del modelo y el acoplamiento metrología-estado, pero sigue siendo una herramienta viable para verificaciones secundarias.

Significancia y Reclamaciones
El artículo afirma que el cierre propuesto es el primer modelo analítico de baja dimensión capaz de capturar la forma en U del calandrado, cerrando la brecha entre las simulaciones de resolución microestructural computacionalmente costosas y los ajustes empíricos simples. Su significancia radica en:

1. Resolución de Proceso: Vincula explícitamente la conductividad térmica con las variables de procesamiento (compresión), yendo más allá de los modelos estáticos basados en la porosidad.
2. Información sobre Mecanismos: Identifica que la forma en U surge de la evolución de la red de contacto (daño y luego entrelazamiento) y, para los ánodos, potencialmente de la reorientación cristalográfica.
3. Utilidad Práctica: La naturaleza diferenciable del cierre permite su uso como un "sensor suave" para el control de calidad en línea (monitoreo de porosidad e interfaz) en la fabricación continua de electrodos, una capacidad que no ofrecen las simulaciones de alta resolución ni las regresiones puramente basadas en datos.
4. Alcance Metodológico: Los autores son modestos respecto a los valores absolutos de los parámetros ajustados, señalando que están "vinculados al método" al objetivo derivado de LFA y cambian significativamente (1.8–2.4×) si se recalibran contra datos de placa caliente guardada. Sus hallazgos robustos son las cantidades agrupadas, las direcciones de tendencia y la capacidad de discriminar mecanismos mediante mediciones ortogonales (XRD, adhesión, SEM).

El artículo concluye que, si bien la división entre contacto y reorientación para los ánodos permanece sin resolver por los datos de conductividad por sí solos, el cierre proporciona un marco falsable. El siguiente paso propuesto es una campaña experimental dedicada que incluya mediciones de gas a dos presiones, análisis de textura XRD y metrología de interfaz emparejada para separar definitivamente estos mecanismos y validar los supuestos físicos del cierre.