智能式断路器灭弧室气流场计算的一种新思路

1、智能式断路器灭弧室气流场计算的一种新思路作者:林莘钟建英转贴自:电源技术应用更新时间:2006-2-14摘要:智能式断路器灭弧室存在着压气室、膨胀室、喷口、动触头、静触头等复杂结构,而且多为不规则形状,使得要准确地对灭弧室的结构建立几何模型是一件比较困难的事情。而几何建模的准确与否关系到能否生成高质量的网格,网格生成技术又是得到精确结果的关键技术之一。所以,针对复杂的几何形状建模困难的情况,该文采用区域扩充法利用PHOENICS计算流体力学软件与自编程相结合的方式对智能式断路器灭弧室内喷口部分的气流场进行了计算,并与其他计算方法取得的结果做了比较,结果证明文中所尝试的方法是有效可行的,为智能式

2、断路器灭弧室气流场的计算提供了一种新思路。关键词:智能式断路器;区域扩充法;计算流体力学软件;气流场计算1 引言智能式断路器灭弧室结构非常复杂,存在着压气室、膨胀室、喷口、动触头、静触头等,而且多为不规则形状,因而所涉及到的气流场的边界条件比较复杂。所以在做几何建模的过程中,要准确地对灭弧室的结构建立几何模型是一件比较困难的事情。而几何建模的准确与否关系到能否生成高质量的网格。网格生成所需要的人力时间占一个全部计算任务时间的60%左右甚至更长,因此网格生成技术是得到精确结果的关键技术之一。在流场计算中,一般情况下采用几种方法来处理复杂区域内的换热和流动,如阶梯型网格、区域扩充法、三角形网格、贴

3、体坐标、以及坐标组合法等。本文认为区域扩充法是流体力学中采用的最多的一种方法,技术比较成熟,特点简单明了,而其它的方法都需要在网格划分上花费很大的精力,甚至还要进行坐标变换,因此可能要引入较多的计算误差。采用区域扩充法在一定程度上会延长计算机计算时间,增加存储量,但是,与用规则的网格来处理任意几何形状的计算区域所带来的好处相比是可以不考虑这一缺点的。所以在本文中利用区域扩充法或者通俗的说是挖洞法来对灭弧室复杂结构进行模拟生成计算网格,从而对其中的气流场进行计算。目前有很多的专门流体计算软件出现,PHOENICS是世界上第一个投放市场的CFD商用软件(1981,由于该软件投放市场较早,所以在工业

4、界得到广泛的应用。例如,空气动力学、电子器件冷却、喷嘴中的流动等等。早期的PHOENICS 在开发时受到基本框架的限制,在人机界面上不很灵活。目前,这个软件虽然在功能与方法上作了较大的改进,但是在建立形状复杂的几何模型的前处理方面仍然有一定的欠缺。另外,智能式断路器的灭弧原理是利用电弧堵塞现象使用电弧自身的能量来熄灭电弧的,所以,对喷口气流场进行分析是一个重要的课题,可以为设计出最佳喷口形状提供理论依据。本文率先采用区域扩充法,利用PHOENICS计算流体力学软件与自编程相结合的方式对空载情况下智能式SF6高压断路器灭弧室内喷口部分的气流场进行了计算,得到了较好的计算结果。2 灭弧室喷口气流场

5、数学模型本文中所用的气流场数学模型是可压缩流体动力学基本方程组。连续方程能量方程式中 Vz 为轴向速度分量;Vr为径向速度分量;为密度;P为气压;T为温度;G p 、分别为定压比热、粘性系数、导热系数;而r、z、t分别为径向坐标、轴向坐标和时间坐标;为修正系数。3 计算方法简介高压断路器在开断的过程中有电弧的产生,所以灭弧室的气流场中涉及到传热的问题。为了充分体现数学模型所表现的能量守恒原理,计算的基础是以单位控制容积为研究对象的。如图1所示,J代表一个典型变量的流量密度。dx,dy及dz为控制容积长、宽、高,J在x 方向的分量Jx是进入面流量密度,在x方向上,控制体流出的净流量为dxdydz

6、。同样可以得到y和z方向的相应流量。S为源项,代表着所有不能归入名义上的扩散项的因子或项;为相对应的扩散系数。式(6就是本计算的理论基础,称为通用微分方程。方程中的四项分别为非稳定项、对流项、扩散项以及源项。因变量可以代表各种不同的物理量,如T,Vz ,Vr,K,等等。式(1(4都可以写成式(6的形式。本文中的区域扩充法就是利用此方程的物理意义和特点对灭弧室中喷口内的气流场进行模拟计算。在计算时,把不规则的计算区域扩充为规则的网格,使得规则网格中的某些控制容积不起作用,让起作用的控制容积构成所需要求解的不规则计算区域。所以,在本文中,计算区域包括流体与固体两部分,并以壁的外表面作为边界。在求解

7、时,取流体的粘度作为流体区域内网格点的扩散系数值,而在固体区域内的网格点的扩散值取一个非常大的数。在给处于固体边界上的网格节点赋值时,源项通常是造成迭代发散的原因,所以需要对源项进行处理,以便使边界上的节点取给定的值。而对源项线性化是可以达到收敛解的关键,因此通常源项可以写成下面的形式式中SC 为与P有关部分的系数。如果在边界上的一个点已知值为B,则利用下列方式就可以把该已知值赋予该点。即设SC =1030B,则SP=-1030。这里1030表示一个大到足以使方程(6离散化后的方程中的其他各项可被忽略的数值。则离散化方程就变为当通用的变量P 是分别代表Vr,Vz,T时,将式(9中的B分别取给定

8、的边界值。4 计算结果与分析智能式SF6断路器灭弧室结构示意图如图2所示。根据上述方法来计算空载情况下智能式SF6高压断路器灭弧室内气流场的情况。在断路器分断操作时,膨胀室内的气体通过喷口流出灭弧室,喷口内流动参数的变化会直接影响到膨胀室内气体状态参数的变化5-8,因而喷口内流动参数的变化对于整个灭弧室内的气流场的情况是十分重要的。本文主要对内外喷口内气流场的情况进行分析。灭弧室内气体基压为0.5MPa,初始温度为300K。场域内各点的气流速度为零。断路器的行程运动曲线如图3所示。计算时膨胀室气流入口速度为气缸运动速度,出口背压为初始压力。由于在开断过程中动触头是运动的,所以涉及到移动边界的问

9、题。本文在整个行程的气流场求解过程中,将动触头的运动看成是非常多的小的行程段组成。在每个小的行程段中,可以认为动静触头位置相对静止,同时将上一个行程段的计算结果作为下一行程段的初值,直至动触头运动结束。在图4和图5中详细地说明利用区域扩充法建立的几何模型,图4为灭弧室喷口部分场域剖分示意图,图中用粗的黑实线表示这些区域内的网格点,其扩散值取一个非常大的数使其成为不活动的区域;图5是局部放大的固体区域处理。被颜色填充的部分为固体区域,其中对于不完全为固体区域的规则网格(如图中的网格1、网格2在编程时乘上一个不大于1的系数来处理,从而可以模拟不规则的形状边界。利用上述方法得到的计算结果如下:由图6

10、可以看出,此时外喷口还处于被静触头阻塞的状态,是一个流动的死区,流速为零,在内喷口中流体在上游区被加速,在喉部区域马赫数达到0.47,在下游喷口截面的扩张区域马赫数迅速降低到0.28。在60%开距时喷口中马赫数的变化以及压力、密度状态图如图7、8、9所示。y表示径向方向,z表示轴向方向,不同的颜色表示不同的数值,由于灭弧室内气流场是关于轴对称的,所以图形取关于轴对称的气流场的一半。由图7可以看出,流体在外喷口的上游区被加速,马赫数在喉部处达到了最大值0.42,同时由于静触头还有阻塞作用,流动截面变小,在此之前的区域流体被加速,此后的区域气流的马赫数迅速降低到0.16;在内喷口喉部区域的马赫数为

11、0.57,下游区降到0.38。同时,由于在内外喷口之间有气流向相反方向吹拂,减弱了气吹作用,所以,在这个位置气流马赫数较低。由图8、9可以看出压力与密度的变化趋势是一致的,随着喷口截面积的减小,压力与密度逐渐减小,在喷口喉部处达到最小值,之后又开始上升,最后下游区的大部分都恢复为最初的状态。这种状态与上面的流速相对应,可以得到喷口中流体状态参数的变化总是与流速的变化相反。而且触头的堵塞引起喷口中流动参数的变化。可以看出,上面的计算结果是符合灭弧室内喷口气流场的变化原理的,说明利用本文提出的方法进行计算是可行的,为了定量说明本文的方法合理正确,本文在相同的条件下利用另外一种常用的几何建模方法¾¾贴体坐标法计算了自能式断路器喷口内气流场的情况9,并把在不同开距下喷口内马赫数的分布情况做了比较,如图10所示。图10是n y=30的内外喷口中马赫数的分布(n y为y方向上的剖分网格数,横坐标表示开距,纵坐标表示马赫数,横坐标中的零点表示内外喷口之间的中心位置,曲线1、3分别表示利用贴体坐标法得到的在开距为20%和喷口完全打开后喷口内的马赫数分布,曲线2、4分别表示利用本文方法得到的在开距为20%和喷口完全打开后喷口内的马赫数分布。从图中可以看到,两种方法得到的结果非常接近,所以从数值计算上来看本文提出的方法也是可行的。5 结论由计算结果可以看出,利用区