气液活塞泵数值计算的修正特征线法.docxVIP

气液活塞泵数值计算的修正特征线 摘要 本文研究了气液活塞泵的运 动机理，提出了一种用于定常振荡流数值计 算的“修正特征线法”.在常规特征线法的 基础上，通过将“时间修正因子”引入运动 方程，不仅统一了变波速多介质段计算节点 的特征线斜率，保证计算在矩形网格上进行， 而且还在一定程度上提高了计算速度.文章 推导了修正特征线法的特征型微分方程及 差分方程：针对单介质管路和多介质管路分 别提出了不同的处理方法； 并在最后给出了 气液活塞泵定常振荡流计算实例，与试验结 果的对比一方面验证了 “修正特征线法” 计算变波速多介质段定常振荡流问题的优 越性，另一方面也提出该方法在精度分析和 敛散性方面尚需进一步研究. 关键词修正特征线法定常振荡流气液活 塞泵 气液活塞泵是一种免维修流体输送机 械，由于气体和液体在泵内作类似固体活塞 的往复运动，故称气液活塞泵.气液活塞运 动属变波速两介质段定常振荡运动， 目前普 遍认为解决定常振荡问题最经济的方法是 传递矩阵法，但对振幅较大的非线性定常振 荡系统，传递矩阵法的线性化假定存在较大 误差.特征线法在这方面有所改进，它对振 幅大小无限制，并可对摩擦项采用二阶估算 形式.因此在处理非线性度较强的瞬变流包 括定常振荡流中，特征线法具有重要的应用 价值.但是，特征线法在处理缓慢瞬变和短 管瞬变时，其计算速度一般较慢.气液活塞 泵内部存在着不同的介质区段， 其密度和波 速均不相同，若采用特征线法直接求解方程， 将出现一种不规则的浮动网格，不能直接获 得特定瞬间和特定截面上的信息.本文针对 上述特点，将常规特征线法进行改进，利用 “时间修正因子” a增加时间步长，提高计 算速度.该法不仅可用于计算单介质、固定 波速的管路非定常流，而且对于变密度、变 波速的气液管路定常振荡运动同样适用， 通 过在不同的时空区段引入不同的 a ,保持矩 形计算网格的稳定，在达到计算精度的前提 下，简化了计算程序.“修正特征线法”拓 展了特征线法在工程中的应用，本文重点讨 论其在气液定常振荡流中的应用^ 1修正特征线法 特征型方程组的推导标准的非定常流特征 型方程 由于a是与计算时步有关的修正系数， 故称“时间修正系数”，其取值的影响因素 十分复杂，目前还未进行严格的误差分析， 本文通过广泛的数值试验， 并与大量现场试 验进行对比发现，当 aW 100时，其压力的 相对误差一般不超过定性结论是，在满足数 值稳定的前提下，a越大计算速度提高越显 着，但精度越差；a越小则反之，a =1时 该法即为常规特征线法.一般情况下，快速 瞬变系统的a较小，采用该法不能节约多 少时间；而短管或缓慢瞬变系统的 a较大， 可显着节约计算时间.本文的做法是根据不 同的课题进行数值试验， 先取大一点的值再 逐步减小，直至获得满意的结果 . 单介质管路定常振荡如上所述，对单介质 管路系统，无论是慢速或快速瞬变，均可根 据需要确定时间步长 △ t,然后用a对运动 方程进行修正.在分析波速变动很大的流体 （如气流或含气液流）中，通过在不同计算时 段采用不同的a，从而使计算保持相等的时 间步长. 作为非定常流的特例我们讨论定常振 荡流，在所有的计算节点中，上、下游边界 点上的变量代表系统的输入和输出，是我们 计算的重点，而对其内部节点的变量究竟如 何传递的问题可不予考虑.假定计算管段可 等分为m段（每段为△ x）,振荡周期T可等 分为n段（每段为At）;在m与n之间建立 某种联系，使之与 a有关，从而建立时段 北t内，上下游边界点之间的变量传递关系 由于特征线C+、C-穿过不止一段（m段）管段, 因此也无须计算管道内节点，只对我们所关 心的边界点进行计算.具体做法是 令： △ t=T/n △ x=L/m a =a % /2i 3 L 1 (6) 式中：n,m——与i有关的正整数；i—— 的计算整数.其值影响数值稳定和精度；计 算时 米用a W 100进行控制. 根据特征线的定义又有 a =aA t/ △ x (7) 由此可以证明式(6)中的整数满足 m/n=1/2i+1 (8) 含时间修正因子a的周期性振荡气流特征 型偏微分方程为 式中：S=(2g △ xmsin 0 )/a2.式(11)、(12) 中下标为计算管段上、下游节点编号，上标 为计算时段编号，其中j表示计算时段所求 参数，j-m表示前m个时段的已知参数.由于 摩擦项采用了二阶估算， 故需用迭代法求解 以下举例说明采用上述方程求解的具 体情况.一个正弦型振荡周期至少需 9个点 描述，即n=8;假设i=0即a =aT/2L可满足 精度要求，则计算管段应等分为 m=4段，计 算网格如图1.图中某边界点经 C砸C-两次 特征线的积分，便可获得对面边界点的信息， 所需时间正好为一个周期(8 At)，即计算