🔬 physics

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

This paper investigates a coupled chemo-mechanical-diffusion problem where an elastic solid transforms into a viscoelastic material via a reaction front, analyzing how the competition between chemical affinity, volume expansion-induced stresses, and stress relaxation governs front propagation kinetics and identifies conditions for front blocking and quasi-equilibrium regimes.

Original authors: Alexander Freidin, Aleksandra Ivanova

Published 2026-07-08

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Alexander Freidin, Aleksandra Ivanova

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: A Chemical Race with a Heavy Backpack

Imagine you are trying to push a heavy, expanding balloon through a narrow hallway. As you push, the balloon grows (chemical reaction), but it also gets stuck against the walls, creating pressure (mechanical stress).

This paper studies a specific type of "chemical race" happening inside a solid material. It looks at what happens when a solid material reacts with a substance that diffuses into it (like oxygen rusting metal or lithium entering a battery).

The unique twist in this study is that the new material created by the reaction isn't just a rigid solid; it's viscoelastic. Think of it like honey mixed with rubber. It has some springiness (elasticity) but also flows slowly over time (viscosity). This "honey-rubber" nature allows the material to slowly relax the pressure built up by the reaction, which changes how fast the reaction can spread.

The Main Characters

The Reaction Front: Imagine a line of marching soldiers moving through a field. The "front" is the leading edge where the old material turns into the new material.
The Expansion: When the soldiers (molecules) change uniforms, they suddenly get bigger. This expansion pushes against the surrounding material, creating stress (pressure).
The Stress: This pressure acts like a traffic jam. If the pressure gets too high, it can stop the soldiers from moving forward. This is called "blocking."
The Viscosity (The Honey): Because the new material is like honey, it doesn't hold that pressure forever. Over time, the honey flows, and the pressure slowly leaks away (relaxes).

The Core Conflict: Chemistry vs. Mechanics

The paper explores a tug-of-war between two forces:

Chemistry: The natural desire of the reaction to happen and spread.
Mechanics: The physical pressure building up that tries to stop the reaction.

The researchers use a concept called the "Chemical Affinity Tensor." Think of this as a "Motivation Meter."

If the meter is high (positive), the reaction wants to move forward.
If the meter hits zero, the reaction stops.
The pressure (stress) from the expanding material lowers this meter. If the pressure gets too high, the meter hits zero, and the reaction front freezes.

Key Discoveries

1. The "Stuck" Front Can't Be Unstuck (Usually)

One of the most interesting findings is about what happens when the reaction front gets blocked.

The Scenario: The front stops because the pressure is too high.
The Question: As the "honey-rubber" material slowly relaxes and the pressure drops, will the front start moving again?
The Answer: No. The paper proves that once the front stops, the relaxation of stress cannot unblock it. It's like a car stuck in deep mud; even if the mud hardens and then softens slightly, the car won't suddenly start moving again on its own. The front stays frozen in place.

(Note: The authors mention a very rare, theoretical exception where the material properties are weird enough that the stress relaxation might jump the front to a different "track," but they say this is outside the scope of their current study.)

2. The "Elastic" vs. "Relaxed" Limits

The researchers looked at two extreme scenarios to understand the range of possibilities:

The Elastic Limit (The Rigid Wall): Imagine the new material is super hard rubber that never flows. The pressure builds up instantly and stays there. The reaction stops very quickly.
The Relaxation Limit (The Flowing Honey): Imagine the material flows instantly. The pressure never builds up enough to stop the reaction.
The Reality (Standard Linear Solid Model): Real materials are somewhere in between. They act like stiff rubber at first, but then slowly start to flow like honey. The paper maps out exactly how the reaction speed changes as the material shifts between these two behaviors.

3. The "Damköhler" Number: Who is Faster?

The authors introduce a special number (similar to a "Damköhler number") to compare two speeds:

Speed A: How fast the chemical reaction wants to happen.
Speed B: How fast the new substance can diffuse (travel) through the material.
If Speed A is much faster: The reaction is limited only by how fast the ingredients can arrive (diffusion). The front moves in a "quasi-equilibrium" mode, staying perfectly balanced.
If Speed A is slower: The reaction is limited by its own chemical speed. The front moves in a "kinetic" mode.

The paper shows that if you make the chemical reaction incredibly fast, the system behaves as if it is always in a perfect, balanced state, even though it is technically moving.

The Takeaway

This paper is a mathematical map of how a chemical reaction behaves when the new material it creates is "squishy" (viscoelastic).

Stress is a brake: The expansion of the new material creates pressure that can stop the reaction.
Viscosity is a slow release valve: It lets that pressure leak away over time.
The Paradox: Even though the pressure leaks away, it doesn't help the reaction start again once it has stopped. The "brake" stays on.
The Outcome: Depending on how "stiff" or "flowy" the new material is, and how much external pressure is applied, the reaction might finish completely, or it might get stuck halfway through.

The authors used the Standard Linear Solid Model (a specific way of mathematically describing that "honey-rubber" behavior) and compared it to the simpler Maxwell Model (just honey) to show that the details of the material's "squishiness" significantly change how the reaction front moves and where it eventually stops.

Technical Summary: On a Coupled Problem of Chemo-mechanics with a Viscoelastic Reaction Product

Problem Statement
The paper investigates a coupled "chemistry–mechanics–diffusion" problem involving a chemical reaction between a solid constituent and a diffusing reactant. The reaction is localized at a sharp interface (the reaction front) and transforms the initial elastic solid into a viscoelastic product. This transformation is accompanied by a volumetric strain (eigenstrain), which generates mechanical stresses. These stresses, in turn, influence the reaction kinetics, potentially accelerating, retarding, or blocking the propagation of the reaction front.

The study focuses on a plane strain problem where a planar reaction front propagates through a plate. Unlike previous studies that assumed zero stress in the direction perpendicular to the front, this work considers external deformations applied in both directions (parallel and perpendicular to the front). The reaction product is modeled using the Standard Linear Solid Model (SLSM), with the Maxwell model treated as a specific limiting case. The primary objective is to analyze the joint kinetics of the reaction front propagation and the relaxation of stresses generated by the transformation and external strains.

Methodology
The theoretical framework is built upon the concept of the chemical affinity tensor, a thermodynamic force driving the reaction. The normal component of this tensor ( $A_{nn}$ ) acts as the conjugate force to the reaction rate at the oriented area element of the front.

Governing Equations:
- Diffusion: Fick's law is applied to the diffusing constituent within the transformed material. A stationary diffusion approximation is adopted, assuming the concentration profile is linear across the transformed layer.
- Mechanics: The problem utilizes small strain theory. The stress-strain state is determined by equilibrium equations and constitutive laws. The viscoelastic behavior of the product is described by the SLSM, decomposing stress and strain into hydrostatic (elastic) and deviatoric (viscoelastic) parts.
- Kinetics: The reaction rate is governed by a kinetic equation where the rate depends on the difference between the actual concentration of the diffusing reactant at the front and an equilibrium concentration ( $c_{eq}$ ). The equilibrium concentration is a function of the chemical affinity, which includes both chemical energy parameters and a mechanical contribution ( $\chi$ ) derived from the stress state.
Analytical Approach:
- The authors derive a system of differential equations describing the evolution of the front position ( $\xi$ ) and the stress ( $\sigma_y$ ) perpendicular to the front.
- To cover the full range of viscosity coefficients, two limiting cases are introduced: the elastic limit ( $\eta \to \infty$ ), representing an instantaneous elastic response, and the relaxation limit ( $\eta \to 0$ ), representing the fully relaxed stress state after the front stops or the transformation completes.
- The study distinguishes between the kinetic mode (driven by non-zero chemical affinity) and a quasi-equilibrium mode (where the affinity remains zero during propagation).

Key Contributions and Results

Stress-Dependent Blocking: The paper demonstrates that mechanical stresses can block the reaction front. Blocking occurs when the normal component of the chemical affinity tensor drops to zero, corresponding to a state where the equilibrium concentration equals the actual concentration. The blocking condition depends on the interplay between chemical energy, transformation strain, elastic moduli, and external strains.
Joint Kinetics of Propagation and Relaxation: The study derives the coupled evolution of the front position and stress. It reveals that stress relaxation can either accelerate or slow down the front velocity depending on the specific material parameters and external loading conditions.
Impossibility of Unblocking: A significant finding is that once a reaction front is blocked, subsequent stress relaxation cannot unblock the front. The analysis shows that the sign of the stress rate at the moment of blocking determines the direction of relaxation, but the stress remains within the "forbidden zone" defined by the blocking conditions immediately after the stop. Unblocking is only theoretically possible in a special case where the fully relaxed stress lies outside the range of blocking stresses (a scenario noted as beyond the scope of the current detailed analysis).
Comparison of Models (SLSM vs. Maxwell): The authors introduce a parameter $\beta$ to quantify the deviation of the SLSM from the Maxwell model. They prove that the relaxation limit stress for the SLSM always lies between the elastic limit stress and the Maxwell relaxation limit stress. Consequently, the front velocity and blocking position for the SLSM fall within the bounds set by the pure elastic and pure Maxwell cases.
Quasi-Equilibrium Mode: The paper introduces a quasi-equilibrium propagation regime where the chemical affinity is zero. A dimensionless number, analogous to the Damköhler number ($Da$), is defined to compare the characteristic times of diffusion/mass supply and chemical reaction. The results show that as the reaction rate parameter increases (high $Da$), the kinetic solution converges to the quasi-equilibrium limit.

Significance
The paper provides a systematic analytical framework for understanding how viscoelasticity influences chemo-mechanical coupling. By moving beyond purely elastic assumptions, the work highlights the complex interplay between stress generation, stress relaxation, and reaction kinetics.

The authors emphasize that the effect of viscosity is not straightforward; it can lead to various regimes including complete transformation, incomplete transformation with blocking, or simultaneous blocking and relaxation completion. The introduction of the chemical affinity tensor allows for a rigorous treatment of the thermodynamic driving forces in the presence of stress. The findings are particularly relevant for systems involving large transformation strains (such as silicon lithiation or oxidation) where viscoelastic relaxation plays a critical role in mitigating or altering stress-induced reaction blocking. The study serves as a continuation of previous work on chemo-mechanics, extending the analysis to two-dimensional strain states and more complex rheological models.