CN109583066A - 一种直流架空线路绝缘子表面污秽沉积模拟方法 - Google Patents
一种直流架空线路绝缘子表面污秽沉积模拟方法 Download PDFInfo
- Publication number
- CN109583066A CN109583066A CN201811396424.9A CN201811396424A CN109583066A CN 109583066 A CN109583066 A CN 109583066A CN 201811396424 A CN201811396424 A CN 201811396424A CN 109583066 A CN109583066 A CN 109583066A
- Authority
- CN
- China
- Prior art keywords
- filthy
- particle
- insulator
- unit
- transmission line
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Insulators (AREA)
Abstract
本发明公开了一种直流架空线路绝缘子表面污秽沉积模拟方法,该方法充分考虑了绝缘子积污过程中的气流场特性、直流电场分布、以及污秽颗粒积聚/出射动态微观过程;在此基础上,本方法利用COMSOL耦合多物理场仿真软件开展了绝缘子表面污秽沉积数值计算,其结果能够较好地再现不同环境参数下的绝缘子表面污层分布情况;此外,通过该方法还可以计算得到绝缘子表面污秽质量密度、污秽不均匀度及积污带电系数。该方法可以做为强有力的工具,为分析输电线路外绝缘污秽问题提供支撑。
Description
技术领域
本发明属于输配电污秽外绝缘技术领域,尤其涉及一种直流架空线路绝缘子表面污秽沉积模拟方法。
背景技术
高压直流输电线路由于其输送距离远、线路造价低、输送容量大等特点而在近几年得到广泛投运。然而由于恒定电场的吸附效应,直流线路集聚的污秽量一般比相同环境下的交流线路高0.5–1倍,使得高海拔、重污秽等复杂环境地区外绝缘配置的可靠稳定运行面临着严峻考验。
工程中往往利用在直流输电线路悬挂参照绝缘子,开展定期污秽度测试,来获得直流线路积污带电系数,从而指导污区划分以及清扫工作的开展。通过悬挂参照绝缘子开展污秽测试的传统方法会消耗大量人力物力,且无法及时反映复杂气象环境下的绝缘子污秽度动态变化。此外,参照绝缘子的布置方式、电场分布与真实运行下的绝缘子情况差异很大,使得目前所积累的积污带电系数取值经验仍缺乏一定的科学性。
基于此,目前国内外科研机构开展了有关绝缘子污秽颗粒运动沉积的建模及仿真研究,旨在替代繁琐的自然积污试验来充分认知绝缘子的积污特性,并且推动自然环境下的绝缘子积污动态预测,从而更科学地指导污区划分、防污及清扫工作的开展。现有的仿真方法虽然对揭示绝缘子积污特性以及污秽度预测提供了重要参考,但是它们普遍得到的是颗粒碰撞系数或颗粒体积分数等无法与积污量直接取得联系的参数,难以通过试验来验证。而且所得仿真结果鲜有直观反映带电积污绝缘子表面污秽沉积及分布情况。
发明内容
针对现有技术存在的问题,本发明充分考虑绝缘子积污过程中的气流场特性、直流电场分布以及污秽颗粒积聚/出射动态微观过程,通过利用COMSOL软件,提供一种直流架空线路绝缘子表面污秽沉积模拟方法。
为实现上述目的,本发明采用的技术方案是:
一种直流架空线路绝缘子表面污秽沉积模拟方法,包括顺序相接的以下步骤:
步骤(1)、按照绝缘子尺寸参数及布置方式,在COMSOL软件内搭建仿真模型,划分计算域;
步骤(2)、采用软件自带静电场及流体力学模块,初始化计算域剖分网格,设置边界条件,迭代计算计算域内的静电场及流场稳态分布;
步骤(3)、采用软件自带流体流动粒子追踪模块,初始化计算域中的污秽颗粒,设置荷电量、电场力、气流曳力及重力,以及绝缘子壁面沉积/出射边界条件;
步骤(4)、开始仿真,软件自行进行网格剖分、迭代计算污秽颗粒的位置及速度,获得绝缘子表面污秽颗粒沉积情况;
步骤(5)、数据后处理,计算得到污秽质量密度、污秽不均匀度、积污带电系数。
优选地,步骤(2)中,计算域内的静电场稳态分布控制方程设为:
其中E为电场强度,单位V/m;D为电位移强度,单位C/m2;U为电位值,单位:伏;ε0为真空绝对介电常数,取8.85×10-12F/m;ε1为介质的相对介电常数;ρe为体电荷密度,单位C/m3;
采用RNG k-ε湍流模型将计算域内的流场稳态分布控制方程设为:
其中k为湍流动能,单位m2·s-2;ε湍流耗散率,单位m2·s-3;ρ为流体密度,单位kg·m-3;Gk为代表由平均速度梯度引起的湍动能项,kg·m-1·s-3;C1ε、C2ε为经验常数;αk为湍动能k的普朗特数,无量纲;αε为耗散率ε的普朗特数,无量纲;μeff为空气粘度和湍流粘度之和,单位Pa·s;ui、uj为平均速度分量;xi、xj为坐标分量。
优选地,步骤(2)中,在设置静电场及流体力学模块计算域边界条件时,绝缘子高压端电势与线路电压等级一致,绝缘子表面设为内壁面,且为糙度表面,其当量砂粒粗糙高度与污秽颗粒粒径一致,入口边界为水平气流速度入口,分别按照经验公式I=0.16(Re)-1/8和L=0.07ld确定气流的湍流强度以及湍流尺度,其中I为湍流强度,L为湍流尺度,ld为水力当量直径,Re为雷诺数;出口边界设置为自由出口;绝缘子表面边界设为无滑移壁面,采用标准壁面函数对近壁区进行处理,从而考虑了壁面边界层内高速度梯度的粘性影响,增加近壁区的求解精确度。
优选地,步骤(3)中,污秽颗粒在计算域均匀释放,带正电荷、带负电荷和中性污秽颗粒的浓度比设为31%、26%、43%,荷电量设置为:
Qp为污秽颗粒荷电量,单位C;E为污秽颗粒所在位置处的电场强度大小,V/m;εp为污秽颗粒的相对介电常数;dp为污秽颗粒粒径,单位μm。计算域中的污秽颗粒初始化时,污秽颗粒在计算域均匀释放,与大气中飞灰荷电情况保持一致。
优选地,步骤(3)中,考虑了污秽颗粒受到的重力、气流曳力、电场力的综合作用,污秽颗粒受力运动控制方程设置为:
其中m为污秽颗粒质量,Vp(t)为污秽颗粒的瞬时速度,Vb为污秽颗粒所处位置的气流速度,Fe、Fd、Fg依次为污秽颗粒所处空间位置的电场力、曳力及重力;E为电场强度,单位V/m;μ为动力粘度,单位1.8×10-5Pa·s;dp为污秽颗粒粒径,单位μm;ρp为污秽颗粒密度,单位kg/m3;g为重力加速度;εp为污秽颗粒的相对介电常数。对于流体流动粒子追踪模块,考虑了污秽颗粒受到的重力、气流曳力、电场力的综合作用,在三维物理场中,任一时间、空间下污秽颗粒受力运动控制方程均可设置为上述污秽颗粒受力运动控制方程。
优选地,步骤(3)中,考虑了污秽颗粒在壁面的沉积和出射动态微观过程:设VpT为污秽颗粒在绝缘子表面的切向速度,VpN为法向速度,eT和eN分别为绝缘子壁面的切向和法向单位矢量,t0表示为污秽颗粒运动至壁面的时刻,则在设置流体流动粒子追踪模块边界条件时,加入沉积/出射判据:
|VpN(t0)|≤VJ,其中
其中e表示颗粒弹性恢复系数,取e=0.5,无量纲;Ec为界面能,kg·m2/s2;dp为颗粒粒径,单位μm;ρp为污秽颗粒密度,单位mg/cm3;
如果污秽颗粒瞬时速度满足不等式,则污秽颗粒沉积,速度赋值为0;
如果未能满足不等式,则污秽颗粒与壁面分离,重新对粒子的切向和法向速度进行赋值:
VpT'=Vp(t0)·eT
其中的VpT’、VpN’依次为粒子与绝缘子壁面碰撞后的切向和法向出射速度。
优选地,通过改变仿真设置中的风速、风向、颗粒浓度以及颗粒粒径来模拟不同的环境参数。
进一步地,步骤(5)中,绝缘子表面污秽质量密度ρm(mg/cm2)计算公式为:
ρm=πdp 3·ρp·ND/6S总
其中,ND为绝缘子表面粘附的污秽颗粒数目;dp为颗粒粒径,单位μm;ρp为污秽颗粒密度,单位mg/cm3;S总为绝缘子表面积,cm2。
进一步地,污秽不均匀度的计算方法为:分区计算绝缘子表面不同位置处(如迎风侧、背风侧或上表面、下表面)的污秽质量密度,求取比值,得到相应的污秽不均匀度。
进一步地,积污带电系数的计算方法为:按照上述步骤依次求得带电、不带电情况下的绝缘子表面污秽质量密度结果,计算两者比值得到积污带电系数。
与现有技术相比,本发明的有益效果是:1)本发明充分考虑绝缘子积污过程中的气流场特性、直流电场分布、以及污秽颗粒积聚/出射动态微观过程,对现有的绝缘子积污仿真方法进行了改进;并利用COMSOL软件开展绝缘子表面污秽沉积仿真,实现了较为接近真实情况的直流架空线路绝缘子表面污秽沉积模拟;2)本发明一方面能够较好地再现不同环境参数下的绝缘子表面污层分布情况,另一方面能够计算得到绝缘子表面污秽质量密度、污秽不均匀度及积污带电系数;3)本发明能够在不消耗人力物力的情况下,为指导污区划分、清扫工作以及外绝缘裕度设计提供有力支撑。
附图说明
图1为本发明方法的流程示意图;
图2为本发明中污秽颗粒沉积/出射判断的流程示意图;
图3为具体实施方式中的绝缘子多物理场仿真模型的示意图;
图4(a)-(c)为具体实施方式中的实施例1的积污仿真结果及实测结果对比示意图;
图5(a)-(c)为具体实施方式中的实施例2的积污仿真结果及实测结果对比示意图;
图6(a)-(b)为具体实施方式中的实施例3的积污仿真结果及实测结果对比示意图;
图7(a)-(b)为具体实施方式中的实施例4的积污仿真结果及实测结果对比示意图;
图8为具体实施方式中的实施例5的绝缘子迎/背风侧分区示意图;
图9为具体实施方式中的实施例6的积污仿真结果及实测结果对比示意图;
图中:1、接地端;2、高压端;3、气流入口;4、细化网格区域;5、静电场计算域;6、流场计算域;7、气流出口。
具体实施方式
下面将结合本发明中的附图,对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动条件下所获得的所有其它实施例,都属于本发明保护的范围。
XP-160绝缘子为普通悬式绝缘子,常被用来作为自然积污和人工污秽试验的样品。故以三片XP-160绝缘子悬垂串为对象,阐述本发明内容。绝缘子结构参数如表1所示。根据目前我国±500kV的直流线路设计经验,绝缘子的耐受梯度需达到70kV/m,换算到XP-160绝缘子则单片至少需承受11kV的直流电压,故,以XP-160瓷绝缘子三片串带电+35kV的情况举例阐述。
表1绝缘子结构参数
本发明提出的直流架空线路绝缘子表面污秽沉积模拟方法的具体流程图如图1,主要步骤为:
(1)按照绝缘子尺寸参数及布置方式,在COMSOL软件内搭建仿真模型(如图3所示),划分计算域;
(2)采用软件自带静电场及流体力学模块,初始化计算域剖分网格,设置边界条件,迭代计算计算域内的静电场及流场稳态分布;
(3)采用软件自带流体流动粒子追踪模块,初始化计算域中的污秽颗粒,设置荷电量、电场力、气流曳力及重力,设置步长、时间、以及绝缘子壁面沉积/出射边界条件;
(4)开始仿真,软件自行进行网格剖分、迭代计算污秽颗粒的位置及速度,获得绝缘子表面污秽颗粒沉积情况。
(5)数据后处理,计算得到污秽质量密度、污秽不均匀度、积污带电系数。
所述步骤(2)中迭代计算计算域内的静电场及流场稳态分布,在具体实施中的步骤为:先采用全耦合求解器、共轭梯度迭代算法计算静电场稳态分布;再采用分离式求解器、GMRES迭代算法计算流场稳态分布;迭代计算时的相对容差设置为0.001;网格剖分时,绝缘子表面设置边界层网格,靠近绝缘子的区域采用自由剖分四面体细化网格,其他部分采用Cooper方法划分锲形/六面体网格。
所述步骤(4)中迭代计算污秽颗粒的位置及速度,在具体实施中的步骤为:剖分网格,绝缘子表面设置边界层网格,靠近绝缘子的区域采用自由剖分四面体细化网格,其他部分采用Cooper方法划分锲形/六面体网格;采用瞬态求解器、PARDISO直接耦合迭代算法;相对容差设置为0.001。
所述步骤(2)中,计算域内的静电场稳态分布控制方程为:
其中E为电场强度,单位V/m;D为电位移强度,单位C/m2;U为电位值,单位:伏;ε0为真空绝对介电常数,取8.85×10-12F/m;ε1为介质的相对介电常数;ρe为体电荷密度,单位C/m3。
流场稳态分布控制方程为:
其中k为湍流动能,单位m2·s-2;ε湍流耗散率,单位m2·s-3;ρ为流体密度,单位kg·m-3;Gk为代表由平均速度梯度引起的湍动能项,单位kg·m-1·s-3;C1ε、C2ε为经验常数;αk为湍动能k的普朗特数,无量纲;αε为耗散率ε的普朗特数,无量纲;μeff为空气粘度和湍流粘度之和,单位Pa·s;ui、uj为平均速度分量;xi、xj为坐标分量。
所述步骤(2)中的静电场及流场分布迭代计算,在设置边界条件时,绝缘子高压端电势与线路电压等级一致,绝缘子表面设为内壁面,且为糙度表面,其当量砂粒粗糙高度与污秽颗粒粒径一致,入口边界为水平气流速度入口,分别按照经验公式I=0.16(Re)-1/8和L=0.07ld确定气流的湍流强度以及湍流尺度,其中I为湍流强度,L为湍流尺度,ld为水力当量直径,Re为雷诺数;出口边界设置为自由出口;绝缘子表面边界设为无滑移壁面,采用标准壁面函数对近壁区进行处理。
所述步骤(3)中的流体流动粒子追踪模块,污秽颗粒受力运动方程设置为:
其中m为污秽颗粒质量,Vp(t)为污秽颗粒的瞬时速度,Vb为污秽颗粒所处位置的气流速度,Fe、Fd、Fg依次为污秽颗粒所处空间位置的电场力、曳力及重力;E为电场强度,单位V/m;μ为动力粘度,单位1.8×10-5Pa·s;dp为污秽颗粒粒径,单位μm;ρp为污秽颗粒密度,单位kg/m3;g为重力加速度;εp为污秽颗粒的相对介电常数。
所述步骤(3)中,计算域污秽颗粒初始化时,设置污秽颗粒在计算域均匀释放,与大气中飞灰荷电情况保持一致,带正电荷、带负电荷和中性污秽颗粒的浓度比设为31%、26%、43%。
所述步骤(3)中的污秽颗粒沉积/出射计算方法如图2所示:设t0时刻污秽颗粒与绝缘子表面发生碰撞,首先提取颗粒碰撞绝缘子表面的瞬时速度;继而提取绝缘子表面的切向、法向单位矢量eT和eN;计算得到污秽颗粒碰撞绝缘子表面时的切向速度VpT和法向速度VpN;判定如下不等式是否满足:
其中e表示颗粒弹性恢复系数,取e=0.5,无量纲;Ec为界面能,kg·m2/s2;dp为颗粒粒径,单位μm;ρp为污秽颗粒密度,单位mg/cm3。
如果满足上述不等式,则污秽颗粒沉积,速度赋值为0;如果未能满足上述不等式,则污秽颗粒与壁面分离,重新对颗粒的切向和法向速度进行赋值:
VpT'=Vp(t0)·eT
其中的VpT’、VpN’依次为粒子与绝缘子壁面碰撞后的切向和法向出射速度。
所述步骤(5)中的数据后处理,其中绝缘子表面污秽质量密度ρm(mg/cm2)计算方法为:
ρm=πdp 3·ρp·ND/6S总
其中,ND为绝缘子表面粘附的污秽颗粒数目,dp为颗粒粒径,单位μm;ρp为污秽颗粒密度,单位mg/cm3;S总为绝缘子表面积,cm2。
污秽不均匀度的计算方法为:分区计算绝缘子表面不同位置处(如迎风侧、背风侧或上表面、下表面)的污秽质量密度,求取比值,得到相应的污秽不均匀度。
积污带电系数的计算方法为:按照上述步骤依次求得带电、不带电情况下的绝缘子表面污秽质量密度结果,计算两者比值得到积污带电系数。
下面分多种情形阐述本发明方法的实施效果。
实施例1
污秽颗粒用二氧化硅模拟,固定风速5m/s、粒径15μm,污秽颗粒浓度为15mg/m3,绝缘子不带电。积污时间取8、16、24h,对比本发明方法在模拟不同积污时间上的效果。图4(a)-(c)表明积污时间分别为8、16、24h时绝缘子表面污层分布情况,其中,图4(a)中的(a-1)、(a-2)表示积污时间为8h时的仿真结果,图(a-1)中的箭头表示风向,(a-3)、(a-4)表示积污时间为8h时的实测结果;图4(b)中的(b-1)、(b-2)表示积污时间为16h时的仿真结果,(b-3)、(b-4)表示积污时间为16h时的实测结果;图4(c)中的(c-1)、(c-2)表示积污时间为24h时的仿真结果,(c-3)、(c-4)表示积污时间为24h时的实测结果。从图中可以看出,不同积污时间下绝缘子表面污秽分布外观没明显差异,均表现为迎风侧积污较轻,背风侧积污较重;随着时间的增加绝缘子迎/背风侧积污量都有所增长;本发明方法很好地体现了绝缘子积污随时间的增长,所得不同积污时间下污秽颗粒分布与实测结果吻合较好。
实施例2
污秽颗粒用二氧化硅模拟,固定粒径15μm,污秽颗粒浓度为15mg/m3,积污时间为16h。风速取1、2、5m/s三种情况,绝缘子带电(+35kV)和不带电两种情况,对比本发明方法在模拟不同风速、带电情况下的效果。图5(a)-(c)表明风速分别取1、2、5m/s时,绝缘子不带电和带电(+35kV)两种情况下仿真结果和实测结果的对比示意图,其中,图5(a)中的(a-1)表示风速1m/s时不带电情况下仿真结果和实测结果的对比,图(a-1)中箭头表示风向s,(a-2)表示风速1m/s时带电(+35kV)情况下仿真结果和实测结果的对比;图5(b)中的(b-1)表示风速2m/s时不带电情况下仿真结果和实测结果的对比,(b-2)表示风速2m/s时带电(+35kV)情况下仿真结果和实测结果的对比;图5(c)中的(c-1)表示风速5m/s时不带电情况下仿真结果和实测结果的对比,(c-2)表示风速5m/s时带电(+35kV)情况下仿真结果和实测结果的对比。图5可以看出,本发明方法可以很好地模拟不同风速下绝缘子表面直流积污现象,体现在:随着风速的增加,本发明方法所得绝缘子表面污秽度与实测结果都呈增加趋势,该结果一致;绝缘子下表面的污秽明显重于上表面,仿真结果和实测结果一致;在带电情况下,绝缘子下表面污层的增量更严重,仿真结果和实测结果一致;;迎/背风侧污秽呈扇状不均匀分布,风速越大,绝缘子背风侧与迎风侧的积污量差异更显著,仿真结果和实测结果也一致。
实施例3
污秽颗粒用二氧化硅模拟,污秽颗粒浓度为15mg/m3,积污时间为16h,固定风速5m/s,污秽颗粒粒径50μm,取绝缘子带电(+35kV)和不带电两种情况,对比本发明方法在模拟不同颗粒粒径、带电情况下的效果。图6(a)、6(b)表明在不同颗粒粒径、带电情况下的对比效果图,其中,图6(a)表示不同颗粒粒径在不带电情况下的对比效果图,图6(b)表示不同颗粒粒径在带电(+35kV)情况下的对比效果图。对比图5和图6可以看出,本发明方法可以很好地体现不同颗粒粒径下绝缘子表面直流积污现象,体现在:粒径增大后绝缘子表面污秽均有减少趋势;带电情况下绝缘子表面污秽显著增加,迎/背风侧积污差异也加重;随着粒径的增加,绝缘子下表面的迎/背风侧污秽分布不均匀性加重;本发明方法得到的污秽沉积仿真结果与实测现象吻合较好。
实施例4
污秽颗粒用二氧化硅模拟,污秽颗粒浓度为15mg/m3,积污时间为16h,取不带电的情况,对比本发明方法在计算绝缘子表面污秽质量密度方面的效果。图7(a)表示在不同风速下绝缘子表面污秽质量密度变化趋势与实测结果及传统方法的对比效果图;图7(b)表示在不同粒径下绝缘子表面污秽质量密度变化趋势与实测结果及传统方法的对比效果图。由图7(a)、7(b)可以看出,本发明方法得到绝缘子表面污秽质量密度计算值与实测值相对误差基本在25%以内,同时得到的不同风速、粒径下绝缘子污秽质量密度变化趋势与实测结果吻合较好;与传统的仿真方法相比,相对误差较小,体现出本发明方法的优越性。
实施例5
污秽颗粒用二氧化硅模拟,污秽颗粒浓度为15mg/m3,积污时间为16h,取绝缘子带电(+35kV)和不带电两种情况,以迎/背风侧不均匀积污为例,背风侧面积比取25%,对比本发明方法在计算绝缘子污秽不均匀度方面的效果。图8为绝缘子迎/背风侧分区示意图。通过数据后处理可以分别得到迎风侧、背风侧的污秽质量密度ρm_迎和ρm_背,然后根据下式求得迎/背风侧污秽不均匀度:
计算得到不同情况下的污秽不均匀度J值,如表2、表3所示:
表2J的计算结果(不带电)
表3 J的计算结果(带电+35kV)
由表中结果可以看出,本发明方法可以有效地计算绝缘子表面污秽不均匀度,并且能够反映出直流电场的影响。
实施例6
污秽颗粒用二氧化硅模拟,污秽颗粒浓度为15mg/m3,积污时间为16h,对比本发明方法在计算绝缘子表面积污带电系数的效果。绝缘子积污带电系数是目前特高压输电线路外绝缘设计及污区划分中需要关注的参量之一。根据带电、不带电情况下的绝缘子污秽质量密度求比值,计算得到绝缘子积污带电系数,如图9所示,并由此可以看出,本发明方法得到的绝缘子积污带电系数与实测结果吻合较好。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。
Claims (7)
1.一种直流架空线路绝缘子表面污秽沉积模拟方法,其特征在于,包括顺序相接的以下步骤:
步骤(1)、按照绝缘子尺寸参数及布置方式,在COMSOL软件内搭建仿真模型,划分计算域;
步骤(2)、采用软件自带静电场及流体力学模块,初始化计算域剖分网格,设置边界条件,迭代计算计算域内的静电场及流场稳态分布;
步骤(3)、采用软件自带流体流动粒子追踪模块,初始化计算域中的污秽颗粒,设置荷电量、电场力、气流曳力及重力,以及绝缘子壁面沉积/出射边界条件;
步骤(4)、开始仿真,软件自行进行网格剖分、迭代计算污秽颗粒的位置及速度,获得绝缘子表面污秽颗粒沉积情况;
步骤(5)、数据后处理,计算得到污秽质量密度、污秽不均匀度、积污带电系数。
2.如权利要求1所述的一种直流架空线路绝缘子表面污秽沉积模拟方法,其特征在于,步骤(2)中,计算域内的静电场稳态分布控制方程设为:
其中E为电场强度,单位V/m;D为电位移强度,单位C/m2;U为电位值,单位:伏;ε0为真空绝对介电常数,取8.85×10-12F/m;ε1为介质的相对介电常数;ρe为体电荷密度,单位C/m3;
采用RNG k-ε湍流模型将计算域内的流场稳态分布控制方程设为:
其中k为湍流动能,单位m2·s-2;ε湍流耗散率,单位m2·s-3;ρ为流体密度,单位kg·m-3;Gk为代表由平均速度梯度引起的湍动能项,单位kg·m-1·s-3;C1ε、C2ε为经验常数;αk为湍动能k的普朗特数,无量纲;αε为耗散率ε的普朗特数,无量纲;μeff为空气粘度和湍流粘度之和,单位Pa·s;ui、uj为平均速度分量;xi、xj为坐标分量。
3.如权利要求1所述的一种直流架空线路绝缘子表面污秽沉积模拟方法,其特征在于,步骤(2)中,在设置静电场及流体力学模块计算域边界条件时,绝缘子高压端电势与线路电压等级一致,绝缘子表面设为内壁面,且为糙度表面,其当量砂粒粗糙高度与污秽颗粒粒径一致,入口边界为水平气流速度入口,分别按照经验公式I=0.16(Re)-1/8和L=0.07ld确定气流的湍流强度以及湍流尺度,其中I为湍流强度,L为湍流尺度,ld为水力当量直径,Re为雷诺数;出口边界设置为自由出口;绝缘子表面边界设为无滑移壁面,采用标准壁面函数对近壁区进行处理,从而考虑了壁面边界层内高速度梯度的粘性影响,增加近壁区的求解精确度。
4.如权利要求1所述的一种直流架空线路绝缘子表面污秽沉积模拟方法,其特征在于,步骤(3)中,污秽颗粒在计算域均匀释放,带正电荷、带负电荷和中性污秽颗粒的浓度比设为31%、26%、43%,荷电量设置为:
其中Qp为污秽颗粒荷电量,单位C;E为污秽颗粒所在位置处的电场强度大小,单位V/m;εp为污秽颗粒的相对介电常数;dp为污秽颗粒粒径,单位μm。
5.如权利要求1所述的一种直流架空线路绝缘子表面污秽沉积模拟方法,其特征在于,步骤(3)中,考虑了污秽颗粒受到的重力、气流曳力、电场力的综合作用,污秽颗粒受力运动控制方程设置为:
其中m为污秽颗粒质量,Vp(t)为污秽颗粒的瞬时速度,Vb为污秽颗粒所处位置的气流速度,Fe、Fd、Fg依次为污秽颗粒所处空间位置的电场力、曳力及重力;E为电场强度,单位V/m;μ为动力粘度,单位1.8×10-5Pa·s;dp为污秽颗粒粒径,单位μm;ρp为污秽颗粒密度,单位kg/m3;g为重力加速度;εp为污秽颗粒的相对介电常数。
6.如权利要求1所述的一种直流架空线路绝缘子表面污秽沉积模拟方法,其特征在于,步骤(3)中,考虑了污秽颗粒在壁面的沉积和出射动态微观过程:设VpT为污秽颗粒在绝缘子表面的切向速度,VpN为法向速度,eT和eN分别为绝缘子壁面的切向和法向单位矢量,t0表示为污秽颗粒运动至壁面的时刻,则在设置流体流动粒子追踪模块边界条件时,加入沉积/出射判据:
|VpN(t0)|≤VJ,其中VpN(t0)=Vp(t0)·eN,
其中e表示颗粒弹性恢复系数,取e=0.5,无量纲;Ec为界面能,单位kg·m2/s2;dp为颗粒粒径,单位μm;ρp为污秽颗粒密度,单位mg/cm3;
如果污秽颗粒瞬时速度满足不等式,则污秽颗粒沉积,速度赋值为0;
如果未能满足不等式,则污秽颗粒与壁面分离,重新对粒子的切向和法向速度进行赋值:
VpT'=Vp(t0)·eT
其中的VpT’、VpN’依次为粒子与绝缘子壁面碰撞后的切向和法向出射速度。
7.如权利要求1所述的一种直流架空线路绝缘子表面污秽沉积模拟方法,其特征在于,通过改变仿真设置中的风速、风向、颗粒浓度以及颗粒粒径来模拟不同的环境参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811396424.9A CN109583066B (zh) | 2018-11-22 | 2018-11-22 | 一种直流架空线路绝缘子表面污秽沉积模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811396424.9A CN109583066B (zh) | 2018-11-22 | 2018-11-22 | 一种直流架空线路绝缘子表面污秽沉积模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109583066A true CN109583066A (zh) | 2019-04-05 |
CN109583066B CN109583066B (zh) | 2023-06-09 |
Family
ID=65923669
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811396424.9A Active CN109583066B (zh) | 2018-11-22 | 2018-11-22 | 一种直流架空线路绝缘子表面污秽沉积模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109583066B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110489864A (zh) * | 2019-08-20 | 2019-11-22 | 国网天津市电力公司电力科学研究院 | 计及风速、污秽、粒径的防湿雪绝缘子防污特性分析方法 |
CN110489867A (zh) * | 2019-08-20 | 2019-11-22 | 国网天津市电力公司电力科学研究院 | 一种防湿雪复合绝缘子的积污特性分析方法 |
CN110516351A (zh) * | 2019-08-27 | 2019-11-29 | 天津大学 | 一种基于绝缘子积污特性的双串绝缘子悬挂角度优化方法 |
CN110728018A (zh) * | 2019-08-27 | 2020-01-24 | 天津大学 | 一种不同空气湿度下大伞裙复合绝缘子积污特性分析方法 |
CN111144012A (zh) * | 2019-12-28 | 2020-05-12 | 中汽研汽车检验中心(天津)有限公司 | 一种用于冷空间内冰粒沉积过程的计算方法 |
CN111241662A (zh) * | 2020-01-02 | 2020-06-05 | 西安交通大学 | 一种仿生波动鳍的推进性能数值计算方法 |
CN111475974A (zh) * | 2020-03-06 | 2020-07-31 | 华北电力大学(保定) | 输电线路中绝缘子串污秽检测模拟方法、装置及存储介质 |
CN112464543A (zh) * | 2021-01-28 | 2021-03-09 | 北京科技大学 | 一种计算vim冶炼过程中的夹杂物运动的方法 |
CN113311300A (zh) * | 2021-05-27 | 2021-08-27 | 华北电力大学(保定) | 一种积污饱和特性试验平台 |
CN113567816A (zh) * | 2021-07-29 | 2021-10-29 | 国家电网有限公司 | 绝缘子污秽测量测量方法 |
CN113947040A (zh) * | 2021-10-12 | 2022-01-18 | 中铁第一勘察设计院集团有限公司 | 基于体积分布系数的绝缘子上下表面不均匀积污分析方法 |
CN115809491A (zh) * | 2023-02-01 | 2023-03-17 | 北京适创科技有限公司 | 一种核电领域颗粒物的沉积分析方法及装置 |
CN116068287A (zh) * | 2023-03-10 | 2023-05-05 | 国网山西省电力公司电力科学研究院 | 一种污秽绝缘子电位测量方法及装置 |
CN116258021A (zh) * | 2023-05-10 | 2023-06-13 | 大唐东北电力试验研究院有限公司 | 基于comsol的光伏板表面积灰模拟分析方法及系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108009130A (zh) * | 2017-11-30 | 2018-05-08 | 国家电网公司 | 基于回归分析的自然环境中绝缘子积污预测方法 |
-
2018
- 2018-11-22 CN CN201811396424.9A patent/CN109583066B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108009130A (zh) * | 2017-11-30 | 2018-05-08 | 国家电网公司 | 基于回归分析的自然环境中绝缘子积污预测方法 |
Non-Patent Citations (1)
Title |
---|
陈杭等: "XHP-160 绝缘子积污特性仿真分析及试验验证", 《电网技术》 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110489867B (zh) * | 2019-08-20 | 2023-05-26 | 国网天津市电力公司电力科学研究院 | 一种防湿雪复合绝缘子的积污特性分析方法 |
CN110489867A (zh) * | 2019-08-20 | 2019-11-22 | 国网天津市电力公司电力科学研究院 | 一种防湿雪复合绝缘子的积污特性分析方法 |
CN110489864A (zh) * | 2019-08-20 | 2019-11-22 | 国网天津市电力公司电力科学研究院 | 计及风速、污秽、粒径的防湿雪绝缘子防污特性分析方法 |
CN110516351A (zh) * | 2019-08-27 | 2019-11-29 | 天津大学 | 一种基于绝缘子积污特性的双串绝缘子悬挂角度优化方法 |
CN110728018A (zh) * | 2019-08-27 | 2020-01-24 | 天津大学 | 一种不同空气湿度下大伞裙复合绝缘子积污特性分析方法 |
CN111144012A (zh) * | 2019-12-28 | 2020-05-12 | 中汽研汽车检验中心(天津)有限公司 | 一种用于冷空间内冰粒沉积过程的计算方法 |
CN111144012B (zh) * | 2019-12-28 | 2023-12-19 | 中汽研汽车检验中心(天津)有限公司 | 一种用于冷空间内冰粒沉积过程的计算方法 |
CN111241662A (zh) * | 2020-01-02 | 2020-06-05 | 西安交通大学 | 一种仿生波动鳍的推进性能数值计算方法 |
CN111475974A (zh) * | 2020-03-06 | 2020-07-31 | 华北电力大学(保定) | 输电线路中绝缘子串污秽检测模拟方法、装置及存储介质 |
CN112464543B (zh) * | 2021-01-28 | 2021-04-06 | 北京科技大学 | 一种计算vim冶炼过程中的夹杂物运动的方法 |
CN112464543A (zh) * | 2021-01-28 | 2021-03-09 | 北京科技大学 | 一种计算vim冶炼过程中的夹杂物运动的方法 |
CN113311300B (zh) * | 2021-05-27 | 2023-10-20 | 华北电力大学(保定) | 一种积污饱和特性试验平台 |
CN113311300A (zh) * | 2021-05-27 | 2021-08-27 | 华北电力大学(保定) | 一种积污饱和特性试验平台 |
CN113567816A (zh) * | 2021-07-29 | 2021-10-29 | 国家电网有限公司 | 绝缘子污秽测量测量方法 |
CN113947040A (zh) * | 2021-10-12 | 2022-01-18 | 中铁第一勘察设计院集团有限公司 | 基于体积分布系数的绝缘子上下表面不均匀积污分析方法 |
CN113947040B (zh) * | 2021-10-12 | 2024-06-25 | 中铁第一勘察设计院集团有限公司 | 基于体积分布系数的绝缘子上下表面不均匀积污分析方法 |
CN115809491A (zh) * | 2023-02-01 | 2023-03-17 | 北京适创科技有限公司 | 一种核电领域颗粒物的沉积分析方法及装置 |
CN115809491B (zh) * | 2023-02-01 | 2023-04-21 | 北京适创科技有限公司 | 一种核电领域颗粒物的沉积分析方法及装置 |
CN116068287A (zh) * | 2023-03-10 | 2023-05-05 | 国网山西省电力公司电力科学研究院 | 一种污秽绝缘子电位测量方法及装置 |
CN116258021B (zh) * | 2023-05-10 | 2023-09-19 | 大唐东北电力试验研究院有限公司 | 基于comsol的光伏板表面积灰模拟分析方法及系统 |
CN116258021A (zh) * | 2023-05-10 | 2023-06-13 | 大唐东北电力试验研究院有限公司 | 基于comsol的光伏板表面积灰模拟分析方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN109583066B (zh) | 2023-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109583066A (zh) | 一种直流架空线路绝缘子表面污秽沉积模拟方法 | |
Zhang et al. | Study on natural contamination performance of typical types of insulators | |
Lv et al. | Simulation study on pollution accumulation characteristics of XP13-160 porcelain suspension disc insulators | |
CN109598065A (zh) | 一种流场环境下绝缘子带电积污分布的获取方法 | |
CN110728018A (zh) | 一种不同空气湿度下大伞裙复合绝缘子积污特性分析方法 | |
CN102540027B (zh) | 一种沙漠干旱区域高压输电线路绝缘子选型方法 | |
Jiang et al. | Dust contamination on surface of transmission line insulators in air-polluted regions in China: statistical characteristics, adhesion mechanism, and environmental impact factors | |
CN110489864A (zh) | 计及风速、污秽、粒径的防湿雪绝缘子防污特性分析方法 | |
CN111665189B (zh) | 车顶绝缘子表面涂层磨蚀试验方法 | |
Huang et al. | Flow field distribution around insulator and contamination uneven characteristic of insulator | |
Lv et al. | Simulation on deposition characteristics of contamination particles in fog-haze environment | |
Lv et al. | Model optimization of pollution deposition criterion of insulators based on energy | |
CN103399980B (zh) | 一种输变电绝缘子的优化选型配置方法 | |
He et al. | Characteristics of dust deposition on suspended insulators during simulated sandstorm | |
CN110489867B (zh) | 一种防湿雪复合绝缘子的积污特性分析方法 | |
Liu et al. | Numerical investigation on collision of pollution particles on outdoor insulators | |
CN108170894A (zh) | 一种输变电设备外绝缘体大气颗粒物干沉降参数化计算方法 | |
Lv et al. | Numerical Simulation of Contamination Accumulation Characteristics of Composite Insulators in Salt Fog Environment | |
CN110096777B (zh) | 静电雾化铣削雾滴输运建模与输运效果评价方法 | |
CN109992890B (zh) | 高速列车-龙卷风的耦合模型及耦合计算方法 | |
CN110516351A (zh) | 一种基于绝缘子积污特性的双串绝缘子悬挂角度优化方法 | |
Othman et al. | Simulation of 33kV String Insulators Using Finite Element Method (FEM) | |
Dong et al. | Simulation Analysis on Flow Field and Collision Characteristics of Contaminated Solid Particles on Outdoor Insulators | |
Huang et al. | An integrated fouling model of insulator surface considering particle deposition and removal | |
Kong et al. | Effects of Air Humidity and Hanging Angle on Accumulation Characteristics of Pollution Particles on Anti-icing Polymer Insulators |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |