气液活塞泵数值计算的修正特征线法Word格式文档下载.docx
- 文档编号:609139
- 上传时间:2023-04-29
- 格式:DOCX
- 页数:5
- 大小:19.09KB
气液活塞泵数值计算的修正特征线法Word格式文档下载.docx
《气液活塞泵数值计算的修正特征线法Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《气液活塞泵数值计算的修正特征线法Word格式文档下载.docx(5页珍藏版)》请在冰点文库上搜索。
1修正特征线法
特征型方程组的推导标准的非定常流特征型方程
由于α是与计算时步有关的修正系数,故称“时间修正系数”,其取值的影响因素十分复杂,目前还未进行严格的误差分析,本文通过广泛的数值试验,并与大量现场试验进行对比发现,当α≤100时,其压力的相对误差一般不超过定性结论是,在满足数值稳定的前提下,α越大计算速度提高越显着,但精度越差;
α越小则反之,α=1时该法即为常规特征线法.一般情况下,快速瞬变系统的α较小,采用该法不能节约多少时间;
而短管或缓慢瞬变系统的α较大,可显着节约计算时间.本文的做法是根据不同的课题进行数值试验,先取大一点的值再逐步减小,直至获得满意的结果.
单介质管路定常振荡如上所述,对单介质管路系统,无论是慢速或快速瞬变,均可根据需要确定时间步长Δt,然后用α对运动方程进行修正.在分析波速变动很大的流体(如气流或含气液流)中,通过在不同计算时段采用不同的α,从而使计算保持相等的时间步长.
作为非定常流的特例我们讨论定常振荡流,在所有的计算节点中,上、下游边界点上的变量代表系统的输入和输出,是我们计算的重点,而对其内部节点的变量究竟如何传递的问题可不予考虑.假定计算管段可等分为m段(每段为Δx),振荡周期T可等分为n段(每段为Δt);
在m与n之间建立某种联系,使之与α有关,从而建立时段mΔt内,上下游边界点之间的变量传递关系.由于特征线C+、C-穿过不止一段(m段)管段,因此也无须计算管道内节点,只对我们所关心的边界点进行计算.具体做法是
令:
Δt=T/nΔx=L/mα=aπ/2iωL≥1(6)
式中:
n,m——与i有关的正整数;
i——α的计算整数.其值影响数值稳定和精度;
计算时
采用α≤100进行控制.
根据特征线的定义又有
α=aΔt/Δx(7)
由此可以证明式(6)中的整数满足
m/n=1/2i+1(8)
含时间修正因子α的周期性振荡气流特征型偏微分方程为
S=(2gΔxmsinθ)/a2.式(11)、(12)中下标为计算管段上、下游节点编号,上标为计算时段编号,其中j表示计算时段所求参数,j-m表示前m个时段的已知参数.由于摩擦项采用了二阶估算,故需用迭代法求解.
以下举例说明采用上述方程求解的具体情况.一个正弦型振荡周期至少需9个点描述,即n=8;
假设i=0即α=aT/2L可满足精度要求,则计算管段应等分为m=4段,计算网格如图1.图中某边界点经C+和C-两次特征线的积分,便可获得对面边界点的信息,所需时间正好为一个周期(8Δt),即计算本周期的上游边界上各时段的变量只需从前一周期相应的时段参数进行两次C+,C-积分即得.
图1特征线法计算网格
修正特征线法可详细模拟系统从静止到定常振荡的变化过程,流场参数的变化仅限于振荡幅值的变化,振荡频率与边界点强迫函数相同,相位关系也固定不变.
多介质段管路系统上述方法是针对单介质系统和固定边界而提出的,对于具有不同介质段的管路系统,例如,气液活塞泵内的脉冲运动,由于气体和液体的密度和声速均不相同,造成计算区域内不同介质段特征线斜率的不同,由于α的引入,通过计算不同管段及不同时段上的介质密度ρ和波速α,利用式(7)不仅可对不同的计算时段调整α,而且在不同的介质段也可采用不同的α,以便在整个计算区域内仍可获得相同的特征线斜率,保持空间步长和时间步长均相等的矩形计算网格(如图1).但由于各网格单元的a、α均不相同,所以不能采取上述单介质系统中,跨越内节点的计算方法,而须逐个对网格的C+、C-进行积分.将式(11)、(12)中的m换成1即为多介质段管路系统的特阵型差分方程.
2方法应用与验证
气液活塞泵简介图2是一种用于特殊场合的气液活塞泵,其主体为气液活塞筒,它与射流泵一起构成了一组利用气液脉冲运动进行传能的装置.该装置具有结构简单,无运动部件,免维修,可靠性高,工作性能好等优点,在不适合近距离操作,例如高温、高压、高放射性、剧毒等工作环境中具有十分重要的实用价值.其工作原理是将脉冲气体作用于活塞筒内液体,经过周期性的压冲排液和反吸充液,以脉冲形式向射流泵输出液体,为之提供动力水源.
(a)压冲排液阶段
(b)反吸充液阶段
图2气液活塞泵传能装置
脉冲运动包括周期性的压冲排液和反吸充液两过程.以一个脉冲周期为例,压冲阶段(td),正压气体将液体从A压至B;
反吸分两阶段(ts=ts1+ts2),第一阶段(ts1)由于液体的惯性,液位继续下移至C,随着活塞泵内负压的形成以及惯性力的减弱进入第二阶段(ts2),泵内液体在吸水箱液位的作用下开始上移至A点,完成一个周期(T=td+ts).应用气液活塞泵作为脉冲发生及传能装置,关键的问题是使液位保持在活塞泵筒内某一固定范围内变化,即保持稳定.其次是系统的优化设计,提高装置效率.进行气液活塞传能机理的研究是装置稳定运行参数和优化结构设计的重要依据,本文采用修正特征线法对活塞泵内气、液流场进行了数值模拟.
计算假定及边界条件
计算假定如图2,将研究的范围定为活塞筒的a-a断面至b-b断面.进行下列假设
假定系统设计可满足稳定振荡要求;
(1)假定流动是等温的一维定常振荡流;
(2)气液活塞泵筒壁的膨胀忽略不计;
(3)气液交界面始终为一平面,在平面上气、液体积流量近似相等;
(4)气液交界面压力传递损失可忽略,界面上气体和液体压力相等;
(5)假设摩阻项与气、液流量及脉冲频率的平方成正比.
边界条件
(1)上游边界条件,本例中脉冲气体发生器是由一对分别产生负压和正压的真空喷射器及压缩喷射器组成,因此上游边界条件在反吸阶段(ts)和压冲阶段(td)应分别满足不同喷射器的特性曲线方程,加上C-方程(12)可求解上游边界的,p.
(2)下游边界条件,下游主要分两个阶段,即充液蓄能阶段(ts1+ts2)和输送液体阶段(td),边界条件由下列方程组给出
A0—液体射流泵喷嘴出口断面面积;
μ—与脉冲液体出流段型式及频率有关的流量系数;
p0—活塞泵出流孔口处的绝对压力,与吸水箱液位有关;
下标b—表示下游边界点.
将下游边界条件(13)与C+方程(11)联立,即可求解下游节点的,p.
(3)内边界条件,在活塞筒内气液交界面是一内边界条件,由假定条件(3)可将内边界条件写为
计算曲线与试验结果的对比如图3,经过试验结果和计算结果的对比,发现采用本方法对定常振荡的振幅、周期及相位的计算是较准确的,尤其是在压冲排液阶段,计算得到波峰和波形与试验几乎完全吻合(图3b).
(a)活塞筒进口气体流量
(b)活塞筒出口液体流量
(c)活塞筒进口气体压力
(d)活塞筒出口液体压力
图3计算与试验对比(——计算--试验)
反吸充液阶段的计算结果不令人十分满意,图3b的反吸流量偏小,且出现时间稍晚.主要原因是假定气液交界面为平面,采用一维处理的计算方法过于简化,没有充分考虑反吸过程初期,倒流液体进入活塞筒突扩断面后,由于瞬时射流和边界层脱流而形成的大量旋滚消耗了做功能量,其次过高估计反向流动的非定常摩阻损失也可能是原因之一,这些因素使实际阻力比预计情况高,因此反吸流量偏小.
活塞筒进、出口压力试验曲线在反吸阶段振荡剧烈,而计算结果由于未考虑二维湍流的影响,没有出现实际上存在的压力脉冲,但作者认为这一因素反而起到了滤波作用,反映出反吸压力的总体变化规律,为我们分析系统提供了方便.
结论
(1)修正特征线法通过引入“时间修正因子”α增加时间步长,在保留了原特征线法优点的同时,可不同程度地提高计算速度;
(2)采用不同的α,可对变波速及多介质段的非定常流计算节点上的特征线斜率进行修正,避免产生特征线交点不确定的自由浮动网格,在处理气液活塞运动时十分灵活方便,成为解决各类非定常流尤其是定常振荡流数值模拟的有效工具.(3)修正特征线法是一种近似的方法,α的取值越大其产生的计算误差也越大,一般α不应超过100,否则易造成发散,且计算精度将受到影响.至于α与计算误差的关系,还需进一步研究.(4)在数值计算中,对管路摩阻特性的模拟还需进一步研究,否则不仅会影响计算结果的精确性,而且有时还会影响问题的收敛.(5)采用本方法对气液活塞式脉冲传能装置进行数值模拟,不仅可对系统的稳定性进行计算,还可充分详细地模拟启动过渡过程,并可通过带入不同的活塞筒几何尺寸、气源压力、脉冲时间等参数,对装置的结构设计和运行参数进行合理性检验,这些对该装置的研究和开发均具有十分重要的意义.
参考文献
[1]怀利EB,斯特里特VL.瞬变流[C].清华大学流体传动与控制教研组译,北京:
水利电力出版社,
索科洛夫ЕЯ,津格尔НМ.喷射器[C].黄秋云译,北京:
科学出版社,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 活塞 数值 计算 修正 特征