双螺杆压缩机流场动力学特性的动态模拟仿真.docxVIP

双螺杆压缩机流场动力学特性的动态模拟仿真 0 数值模拟研究 双螺杆泵具有可靠性、使用寿命长、结构简单、技术素质等优点。广泛应用于工业生产。国内外科学家对双螺杆泵的研究日益深入。Wu等建立了双螺杆压缩机的新颖数学模型, 通过实验研究验证了其合理性。Kovacevic等将数值模拟的方法运用到螺杆机械的研究之中。Cao等对双螺杆压缩机的内部压力分布进行了实验研究, 阐述了进气口与排气口处压力的分布规律。张小军等结合试验研究了双螺杆压缩机排气压力的脉动, 发现产生脉动的原因是因为实际排气压力与设计值偏离。本文对双螺杆压缩机建立了合理有效的流体模型, 运用动网格技术和Fluent软件对双螺杆压缩机的性能进行动态模拟, 并结合实验测试验证了动态模拟仿真的结果的正确性, 为双螺杆压缩机的改进设计提供一种有效的手段。 1 双螺杆压缩机的工作过程 双螺杆压缩机的结构如图1所示, 在双螺杆压缩机的机体中, 平行放置着一对相互啮合的螺旋形转子。通常把节圆外具有凸齿的转子称为阳转子或阳螺杆, 把节圆内具有凹齿的转子称为阴转子或阴螺杆。一般来说, 阳转子与原动机相连接, 由阳转子与阴转子之间的同步齿轮实现阳转子带动阴转子转动。在双螺杆压缩机机体的两端, 分别开设一定形状和大小的孔口。一个供吸气用, 称为进气口;另一个供排气用, 称为排气口。 双螺杆压缩机的工作循环可分为进气、压缩和排气三个过程。随着转子旋转, 每对相互啮合的齿相继完成相同的工作循环。具体的工作过程如下:转子转到某一角度后, 两个孤立的齿间基元容积由于阳转子的凸齿侵入阴转子的凹齿, 基元容积同时开始缩小, 实现气体的压缩过程, 直到一对基元容积与排气孔口相连通的瞬间停止;基元容积和排气孔口相连通后, 排气过程开始, 排气过程一直持续到两个齿完全啮合, 即两个基元容积因两个转子完全啮合而等于零时。 2 气体等熵流动 本文中双螺杆压缩机动态数学模型的建立基于以下假设: (1) 压缩气体为理想气体; (2) 工作介质没有向外界环境发生泄漏; (3) 润滑油对工作介质的性质没有改变; (4) 通过进气口、排气口的气体为等熵流动; (5) 双螺杆压缩机的工作循环为一个理论循环过程。 2.1 压缩机边界条件 双螺杆压缩机的工作过程比较复杂, 为了表现每一次工作循环中的进气、压缩和排气过程, 必须实现基元容积与进气口和排气口之间的连通和关闭, 需要对螺杆压缩机的流体模型分割成三个部分:进气口的流体、基元容积内的流体、排气口的流体。通过设置边界条件使其控制基元容积与进气口与排气口之间的流通与封闭。在建立基元容积内的流体时, 根据动网格计算要求运动部件之间要留有一定的间隙, 设置阴阳转子与腔壁之间的距离为0.5mm, 两转子之间的中心距增加0.5mm。本文的计算域模型如图2所示。 2.2 双螺杆压缩机的流场模型 双螺杆压缩机内部流体是一种高度复杂的三维非稳态、旋转的不规则流动, 是典型的湍流。修正的RNG/ (k-ε) 模型在强流线弯曲、漩涡和旋转流模型中有很好的效果, 利用经压缩性修正的RNG/ (k-ε) 模型对腔体内流动的模拟可以避免标准的k-ε模型在用于强旋流或带有弯曲壁面流动时出现的失真, 更接近于实际状况。因此, 本文采用RNG/ (k-ε) 模型来模拟双螺杆压缩机的流体流动。 RNG/ (k-ε) 模型中, 通过在大尺度运动和修正后的黏度项体现小尺度的影响, 而使这些小尺度运动系统地从控制方程中去除。由此, 湍流脉动动能方程 (k方程) 和湍流耗散率方程 (ε方程) 分别为 式中, p为该场压力;μ为动力黏性系数;μi、μj分别为气流速度在i、j方向上的分量;μt为湍流黏性系数;ρ为密度;Cμ为常数;δε为湍流脉动动能耗散率ε的湍流Prandtl数;δk湍流脉动动能k的湍流Prandtl数;sij、sil均为流体变形率;gi为i方向的重力加速度。 式 (2) 中Cμ、δk、δε、Cε1、Cε2、Cε3、Cε4、η0、β等均为经验系数, 其取值如表1所示。 2.3 流体边界条件设置 采用Fluent软件对双螺杆压缩机的流场分布规律进行动态仿真, 需要首先确定模型的边界条件及相关的设置。在双螺杆压缩机的流体模型中, 对基元容积处的流体与进气口和排气口处的流体的边界设置是影响模拟进气、压缩、排气过程与实际工况是否符合的关键。 流体边界条件类型设置见表2。在Fluent软件中通过Grid Interface设置交界面1-1与1-2、2-1与-2-2、3-1与3-2-1和3-2-2使其配对, 保证重合边界的属性为Interior, 实现流体从进口到出口的连通。进气口、排气口的初始边界条件设置见表3。 2.4 流体网格划分 双螺杆压缩机的流体模型的原始网格利用Gambit软件生成, 三维流体模型的网格