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

Resumo Técnico: Capturando a Forma em U de Calandragem na Condutividade Térmica de Eletrodos de Íon-Lítio

Definição do Problema
O gerenciamento térmico em baterias de íon-lítio depende de modelos precisos da condutividade térmica efetiva através do plano ($\lambda_{eff}$) para cada camada porosa. Existe uma lacuna crítica na modelagem da dependência da $\lambda_{eff}$ em relação ao processo de calandragem (compactação por rolo). Dados experimentais revelam uma dependência não monotônica, em formato de "u", da condutividade em relação à porosidade durante a calandragem: a condutividade cai inicialmente antes de subir novamente. Modelos analíticos existentes, como modelos baseados apenas em porosidade ou teorias de meio efetivo de contato estático (ex: Bruggeman, Maxwell-Eucken, Zehner–Bauer–Schlünder padrão), falham em reproduzir essa forma em U, prevendo apenas um aumento monotônico. Além disso, os modelos padrão frequentemente ignoram o efeito Knudsen (condução de gás rarefeito em poros micrométricos) e a evolução das redes de contato de partículas ou a reorientação cristalográfica durante o processamento. Essa falta de um fechamento analítico que explique a forma em U introduz erros dependentes do processo que podem chegar a uma ordem de magnitude.

Metodologia
Os autores propõem um fechamento consciente da calandragem baseado na teoria de meio efetivo Zehner–Bauer–Schlünder (ZBS), aprimorado com três capacidades específicas:

1. Base Corrigida por Knudsen: O modelo incorpora a correção de Smoluchowski/Knudsen para a condutividade do gás nos poros, o que é significativo em poros de escala micrométrica (especialmente sob purga de hélio usada em Análise de Flash de Laser, LFA).
2. Termo de Contato Consciente da Calandragem: Para abordar a forma em U, os autores introduzem um termo de contato indexado pela compressão, $\phi(\Pi)$, onde $\Pi$ é a taxa de compressão. Este termo modela a evolução da rede de pontes sólido-sólido. Ele é formulado como uma função quadrática, $\phi(\Pi) = \max(0, \phi_0 + a\Pi + b\Pi^2)$, para capturar o dano inicial à rede de aglutinante (cisalhamento, $a < 0$) seguido pela recuperação via intertravamento de partículas ($b \geq 0$) em compressões mais altas.
3. Reorientação Anisotrópica (Ânodo de Grafite): Para ânodos de grafite, o modelo inclui um mecanismo de reorientação limitado. As escamas de grafite são altamente anisotrópicas ($\lambda_a \gg \lambda_c$). O modelo postula que a calandragem alinha as escamas no plano, forçando o caminho de calor através do plano a cruzar o fraco eixo $c$. Isso é modelado via um fator de orientação de Hermans $S(\Pi)$ que evolui com a compressão, reduzindo a condutividade sólida efetiva $\lambda_s^\perp$ em direção ao limite do eixo $c$.

O modelo é implementado como um fechamento diferenciável de baixa dimensão usando JAX, permitindo tanto a predição direta quanto a inferência inversa (ex: estimar a porosidade a partir da condutividade). A calibração visa a condutividade aparente do revestimento através do plano derivada de medições LFA após subtrair a contribuição do coletor de corrente.

Principais Contribuições

• Contribuição Primária: Um fechamento consciente da calandragem que reproduz com sucesso a dependência em forma de U da $\lambda_{eff}$ na porosidade em quatro famílias de eletrodos (dois ânodos de grafite, duas cátodos de NMC). Este é o primeiro variante de baixa dimensão encontrado que captura esse comportamento, reduzindo o Erro Médio Percentual Absoluto (MAPE) intra-amostra de 31,1% (apenas porosidade) para 4,5%.
• Contribuição Secundária (Específica para o Ânodo): Uma conta de reorientação limitada para ânodos de grafite que explica a condutividade sólida ajustada anormalmente baixa ao alinhar as escamas em direção ao eixo $c$. Embora os dados de condutividade sozinhos não possam distinguir unicamente isso de um ajuste apenas de contato, o modelo de reorientação é mais parcimonioso e fisicamente limitado.
• Inverso Diferenciável: A implementação permite a inversão do fechamento para inferir a porosidade do revestimento ou a qualidade da interface a partir de medições térmicas, suportando o controle de qualidade em linha.
• Análise de Identificabilidade: O artigo aborda rigorosamente a identificabilidade dos parâmetros dado o pequeno conjunto de dados (6–8 estados por família). Ele reporta quantidades agrupadas (ex: condutância da ponte $G_0 = \phi_0 \lambda_b$) e vales de verossimilhança de perfil em vez de reivindicar estimativas de ponto únicas para parâmetros confundidos como $\lambda_s$ e $\phi_0$.

Resultados

• Precisão: O fechamento calibrado atinge um MAPE intra-amostra de 4,5% em 27 estados de calandragem, correspondendo à incerteza estimada do alvo LFA reconstruído. O erro de "deixar um estado de fora" (held-out) é de 7,8%.
• Discriminação de Mecanismo:
  • Cátodos: A forma em U é impulsionada inteiramente pela evolução da rede de contato (dano e recuperação). A reorientação é inerte para partículas de NMC quase isotrópicas.
  • Ânodos: Tanto um modelo de apenas contato quanto um modelo acoplado à reorientação se ajustam bem aos dados (MAPE ~1,8–5,4%). O modelo de reorientação é favorecido pela parcimônia e por sua capacidade de explicar o baixo $\lambda_s$ via alinhamento para o eixo $c$, mas permanece degenerado com o modelo de contato até que a textura de folha de XRD seja medida.
• Consistência Externa: O modelo passa em vários testes de consistência:
  • Prevê o sinal correto do coeficiente de temperatura (negativo, indicando controle pelo sólido) sem ajustar dados de temperatura.
  • Prevê a condutividade no plano ($\lambda_{\parallel}$) consistente com medições fototérmicas independentes, apesar de ter sido ajustado apenas para dados através do plano.
  • Extrapola para um stack de célula 18650 comercial com erro de ~15% quando a resistência de interface é inferida.
• Viabilidade de Inversão: Em dados sintéticos, o inverso diferenciável recupera a porosidade com um RMSE de 0,008. Em dados reais, o erro é mais grosseiro (0,05) devido a resíduos de forma de modelo e acoplamento de metrologia-estado, mas permanece uma ferramenta viável para verificações secundárias.

Significância e Alegações
O artigo afirma que o fechamento proposto é o primeiro modelo analítico de baixa dimensão capaz de capturar a forma em U da calandragem, preenchendo a lacuna entre simulações resolvidas por microestrutura computacionalmente caras e ajustes empíricos simples. Sua significância reside em:

1. Resolução de Processo: Ele vincula explicitamente a condutividade térmica às variáveis de processamento (compressão), indo além dos modelos estáticos baseados em porosidade.
2. Insight de Mecanismo: Identifica que a forma em U surge da evolução da rede de contato (dano seguido de intertravamento) e, para ânodos, potencialmente da reorientação cristalográfica.
3. Utilidade Prática: A natureza diferenciável do fechamento permite seu uso como um "sensor suave" para controle de qualidade em linha para porosidade e monitoramento de interface em manufatura contínua de eletrodos, uma capacidade não oferecida por simulações resolvidas ou regressões puramente baseadas em dados.
4. Escopo Metodológico: Os autores são modestos quanto aos valores absolutos dos parâmetros ajustados, observando que eles são "dentro do escopo do método" para o alvo derivado de LFA e mudam significamente (1,8–2,4×) se recalibrados contra dados de placa quente guardada. As descobertas robustas são os parâmetros agrupados, as direções de tendência e a capacidade de discriminar mecanismos via medições ortogonais (XRD, adesão, SEM).

O artigo conclui que, embora a divisão entre contato versus reorientação para ânodos permaneça não resolvida pelos dados de condutividade sozinhos, o fechamento fornece um framework falsificável. O próximo passo proposto é uma campanha experimental dedicada envolvendo medições de gás de duas pressões, análise de textura de XRD e metrologia de interface pareada para separar definitivamente esses mecanismos e validar as suposições físicas do fechamento.