On a coupled problem of chemo-mechanics with a viscoelastic reaction product

Resumen Técnico: Sobre un problema acoplado de quimio-mecánica con un producto viscoelástico

Planteamiento del problema
El artículo investiga un problema acoplado de "química–mecánica–difusión" que involucra una reacción química entre un constituinte sólido y un reactivo en difusión. La reacción se localiza en una interfaz nítida (el frente de reacción) y transforma el sólido elástico inicial en un producto viscoelástico. Esta transformación viene acompañada de una deformación volumétrica (deformación propia o eigenstrain), la cual genera tensiones mecánicas. Estas tensiones, a su vez, influyen en la cinética de la reacción, pudiendo acelerar, retardar o bloquear la propagación del frente de reacción.

El estudio se centra en un problema de deformación plana, donde un frente de reacción plano se propaga a través de una placa. A diferencia de estudios previos que asumían una tensión nula en la dirección perpendicular al frente, este trabajo considera deformaciones externas aplicadas en ambas direcciones (paralela y perpendicular al frente). El producto de la reacción se modela mediante el Modelo de Sólido Lineal Estándar (SLSM), tratando el modelo de Maxwell como un caso límite específico. El objetivo principal es analizar la cinética conjunta de la propagación del frente de reacción y la relajación de las tensiones generadas por la transformación y las deformaciones externas.

Metodología
El marco teórico se construye sobre el concepto de tensor de afinidad química, una fuerza termodinámica que impulsa la reacción. La componente normal de este tensor ($A_{nn}$) actúa como la fuerza conjugada a la velocidad de reacción en el elemento de área orientado del frente.

1. Ecuaciones de gobierno:

   • Difusión: Se aplica la ley de Fick al constituinte en difusión dentro del material transformado. Se adopta una aproximación de difusión estacionaria, asumiendo que el perfil de concentración es lineal a través de la capa transformada.
   • Mecánica: El problema utiliza la teoría de pequeñas deformaciones. El estado de tensión-deformación se determina mediante ecuaciones de equilibrio y leyes constitutivas. El comportamiento viscoelástico del producto se describe mediante el SLSM, descomponiendo la tensión y la deformación en partes hidrostáticas (elásticas) y desviadoras (viscoelásticas).
   • Cinética: La velocidad de reacción está gobernada por una ecuación cinética donde la velocidad depende de la diferencia entre la concentración real del reactivo en difusión en el frente y una concentración de equilibrio ($c_{eq}$). La concentración de equilibrio es una función de la afinidad química, la cual incluye parámetros de energía química y una contribución mecánica ($\chi$) derivada del estado de tensión.

2. Enfoque analítico:

   • Los autores derivan un sistema de ecuaciones diferenciales que describen la evolución de la posición del frente ($\xi$) y la tensión ($\sigma_y$) perpendicular al frente.
   • Para cubrir todo el rango de coeficientes de viscosidad, se introducen dos casos límite: el límite elástico ($\eta \to \infty$), que representa una respuesta elástica instantánea, y el límite de relajación ($\eta \to 0$), que representa el estado de tensión totalmente relajado después de que el frente se detiene o la transformación se completa.
   • El estudio distingue entre el modo cinético (impulsado por una afinidad química no nula) y un modo de cuasi-equilibrio (donde la afinidad permanece en cero durante la propagación).

Contribuciones clave y resultados

• Bloqueo dependiente de la tensión: El artículo demuestra que las tensiones mecánicas pueden bloquear el frente de reacción. El bloqueo ocurre cuando la componente normal del tensor de afinidad química cae a cero, lo que corresponde a un estado donde la concentración de equilibrio es igual a la concentración real. La condición de bloqueo depende de la interacción entre la energía química, la deformación de transformación, los módulos elásticos y las deformaciones externas.
• Cinética conjunta de propagación y relajación: El estudio deriva la evolución acoplada de la posición del frente y la tensión. Revela que la relajación de la tensión puede acelerar o ralentizar la velocidad del frente dependiendo de los parámetros específicos del material y las condiciones de carga externa.
• Imposibilidad de desbloqueo: Un hallazgo significativo es que, una vez que un frente de reacción es bloqueado, la relajación de la tensión subsecuente no puede desbloquear el frente. El análisis muestra que el signo de la tasa de tensión en el momento del bloqueo determina la dirección de la relajación, pero la tensión permanece dentro de la "zona prohibida" definida por las condiciones de bloqueo inmediatamente después de la detención. El desbloqueo es solo teóricamente posible en un caso especial donde la tensión totalmente relajada se encuentra fuera del rango de las tensiones de bloqueo (un escenario señalado como más allá del alcance del análisis detallado actual).
• Comparación de modelos (SLSM vs. Maxwell): Los autores introducen un parámetro $\beta$ para cuantificar la desviación del SLSM respecto al modelo de Maxwell. Demuestran que la tensión del límite de relajación para el SLSM siempre se sitúa entre la tensión del límite elástico y la tensión del límite de relajación de Maxwell. En consecuencia, la velocidad del frente y la posición de bloqueo para el SLSM caen dentro de los límites establecidos por los casos de elástico puro y de Maxwell puro.
• Modo de cuasi-equilibrio: El artículo introduce un régimen de propagación de cuasi-equilibrio donde la afinidad química es cero. Se define un número adimensional, análogo al número de Damköhler ($Da$), para comparar los tiempos característicos de difusión/suministro de masa y la reacción química. Los resultados muestran que a medida que aumenta el parámetro de velocidad de reacción (alto $Da$), la solución cinética converge al límite de cuasi-equilibrio.

Significancia
El artículo proporciona un marco analítico sistemático para comprender cómo la viscoelasticidad influye en el acoplamiento quimio-mecánico. Al ir más allá de las suposiciones puramente elásticas, el trabajo destaca la compleja interacción entre la generación de tensión, la relajación de la tensión y la cinética de la reacción.

Los autores enfatizan que el efecto de la viscosidad no es sencillo; puede conducir a varios regímenes, incluyendo la transformación completa, la transformación incompleta con bloqueo, o el bloqueo simultáneo y la finalización de la relajación. La introducción del tensor de afinidad química permite un tratamiento riguroso de las fuerzas impulsoras termodinámicas en presencia de tensión. Los hallazgos son particularmente relevantes para sistemas que involucran grandes deformaciones de transformación (como la litación de silicio o la oxidación) donde la relajación viscoelástica juega un papel crítico en la mitigación o alteración del bloqueo de la reacción inducido por la tensión. El estudio sirve como una continuación del trabajo previo sobre quimio-mecánica, extendiendo el análisis a estados de deformación bidimensionales y modelos reológicos más complejos.