图1 实验平台
Fig.1 Experimental platform
摘要 聚酰亚胺（PI）薄膜被广泛应用在固态变压器（SST）的绝缘设计中，而沿面放电是导致SST绝缘失效的重要原因之一。该文使用COMSOL Multiphysics软件为针电极实验装置建立了二维自洽等离子体仿真模型，通过耦合等离子体化学反应、粒子输运方程和泊松方程来分析高频电应力下的沿面放电发展过程；并建立了一个简化模型，用于验证等离子体模型中的空间电荷密度。实验结果表明，沿面放电会使电介质表面碳化并增加PI表面电导率。当电压频率为50Hz时，表面放电呈树枝状，放电出现在该阶段正半周期的上升和下降沿，最大放电幅度为0.013V，放电次数为2 180。在高频下，随着绝缘老化，放电幅度从0.009V增加到0.015V，最后阶段记录的放电总数为3 610，高频电应力下的沿面放电造成的损伤比工频更为显著。通过仿真模拟，该文获取了空间电荷、电子密度、反应速率和电子温度等微观量。高频电应力下PI薄膜沿面放电轨迹近似线性，在176×10-6s时电子密度达到1.1×1019m-3，随着电场强度的增加，电子密度和电子温度分布呈现上升趋势；正负离子密度分布范围与电子密度分布有很大的相似性。对电子密度和反应速率的分析表明，背景气体中的化学反应分析是正确的。
关键词：聚酰亚胺 沿面放电 空间电荷 电子密度 电子温度 反应速率
随着全球清洁电能的生产越来越多，未来几年中，可再生能源将在现有电力基础设施中使用得更加普遍。高频电力变压器（High-Frequency Power Transformer, HFPT），也称为固态变压器（Solid-State Transformer, SST），有望在未来的智能电网中发挥重要作用。与传统变压器相比，SST可以提供更优的电能质量，其在输出电压稳定性、无功补偿、电压调节和功率因数校正等方面具有可控性优势[1-5]。为保证设备的紧凑性和高效率，工程中通常采用气固干式绝缘系统[6]。聚酰亚胺（Polyimide, PI）因其优异的耐热性、化学和物理稳定性及电气绝缘性而被广泛用作干式电机（包括变压器）的绝缘膜[7-8]。在高频电应力下，普通PI薄膜的绝缘能力大大降低，使得高压大容量HFPT的发展受到限制[9]。目前许多学者对于PI薄膜的绝缘研究都集中在其放电耐久性上，而对于热稳定性及其他理想特性的研究较少[10-11]。由于高频变压器的高功率密度、高效率和快速动态性，会出现一些相应的现象，例如施加电压脉冲的高重复频率[1-4]、快速极性反转（dV/dt）[12-13]及相当高的热应力[2,5]。高频电力变压器绝缘系统在极高的频率和温度下运行，具有很大的安全风险[14]。
沿面放电会导致绝缘的劣化和过早失效[15]，因此研究沿面放电对于构建SST的高频绝缘结构至关重要。介电材料易在电气设备中引起放电，电子崩和辉光放电可能发生在绝缘体周围的应力区，这些放电最终可能导致绝缘体表面发生闪络或电气击穿。早期对于沿面放电的研究主要集中在对闪络电压的测量上[16]，目前研究人员对沿面放电和闪络已使用高速相机进行了更深入的研究[17-18]。
许多学者研究了PI击穿后的表面电荷效应[19-23]，文献[24]进行了表面介质阻挡放电实验，发现了流注速度和流注电场间的规律。然而，由于需要具有微米空间分辨率和纳秒时间分辨率的非侵入式诊断方法，对大气压下等离子体表面的相互作用进行微观分析极其困难[25]。文献[26]验证了其正弦交流电压下的沿面放电特征，并探索了影响PI绝缘沿面放电的因素。
为了模拟沿面放电特性，一些研究人员采用了分形几何和蒙特卡罗方法[21-24]。随着计算机技术的进步，人们通过表面传播的动力学对流注进行了许多数值模拟，揭示了早期实验未知的微观机制。大多数值模拟都采用流体动力学模型，文献[27]使用脉冲电应力下的空气-聚酰亚胺的非平衡等离子体模型来模拟沿面放电，发现在重复和单脉冲应力下，从针电极到接地电极间会产生沿面放电。文献[28]使用电流体动力学方法，模拟了具有针电极配置的固液绝缘系统的沿面放电，并确定了放电传播速度。文献[29]使用流体动力学模型模拟了介电势垒内的电荷传输过程，并探究了固体介电势垒对针电极和板电极间空气放电形成的影响。在模拟放电时，研究人员将流体动力学理论与等离子体反应相结合，但对于高频电场下的沿面放电机理及其微观研究还没有准确的描述。
本文创建了二维等离子体自洽仿真模型，该简化模型采用自适应网格细化进行简单高效的仿真。将针电极结构视为一种几何结构[30]，流注沿着电介质表面由该几何结构进入气体中。本文通过实验和计算机模拟相结合的方法揭示了聚酰亚胺在高频电应力下的表面闪络机理，为SST绝缘结构设计提供了理论基础和工程建议。
高频沿面闪络实验平台如图1所示。高频高压电源可输出频率为5～50kHz、峰值为0～15kV的持续可调的正弦波电压。本文中的电极结构为针电极，其长度为28mm，直径为4mm，曲率半径为50mm，距离环氧板10mm，与样品夹角为15°。使用带宽300kHz～100MHz的ETS-93686高频脉冲电流传感器测量放电信号，使用横河DL6154四通道示波器收集外加电压和放电信号数据，实时采集的数据通过USB串行总线和LabVIEW的数据采集存储系统发送到计算机进行存储，并使用针电极上方的高速摄像机对沿面放电的发展过程进行连续拍照和记录。
图1 实验平台
Fig.1 Experimental platform
在室温下，施加30kHz的频率，使用200V/s的短期快速升压方法，直到发生沿面闪络。保持电压为10kV、频率为30kHz不变，采用恒定电压法记录五组放电起始电压和沿PI表面的闪络电压。在该频率和电压下，沿面放电可以继续发展而速度不会太快，确保可以收集到足够多的数据进行分析。放电信号通过数据采集系统进行实时存储，摄像头捕捉样品上方沿面放电的发展过程。
本文使用COMSOL Multiphysics仿真软件来揭示时间-空间参数，通过将流体模型与动力学模型[26]耦合，构建了流体-动力学数值模型。此外，该模型在逻辑上简化了复杂的化学反应，以解决数值振荡问题。本文研究了大气压下的气固绝缘系统的沿面放电机制，采用自洽等离子体模型和简化模型解决实验和模拟的有效性问题。
等离子体模型在对流-扩散近似中求解与泊松方程自洽耦合的电子、离子连续性和动量方程。用自洽等离子体模型对沿面放电进行建模，计算量较大，其中给定电场中的针状电极负责维持放电。等离子体模型及其网格划分如图2所示。
图2 等离子体模型及其网格划分
Fig.2 Plasma model and its meshing for solution
简化模型的建立基于电场和流经载流子上的电流守恒。该模型是非自洽的，以带电载流子流经的维持电流一定为前提，其几何形状与图2所示的等离子体模型一致。
对于电子密度的计算，使用对流扩散方程，即
（1）
式中，Re为粒子反应决定的电子源项；μe和De分别为电子迁移率和扩散率；ne为电子密度；E为电场强度。使用式（2）从电子迁移率计算电子扩散率，即
（2）
式中，Te为电子温度；电子迁移率μe=3.74´ 1024(1021E)-0.22 (V·m·s)-1。等离子体化学中使用速率系数确定上述等式中的电子源项。假设反应 M（M=M1+M2）有助于电子密度增长或衰减，且发生了P（P=P1+P2）次电子-中性粒子的非弹性碰撞。电子源项在速率系数的情况下给出，即
（3）
（4）
式中，xi和xj分别为参与反应i和j的目标物质的摩尔分数；ki和αj分别为反应i和j的速率系数和汤森系数；Nn为中性粒子的数量密度；
和
分别为反应i和反应j的能量损失；|Γe|为电子密度通量的模；
为由于P次非弹性碰撞引起的能量损失或增益。对于非电子物质，使用式（5）求解每种物质的质量分数，即
（5）
式中，wk为粒子k的质量分数，k= 1, 2,
, N；ρ为混合物密度；u为平均流体速度；Rk为粒子k的源项；Sp为光致电离；jk为扩散通量。仿真计算时，离子迁移系数μi设置为6×1021(V·m·s)-1。静电场的计算公式为
（6）
（7）
式中，ρv为空间电荷密度；
为真空介电常数；
为相对介电常数；Zk为粒子k携带的电荷数；nk为粒子k的密度；q为基本电荷量；V为电动势。
简化模型的沿面放电通过电荷守恒和泊松方程来求解电荷载体的传输过程，即
（8）
（9）
（10）
式中，z为电荷数；m为空间电荷的迁移率；rs为空间电荷数密度；J为电流密度；S为电流源项。通过求解此方程组可以得到传输方程
（11）
当电流源项被忽略时，该方程采取式（12）所示的常用形式来模拟空间电荷密度，即
（12）
电子密度通量的边界条件为
（13）
式中，γ为二次发射系数，是由于二次发射效应而产生的电子增益；re为反射系数，在反射中被设定为0；ve,th为电子的热速度；
为热发射通量；
为二次发射通量。关于气隙放电模拟，设定阴极二次电子发射系数为0.004，PI表面的二次电子发射系数为0.05。设定气压为一个标准大气压（1.013×105Pa），环境温度为20℃。正负离子的初始密度为1017m-3，电子的初始密度为1010m-3。PI的相对介电常数为3.7[27]。由于表面反应，离子未越过壁垒。
（14）
（15）
式中，
为表面反应的速率常数；下标k为重粒子的类型，k = 1, 2, , K；Γtot为所有粒子的表面总浓度；γf为正向粘附系数；R为气体常数；T为气体温度；σk为势长；vk为化学计量矩阵；Mk、ck、µk和Zk分别为粒子k的摩尔质量、物质的量浓度、迁移率和电荷数；[8RT/(πMk)]0.5为粒子k的热运动速度。表面电荷积累边界条件为[31]
（16）
式中，D1和D2分别为介质绝缘表面两侧的电位移矢量；σs为表面电荷密度，可通过积分流动方程计算，即
（17）
式中，Ji和Je分别为总离子流密度和总电子流密度，它们在壁面的正态分量为
和
。
泊松方程使用针电极上的电场强度E的法向分量作为边界条件。
（18）
式中，
为击穿电场强度。
针电极和PI表面的电场必须足够接近真实的电场，才可以产生可预测的物理结果。本文采用皮克定律，即
（19）
式中，d为归一化为1.013×105Pa（一个标准大气压）和293.15K时的气体数密度；
为电晕电极的半径。使用二维等离子体仿真模型的非均匀网格划分如图2所示，网格被划分为221 287个单元。
在气-固绝缘系统中的等离子体化学反应较为复杂，对关键的激发态进行详细分析可能会包含数百个反应[32-33]。本文所建立的模型利用一组简化的反应，准确地描述了在气体背景下带电物质的形成和破坏过程。
该过程中所考虑的碰撞和化学反应见表1，其中氮和氧没有作为流体方程中的详细化学反应被单独考虑。背景气体采用通用物质A表示，A可以被电离，产生正离子p；也可以附着电子，产生负离子n。表1中，kf为正向速率系数，
为反应损失的能量。
表1 碰撞和反应
Tab.1 Collisions and reactions
反应公式类型/eVkf 1e+A→p+2e电离15— 2e+A→n附着—— 3e+2A→n+A附着—— 4e+p→A反应—5×10-14 5n+p→2A反应—5×10-12
电离系数、附着系数和三体附着的速率常数描述了体积中电子的形成和破坏过程。模型中在针电极边界和PI表面设置二次发射系数分别为0.004和0.05，使用80%的氮气和20%的氧气的混合物，并应用玻耳兹曼求解器计算了电子-离子和离子-离子重组的速率系数[34]，得出的平均电子能量和简化电场之间的关系如图3所示。
图3 平均电子能量和简化电场之间的关系
Fig.3 Relationship between average electron energy and simplified electric field
等离子体中的主要带电物质是离子，因此初始条件中正负离子的密度相等（均为1017m-3），电子的密度很小（1010m-3）。可以认为，当离子撞击壁垒时，它们将恢复为中性原子[35-36]。
为了了解沿面放电的发展过程，必须分析电子密度和每个阶段的发射数量。由于针电极、空气和绝缘表面三点交界处的电场强度较大，会使空气发生电离并产生有效的初始电子。随着电子崩的发展和电离程度的提高，沿面放电进入流注发展阶段。大量电荷聚集于电子崩的头部，使得崩头处的电场增强。流注头部的电荷不断得到补充，流注继续向前移动。当等离子体区域延伸到接地电极时，通过电极形成一个高导电性的路径，在表面上引起火花或电弧。
仿真模拟、高速相机所拍摄的沿面放电发展过程及沿面放电实验的局部放电相位分布（Phase Resolved Partial Discharge, PRPD）谱图如图4所示，数值模拟和实验获得的沿面放电统计数据见表2。
图4 沿面放电发展阶段
Fig.4 Surface discharge development stages
表2 沿面放电数据
Tab.2 Surface discharge data
放电阶段持续时间/μs电子密度/m-3放电次数最大放电幅值/V 11514.1×10132 5300.009 21673.5×10153 0400.012 31723.9×10173 3500.013 41761.1×10193 6100.015
根据PRPD统计数据和模拟结果，沿面放电具有明显的线性特征。脉冲放电电流在正弦波峰值电压和过零点附近较强。此外，施加电压的正半周期的放电幅度比负半周期更大。由于峰值附近的瞬时电压相对较高，放电区域的总电场强度非常大。由于空间电荷的滞后效应，当外加电压的极性改变时，空间电荷所形成的电场在电压的过零点不能迅速改变，导致整体电场强度在短时间内较大。在针电极结构下，材料表面的切向电场分布相对集中。本研究将沿面闪络发展分为四个阶段进行分析。
第一阶段，沿面放电初期，针电极附近的强电场会导致电晕放电，空气电离产生大量带电粒子。这些高能粒子会不断轰击样品表面，针尖附近出现白点状放电区域。在此阶段，电子密度从最初的1010m-3增加到4.1×1013m-3，放电次数达到2 530次。
第二阶段，随着加压时间的延长，放电通道逐渐向接地电极发展，白点状放电区域面积逐渐扩大。在放电过程中，等离子体通道随着放电的发展延伸到接地电极。白点区域的半径也随之扩大，并呈现出棒状。此阶段电子密度在167μs时达到3.5×1015m-3，放电次数达到3 040次。
第三阶段，在放电的后期阶段，白点区域呈鱼尾状不断向接地极移动，颜色继续加深。当等离子体通道和接地电极之间的距离减小到一定程度时，表面闪络会瞬间发生。带电粒子继续侵蚀PI薄膜的绝缘表面。电子密度在172μs时达到3.9×1017m-3，放电次数达到3 350次。
第四阶段，在高频电压下绝缘劣化的后期，放电迅速从等离子体放电通道转换为沿表面的闪络，对绝缘材料产生明显的破坏，甚至在材料表面发生燃烧。试样表面的等离子体通道逐渐向接地电极扩展，最终穿透两个电极并突破表面，产生闪络。此阶段电子密度在176μs时达到1.1×1019m-3，放电次数达到3 610次。
高频沿面放电通道与工频不同，两者都会在针电极到板电极间形成明亮的电弧，但形状不同，实验结果如图5所示。由图5a可知，工频下的沿面放电呈树枝状。在电压频率为50Hz时，最大放电幅值为0.013V，放电仅出现在正半周期的上升沿和下降沿，放电次数为2 180，在负半周期内没有出现放电信号。由图5b可以看出，频率增加至30kHz时，高频下PI沿面放电轨迹近似为线性。从第一阶段到最
图5 工频与高频电压下的沿面放电
Fig.5 Spark discharge under power frequency voltage and high-frequency voltage
后阶段，放电幅值从0.009V增加到0.015V，最后阶段记录的放电总数为3 610。高频电压下，正半周期的下降沿没有放电现象，而在负半周期的上升沿和下降沿都存在放电。
高频与工频损伤的PI材料如图6所示。在工频电应力作用下，由于火花放电的多次冲击，绝缘材料表面受到破坏。尽管如此，其损伤仍远小于高频电应力造成的损伤。
图6 高频与工频损伤的PI材料
Fig.6 PI material damaged by high-frequency and power frequency
对于沿面放电的动态演变，研究其每个阶段的电子密度、电场强度和电子温度是非常必要的。本文在10kV和30kHz的高频正弦电压下分析了针电极结构下PI薄膜的沿面放电特性。放电区域内的电场和电子密度的分布规律可以清楚地描述沿面放电的发展过程，这对于其机理研究至关重要。沿面放电的动态演变过程如图7所示，电子密度、电场强度和电子温度之间具有极强的关联性。随着电晕放电强度的增加，电子密度和温度随着电场强度的增加而增加。在放电过程中，针电极附近的电场强度最大，其中针电极尖端处的放电起始场强约为3MV/m。在强电场下，针电极上会发生许多碰撞电离，产生大量的电子。
如图7所示，电子密度呈现出稳健的线性化特征。在整个沿面放电过程中，电子密度保持在1013～1019m-3数量级。最大的电子密度发生于电场强度最大的地方。在沿面放电过程中，电子温度保持在1～10eV。
电子在电场中吸收能量，并通过与离子和中性粒子的碰撞传递能量。在沿面放电的早期阶段，针电极上的电子能量较弱。电子温度是其能量的量度，并被认为是沿面放电中能量传递的基本媒介。电场分布、电子密度和温度对于研究沿面放电过程和能量传输速率至关重要。
图7 沿面放电动态演变
Fig.7 Dynamic evolution of surface discharge
使用自洽等离子体模型研究电子、正负离子密度分布。在沿面放电过程中，电子密度分布对放电机制的理解具有重要意义。图8显示了从针电极到板电极之间的电子密度分布与时间尺度上的电子密度分布。在放电的初始阶段，针电极附近的高场强使空气发生电离，产生大量电子、负离子和正离子。空间电荷强烈扭曲了外部电场，电场强度最大的区域向前发展，这促使绝缘表面上的电子密度分布逐渐延伸到接地电极。由于针电极结构特点，最大电场强度存在于针电极附近，空间电荷进一步增加了针电极附近的电场强度。流注头部电子密度高的原因是头部前方的电场强度持续增加，导致发生强烈电离，电子密度增加。
沿绝缘表面的正负离子密度分布分别如图9、图10所示。因为负离子主要是通过电子与中性分子的相互作用形成的，所以负离子密度分布规律与电子密度差别不大。在放电最后阶段，正离子密度显著高于负离子密度，这是因为沿面放电中的冲击电离会产生许多正离子，且在强电场下，由于电子的快速移动，中性粒子难以捕捉电子并形成负离子。随着沿面放电的发展，绝缘材料表面的电子密度、正离子密度和负离子密度在整个放电区域持续增加。电场强度越高，电离越强烈，电离产生的电子也越多，所以最大电子密度值总是出现在电场强度的最大值附近。
图8 电子密度分布
Fig.8 Electron density distribution
图9 正离子分布
Fig.9 Positive ion distribution
图10 负离子分布
Fig.10 Negative ion distribution
使用简化模型与自洽等离子体模型研究空间电荷分布。对于自洽等离子体模型，取与绝缘介质表面平行，并位于表面上方0.2mm处的截面线，得到放电区域的空间电荷密度分布如图11所示。在放电发展的各个阶段，空间电荷都主要分布在放电通道的前端。碰撞电离产生的电子运动速度较快，在外部电场的作用下，以较高的速度向板电极移动；正电荷运动速度较慢，在外部电场的作用下，缓慢地向针电极移动。此外，针电极尖端附近的强电场导致许多正离子被电离，空间电荷呈现正极性，流注产生前电离产生的正离子也会在电场的作用下移动到针电极上。在第一阶段，针尖附近的正电荷进一步增加，这主要是因为电场增强，电离作用更加强烈。此后，针电极附近的正电荷不再增加，反而开始减少。随着流注的发展，它离针电极越来越远，这种影响也越来越弱。第一阶段后，针尖附近的空间电荷将逐渐减少。在最后阶段，离针电极较远处的空间电荷密度值明显增加，并且当等离子体通道穿透两极时，绝缘表面附近的空间电荷密度达到最大。
图11 自洽等离子体模型的空间电荷分布
Fig.11 Space charge distribution of plasma model
简化模型得到的空间电荷密度空间分布如图12所示。简化模型中电极处的电场计算采用皮克定律，不包括源项。由图11和图12可知，简化模型与自洽等离子体模型的空间电荷密度分布高度相似，但皮克定律计算所得的电场强度较小，这是两种模型所得空间电荷分布有差异的主要原因。
图12 简化模型的空间电荷分布
Fig.12 Space charge distribution of simplified model
反应速率分析可以对模拟结果进行验证，在数值模拟工作中具有重要意义。图13显示了使用自洽等离子体模型分析所得的所有带电物质反应速率的逐级总和的空间分布。由图13可知，反应速率分布与电子密度分布基本相同，这表明背景气体中的化学反应分析是正确的。反应速率受电场的影响，随着沿面放电的发展，材料绝缘表面的反应速率在整个放电区域内继续增加。最后阶段放电期间所有带电物质反应速率的总和约为1 200mol/(m3·s)。图14和图15显示了使用自洽等离子体模型确定的正离子和负离子的反应速率。通过电子碰撞电离产生正离子，电离区和传输区的电子附着产生负离子。正离子和负离子的速率表达与第4.1节的结果表现出良好的一致性。
图13 反应速率总和
Fig.13 Sum of all reactions rate
图14 正离子反应速率
Fig.14 Positive ion reaction rate
图15 负离子反应速率
Fig.15 Negative ion reaction rate
本文介绍了一种二维流体动力学等离子体模拟模型，用于研究高频正弦电压引起的聚酰亚胺（PI）沿面放电，并通过寻找实验上无法获得的微观参数来揭示沿面放电的演变过程，得到如下结论：
1）高频电应力下沿面放电发展过程可以划分为四个阶段。在每个阶段，电子密度都表现出线性化特征，在176×10-6s时电子密度达到1.1×1019m-3，放电次数达到3 610次。本研究通过PRPD分析，发现在正半周期内放电次数和放电幅度较大。工频电应力下无火花放电，高频电应力下的表面闪络损伤比工频更严重。
2）在高频电应力下，PI薄膜的沿面放电非常严重。电子密度与温度随着外加电压幅度和频率的升高而升高，正负离子密度与电子密度分布趋势相似。
3）可以使用简化模型获得空间电荷密度，所得结果与自洽等离子体模型所得空间电荷密度基本一致。
4）随着沿面放电过程的发展，反应速率受电场的影响逐渐增加，在放电最后阶段，反应速率总和可达1 200mol/(m3·s)，其分布趋势与电子密度分布相似，表明了背景气体中化学反应分析的正确性。
5）使用简化模型进行分析较为灵活方便，可用于工程应用，该方法的发展前景较好，可以为表面闪络的研究开辟新的途径。
参考文献
[1] 高范强, 李子欣, 李耀华, 等. 面向交直流混合配电应用的10kV-3MV·A四端口电力电子变压器[J]. 电工技术学报, 2021, 36(16): 3331-3341. Gao Fanqiang, Li Zixin, Li Yaohua, et al. 10kV-3MV·A four-port power electronic transformer for AC-DC hybrid power distribution applications[J]. Transactions of China Electrotechnical Society, 2021, 36(16): 3331-3341.
[2] She Xu, Huang A Q, Burgos R. Review of solid-state transformer technologies and their application in power distribution systems[J]. IEEE Journal of Emerging and Selected Topics in Power Electronics, 2013, 1(3): 186-198.
[3] 涂春鸣, 黄红, 兰征, 等. 微电网中电力电子变压器与储能的协调控制策略[J]. 电工技术学报, 2019, 34(12): 2627-2636. Tu Chunming, Huang Hong, Lan Zheng, et al. Coordinated control strategy of power electronic transformer and energy storage in microgrid[J]. Transactions of China Electrotechnical Society, 2019, 34(12): 2627-2636.
[4] Falcones S, Ayyanar R, Mao Xiaolin. A DC-DC multiport-converter-based solid-state transformer integrating distributed generation and storage[J]. IEEE Transactions on Power Electronics, 2013, 28(5): 2192-2203.
[5] Leibl M, Ortiz G, Kolar J W. Design and experimental analysis of a medium-frequency transformer for solid-state transformer applications[J]. IEEE Journal of Emerging and Selected Topics in Power Electronics, 2017, 5(1): 110-123.
[6] Qin Hengsi, Kimball J W. Solid-state transformer architecture using AC-AC dual-active-bridge converter[J]. IEEE Transactions on Industrial Electronics, 2013, 60(9): 3720-3730.
[7] Chen Hao, Prasai A, Moghe R, et al. A 50-kVA three-phase solid-state transformer based on the minimal topology: dyna-C[J]. IEEE Transactions on Power Electronics, 2016, 31(12): 8126-8137.
[8] 周力任, 吴广宁, 高波, 等. 聚酰亚胺薄膜中电荷输运机理和空间电荷特性[J]. 电工技术学报, 2009, 24(12): 6-11. Zhou Liren, Wu Guangning, Gao Bo, et al. Charge transport mechanism and space charge characteristic in polyimide film[J]. Transactions of China Electrotechnical Society, 2009, 24(12): 6-11.
[9] Liu Tao, Li Qingmin, Huang Xuwei, et al. Partial discharge behavior and ground insulation life expectancy under different voltage frequencies[J]. IEEE Transactions on Dielectrics and Electrical Insulation, 2018, 25(2): 603-613.
[10] Montanari G C, Seri P, Stöck M, et al. Investigating electrical performance of polyimide insulation tapes for automotive application[C]//2018 IEEE Electrical Insulation Conference, San Antonio, TX, USA, 2018: 369-372.
[11] Zhang Peihong, Zhang Weiguo, Liu Yan, et al. Study on corona-resistance of polyimide-nano inorganic composites[C]//Proceedings of the 7th International Conference on Properties and Applications of Dielectric Materials (Cat. No.03CH37417), Nagoya, Japan, 2003: 1138-1141.
[12] Kolar J, Ortiz G. Solid-state-transformers: key components of future traction and smart grid systems[C]//the 2014 International Power Electronics Conference, Hiroshima, Japan, 2014: 22-35.
[13] Huber J E, Kolar J W. Optimum number of cascaded cells for high-power medium-voltage AC-DC converters[J]. IEEE Journal of Emerging and Selected Topics in Power Electronics, 2017, 5(1): 213-232.
[14] Guillod T, Färber R, Krismer F, et al. Computation and analysis of dielectric losses in MV power electronic converter insulation[C]//2016 IEEE Energy Conversion Congress and Exposition, Milwaukee, WI, USA, 2016: 1-8.
[15] 张冠军, 赵文彬, 郑楠, 等. 真空中固体绝缘沿面闪络现象的研究进展[J]. 高电压技术, 2007, 33(7): 30-35. Zhang Guanjun, Zhao Wenbin, Zheng Nan, et al. Research progress on surface flashover phenomena across solid insulation in vacuum[J]. High Voltage Engineering, 2007, 33(7): 30-35.
[16] Sudarshan T S, Dougal R A. Mechanisms of surface flashover along solid dielectrics in compressed gases: a review[J]. IEEE Transactions on Electrical Insulation, 1986, EI-21(5): 727-746.
[17] Trienekens D J M, Nijdam S, Ebert U. Stroboscopic images of streamers through air and over dielectric surfaces[J]. IEEE Transactions on Plasma Science, 2014, 42(10): 2400-2401.
[18] Sobota A, van Veldhuizen E M, Stoffels W W. Discharge ignition near a dielectric[J]. IEEE Transactions on Plasma Science, 2008, 36(4): 912-913.
[19] Chen Geng, Tu Youping, Wu Shaocong, et al. Intrinsic hetero-polar surface charge phenomenon in environmental friendly C3F7CN/CO2 gas mixture[J]. Journal of Physics D: Applied Physics, 2020, 53(18): 18LT03.
[20] Kumara S, Alam S, Hoque I R, et al. DC flashover characteristics of a polymeric insulator in presence of surface charges[J]. IEEE Transactions on Dielectrics and Electrical Insulation, 2012, 19(3): 1084-1090.
[21] Li Chuanyang, Hu Jun, Lin Chuanjie, et al. Surface charge migration and dc surface flashover of surface-modified epoxy-based insulators[J]. Journal of Physics D: Applied Physics, 2017, 50(6): 065301.
[22] Shao Tao, Kong Fei, Lin Haofan, et al. Correlation between surface charge and DC surface flashover of plasma treated epoxy resin[J]. IEEE Transactions on Dielectrics and Electrical Insulation, 2018, 25(4): 1267-1274.
[23] Liu Yanqin, Wu Guangning, Gao Guoqiang, et al. Surface charge accumulation behavior and its influence on surface flashover performance of Al2O3-filled epoxy resin insulators under DC voltages[J]. Plasma Science and Technology, 2019, 21(5): 055501.
[24] Huang Bangdou, Zhang Cheng, Qiu Jintao, et al. The dynamics of discharge propagation and X-ray generation in nanosecond pulsed fast ionisation wave in 5 mbar nitrogen[J]. Plasma Sources Science and Technology, 2019, 28(9): 095001.
[25] Huang Bangdou, Zhang Cheng, Adamovich I, et al. Surface ionization wave propagation in the nanosecond pulsed surface dielectric barrier discharge: the influence of dielectric material and pulse repetition rate[J]. Plasma Sources Science and Technology, 2020, 29(4): 044001.
[26] Ayubi B I, Zhang Li, Xu Huangkuan, et al. Surface discharge characteristics and influencing factors of polyimide insulations under high-frequency sinusoidal voltages[J]. IET Generation, Transmission & Distribution, 2022, 16(5): 1000-1012.
[27] Dong Guojing, Li Qingmin, Liu Tao, et al. Finite-element analysis for surface discharge on polyimide insulation in air at atmospheric pressure under pulsed electrical stress[J]. High Voltage, 2020, 5(2): 166-175.
[28] Lee H Y, Kang I M, Lee S H. Fully coupled finite-element analysis for surface discharge on solid insulation in dielectric liquid with experimental validation[J]. IEEE Transactions on Magnetics, 2016, 52(3): 1-4.
[29] Singh S, Serdyuk Y V, Summer R. Streamer propagation in hybrid gas-solid insulation[C]//2015 IEEE Conference on Electrical Insulation and Dielectric Phenomena, Ann Arbor, MI, USA, 2015: 387-390.
[30] Firuzi K, Vakilian M, Phung B T, et al. Transformer components impact on compatibility of measured PDs: comparison of IEC60270 and RF methods[J]. High Voltage, 2019, 4(1): 33-40.
[31] 王泽忠, 全玉生, 卢斌先. 工程电磁场[M]. 2版. 北京: 清华大学出版社, 2011.
[32] Shibkov V, Konstantinovskij R. Kinetic model of ignition of hydrogen-oxygen mixture under conditions of non-equilibrium plasma of the gas discharge[C]// 43rd AIAA Aerospace Sciences Meeting and Exhibit, Reno, NV, 2005: 1-10.
[33] Tsai J H, Hsu C M, Hsu C C. Numerical simulation of downstream kinetics of an atmospheric pressure nitrogen plasma jet using laminar, modified laminar, and turbulent models[J]. Plasma Chemistry and Plasma Processing, 2013, 33(6): 1121-1135.
[34] Tejero-del-Caz A, Guerra V, Gonçalves D, et al. The LisbOn KInetics Boltzmann solver[J]. Plasma Sources Science and Technology, 2019, 28(4): 043001.
[35] Du Boxue, Liang Hucheng, Li Jin, et al. Electrical field distribution along SF6/N2 filled DC-GIS/GIL epoxy spacer[J]. IEEE Transactions on Dielectrics and Electrical Insulation, 2018, 25(4): 1202-1210.
[36] Sima Wenxia, Liu Chunxiang, Yang Ming, et al. Study of a nonequilibrium plasma model of surface discharge and the influencing factors[J]. IEEE Transactions on Plasma Science, 2017, 45(6): 906-912.
Abstract In the insulation design of solid-state transformers (SSTs), polyimide (PI) film has been widely used. The most significant reason for the insulation failure of SSTs is surface discharge. This paper proposes a two-dimensional fluid-kinetic plasma simulation model for a needle plate electrode experimental setup using COMSOL Multiphysics software to address this issue. We analyzed the surface discharge development process under high-frequency electric stress by coupling plasma chemical reaction, particle transport equation, and Poisson equation. We developed a simplified model that validates the space charge density from the fluid-kinetic plasma model.
Firstly, at room temperature, applied a frequency of 30kHz and 10kV with a 200V/s short-term rapid boosting method. Record the discharge starting voltage and the flashover voltage along the surface of PI. The discharge signal was measured using the ETS-93686 high-frequency pulsed current sensor. The applied voltage and discharge signal data were collected using a Yokogawa DL6154 four-channel oscilloscope. The real-time collected data is sent to the computer for storage through the USB serial bus and LabVIEW's data acquisition storage system. The development of surface discharge is continuously photographed and recorded using a high-speed camera above the needle electrode.
The experimental results show that the surface discharge carbonizes the dielectric surface and increases the PI surface conductivity. When the voltage frequency is 50Hz, the surface discharge is dendritic, and the discharges appear on both rising and falling edges during a positive half-cycle at this stage, the maximum discharge amplitude is 0.013V, and the number of discharges is 2 180. At high-frequency with the insulation's ageing, there is an increase in the discharge amplitude from 0.009V to 0.015V, and the total number of discharges is 3 610 recorded in the final stage.
To reveal temporal-spatial parameters a fluid-kinetic plasma simulation model was used. This study constructed a fluid-kinetic numerical simulation model by coupling the fluid model with the kinetic model. In addition, the model solves the complex chemical reactions logically to tackle the numerical oscillation problem. The simplified model is predicated on the maintained current carried by charged carriers. Microscopic quantities such as space charge, electron density, reaction rates, and electron temperature are obtained by simulation. The PI surface discharge trajectory at high frequency is linear, and the electron density reaches up to 1.1×1019 m-3 at 176×10-6 s.
The following conclusions can be drawn：(1)The development process of surface discharge can be divided into four stages. At each stage, the electron density exhibits a linearized characteristic. In this study, through phase resolved partial discharge (PRPD) analysis, it was found that the number of discharges and the discharge amplitude were larger in the positive half cycle. (2) The electron density and temperature increase with the applied voltage amplitude and frequency increase, and the distribution of positive and negative ion density and electron density is similar. (3) The space charge density distribution of surface discharge can be obtained using the simplified model, and the results are consistent with the self-consistent plasma model. (4) With the development of the surface discharge process, the electric field's influence gradually increases the reaction rate. In the final stage of the discharge, the total reaction rate can reach 1 200mol/(m3·s), and its distribution trend is similar to the electron distribution which indicates the correct chemical reaction in the background gas.
Keywords：Polyimide, surface discharge, space charge, electron density, electron temperature, the reaction rate
国家自然科学基金资助项目（52077127, 51929701）。
收稿日期 2022-06-30
DOI: 10.19595/j.cnki.1000-6753.tces.221295
中图分类号：TM855
Bilal Iqbal Ayubi 男，1992年生，博士研究生，研究方向为高频高压绝缘放电。E-mail：bilal.ayubi786@mail.sdu.edu.cn
张 黎 男，1979年生，副教授，博士生导师，研究方向为高电压与绝缘技术、电力系统电磁兼容。E-mail：zhleee@sdu.edu.cn（通信作者）
改稿日期 2022-07-26
（编辑 李冰）