4. 多物理场耦合框架

在锂离子电池的相场模拟中，单一物理场只能描述局部行为，而真实电极颗粒内部的演化几乎总是多场同时作用的结果。以三元正极颗粒为例：锂离子嵌入/脱出产生浓度梯度，浓度梯度通过扩散引起化学势变化，化学势与电势共同决定离子迁移驱动力；同时，晶格参数随锂含量变化产生本征应变，本征应变在颗粒内部形成应力场，应力又反过来改变化学势，影响扩散路径；此外，充放电过程中的焦耳热与反应热形成温度场，温度通过Arrhenius关系影响扩散系数与反应速率。四个场之间并非简单叠加，而是通过源项、系数和本构关系双向耦合。

4.1 四场耦合关系与物理机制

浓度场 c(r,t) 是相场模拟的核心变量，通常用Cahn-Hilliard方程描述其守恒演化。电势场 φ(r,t) 由Poisson方程与电荷守恒约束，在电极颗粒内部表现为离子迁移与电子传导的耦合。应力场 σ(r,t) 由力学平衡方程给出，其本征应变 ε₀(c) 直接依赖浓度。温度场 T(r,t) 由热传导方程控制，热源包括不可逆焦耳热和可逆熵变热。

耦合路径可以概括为：浓度→本征应变→应力；浓度→化学势→扩散通量；电势→迁移通量→浓度；温度→扩散系数/反应速率→浓度与电势；应力→化学势修正→浓度。其中应力对化学势的修正项为 −Ωσ_h，Ω为偏摩尔体积，σ_h为静水应力。这一项在硅负极等大体积变化材料中尤为关键，可导致扩散势垒改变甚至相分离路径偏移。

4.2 控制方程与弱形式推导

浓度场采用Cahn-Hilliard方程，其自由能泛函为：

F = ∫[f(c) + κ|∇c|² + Ωσ_h·c] dV

其中 f(c) 为均匀自由能密度，κ为梯度能系数，Ωσ_h·c 为应力贡献项。化学势 μ = δF/δc = f'(c) − κ∇²c + Ωσ_h。守恒方程 ∂c/∂t = ∇·(M∇μ) 中，M为迁移率，通常取 M = D·c(1−c)/RT。

电势场满足 ∇·(σ_e∇φ) = 0，其中 σ_e 为电子电导率，在电极颗粒内部与浓度相关。应力场满足 ∇·σ = 0，本构关系 σ = C:(ε − ε₀(c))，ε₀(c) = β(c − c_ref)I，β为膨胀系数。温度场满足 ρC_p∂T/∂t = ∇·(k∇T) + Q，Q = I²R + TΔS·∂c/∂t。

弱形式推导以浓度场为例。将Cahn-Hilliard方程两边乘以测试函数 η，在域Ω上积分并分部积分：

∫(∂c/∂t)·η dV = −∫M∇μ·∇η dV + ∫M(∂μ/∂n)·η dS

取自然边界条件 ∂μ/∂n = 0，得到弱形式：

∫(∂c/∂t)·η dV + ∫M∇μ·∇η dV = 0, ∀η ∈ H¹(Ω)

同理，电势场弱形式为 ∫σ_e∇φ·∇ψ dV = 0，应力场弱形式为 ∫σ:∇_sξ dV = ∫t·ξ dS，温度场弱形式为 ∫ρC_p(∂T/∂t)θ dV + ∫k∇T·∇θ dV = ∫Qθ dV。四个弱形式通过耦合项（μ中的Ωσ_h、σ_e中的c、Q中的∂c/∂t）联立，形成非线性方程组。

4.3 耦合求解策略与示例参数

实际求解时，通常采用交错迭代或全耦合Newton法。交错迭代将四个场依次求解，每步更新耦合项，适合弱耦合场景；全耦合Newton法将四个场变量同时求解，收敛快但雅可比矩阵规模大。以下给出一个教学示例参数集（属于合理假设，非实测数据）：

在COMSOL或MOOSE中实现时，可将四个弱形式分别写入四个物理场接口，通过耦合变量传递。例如在MOOSE中，浓度场用CahnHilliard内核，电势场用Poisson内核，应力场用TensorMechanics内核，温度场用HeatConduction内核，耦合项通过Material属性传递。求解器建议采用Newton法，相对容差1×10⁻⁶，绝对容差1×10⁻⁸。

4.4 结果判读与耦合强度判断

判断耦合是否显著，可比较各场特征时间尺度：扩散时间 τ_d = L²/D，热扩散时间 τ_t = L²/α，力学松弛时间 τ_m = L/c_s（c_s为声速）。若 τ_d ≫ τ_t，温度场可视为准静态；若 τ_m ≪ τ_d，力学平衡可瞬时满足。对于典型三元正极颗粒（L=5 μm，D=1×10⁻¹⁴ cm²/s），τ_d ≈ 2500 s，而τ_t ≈ 0.1 s，τ_m ≈ 1×10⁻⁹ s，因此温度场和应力场均可采用准静态近似，仅浓度场需动态求解。

观察结果时，重点关注：浓度分布是否出现相分离（Cahn-Hilliard特征），应力集中区域是否与浓度梯度最大处重合，电势降是否主要发生在颗粒表面，温度升高是否超过安全阈值（如60°C）。若应力修正项 −Ωσ_h 导致化学势变化超过RT（约2.5 kJ/mol），则应力耦合不可忽略。