氧化铜颗粒介质中Branly效应的多尺度电热建模

《Materials Today Communications》：Multiscale electro-thermal modeling of the Branly effect in oxidized copper granular media

　　

该研究探讨直流（DC）应力下覆盖本征半导体氧化物层（Cu2O）的铜颗粒介质的介电击穿现象。这种绝缘体向导体的转变被称为Branly效应，通过破坏本征氧化物扩散势垒引发火花等离子烧结（Spark Plasma Sintering, SPS）中的粉末固结。该现象通过结合离散元法（Discrete Element Method, DEM）、有限元法（Finite Element Method, FEM）与分析建模的多尺度方法进行剖析。在局部尺度，三维电热FEM模型验证了基于热失控的解析击穿准则。在宏观尺度，不同于聚焦于击穿后金属渗流（percolation）的模型，三维DEM框架整合了氧化物的热降解以主动驱动该转变。结果揭示了一种两阶段宏观击穿：由功率驱动的成丝（filamentation）稳定电场的“软”机制，继而是微小接触分数（<0.1%）的金属失效建立渗流路径的“硬”机制。这种数值方法深化了对击穿机制及诱导电流渗流的理解，为优化SPS早期活化与致密化阶段提供了物理见解。

论文发表在《Materials Today Communications》。研究背景方面，电流辅助烧结（Electric Current-Assisted Sintering, ECAS）技术尤其是放电等离子烧结（Spark Plasma Sintering, SPS）或场辅助烧结技术（Field Assisted Sintering Technology, FAST）已成为粉末冶金的关键技术，能够在较低温度与较短时间内实现全致密化并保留细微微观结构。不同于常规辐射加热，SPS依赖脉冲直流电（Direct Current, DC）流经工装与导电粉末压坯产生焦耳热（Joule heating），能量在作为高阻收缩区的颗粒间接触处局部内部生成，是高升温速率与高能效的主要驱动力。然而，松散粉末床中电流流动的启动受表面物理而非块体性质支配。可氧化金属粉末如雾化铜在空气中处理时覆盖本征氧化物层，铜的表面层为亚铜氧化物（cuprous oxide, Cu₂O），一种p型半导体，电导率依赖热激活与晶格缺陷，电阻率比金属核高一个数量级。因此初始压坯表现为电绝缘体，其活化与这些接触界面的局域热力学状态相关。允许电流流动的转变为跨越绝缘氧化物层的介电击穿，这种现象在现象学上被识别为Branly效应。颗粒物理研究表明电热机制主导Branly效应，电压饱和是电热耦合的特征，可通过Holm的“烧结（fritting）”理论解释，即介电膜击穿及后续局域熔化与电迁移导致的金属微桥形成。此外，电击穿不仅启用电流流动，还触发物理烧结，通过破坏与消除作为界面扩散势垒的本征氧化物层，电热转变通过解锁金属质量传输启动粉末固结与致密化。由于局域击穿决定烧结起始，理解其宏观动力学是控制整体SPS工艺的前提，这些局域击穿事件的空间分布与序列决定电流网络均匀性与后续温度场，空间异质击穿可导致局域热点、过早熔化与异常晶粒生长，损害最终烧结部件性能，因此理解该转变对控制SPS中氧化金属粉末微观结构发展是必要的。现有离散元法（Discrete Element Method, DEM）模型关注击穿后渗流并假设金属接触性质，依赖简化经验电接触定律，即便先进耦合热电模型也假设欧姆接触电阻受颈部几何与块体性质支配，从而忽略表面氧化物层的复杂物理。本研究引入首个基于物理的多尺度框架致力于电转变起始，填补亚微米氧化物层物理与三维颗粒堆积宏观介电击穿之间的空白，将颗粒响应锚定于本征氧化物的固态热力学，展示其热激活与结构降解如何驱动绝缘体向导体转变。研究人员开展了多尺度电热建模研究，结合三维有限元法（Finite Element Method, FEM）与三维离散元法（Discrete Element Method, DEM）及解析建模，在局部尺度验证氧化物层电热击穿物理，在宏观尺度模拟颗粒压坯行为，得出击穿为两阶段过程等结论，意义在于深化对击穿机制与电流渗流的理解，为优化SPS早期阶段提供物理见解。

关键技术方法方面，研究人员采用三维FEM在单颗粒双体体系上验证局部击穿物理与解析热失控准则，使用COMSOL Multiphysics构建含Cu₂O层的铜半球几何并赋予温度依赖材料属性进行准静态电热耦合求解；采用开源C++离散元库GranOO构建三维DEM宏观模型，生成匹配实验粒径分布与相对密度的颗粒堆积，耦合机械Hertz–Mindlin接触定律与基于图论的电路求解器及局域焦耳热与散热算法，引入热失衡判据区分稳定欧姆态、软击穿与硬击穿三种接触模式并更新电阻网络迭代计算。

研究结果部分，2 理论与解析建模及电热击穿解析模型小节中，研究人员通过计算毕渥数（Biot number, Bi）验证均匀温度假设，得出铜颗粒内部导热阻力远小于表面散热阻力、氧化物层与间隙介质集中热梯度、铜核为准均匀温度的结论；基于Arrhenius型电导率定律建立氧化物层焦耳产热与牛顿冷却耗散平衡，通过临界切点条件推导热失控临界温升与临界击穿电场解析表达式，指出击穿场强与接触半径成反比、与半导体电参数平方根成反比、与传热系数平方根成正比。3 多尺度数值实现小节中，3.1 接触尺度三维FEM部分研究人员通过FEM获取接触级空间场分布，验证电位降集中于氧化物层、电流收缩效应最大于接触周边、铜核为热沉的假设，并对比FEM与解析模型得出解析模型因忽略边缘场（fringing field）效应轻微高估击穿场强但可提供介电强度上限；3.2 宏观尺度DEM部分研究人员通过落雨法生成相对密度54.7%的圆柱颗粒堆积，在固定几何上运行耦合电热算法，将热失衡发散分为低能软击穿（指数损伤律受限至饱和比（Saturation Ratio, S_ratio））与高能硬击穿（设为金属收缩电阻约10^-3 Ω），更新电阻网络迭代。4 结果与讨论小节中，4.1 局部尺度FEM结果与解析模型验证部分研究人员通过FEM与解析对比确认击穿场强随接触半径增大而非线性衰减、随颗粒直径增大而增加的趋势，FEM因捕捉电流侧向扩展导致略早热失控；4.2 中尺度FEM与DEM比较部分研究人员在单双体与三颗粒链对比中得出DEM因一维圆柱约束低估击穿电流密度但与FEM在欧姆区至失稳点对齐，验证DEM电热耦合可靠性；4.3 宏观尺度DEM与实验部分研究人员采用分层校准策略约束预指数电导率与热敏参数γ及饱和比S_ratio、近绝热传热系数h与损伤动力学参数C_Power与β，对比实验电压扫描数据得出完整电热损伤模型复现初始欧姆区、软击穿非线性电阻下降与硬击穿突变，纯热模型仅Arrhenius加热则过早发散；通过追踪接触状态数量得出软击穿分两阶段，初期仅约12个接触沿局域通道活化且电场稳定，后期约8000个接触成丝分支使电阻分布式下降，硬击穿在临界场约50 kV/m时仅约8个接触发生金属失效形成连续渗流路径，占比<0.1%，为动态雪崩而非随机渗流；讨论部分研究人员指出DEM因一维接触忽略三维边缘场致稍低估击穿电流密度，理想球忽略表面粗糙度局域场增强与氧化物厚度随机性，刚性颗粒忽略低压下机械变形但与1 MPa实验条件匹配，直流稳态模型为脉冲SPS的热积累极限情形，排除瞬时介电击穿、微等离子体与非击穿期电迁移主导说，强调热失控为核心机制。

总结讨论部分，研究人员指出解析模型给出击穿场强与接触半径反比的理论依据，FEM验证铜核等温与产热集中于氧化物层，DEM揭示Branly效应为软击穿成丝与硬击穿微量接触金属失效的两阶段过程，软击穿受几何收缩限域维持电压饱和，硬击穿为渗流事件。对比实验存在数值微小偏差源于DEM一维简化与理想几何，但提供保守阈值与内部不可观测的渗流动力学。框架基于物理氧化物动力学而非经验定律，可预测初始参数如氧化层厚度影响，优化SPS活化阶段压力与电压速率以防过热。未来需集成弹塑接触力学以扩高压范围及瞬态脉冲热方程以近工业SPS条件。翻译研究结论部分为：本文提出了氧化铜颗粒介质中介电击穿的多尺度分析，架起了氧化物层局域电热固态物理与压坯宏观活化之间的桥梁。现有文献中DEM模型在击穿事件后模拟电流渗流，假设颗粒核为金属性质同时忽略薄氧化物层，为填补此空白，本框架整合了该本征氧化物的物理行为，展示其热激活与结构降解如何驱动绝缘体向导体转变。研究人员基于热失控理论解析推导了介电击穿的物理准则，该模型强调了击穿场强对接触半径的反比依赖，为理解SPS工艺中晶粒尺寸与初始压实压力的影响提供了理论基础。数值上，三维FEM模拟证实了这些解析假设的有效性，表明产热局限于阻性半导体层而铜核充当等温热沉。在宏观层面，耦合电热三维DEM模拟为金属粉末的Branly效应提供了物理解释。绝缘体向导体转变被描述为受不同物理机制支配的多阶段过程：起始于“软”击穿相，数千接触因局域焦耳热发展导电细丝，导致宏观电阻率分布式渐进下降；该成丝相受氧化物内导电路径几何收缩自限制，使颗粒介质在维持电压饱和的同时容纳增长电流；此先驱相后续为“硬”击穿，被识别为渗流事件，由微小接触分数（<0.1%活跃网络）的金属失效建立首条桥接电极的连续导电路径所驱动。与实验数据对比揭示了烧结早期阶段的见解，虽通过受宏观数据集与微观物理观测约束的多尺度校准复现了基线宏观击穿电压阈值，但对硬击穿起始电流密度的数值低估凸显了数值与物理因素：数值上DEM一维接触简化忽略三维边缘场效应，此几何约束高估局域电阻与温升从而加速热失控，在网络尺度加速接触顺序失效，使模型在低于三维物理介质的宏观电流密度下预测宏观硬渗流（连续导电路径形成）以提供保守阈值估计；物理上局域加热引起微凸体软化与早期颈长大（自适应机制）在宏观熔化前动态增大有效接触面积并降低收缩电阻，此现象强调了电活化与机械致密化的关联。尽管有这些简化，数值框架的预测能力在于解析不可及的内部渗流动力学，揭示软击穿网络的时空分支与支配转变的关键接触分数；通过显式建模氧化物层动力学而非依赖现象学击穿后接触定律，该框架为绝缘体向导体转变起始提供了基础性理解；一旦在参考粉末状态校准，模型可模拟初始参数变化如初始氧化物厚度的影响。这些能力为优化SPS工艺初始“活化”阶段提供了见解，特别是关于控制压力与电压扫描速率以确保电流均匀与防止局域过热的策略。该框架未来扩展旨在模拟施加机械压力的影响，集成弹塑性接触力学与氧化物层破裂将丰富DEM模型并将其应用范围拓宽至标准SPS高压条件模拟；此外，虽本研究捕获DC应力下击穿机制，工业SPS使用脉冲电流，集成瞬态热方程与脉冲电边界条件构成必要未来发展以探究脉冲频率与占空比如何在致密化最早阶段调制热积累、渗流动力学与微观结构均匀性。