1. 电池材料相场模拟导论

锂离子电池的性能退化、枝晶生长、相分离与裂纹演化，本质上都是介观尺度上多物理场耦合驱动的微结构演化问题。相场方法（Phase-Field Method）通过引入连续变化的序参量场，把尖锐界面问题转化为体相偏微分方程求解，天然适合描述这类"界面在演化、拓扑在变化"的过程。本章回答三个问题：相场在电池材料研究中处于什么位置、为什么必须耦合多物理场、以及本课程将沿着怎样的路线展开。

一、相场方法在电池材料中的定位

电池材料模拟通常按空间尺度分层：原子尺度用第一性原理（DFT）计算形成能、迁移势垒；微观尺度用分子动力学（MD）观察局域结构；介观尺度（约1 nm～10 μm）用相场描述颗粒、枝晶、孔隙、裂纹的形貌演化；宏观尺度用等效电路或P2D模型描述电芯行为。相场恰好填补了"原子尺度算得准但算不大、宏观模型算得大但看不见形貌"之间的空白。

相场的核心思想是：不再显式追踪界面位置，而是用一个在空间上连续过渡的序参量 φ 表示相态。例如 φ=1 表示固相、φ=0 表示液相，界面处 φ 在几个网格宽度内平滑过渡。系统的自由能泛函写为

F = ∫ [ f_local(φ, c) + (κ/2)|∇φ|² ] dV

其中 f_local 是局部自由能密度（决定相是否稳定），κ 是梯度能系数（决定界面厚度与界面能）。序参量随时间的演化由 Allen–Cahn 方程（非保守场）或 Cahn–Hilliard 方程（保守场）给出：

∂φ/∂t = -L ( δF/δφ )          (Allen–Cahn)
∂c/∂t = ∇·( M ∇( δF/δc ) )     (Cahn–Hilliard)

L 是界面迁移率，M 是原子迁移率。这两个方程是后续所有电池相场模型的骨架：枝晶生长改的是 L 与界面各向异性，相分离改的是 f_local 的形式，而浓度场 c 与电势场 φ_e 的耦合则通过电化学驱动力进入自由能或动力学系数。

二、为什么必须多物理场耦合

纯相场只能描述"形貌怎么变"，却回答不了"为什么这样变"。电池材料中的相变几乎都由外部场驱动：锂离子浓度梯度提供化学驱动力，电势分布决定反应速率，应力场抑制或促进相变，温度场改变扩散系数。因此实际模型至少需要耦合四类场：

• 浓度场 c：由 Cahn–Hilliard 或扩散方程演化，决定局部锂含量与相稳定性；
• 电势场 φ_e：由 Poisson 方程 ∇·(σ∇φ_e)=0 求解，给出电解液与电极内的电位分布；
• 应力/位移场 u：由力学平衡 ∇·σ=0 求解，本征应变来自浓度变化（Vegard 定律 ε₀=Ω(c−c₀)/3）；
• 温度场 T：由热传导方程求解，影响扩散系数与反应速率。

耦合的关键在于双向反馈：浓度变化产生本征应变，应力反过来改变化学势（μ=μ₀−Ωσ_h），从而改变相变驱动力；电势决定界面反应速率，反应速率又消耗锂离子改变浓度。若只做单向耦合，枝晶形貌会明显失真——这正是许多早期模型预测枝晶"长得太规则"的原因。

三、典型应用场景

相场在电池材料中的典型应用可归纳为四类，它们对应的物理过程与关键参数各不相同：

以锂枝晶为例：在电解液中施加过电位 η，锂离子在负极表面还原沉积。若界面能各向同性，沉积界面会保持平坦；一旦引入各向异性（如 κ 随晶向变化），尖端处曲率大、局部过电位高，沉积速率加快，形成正反馈，最终发展为枝晶。相场模型能定量给出"临界电流密度"——超过该值枝晶必然失稳，这与实验观察到的 Sand's time 行为一致。

四、课程路线图

本课程共30章，按"基础→单场→耦合→应用→验证"的逻辑推进。下图给出课程的整体知识结构：从相场数学基础出发，依次建立浓度场、电势场、应力场、温度场的求解能力，再通过双向耦合形成完整的多物理场模型，最终落到枝晶、相分离、裂纹、固态界面四类典型应用，并以实验对比与参数标定收尾。

具体而言，第2～5章建立相场数学与数值基础（自由能泛函、变分推导、离散格式、边界条件）；第6～12章分别求解浓度场、电势场、应力场与温度场，并完成两两耦合；第13～20章进入双向耦合与求解器实现，包括非线性迭代、时间步进、并行加速；第21～27章针对枝晶、相分离、裂纹、固态界面四类场景做完整案例；第28～30章讨论参数标定、实验对比与模型可信度评估。

学习本课程不需要预先掌握相场理论，但需要具备偏微分方程、有限差分/有限元基础，以及 Python 或 C++ 编程能力。课程中所有示例参数均为教学假设值，用于说明模型行为，不代表特定材料的实测数据；实际研究时应以文献或实验标定值为准。