基于FEM - FCT法洞察SF6N2混合气体流注电晕放电特性与机制

基于FEM-FCT法洞察SF6N2混合气体流注电晕放电特性与机制一、引言1.1研究背景与意义1.1.1SF6气体的应用现状与问题在现代电力系统中，六氟化硫（SF6）气体凭借其优异的电气性能，如高绝缘强度、良好的电弧熄灭能力和稳定的化学性质，被广泛应用于各类高压电气设备，如气体绝缘金属封闭开关设备（GIS）、变压器、高压断路器等。其绝缘强度约为空气的2.5倍以上，在高压环境下能够有效抑制放电现象，保障设备的安全稳定运行；良好的电弧熄灭性能使得它在切断电路时能迅速熄灭电弧，大大提高了设备的开断能力和可靠性。然而，随着全球对环境保护的关注度不断提高，SF6气体的高全球变暖潜势（GWP）问题日益凸显。SF6的GWP值是二氧化碳的23500倍，在大气中可存在长达3200年。这使得它成为《京都议定书》中明确限制排放的六种温室气体之一。据统计，全球每年排放的SF6总量相当于1.25亿吨CO2气体，且还在以每年近10%的速率继续增长。随着各国对温室气体排放限制的法规政策愈发严格，如美国加州提出从2020年开始逐年降低电气领域SF6使用量，2025年38kV以下电压等级设备将禁止使用SF6；欧盟计划在2030年将SF6排放量缩减到2014年的2/3，部分国家开始征收SF6碳排放税，SF6气体在电力设备中的应用面临着严峻的挑战。寻找一种既能满足电气性能要求，又具有良好环境友好性的替代气体或混合气体，成为了电力行业亟待解决的关键问题。1.1.2SF6/N2混合气体的研究价值SF6/N2混合气体作为潜在的SF6替代方案，近年来受到了广泛的关注。氮气（N2）是地球上含量最为丰富的气体，约占大气的78%，具有环境友好、成本低廉的显著优势。将SF6与N2混合，能够在一定程度上平衡两者的优缺点。一方面，通过添加N2，可以减少SF6的使用量，从而降低混合气体的全球变暖潜势，减轻对环境的压力；另一方面，SF6的高绝缘强度和良好的灭弧性能能够弥补N2电气性能的不足。研究表明，SF6/N2混合气体的绝缘性能和灭弧性能在一定混合比例范围内能够满足电力设备的基本要求，例如在SF6含量为20%-80%的混合气体中，能够实现良好的电气性能，在高压开关设备、变压器等电力设备中具有潜在的应用前景。深入研究SF6/N2混合气体的流注电晕特性具有重要意义。流注电晕是气体放电过程中的一个重要阶段，对混合气体的绝缘性能和电气设备的安全运行有着直接影响。通过研究流注电晕的发展过程、影响因素以及与混合气体成分、压力、电场强度等参数的关系，可以为SF6/N2混合气体在电力设备中的应用提供更坚实的理论基础，有助于优化电气设备的绝缘设计，提高设备的可靠性和稳定性，推动电力行业向绿色、可持续方向发展。1.2国内外研究现状1.2.1SF6/N2混合气体相关研究进展国内外学者对SF6/N2混合气体的电气性能进行了大量研究。在绝缘性能方面，伊斯坦布尔科技大学的E.Onal等人采用杆-面电极、间隙可调的试验装置，对非均匀电场下SF6、N2及SF6/N2混合气体的交流击穿强度进行了实验研究。结果表明，N2的击穿特性随压强近似线性变化，而SF6和SF6/N2混合气体均有最大最小值，尤其在间隙大于15mm时表现明显；N2+0.125%SF6的击穿电压在间隙5-20mm、压强100-500kPa时均高于N2，而较SF6低，仅在400-500kPa、间隙为25mm时击穿电压跌落到N2之下。西安交通大学邱毓昌等人通过不同混合比、不同压力下SF6/N2混合气体的工频耐压试验和操作冲击耐压试验，对其绝缘特性进行了深入研究。发现均匀场中，当混合气体各自组分的电子能量分布相同，才能用Wieland插值计算其电子碰撞电离系数及电子附着系数；在极不均匀场下，混合气体中的电晕稳化作用比纯SF6气体中的要强；在气压较高时，混合气体在较高频率VFT0作用下的击穿电压比较低频率VFT0作用下的要低，且频率越高击穿电压降低得越多；均匀场中混合气体的击穿电压与气体的压强基本呈线性关系。在灭弧特性研究方面，三菱公司的Kokura等人研究了在隔离开关中采用SF6/N2混合气体的开断性能，并采用永磁体产生的磁场驱动电弧运动，通过燃弧时间的长短来衡量其开断能力。研究表明，加入的驱动磁场强度对于纯SF6的开断性能几乎没有影响，但当采用SF6/N2混合气体时，其燃弧时间随着磁场的增强而显著缩短。当磁感应强度为0.59T，开断电流为0.8kA(rms)时，含10%SF6的混合气体具有和纯SF6相同的燃弧时间。韩国晓星公司的Song等学者利用实际的170kV/50kA的压气式GCB进行实验，研究了SF6/N2混合气体的开断性能，通过弧隙所能承受的最大的恢复电压上升率（RRRV）来评估。结果表明，对于SF6/N2混合气体，当N2含量为25%时和纯SF6有相同的开断性能，而50%的N2时开断性能较差。在混合比例优化方面，众多研究致力于寻找在满足电气性能要求的前提下，能够最大程度减少SF6使用量的最佳混合比例。目前的研究成果表明，不同应用场景和电气设备对混合比例的要求存在差异，需要综合考虑设备类型、运行条件等因素进行优化。1.2.2FEM-FCT法在气体放电仿真中的应用情况有限元-通量校正传输法（FEM-FCT）在气体放电仿真领域得到了一定的应用。20世纪80年代，国外科学家基于有限差分-通量校正传输法（FDM-FCT）分析了均匀场中空气、N2和SF6等气体的放电机理。21世纪初，剑桥大学首次采用FEM-FCT法对空气间隙中流注放电机理进行了二维模拟。国内部分高校也分别采用FDM-FCT法和FEM-FCT法分析了均匀场中SF6二维流注放电机理和极不均匀场中短空气间隙的流注放电机理。在SF6/N2混合气体流注电晕仿真研究方面，虽然已有一些相关工作，但仍存在不足。现有研究对于混合气体中复杂的物理过程，如不同气体成分之间的相互作用、电子碰撞电离和附着过程在混合气体中的特性等，考虑还不够全面和深入。此外，对于不同电场条件下，特别是在实际电气设备中存在的非均匀、复杂电场环境下，SF6/N2混合气体流注电晕的仿真研究还相对较少，缺乏系统性和针对性。这限制了对混合气体在实际应用中放电特性的准确理解和预测，需要进一步深入研究以完善相关理论和仿真方法。1.3研究内容与方法1.3.1研究内容概述本研究基于FEM-FCT法对SF6/N2混合气体流注电晕进行全面深入的仿真研究。首先，建立适用于SF6/N2混合气体流注电晕放电的二维或三维物理模型，综合考虑混合气体中电子碰撞电离、附着、复合等基本物理过程，以及空间电荷对电场的畸变作用、空间光电离与阴极光发射等因素。其次，利用FEM-FCT法对所建立的模型进行数值求解，模拟不同混合比例、压力、电场强度等条件下SF6/N2混合气体流注电晕的发展过程，分析流注电晕的起始条件、发展速度、形态特征以及对电场分布的影响。通过仿真结果，研究混合气体参数与流注电晕特性之间的内在关系，揭示其放电机理。再者，对仿真结果进行详细的参数分析，包括混合气体中SF6含量、N2含量的变化对流注电晕特性的影响，不同压力和温度条件下的放电特性差异，以及电场均匀度、电极形状等因素对流注电晕发展的作用。通过参数分析，确定影响SF6/N2混合气体流注电晕特性的关键因素，为优化混合气体绝缘性能和电气设备设计提供理论依据。最后，将仿真结果与相关实验数据或已有研究成果进行对比验证，评估仿真模型和方法的准确性和可靠性。针对仿真结果与实验数据之间的差异进行分析，进一步完善仿真模型和参数设置，提高仿真的精度和可信度。1.3.2研究方法阐述本研究主要采用FEM-FCT法结合相关物理模型进行仿真研究。FEM-FCT法将有限元方法（FEM）的高精度和通量校正传输法（FCT）处理非线性问题的优势相结合，能够有效地求解气体放电过程中的复杂物理方程。在建立物理模型时，依据气体放电的基本理论，考虑SF6/N2混合气体的特性，构建包含电子、离子、中性粒子等多组分的输运方程，以及电场的泊松方程，全面描述混合气体流注电晕放电过程中的物理现象。为了验证仿真结果的准确性，将结合相关实验数据进行对比分析。若有条件，将设计并开展针对SF6/N2混合气体流注电晕的实验研究，获取不同条件下的放电数据，如起始放电电压、流注发展速度、放电电流等。通过将实验数据与仿真结果进行对比，评估仿真模型的可靠性，对模型进行修正和完善。此外，还将采用对比分析的方法，研究不同参数条件下SF6/N2混合气体流注电晕特性的差异。通过改变混合气体的比例、压力、电场强度等参数，对比分析仿真结果，深入探讨各参数对放电特性的影响规律，从而为混合气体在电力设备中的应用提供优化建议。二、理论基础2.1SF6/N2混合气体特性2.1.1SF6气体特性SF6气体作为一种人工合成的惰性气体，在电力领域具有独特且卓越的性能优势。从电气性能角度来看，其绝缘强度表现极为突出，在均匀电场条件下，它的绝缘强度约为空气的2.5倍以上。这使得它能够在高电压环境中有效抑制气体放电现象，保障电气设备的可靠运行。例如在气体绝缘金属封闭开关设备（GIS）中，SF6气体填充于设备内部，为高压带电部分提供可靠的绝缘屏障，确保设备在运行过程中不会因绝缘问题而发生故障。在灭弧性能方面，SF6气体更是表现出色。当电气设备中的电弧产生时，SF6气体能够迅速与电弧发生相互作用。其分子在电弧高温作用下会分解产生活性物质，这些活性物质能够与电弧中的带电粒子发生反应，从而有效减弱电弧的能量。同时，SF6气体具有较高的热导率，能够快速吸收电弧的热量，使电弧温度迅速降低，进而实现快速灭弧。以高压断路器为例，在切断电路时，强大的电弧会瞬间产生，SF6气体凭借其优异的灭弧性能，能够在极短的时间内熄灭电弧，大大提高了断路器的开断能力和可靠性，保障了电力系统的安全稳定运行。从化学性质方面分析，SF6气体具有高度的稳定性。其分子结构呈八面体排布，键合距离小、键合能高，这使得在温度不超过180℃时，它与电气结构材料具有良好的相容性，如同氮气一样稳定，不会与设备中的金属、绝缘材料等发生化学反应，保证了设备的长期稳定运行。然而，SF6气体也存在着严重的缺点，其中最为突出的是其高全球变暖潜势（GWP）。根据相关研究，SF6的GWP值是二氧化碳的23500倍，并且在大气中可稳定存在长达3200年。这意味着即使是极少量的SF6气体排放到大气中，也会对全球气候变暖产生巨大的影响。随着全球对环境保护的重视程度不断提高，SF6气体的排放受到了越来越严格的限制。在《京都议定书》中，SF6被明确列为需要限制排放的六种温室气体之一。此外，SF6气体在电气设备中使用时，还存在其他一些问题。当它在强烈的电弧放电作用下，会分解产生一些含硫的低氟化物，这些物质具有腐蚀性，可能会对电气设备中的某些绝缘材料和金属部件造成损害，影响设备的使用寿命和性能。而且，由于电气设备制造和安装过程中存在的差异性，以及设备长时间运行后材料老化等因素，可能会导致SF6气体泄漏。一旦发生泄漏，SF6气体及其分解物会在室内低层空间积聚，造成局部缺氧和带毒的危险环境，对进入室内的工作人员的生命安全构成严重威胁。2.1.2N2气体特性N2是地球大气中含量最为丰富的气体，约占大气总体积的78%，具有诸多独特的性质，使其在多个领域得到广泛应用，在电力领域中作为潜在的绝缘气体替代方案也备受关注。从绝缘性能方面来看，N2的绝缘强度相对较低，大约为空气的1.0-1.2倍，这使得它在单独作为绝缘气体应用于高电压电气设备时存在一定的局限性。然而，在一些对绝缘要求相对不那么苛刻的场合，或者与其他绝缘性能优异的气体混合使用时，N2的绝缘性能仍能发挥一定的作用。N2具有化学惰性，其分子中的氮原子之间形成了三键，这种化学键的键能非常高，使得N2在常温常压下几乎不与其他物质发生化学反应。在焊接和切割过程中，使用含氮气体的混合物可以营造一个惰性气氛，防止金属被氧化，保证焊接和切割的质量。在电力设备中，这种化学稳定性也有助于维持设备内部环境的稳定，减少因化学反应导致的设备故障风险。N2还具有显著的环保优势。它是大气的自然组成部分，不会对环境造成额外的污染，也不存在温室效应等环境问题。与其他一些可能对环境产生负面影响的气体相比，N2的广泛使用不会增加全球变暖的压力，符合当今社会对可持续发展和环境保护的要求。在资源方面，由于N2在大气中含量丰富，获取成本相对较低，这使得它在大规模应用时具有经济可行性。在工业生产中，通过空气分离技术可以较为容易地制取高纯度的N2，为其在电力领域等多个行业的应用提供了充足的资源保障。2.1.3SF6/N2混合气体特性SF6/N2混合气体的电气性能会随着混合比例的变化而呈现出显著的改变。当SF6含量较高时，混合气体的绝缘强度和灭弧性能更接近纯SF6气体。研究表明，在均匀电场中，随着SF6含量的增加，混合气体的击穿电压逐渐升高。当SF6含量达到一定比例时，混合气体的击穿电压可与纯SF6气体相当，这使得在一些对绝缘性能要求较高的电气设备中，能够通过调整混合比例来满足绝缘需求。在高压断路器的灭弧应用中，较高含量的SF6可以确保在开断电路时，混合气体能够迅速熄灭电弧，保障设备的正常运行。随着N2含量的增加，混合气体的绝缘强度和灭弧性能会逐渐下降，但同时其全球变暖潜势也会显著降低。当N2含量较高时，虽然混合气体的电气性能有所减弱，但在一些对电气性能要求不是特别严格，同时对环保要求较高的场合，仍具有应用价值。在一些中低压电气设备或对环境要求较高的特殊应用场景中，可以适当增加N2的比例，以减少对环境的影响，同时通过优化设备设计等方式来弥补电气性能的不足。混合气体的绝缘性能还受到温度和压力等因素的影响。在高温环境下，分子的热运动加剧，电子的能量增加，这会使得混合气体中的电子更容易发生碰撞电离，从而降低混合气体的绝缘强度。当温度升高时，混合气体的击穿电压会下降，设备的绝缘可靠性面临挑战。在高压条件下，气体分子的密度增加，电子与分子之间的碰撞频率提高，这也会对混合气体的电气性能产生影响。压力的变化可能会导致混合气体的击穿电压发生改变，需要在实际应用中进行充分的考虑和研究。在设计和运行使用SF6/N2混合气体的电气设备时，需要综合考虑这些因素，通过合理选择混合比例、优化设备结构和运行条件等方式，确保设备在各种工况下都能安全可靠地运行。2.2流注电晕理论2.2.1流注放电基本原理流注放电是气体放电过程中的一个关键阶段，其从初始电子崩发展为流注的过程涉及一系列复杂的物理机制。在气体中，当存在足够强的电场时，自由电子会在电场力的作用下加速运动，获得足够的能量后与气体分子发生碰撞。这种碰撞会使气体分子发生电离，产生新的电子和正离子，从而形成电子崩。初始电子崩中的电子在向阳极运动的过程中，不断与气体分子碰撞，导致电离过程持续进行，电子崩中的电子数量呈指数增长。当电子崩头部的电子数达到一定数量时，电子崩头部的空间电荷会对电场产生显著的畸变作用。电子崩头部集中着大量的正离子和几乎全部电子，原有均匀场强在电子崩前方和尾部处都得到增强，而崩头内正、负电荷区域之间的电场则被削弱。随着电子崩的发展，崩头前后电场的增强有利于发生分子和离子的激励现象。当这些分子和离子从激励状态回复到正常状态时，会放射出光子。同时，崩头内部正、负电荷区域之间电场的削弱有助于发生复合过程，同样也会发射出大量光子。这些光子会引发空间光电离，产生新的光电子。光电子在电场作用下向正空间电荷区运动，进而引发新的碰撞游离，形成二次电子崩。二次电子崩的电子迅速移向阳极方向，在正电荷区域内形成正负带电粒子混合通道，这个电离通道即为流注。流注端部又有二次电子崩留下的正电荷，进一步加强了电场，促使更多的新电子崩相继产生并与之汇合，使流注不断向前发展。最后，流注发展到阴极，将两极接通，导致气隙被击穿，完成气体放电过程。在流注放电过程中，电子的迁移速度比正离子的要大两个数量级，这导致在电子崩发展过程中出现了大量的空间电荷，崩头最前面集中着电子，其后直到尾部则是正离子，而其外形则好似球头的锥体。这种空间电荷分布和电子、离子的运动特性对流注的形成和发展起到了关键作用。2.2.2电晕放电原理电晕放电是一种在强电场作用下，气体在电极表面局部区域发生的自持放电现象。其产生需要满足一定的条件，其中最关键的是电场强度要达到或超过气体的起始电晕场强。在实际的电气设备中，当电极表面存在曲率半径较小的部位，或者电极之间的距离较大且施加的电压较高时，就容易在这些部位形成强电场，从而引发电晕放电。在电晕放电起始阶段，电极表面附近的气体分子首先发生电离，产生少量的电子和正离子。这些电子在强电场的作用下获得能量，与周围的气体分子发生碰撞，导致更多的气体分子电离，形成电子崩。随着电子崩的发展，空间电荷逐渐积累，电场分布发生畸变。在电极附近，正离子聚集形成正空间电荷区，使电场强度进一步增强，而在远离电极的区域，电场强度则相对减弱。电晕放电具有一些明显的特征。在外观上，电晕放电通常会在电极表面出现微弱的发光现象，呈现出蓝紫色的光晕，这是由于气体分子在电离和复合过程中发射出光子所致。电晕放电还会伴随着声音的产生，通常是一种轻微的“嘶嘶”声，这是由于气体分子的电离和激发导致气体振动所引起的。电晕放电会产生一定的能量损耗，这会导致电气设备的效率降低，同时还可能对周围的电磁环境产生干扰，影响其他设备的正常运行。2.2.3SF6/N2混合气体流注电晕特性在SF6/N2混合气体中，流注电晕的发展具有独特的特点，并且与纯气体存在明显的差异。由于SF6和N2的物理性质不同，混合气体中电子与不同气体分子的碰撞电离和附着过程也会有所变化。SF6是强电负性气体，其分子极易吸附自由电子而形成质量大的负离子，这会削弱气体中碰撞电离过程。在混合气体中，随着SF6含量的增加，电子更容易被SF6分子捕获，从而减少了参与碰撞电离的自由电子数量，使得流注电晕的发展速度相对减慢。当SF6含量较高时，混合气体中的电子崩发展相对较为缓慢，流注的起始和发展需要更高的电场强度。相比之下，N2的绝缘强度较低，电子在N2中的碰撞电离相对容易发生。当N2含量增加时，混合气体中的碰撞电离过程会相对增强，这可能会使流注电晕的发展速度加快，但同时也会导致混合气体的绝缘性能下降。在一定混合比例下，N2含量的增加可能会使混合气体在较低电场强度下就出现流注电晕现象，但此时混合气体的击穿电压也会相应降低。混合气体中不同气体成分之间的相互作用也会对流注电晕特性产生影响。在放电过程中，SF6和N2分子之间可能会发生能量转移和化学反应，这些过程会改变电子的能量分布和气体分子的状态，进而影响流注电晕的发展。由于这些复杂的相互作用，混合气体的流注电晕特性并非简单地是纯SF6和纯N2特性的线性叠加，而是需要综合考虑多种因素来进行分析和研究。2.3FEM-FCT法原理2.3.1有限元法（FEM）基础有限元法是一种借助电子计算机解决工程问题的强大近似方法，其基本思想可以概括为“化整为零，集零为整”。在处理实际问题时，有限元法将原本连续的物理场或结构体离散化为一系列相互连接、形状规则或不规则的有限大小单元，这些单元通过节点相互连接，构建出一个逼近真实物理模型的数字框架。在对一个复杂的结构进行力学分析时，首先将该结构假想地分割成由有限个三角形、四边形等单元组成的集合体，这个过程称为“离散化”。离散化的关键在于合理地划分单元和确定节点位置，以确保能够准确地描述结构的几何形状和物理特性。对于一个具有复杂外形的机械零件，需要根据其形状特点和受力情况，选择合适的单元类型和划分方式，使得离散后的模型能够尽可能地接近真实结构。在完成离散化后，需要对每个单元进行力学特征分析。通过基于变分原理或加权余量法，建立单元节点力和节点位移之间的关系。这个过程涉及到对单元内部物理行为的数学描述，根据单元的材料属性、几何形状以及所受的载荷，推导出相应的数学方程，以描述单元在受力情况下的响应。对于一个弹性力学问题，需要根据弹性力学的基本原理，建立单元的应力-应变关系，从而得到单元节点力和节点位移之间的关系式。把所有单元的这种关系式集合起来，形成整个结构的力学特性关系，即得到一组以节点位移为未知量的代数方程组。通过对这些代数方程组进行处理，利用计算机的强大计算能力进行求解，就可以得到结点的位移。在得到节点位移后，进一步根据相关的物理公式和原理，计算出结构的应变和应力，从而对结构的力学性能进行评估和分析。有限元法是一种近似方法，它通过用离散单元的组合体来逼近原始结构，在几何形态上实现了对真实结构的逼近，同时在数学层面采用预设的近似函数来逼近未知变量的真实解。通过巧妙的插值技术，确保了近似解在单元边界上的连续性和协调性，从而在整体上达到了较高的精度，能够有效地解决各种复杂的工程问题。2.3.2通量校正传输法（FCT）原理通量校正传输法（FCT）主要用于解决数值计算中出现的数值震荡问题，以保证计算精度。在气体放电等复杂物理过程的数值模拟中，由于物理量的变化剧烈，传统的数值方法往往会导致数值解出现非物理的震荡，这会严重影响计算结果的准确性和可靠性。FCT法的核心原理是基于通量校正的思想。在数值计算中，物理量的传输是通过通量来描述的。FCT法通过对数值计算中得到的通量进行校正，来消除数值震荡。它首先计算出一个初步的数值解，这个解可能包含了数值震荡。然后，通过分析数值解的特性，确定需要校正的通量部分。在对流扩散问题的数值模拟中，传统的有限差分方法可能会导致解出现震荡。FCT法会根据解的变化情况，判断哪些区域的通量需要校正。它会计算出一个校正通量，这个校正通量是根据数值解的局部特性和物理过程的特点来确定的。通过将校正通量添加到原始的通量中，得到校正后的通量，再用校正后的通量重新计算数值解，从而有效地抑制了数值震荡，使计算结果更加接近真实的物理过程。FCT法通常采用两步法来实现通量校正。第一步是进行预测步骤，使用一种相对简单的数值方法计算出一个初步的数值解和通量分布。第二步是校正步骤，根据预测步骤得到的结果，计算出校正通量，并对初步的数值解进行修正。通过这种方式，FCT法能够在保持计算效率的同时，有效地提高计算精度，使数值模拟结果更加准确地反映物理过程的真实特性。2.3.3FEM-FCT法的耦合机制FEM-FCT法的耦合是将有限元法（FEM）和通量校正传输法（FCT）的优势相结合，以实现对气体放电过程的精确仿真。有限元法在处理复杂几何形状和边界条件方面具有独特的优势，能够准确地离散化物理场域，建立精确的数学模型；而通量校正传输法能够有效地解决数值计算中的震荡问题，保证计算结果的精度和稳定性。在耦合过程中，首先利用有限元法对气体放电区域进行离散化，将连续的场域划分为有限个单元，并建立每个单元的控制方程。在处理一个包含复杂电极结构的气体放电问题时，有限元法可以根据电极的形状和位置，合理地划分单元，准确地描述电场的分布情况。通过有限元法得到的数值解，虽然能够反映物理场的大致分布，但可能会存在数值震荡，影响计算结果的准确性。此时，引入通量校正传输法对有限元法得到的数值解进行处理。根据有限元法计算得到的物理量分布，如电场强度、电子密度等，FCT法计算出相应的通量，并对通量进行校正。在计算电子在电场中的传输时，FCT法会根据有限元法得到的电场分布和电子密度分布，计算出电子的通量。由于数值计算的误差，这个通量可能会导致数值解出现震荡。FCT法通过分析通量的特性，计算出校正通量，对原始通量进行修正，从而消除数值震荡。通过这种耦合机制，FEM-FCT法既能够充分利用有限元法处理复杂几何和边界条件的能力，又能够借助通量校正传输法保证计算结果的精度和稳定性，实现对气体放电过程中复杂物理现象的准确模拟和分析，为研究SF6/N2混合气体流注电晕特性提供了有力的工具。三、仿真模型建立3.1物理模型构建3.1.1几何模型建立以典型的棒-板电极结构为例构建用于仿真的几何模型，这种电极结构在研究气体放电特性时被广泛应用，能够产生较为典型的非均匀电场，对于研究SF6/N2混合气体流注电晕特性具有重要意义。在实际的高压电气设备中，如气体绝缘金属封闭开关设备（GIS）中的母线连接部位、变压器的绕组端部等，都存在类似的非均匀电场结构，因此选择棒-板电极结构进行研究具有很强的实际应用背景。设定棒电极采用圆柱形，半径为r1=5mm，长度为L1=50mm；板电极采用圆形平板，半径为r2=100mm，厚度为t=10mm。棒电极与板电极之间的距离d=20mm，这个距离在实际电气设备中是一个常见的绝缘间隙尺寸，通过研究该间隙下的放电特性，可以为实际设备的绝缘设计提供参考。整个几何模型置于一个边长为a=200mm的正方体计算域内，计算域的边界条件设置为：棒电极施加直流高电压U，电压值根据具体仿真需求在0-100kV范围内变化，以模拟不同电场强度下的放电情况；板电极接地，电势为0V；计算域的其他边界设置为绝缘边界条件，即垂直于边界的电场强度分量为0，以确保计算域内的电场分布不受外部干扰。利用专业的建模软件，如COMSOLMultiphysics自带的几何建模模块，首先创建棒电极的圆柱体模型，通过输入半径和长度参数确定其几何形状。接着创建板电极的圆形平板模型，设置好半径和厚度。将棒电极和板电极按照设定的距离进行定位，确保两者的相对位置准确。然后创建正方体计算域模型，并将棒-板电极结构放置在计算域的中心位置。在建模过程中，对电极边缘和计算域边界进行适当的平滑处理，以避免因几何形状的突变而导致电场计算出现奇异性，影响仿真结果的准确性。通过这种方式，构建出一个准确、合理的几何模型，为后续的物理过程仿真提供了可靠的基础。3.1.2物理过程描述在SF6/N2混合气体流注电晕放电过程中，涉及一系列复杂的物理过程，这些过程相互作用，共同决定了放电的特性和发展。碰撞电离是其中的关键过程之一。当混合气体中的电子在电场作用下获得足够的能量时，会与气体分子（包括SF6和N2分子）发生碰撞，使分子电离，产生新的电子和正离子。在强电场区域，电子的加速运动使其能量不断增加，当电子能量达到或超过气体分子的电离能时，就会引发碰撞电离。对于SF6分子，其电离能约为15.6eV，N2分子的电离能约为15.58eV。电子与SF6分子碰撞时，可能会使SF6分子电离为SF6+和一个自由电子；与N2分子碰撞时，会产生N2+和自由电子。这个过程导致电子和正离子的数量不断增加，为流注电晕的发展提供了初始的带电粒子源。附着过程是电负性气体SF6特有的重要物理过程。由于SF6是强电负性气体，其分子极易吸附自由电子而形成质量大的负离子。当电子与SF6分子碰撞时，除了可能发生碰撞电离外，还很容易被SF6分子捕获，形成SF6-负离子。这种附着过程会减少参与碰撞电离的自由电子数量，对放电的发展起到抑制作用。在混合气体中，随着SF6含量的增加，附着过程的影响会更加显著，使得流注电晕的起始和发展需要更高的电场强度。扩散过程在混合气体放电中也起着重要作用。带电粒子（电子、离子）由于热运动，会从高浓度区域向低浓度区域扩散。在流注电晕发展过程中，电子和离子在空间中的分布并不均匀，扩散过程使得带电粒子能够在空间中重新分布，影响电场的分布和放电的发展。电子的扩散速度相对较快，这是因为电子质量小，热运动速度大。在电场和扩散的共同作用下，电子和离子的运动轨迹变得复杂，进一步影响了放电的特性。空间电荷对电场的畸变作用也是不可忽视的物理过程。在放电过程中，产生的大量电子和离子形成空间电荷，这些空间电荷会改变原有的电场分布。在电子崩头部，集中着大量的正离子和电子，使得电子崩前方和尾部的电场得到增强，而崩头内正、负电荷区域之间的电场则被削弱。这种电场畸变会影响后续电子的运动和碰撞电离过程，对流注的形成和发展产生重要影响。当电场畸变到一定程度时，会引发新的电子崩和流注的产生，促使放电进一步发展。空间光电离与阴极光发射也是流注电晕放电过程中的重要环节。在放电过程中，气体分子和离子的激励和复合会发射出光子。这些光子在空间中传播，当遇到其他气体分子时，可能会引发空间光电离，产生新的光电子。光电子在电场作用下向正空间电荷区运动，进而引发新的碰撞游离，形成二次电子崩，加速流注的发展。阴极表面在受到光子照射时，也可能会发射出光电子，这些光电子同样会参与到放电过程中，为流注的发展提供更多的带电粒子。3.2数学模型建立3.2.1连续性方程连续性方程用于描述电子、离子等带电粒子浓度随时间和空间的变化情况，是研究气体放电过程的重要基础方程之一。对于电子浓度ne，其连续性方程可以通过质量守恒定律推导得到。在一个微小的空间体积元内，电子浓度的变化率由三部分组成：电子的产生率、电子的消失率以及电子的净通量。假设电子的产生主要来源于碰撞电离过程，碰撞电离率为αne，其中α为电子碰撞电离系数，与电场强度、气体种类和温度等因素有关；电子的消失主要是由于附着过程，附着率为ηne，其中η为电子附着系数，对于SF6/N2混合气体，其值与混合气体的比例有关；电子的净通量由电子的漂移速度ve和扩散系数De决定，根据菲克定律，电子的扩散通量为-De∇ne，电子在电场E作用下的漂移通量为vene。根据质量守恒定律，在单位时间内，进入和离开体积元的电子净通量等于体积元内电子浓度的变化率，因此可以得到电子连续性方程：\frac{\partialn_e}{\partialt}=\alphan_e-\etan_e-\nabla\cdot(v_en_e-D_e\nablan_e)对于正离子浓度ni+和负离子浓度ni-，其连续性方程的推导过程与电子类似。正离子的产生主要来源于碰撞电离，消失主要是由于复合过程，复合率为βni+ne，其中β为复合系数；正离子的净通量同样由漂移速度vi+和扩散系数Di+决定。负离子的产生主要是由于电子附着，消失主要是由于复合过程，其净通量由漂移速度vi-和扩散系数Di-决定。因此，正离子和负离子的连续性方程分别为：\frac{\partialn_{i^+}}{\partialt}=\alphan_e-\betan_{i^+}n_e-\nabla\cdot(v_{i^+}n_{i^+}-D_{i^+}\nablan_{i^+})\frac{\partialn_{i^-}}{\partialt}=\etan_e-\betan_{i^-}n_e-\nabla\cdot(v_{i^-}n_{i^-}-D_{i^-}\nablan_{i^-})这些连续性方程全面地描述了电子、离子在混合气体中的产生、消失和输运过程，为深入研究SF6/N2混合气体流注电晕放电特性提供了重要的数学依据。3.2.2泊松方程泊松方程在研究气体放电过程中起着关键作用，它主要用于计算空间电场分布。在气体放电的空间中，电场强度E与电荷密度ρ之间存在着紧密的联系，这种联系通过泊松方程得以体现。泊松方程的数学表达式为：\nabla^2\varphi=-\frac{\rho}{\epsilon_0}其中，\varphi表示电势，\epsilon_0是真空介电常数，其值约为8.854×10^{-12}F/m，\rho为空间电荷密度，它是由电子、正离子和负离子的电荷密度共同构成，即\rho=e(n_{i^+}-n_{i^-}-n_e)，这里e为元电荷，其值约为1.602×10^{-19}C。在实际应用中，泊松方程与连续性方程相互耦合。连续性方程描述了电子、离子等带电粒子浓度的变化，而这些带电粒子的分布又会影响空间电荷密度，进而通过泊松方程影响电场分布。电场分布的变化反过来又会作用于带电粒子的运动，影响碰撞电离、附着等物理过程，从而改变带电粒子的浓度分布，形成一个相互影响、相互作用的复杂体系。在求解泊松方程时，需要根据具体的边界条件进行数值计算。对于前面构建的棒-板电极几何模型，棒电极施加直流高电压U，板电极接地，电势为0V，计算域的其他边界设置为绝缘边界条件，即垂直于边界的电场强度分量为0。通过合理地设置这些边界条件，可以准确地求解出空间电势分布，进而得到电场强度分布，为研究SF6/N2混合气体流注电晕放电过程提供准确的电场信息。3.2.3其他相关方程在研究SF6/N2混合气体放电过程中，除了连续性方程和泊松方程外，反应速率方程也是非常重要的辅助方程。反应速率方程用于描述混合气体中各种物理过程的反应速率，它与混合气体的成分、温度、电场强度等因素密切相关。对于碰撞电离过程，其反应速率方程可以表示为：R_{\alpha}=\alphan_e其中，R_{\alpha}表示碰撞电离反应速率，\alpha为电子碰撞电离系数，它是电场强度和气体温度的函数。在不同的电场强度和温度条件下，电子获得的能量不同，与气体分子碰撞时发生电离的概率也不同，因此\alpha的值会发生变化。通过实验测量和理论计算，可以得到不同条件下的\alpha值，从而准确地描述碰撞电离过程的反应速率。对于附着过程，其反应速率方程为：R_{\eta}=\etan_e其中，R_{\eta}表示附着反应速率，\eta为电子附着系数，它主要取决于混合气体中SF6的含量。由于SF6是强电负性气体，其含量越高，电子被捕获形成负离子的概率就越大，\eta的值也就越大。通过研究不同混合比例下的附着系数，可以深入了解附着过程对混合气体放电特性的影响。复合过程的反应速率方程为：R_{\beta}=\betan_{i^+}n_e=\betan_{i^-}n_e其中，R_{\beta}表示复合反应速率，\beta为复合系数，它与气体温度、压力等因素有关。在不同的温度和压力条件下，离子和电子的热运动速度和碰撞频率不同，复合的概率也会发生变化，因此\beta的值会相应改变。通过实验和理论研究，确定不同条件下的复合系数，有助于准确地描述复合过程的反应速率。这些反应速率方程与连续性方程、泊松方程相互配合，全面地描述了混合气体放电过程中的物理现象，为深入研究SF6/N2混合气体流注电晕特性提供了完整的数学模型。3.3FEM-FCT法求解过程3.3.1离散化处理在运用FEM-FCT法求解之前，需要对建立的数学模型进行离散化处理，将连续的数学模型转化为适合数值计算的离散形式。对于有限元法（FEM）部分，首先对计算域进行网格划分。采用三角形或四边形单元对前面构建的棒-板电极几何模型及其所在的正方体计算域进行离散，在电极附近和电场变化剧烈的区域，适当加密网格，以提高计算精度。在棒电极和板电极的边缘，以及流注电晕可能发展的区域，网格尺寸设置得较小，例如在电极边缘区域，网格边长可以设置为0.1mm，而在远离电极的区域，网格边长可以适当增大，如0.5mm。这样的网格划分策略能够在保证计算精度的同时，控制计算量，提高计算效率。在每个单元上，定义插值函数来近似表示物理量的分布。对于电势\varphi，采用线性插值函数，即在三角形单元中，设三个顶点的电势分别为\varphi_1、\varphi_2、\varphi_3，则单元内任意一点的电势\varphi(x,y)可以表示为：\varphi(x,y)=N_1(x,y)\varphi_1+N_2(x,y)\varphi_2+N_3(x,y)\varphi_3其中，N_1(x,y)、N_2(x,y)、N_3(x,y)为三角形单元的形状函数，它们是关于坐标(x,y)的线性函数，满足在顶点处的值为1，在其他顶点处的值为0的条件。通过这种方式，将连续的电势分布用离散的节点电势来近似表示。对于通量校正传输法（FCT）部分，对连续性方程中的通量项进行离散处理。以电子通量为例，将电子的漂移通量v_en_e和扩散通量-D_e\nablan_e在每个单元上进行离散计算。在离散过程中，考虑相邻单元之间的通量传递，采用合适的数值格式来计算通量，如中心差分格式或迎风格式。对于电子的漂移通量，在每个单元的边界上，根据单元内的电子浓度和漂移速度，利用相应的数值格式计算出通过边界的通量；对于扩散通量，同样根据单元内的电子浓度梯度和扩散系数，在边界上计算出扩散通量。通过这种离散化处理，将连续性方程转化为适合FCT法求解的形式。3.3.2迭代求解步骤FEM-FCT法的迭代求解过程是一个逐步逼近真实解的过程，主要包括以下几个关键步骤：预测步骤：利用有限元法对泊松方程进行求解，得到电势分布\varphi的初步解。在这个过程中，根据前面离散化处理得到的有限元方程，采用迭代求解算法，如共轭梯度法或高斯-赛德尔迭代法，求解节点电势。通过求解泊松方程得到的电势分布，进而计算出电场强度E。根据电场强度E和相关物理参数，利用连续性方程计算出电子、离子浓度的预测值。在计算电子浓度的预测值时，根据前面推导的电子连续性方程，将电场强度代入碰撞电离系数和漂移速度的计算公式中，结合初始的电子浓度分布，计算出下一时间步的电子浓度预测值。同样的方法可以计算出正离子和负离子浓度的预测值。通量计算：根据预测步骤得到的电子、离子浓度和电场强度，计算出各个单元的通量。对于电子通量，计算漂移通量v_en_e和扩散通量-D_e\nablan_e；对于正离子和负离子通量，也进行类似的计算。在计算通量时，采用前面离散化处理中确定的数值格式，确保通量计算的准确性。在计算电子的漂移通量时，根据电场强度和电子迁移率计算出电子的漂移速度，再结合电子浓度计算出漂移通量；对于扩散通量，根据电子浓度梯度和扩散系数计算出扩散通量。校正步骤：利用通量校正传输法对预测步骤得到的通量进行校正。通过比较预测通量和物理上合理的通量范围，计算出校正通量。具体来说，根据FCT法的原理，确定一个通量限制器，通过通量限制器对预测通量进行调整，得到校正后的通量。通量限制器的选择需要考虑计算的稳定性和精度，常用的通量限制器有minmod限制器、superbee限制器等。以minmod限制器为例，它根据相邻单元之间的通量变化情况，对预测通量进行限制，使得校正后的通量更加符合物理实际。更新步骤：根据校正后的通量，更新电子、离子浓度。将校正后的通量代入连续性方程，计算出新的电子、离子浓度值。在更新电子浓度时，根据电子连续性方程，将校正后的通量代入相应的项中，计算出下一时间步的电子浓度。同样的方法可以更新正离子和负离子浓度。然后，将更新后的电子、离子浓度作为下一次迭代的初始值，重复上述预测、通量计算、校正和更新步骤，直到计算结果收敛。3.3.3计算收敛条件为了确保仿真计算结果的准确性和可靠性，需要确定合理的计算收敛条件。在FEM-FCT法的迭代求解过程中，通常采用以下两种方式来判断计算是否收敛：残差收敛准则：计算每次迭代过程中方程的残差，当残差小于设定的阈值时，认为计算收敛。以泊松方程为例，定义残差R_{\varphi}为当前迭代步计算得到的电势与上一次迭代步电势的差值的范数，即：[R_{\varphi}四、仿真结果与分析4.1流注电晕发展过程分析4.1.1不同时刻流注形态变化通过仿真得到了不同时刻SF6/N2混合气体中流注电晕的形态演变情况，结果如图1所示。在t=0.1ns时，棒电极附近开始出现少量电子崩，此时电子崩的范围较小，呈分散状分布。这是因为在初始阶段，电场强度相对较低，只有少数电子获得足够能量引发碰撞电离，形成初始电子崩。随着时间推移，在t=0.5ns时，电子崩逐渐发展并开始合并，形成了较为明显的流注通道，流注头部呈现出锥形，且有向板电极发展的趋势。这是由于电子崩中的空间电荷对电场产生畸变作用，使得电场在电子崩头部和尾部增强，促进了新的电子崩产生并与原电子崩合并，从而形成流注通道。到t=1.0ns时，流注通道进一步发展，长度增加，且分支增多，流注头部的电场强度进一步增强，使得流注能够继续向板电极快速推进。在这个阶段，空间光电离和阴极光发射产生的新电子不断加入流注发展过程，导致流注通道不断扩展和分支。从图中可以清晰地看出流注电晕从初始阶段到发展阶段的动态变化过程。在初始阶段，电子崩的产生和发展较为缓慢，分布较为分散；随着时间的增加，电子崩逐渐合并形成流注，流注的发展速度逐渐加快，形态变得更加复杂，分支增多，这与流注放电的理论预期相符。流注电晕的发展过程受到电场强度、混合气体成分以及空间电荷等多种因素的共同作用，这些因素相互影响，导致流注形态不断变化。4.1.2带电粒子浓度分布变化在流注电晕发展过程中，电子、离子等带电粒子浓度在空间的分布变化情况对于理解放电机制至关重要。仿真结果显示，在t=0.1ns时，电子浓度主要集中在棒电极附近，这是因为棒电极附近电场强度较高，电子更容易获得能量发生碰撞电离。随着时间的推移，到t=0.5ns时，电子浓度在流注通道内逐渐增加，且在流注头部达到最大值，这是由于流注头部的电场强度最强，电子的碰撞电离过程最为剧烈。正离子浓度在流注通道内也有明显分布，且与电子浓度分布呈现一定的相关性，在电子浓度较高的区域，正离子浓度也相对较高，这是因为电子与气体分子碰撞电离产生电子的同时也产生了正离子。在t=1.0ns时，电子浓度和正离子浓度在流注通道内进一步增加，流注通道内的带电粒子浓度分布更加不均匀，流注头部的带电粒子浓度远远高于其他区域。负离子浓度在整个放电过程中相对较低，这是因为SF6虽然是强电负性气体，但在本仿真条件下，电子附着形成负离子的过程相对较弱。随着流注的发展，带电粒子浓度分布的变化反映了流注电晕中物理过程的动态变化，带电粒子浓度的不均匀分布会进一步影响电场分布，从而对流注的发展产生重要影响。4.1.3电场分布变化电场强度和电位分布随流注电晕发展发生了显著改变。在t=0.1ns时，电场强度在棒电极附近最高，随着距离棒电极距离的增加而逐渐降低，呈现出典型的非均匀电场分布特征。这是由于棒-板电极结构本身会产生非均匀电场，棒电极曲率半径小，电场集中在其附近。当流注电晕开始发展后，在t=0.5ns时，流注头部的电场强度明显增强，这是因为流注头部的空间电荷对电场产生了畸变作用，使得电场在流注头部得到增强。同时，在流注通道内，电场强度也有所变化，与初始电场分布相比，流注通道内的电场强度分布变得更加不均匀。到t=1.0ns时，流注头部的电场强度进一步增强，且流注通道内的电场强度分布更加复杂。在流注分支处，电场强度也出现了局部增强的现象。电位分布也随着流注电晕的发展而发生改变，在棒电极附近电位最高，随着流注的发展，流注通道内的电位逐渐降低，到板电极处电位为0V。电场分布的变化对流注电晕的发展起到了关键作用，电场强度的增强会促进电子的碰撞电离和流注的发展，而电场分布的不均匀性会导致流注的形态和发展方向发生变化。4.2混合气体比例对流注电晕的影响4.2.1不同混合比例下的起始放电电压通过仿真对比了不同SF6/N2混合比例时的起始放电电压，结果如图2所示。可以看出，随着SF6含量的增加，起始放电电压呈现先升高后降低的趋势。当SF6含量较低时，由于N2的绝缘强度相对较低，混合气体中电子的碰撞电离相对容易发生，起始放电电压较低。随着SF6含量的增加，SF6的强电负性作用逐渐显现，电子更容易被SF6分子捕获，抑制了碰撞电离过程，使得起始放电电压升高。当SF6含量超过一定比例后，由于混合气体中SF6分子的浓度过高，电子与SF6分子的碰撞概率增大，更容易发生电离，导致起始放电电压又开始降低。在SF6含量为30%-50%时，起始放电电压相对较高，这表明在这个混合比例范围内，混合气体的绝缘性能较好。这是因为在这个比例下，SF6和N2的优势得到了较好的结合，SF6的强电负性抑制了电子的碰撞电离，而N2的存在又在一定程度上保证了混合气体的基本绝缘性能。不同混合比例下起始放电电压的变化规律对于选择合适的混合气体比例用于实际电气设备具有重要指导意义，在实际应用中，可以根据设备的电压等级和绝缘要求，选择起始放电电压满足要求的混合气体比例，以确保设备的安全运行。4.2.2流注发展速度与混合比例的关系研究发现混合气体比例对流注发展速度有着显著影响。随着SF6含量的增加，流注发展速度先降低后升高。当SF6含量较低时，混合气体中主要以N2为主，N2的绝缘强度较低，电子在其中的碰撞电离相对容易，流注发展速度较快。随着SF6含量的增加，由于SF6的强电负性，电子更容易被SF6分子捕获形成负离子，参与碰撞电离的自由电子数量减少，这使得流注发展速度降低。当SF6含量继续增加到一定程度后，虽然电子被捕获的概率增大，但由于SF6分子的存在使得混合气体的整体电离能力增强，在高电场区域，电子仍能获得足够能量引发碰撞电离，且电离产生的带电粒子在强电场作用下加速运动，导致流注发展速度又开始升高。在SF6含量为40%左右时，流注发展速度相对较慢，这说明此时混合气体的绝缘性能相对较好，能够有效地抑制流注的发展。这是因为在这个比例下，SF6对电子的捕获作用和混合气体的整体电离特性达到了一个相对平衡的状态，使得流注的发展受到了较好的抑制。混合气体比例对流注发展速度的影响机制较为复杂，涉及到电子与不同气体分子的碰撞电离、附着等多种物理过程，这些过程相互作用，共同决定了流注的发展速度。在实际应用中，了解这种关系可以帮助我们通过调整混合气体比例来控制流注的发展，从而提高电气设备的绝缘性能和可靠性。4.2.3放电稳定性分析不同混合比例下放电的稳定性存在差异。当SF6含量较低时，由于N2的绝缘性能相对较差，混合气体中的放电容易出现不稳定现象，如流注的发展呈现出不规则的形态，分支较多且杂乱无章。这是因为在低SF6含量下，电子的碰撞电离较为容易，放电过程中产生的带电粒子数量较多且分布不均匀，导致电场分布不稳定，从而使得流注的发展不稳定。随着SF6含量的增加，放电稳定性逐渐提高，流注的发展形态更加规则，分支相对较少。这是因为SF6的强电负性抑制了电子的碰撞电离，使得放电过程中的带电粒子数量相对较少且分布较为均匀，电场分布也更加稳定，从而提高了放电的稳定性。当SF6含量过高时，虽然放电稳定性仍然较好，但由于混合气体的整体电离能力增强，可能会导致放电能量增大，对电气设备造成潜在的危害。在选择混合气体比例时，需要综合考虑放电稳定性和其他性能指标，以确保电气设备的安全可靠运行。在一些对放电稳定性要求较高的电气设备中，如高压断路器，应选择能够保证放电稳定的混合气体比例，同时也要兼顾设备的其他性能要求，如绝缘强度、灭弧能力等。4.3外施电压对流注电晕的影响4.3.1电压幅值对放电特性的影响外施电压幅值变化时，流注电晕的发展速度、形态和放电稳定性等均发生了显著变化。随着电压幅值的增加，流注发展速度明显加快。这是因为更高的电压幅值会产生更强的电场强度，电子在电场中获得的能量更大，碰撞电离的概率增加，从而使得流注能够更快地发展。在低电压幅值下，流注的发展较为缓慢，形态相对规则，分支较少。随着电压幅值的升高，流注的形态变得更加复杂，分支增多，且流注头部的电场强度进一步增强，导致流注向各个方向发展的趋势更加明显。电压幅值的增加还会影响放电的稳定性。当电压幅值较低时，放电相对稳定，流注的发展过程较为有序。随着电压幅值的不断增大，放电的稳定性逐渐降低，容易出现不稳定放电现象，如流注的突然分叉、熄灭和重新发展等。这是因为高电压幅值下，放电过程中的能量较大，电场分布更加不均匀，容易导致放电过程的失控。在实际电气设备的运行中，需要合理控制外施电压幅值，以确保流注电晕的发展在安全稳定的范围内，避免因电压过高导致设备绝缘损坏或其他故障的发生。4.3.2电压频率对放电特性的影响外施电压频率对流注电晕放电特性也有着重要作用。研究发现，随着电压频率的增加，起始放电时间呈现出先缩短后增加的趋势。在较低频率下，随着频率的增加，电场的变化速度加快，电子在电场中获得能量的机会增多，更容易引发碰撞电离，从而使得起始放电时间缩短。当频率增加到一定程度后，由于电场的快速变化，电子在电场中来不及获得足够的能量引发碰撞电离，反而会受到电场变化的干扰，导致起始放电时间增加。电压频率的变化还会影响流注的发展形态和放电稳定性。在低频时，流注的发展形态相对较为稳定，分支较少；随着频率的升高，流注的发展形态变得更加复杂，分支增多，放电稳定性降低。这是因为高频电场的快速变化会导致电场分布的不稳定，使得流注的发展受到更多的干扰。在实际应用中，对于不同频率的外施电压，需要根据其对放电特性的影响，合理设计电气设备的绝缘结构和参数，以保证设备在不同频率下都能安全可靠地运行。五、实验验证与对比5.1实验设计与搭建5.1.1实验目的与方案本实验的核心目的在于对基于FEM-FCT法的SF6/N2混合气体流注电晕仿真结果进行验证，通过实际测量与仿真数据的对比，评估仿真模型的准确性和可靠性。实验采用的SF6/N2混合气体设定了多种不同的比例，包括20%SF6/80%N2、30%SF6/70%N2、40%SF6/60%N2、50%SF6/50%N2，以此全面研究混合气体比例对放电特性的影响。这些比例涵盖了在实际应用中可能使用的范围，通过对不同比例下放电特性的研究，可以为实际电气设备中混合气体的选择提供更全面的依据。电极结构采用与仿真模型一致的棒-板电极结构，以确保实验条件与仿真环境的一致性。棒电极采用不锈钢材质，直径为5mm，长度为50mm；板电极采用铜板，直径为100mm，厚度为10mm，棒电极与板电极之间的距离固定为20mm。这种电极结构能够产生典型的非均匀电场，有利于研究SF6/N2混合气体在非均匀电场下的流注电晕特性。测量方法方面，使用高电压测量仪测量起始放电电压，通过观察放电现象，记录首次出现稳定放电时的电压值。对于流注形态，采用高速摄像机进行拍摄记录，高速摄像机的拍摄帧率设置为10000帧/秒，能够清晰捕捉流注发展的瞬间形态。放电电流则通过电流传感器进行测量，选用的电流传感器精度为0.1A，能够准确测量放电过程中的电流变化。通过这些测量方法，获取实验条件下的关键放电参数，以便与仿真结果进行对比分析。5.1.2实验设备与仪器实验所需的高压电源采用直流高压电源，其输出电压范围为0-100kV，电压稳定度优于±1%，能够满足不同电场强度下的实验需求，确保实验过程中电压的稳定输出，为研究混合气体在不同电压条件下的放电特性提供可靠的电源保障。电极装置即上述的棒-板电极结构，棒电极和板电极经过精密加工，表面粗糙度小于0.1μm，以保证电极表面的光滑度，减少因电极表面缺陷对放电特性的影响。电极安装在一个绝缘支架上，绝缘支架采用聚四氟乙烯材料，具有良好的绝缘性能和机械强度，能够确保电极在实验过程中的稳定性和安全性。测量仪器包括高电压测量仪，其测量精度为±0.5%，能够准确测量起始放电电压；高速摄像机，分辨率为1920×1080，能够清晰捕捉流注发展的瞬间形态；电流传感器，测量精度为0.1A，能够实时监测放电电流的变化。还配备了数据采集系统，能够对测量仪器采集到的数据进行实时采集和存储，数据采集频率为10kHz，确保能够准确记录放电过程中的数据变化。5.1.3实验步骤与流程实验开始前，首先将实验装置进行清洁和检查，确保电极表面无杂质、无损伤，各连接部位牢固可靠。然后，将SF6和N2气体按照设定的比例通过气体混合装置进行混合，气体混合装置采用动态混合方式，能够确保混合气体的均匀性。混合好的气体充入放电腔室，放电腔室采用不锈钢材质，内部体积为10L，能够提供稳定的实验环境。充气过程中，使用气体纯度分析仪监测混合气体的纯度，确保混合气体的比例符合实验要求。连接好高压电源、测量仪器和数据采集系统，设置好各仪器的参数。将高压电源的输出电压设置为0，逐渐升高电压，升压速率为1kV/s，同时观察高速摄像机和电流传感器的数据。当观察到放电现象时，记录此时的电压值作为起始放电电压，同时通过高速摄像机拍摄流注形态，电流传感器记录放电电流。保持电压稳定，持续记录一段时间内的放电电流和流注形态变化。每次实验结束后，将放电腔室内的气体排空，并用干燥氮气冲洗3次，以避免残留气体对下一次实验的影响。按照上述步骤，对不同比例的SF6/N2混合气体进行多次实验，每组混合气体重复实验5次，取平均值作为实验结果，以提高实验数据的准确性和可靠性。5.2实验结果与仿真对比5.2.1起始放电电压对比实验测得的不同混合比例下的起始放电电压与仿真结果对比如表1所示。从表中数据可以看出，在20%SF6/80%N2混合比例下，实验测得的起始放电电压为32.5kV，仿真结果为30.8kV，相对误差为5.2%；在30%SF6/70%N2混合比例下，实验值为38.6kV，仿真值为36.2kV，相对误差为6.2%；在40%SF6/60%N2混合比例下，实验值为45.3kV，仿真值为42.8kV，相对误差为5.5%；在50%SF6/50%N2混合比例下，实验值为48.1kV，仿真值为45.6kV，相对误差为5.2%。整体来看，实验值与仿真值较为接近，相对误差均在7%以内，说明仿真模型在预测起始放电电压方面具有较高的准确性。起始放电电压的仿真值与实验值存在一定差异，这可能是由于实验过程中存在一些难以精确控制的因素，如气体的均匀性、电极表面的微观状态等。实验环境中的微小干扰也可能对起始放电电压产生影响。仿真模型虽然考虑了主要的物理过程，但在一些细节上可能与实际情况存在差异，如模型中对气体分子间相互作用的简化处理等。尽管存在这些差异，但仿真结果与实验结果的一致性表明，基于FEM-FCT法的仿真模型能够较好地反映SF6/N2混合气体起始放电电压随混合比例的变化趋势，为实际应用提供了有价值的参考。5.2.2流注形态对比实验观测到的流注形态与仿真模拟的流注形态对比如图3所示。从图中可以看出，在实验中，流注从棒电极开始发展，呈现出树枝状的形态，向板电极延伸。在不同混合比例下，流注的分支数量和长度有所不同，随着SF6含量的增加，流注的分支数量相对减少，长度相对增加。在20%SF6/80%N2混合比例下，流注分支较多且较为杂乱；而在50%SF6/50%N2混合比例下，流注分支相对较少，形态更加规则。仿真模拟的流注形态与实验观测结果具有相似性，同样呈现出从棒电极向板电极发展的树枝状形态。在分支数量和长度的变化趋势上，仿真结果与实验结果基本一致。在一些细节上，两者仍存在差异，实验中流注的边缘可能更加模糊，这可能是由于实验中存在一些微观的物理过程，如气体分子的热运动、杂质的影响等，这些因素在仿真模型中难以完全精确地考虑。仿真模型中对流注发展过程中的一些物理参数，如电子碰撞电离系数、附着系数等，采用了近似的计算方法，也可能导致流注形态的模拟与实际情况存在一定偏差。总体而言，仿真模型能够较好地模拟流注的宏观形态和发展趋势，为进一步研究流注电晕的物理机制提供了有效的工具。5.2.3放电特性参数对比除了起始放电电压和流注形态，还对其他放电特性参数，如放电电流、放电时间等进行了实验与仿真对比分析。在放电电流方面，实验测得的不同混合比例下的放电电流峰值与仿真结果对比如表2所示。在20%SF6/80%N2混合比例下，实验测得的放电电流峰值为0.85A，仿真结果为0.80A，相对误差为5.9%；在30%SF6/70%N2混合比例下，实验值为1.02A，仿真值为0.96A，相对误差为5.9%；在40%SF6/60%N2混合比例下，实验值为1.20A，仿真值为1.13A，相对误差为5.8%；在50%SF6/50%N2混合比例下，实验值为1.35A，仿真值为1.28A，相对误差为5.2%。可以看出，放电电流峰值的实验值与仿真值较为接近，相对误差在6%以内。在放电时间方面，实验测得的从起始放电到稳定放电的时间与仿真结果对比如表3所示。在20%SF6/80%N2混合比例下，实验测得的放电时间为5.5μs，仿真结果为5.2μs，相对误差为5.5%；在30%SF6/70%N2混合比例下，实验值为4.8μs，仿真值为4.5μs，相对误差为6.3%；在40%SF6/60%N2混合比例下，实验值为4.2μs，仿真值为3.9μs，相对误差为7.1%；在50%SF6/50%N2混合比例下，实验值为3.8μs，仿真值为3.5μs，相对误差为7.9%。随着混合气体中SF6含量的增加，放电时间逐渐缩短，仿真结果与实验结果在变化趋势上基本一致，但在具体数值上存在一定的相对误差，最大相对误差为7.9%。这些放电特性参数的实验值与仿真值存在差异的原因，一方面与实验过程中的测量误差、环境干扰等因素有关；另一方面，仿真模型中的一些假设和简化处理，如对气体放电过程中复杂化学反应的简化、对电极表面微观结构的忽略等，也可能导致仿真结果与实际情况存在偏差。尽管存在这些差异，仿真结果与实验结果在放电特性参数的变化趋势上的一致性，表明基于FEM-FCT法的仿真模型能够在一定程度上准确地预测SF6/N2混合气体的放电特性，为电气设备的设计和优化提供了重要的参考依据。5.3误差分析与改进措施5.3.1误差来源分析实验与仿真结果存在差异的原因是多方面的。从模型简化角度来看，仿真模型在建立过程中对一些复杂的物理过程进行了简化处理。在考虑电子碰撞电离和附着过程时，虽然采用了相关的系数来描述这些过程，但实际的气体放电过程中，电子与气体分子的碰撞行为受到多种因素的影响，如气体分子的振动、转动能级，以及分子间的相互作用等，这些因素在模型中难以完全精确地考虑。模型中对空间电荷的分布和运动的描述也可能存在一定的近似，实际的空间电荷分布可能更加复杂，其运动受到电场、磁场以及其他带电粒子的相互作用，而仿真模型可能无法完全捕捉这些复杂的相互作用。测量误差也是导致实验与仿真结果差异的重要原因之一。在实验过程中，测量仪器的精度限制会引入误差。高电压测量仪虽然精度为±0.5%，但在实际测量中，由于电压的波动、测量仪器的稳定性等因素，可能会导致测量结果存在一定的偏差。高速摄像机在拍摄流注形态时，由于拍摄角度、光线条件等因素的影响，可能无法完全准确地捕捉流注的真实形态。实验环境中的干扰因素，如环境温度、湿度的变化，以及外界电磁场的干扰等，也可能对实验结果产生影响。在不同的环境温度下，混合气体的物理性质可能会发生变化，从而影响放电特性，而实验过程中难以完全消除这些环境因素的干扰。5.3.2改进措施探讨针对误差来源，可以采取一系列改进方法来提高仿真和实验的准确性。在优化模型方面，需要进一步完善物理模型，考虑更多的物理过程和因素。对于电