约束反应性固体混合物的化学-热-力学耦合理论框架

《European Journal of Mechanics - A/Solids》：On the chemo-thermo-mechanics of constrained reactive mixtures of solids

　　

本手稿提出了一个严谨的理论框架，用于描述固体混合物中可能伴随体积变化的化学反应过程。研究重点关注参考构型（reference configuration）的概念及其一致形式化定义所需假设，揭示了该定义对所得控制方程的重要影响。尽管该公式是在一般性框架下建立的，但其涵盖广泛的应用范围。代表性实例包括：具有转化型电极（conversion-type electrodes）的离子电池充放电过程；机械生物学中聚合驱动的现象，如血凝块形成、肿瘤转移过程中的细胞运动及血管生成（angiogenesis）；金属的热化学处理，包括渗碳（carburizing）和渗氮（nitriding）；以及固体推进剂的行为。

论文解读：约束反应性固体混合物的化学-热-力学耦合理论

一、研究背景与问题

在过去的几十年中，大量科学与工程应用中出现了一系列紧密耦合的物理化学过程——包括质量输运、力学变形、化学反应、热分解、相变及产气反应——这些过程共同主导着固体混合物的行为。这些应用涵盖机械生物学现象（如血管生成和肿瘤转移过程中的细胞运动）、胚胎发育、软骨演化、生物膜形成及血液凝固，也涉及冶金学中的渗碳和渗氮工艺、储能领域中的可充电电池运行，以及固体火箭推进剂的燃烧等。在这些异质应用中，质量传递往往伴随显著的体积变化，这可能深刻影响固体相的力学响应及耦合多物理场行为。然而，现有文献尚未对这些效应给予应有的深入考察。填补这一空白正是本研究的核心目标。

混合物理论（mixture theory）是处理多连续介质相互作用复杂系统的有力框架。现有理论可划分为连续介质混合物理论和多相理论两大框架。根据平衡定律的复杂度，混合物可分为四类：第I类包含N个质量平衡、单一动量平衡及混合物整体能量和熵平衡；第II类考虑N个质量和动量平衡但仅单一能量和熵平衡；第III类包含N个质量、动量和能量平衡但单一熵平衡；第IV类则为每个组分单独建立熵平衡。Larché和Cahn在1973年和1978年的开创性工作中发展了用于多组分固体热力学平衡的连续介质混合物理论，并将其与有限应变几何演化耦合，成为材料化学-力学领域的里程碑。然而，这些理论在处理反应性固体混合物时存在关键局限：参考构型的定义、体积分数在参考构型与当前构型之间的变换关系，以及化学反应驱动的相变与体积变化之间的耦合机制，均未得到严格处理。

二、研究内容与结论

本研究建立了一个有限应变、化学-力学连续介质公式，用于描述固体混合物随时间演化的非平衡过程。该框架建立在以下基本假设之上：所有固相共享同一速度场，因此不存在相对运动。通过假设单一运动控制整个物体的运动学，自然得到一个唯一的参考构型。这一视角将该框架与允许各组分独立运动的无约束混合物理论，以及仅在当前构型中共同运动的约束混合物理论区分开来。

研究揭示了一个与无反应混合物相比的根本性区别：尽管物体的运动是唯一定义的，但各相在参考构型中的体积分数通常不会在当前构型中保持不变。这一结果源于参考构型中所有相具有任意但相同的参考摩尔体积（referential molar volume）这一假设，使得参考构型成为抽象的理想化构造而非初始时刻的实际物理构型。研究还表明，由于参考摩尔体积是任意选择的量，各相对应的参考密度也必须被视为任意的。此外，经典的Larché-Cahn网络理论及其后续扩展——包括质量输运和化学反应——被证明是本研究提出的更一般框架的特例。

三、关键技术方法

研究人员主要采用了以下关键技术方法：一是基于乘法分解（multiplicative decomposition）的运动学框架，将变形梯度分解为热力学张量和化学-热膨胀张量的乘积；二是通过代表性体积单元（Representative Volume Element, RVE）视角将物体视为固体混合物，定义各相体积分数；三是建立拉格朗日形式的控制方程，基于反应动力学、质量和动量平衡定律以及输运、化学反应和力学耦合的热力学分析；四是将化学与温度耦合，假设温度变化适度且不触发相变，引入化学-热膨胀因子。

四、研究结果

2.1 相、物种、构型与体积分数：研究明确了相（phase）与物种（species）的经典热力学概念区分，以锡纳米颗粒电化学钠化过程中的微观结构演化为示例，展示了多相共存和物种迁移的物理图景。定义了参考构型和当前构型中的体积分数，并指出二者在化学反应过程中通常不相等。

2.2 中间构型：引入膨胀因子（swelling factor）概念，建立了参考构型与中间构型之间体积变化的关系式，推导出化学膨胀因子J^c的表达式，这是本研究的主要成果之一，扩展了文献中的经典表达式。浓度和体积分数通过J^c在构型之间变换，但变换规律不同。化学反应导致的体积变化源于各反应相摩尔体积的差异，尽管质量在所有构型中守恒，体积通常不守恒。网络模型被证明是本框架的特例，化学-温度耦合通过引入温度依赖的摩尔体积关系实现。

五、讨论与结论

本研究的讨论部分强调了参考构型定义的重要性及其对控制方程的影响。尽管公式在一般性框架下建立，但涵盖广泛的应用范围，包括转化型电极离子电池的充放电过程、机械生物学中的聚合驱动现象（如血凝块形成、肿瘤转移中的细胞运动和血管生成）、金属的热化学处理（渗碳和渗氮）以及固体推进剂的行为。研究结论指出，所提出的框架为反应性固体混合物的化学-热-力学耦合行为提供了严格的理论基础，其膨胀响应表达式区别于文献中的标准处理方法，对所得控制方程具有直接影响，是本研究的主要贡献之一。