新型油气施工流程泵构件材料内场强数控运算.docVIP

新型油气施工流程泵构件材料内场强数控运算 　　1模型的建立及网格划分 　　该模型由转轮、蜗壳、过渡流道3部分组成。转轮为双吸叶轮，叶片数为5，蜗壳为180°对称双蜗壳。 　　根据叶轮和双蜗壳流道的结构参数、水力图以及叶片木模图，利用Pro /E三维绘图软件进行整体造型，整体流道模型。由于几何结构比较复杂，本模型采用ANSYS公司的前处理软件ICEM划分网格。ICEM软件最大的优点是根据模型各面的特征，将其分成不同的part，根据需要定义不同的尺寸，可大大提高网格质量，优化网格结构，为下一步的流场计算打下很好的基础。笔者采用了对复杂边界适应性强的非结构化四面体网格对计算域进行网格划分。计算时，逐渐加密网格，通过考察相对速度的最大、最小值及等值线分布情况来判别是否获得与网格无关的数值解。最后针对本模型，同时考虑计算资源，选用网格数为249 205，最终基本得到与网格无关的数值解。 　　2控制方程 　　将不可压流体的N - S方程进行时间平均可得到雷诺平均N - S方程，其在直角坐标系下的张量为 　　5ρui5xi=0 　　式中u i（ i = 1，2，3）为时均速度；ρ为密度。 　　由Boussinesq假设，雷诺应力模型μk 2ε（4）式中μt为涡粘系数； k为湍动能； Cμ为常数；ε为湍动能耗散率。 　　C 1ε， C 2ε为常数；σk，σε分别是k和ε的普朗特数； G k为平均速度梯度引起的湍动能生成项。以上方程中的模型常数。 　　以上控制方程可直接应用于静止部件（蜗壳与过渡段）内流场的计算，转动部件（转轮）内流场计算在转动参照系下进行，其控制方程与静止部件内流场的控制方程的形式相同，但速度采用转动参照系下的相对速度。 　　3数值方法 　　3. 1方程的离散 　　动量、湍动能和湍动能耗散率方程的对流项均采用一阶迎风格式。压力速度耦合采用了SIMPLE算法。 　　3. 2边界条件给定 　　（1）进口边界条件：进口位于叶轮进口断面，给定其速度大小和方向，湍流参数通过给定湍流度和进口特征尺度给出。双吸叶轮有两个进口，设置完全一样。 　　（2）出口边界条件：出口位于过渡段的末端，采用充分发展湍流假设，即变量沿出口截面的法向导数为零。双蜗壳结构两个出口，定义完全一样，均为质量出口。 　　（3）固壁边界条件：固壁边界采用无滑移边界条件，即流体速度等于边界速度。动静接触边界处理和多块网格间的内部界面一样，只是在界面两侧分别采用静止和旋转参照系下的控制方程。紊动能、紊动能耗散率由壁面函数决定。在计算中要用到动静参照系下的相对速度和绝对速度之间的变换公式。这样的动静接触边界处理可以允许蜗壳进口的周向不均匀分布。 　　4计算结果分析 　　2为双蜗壳流道的静压分布图，由此图可以看出，双蜗壳内的压力成180°对称分布，从隔舌到扩散段出口，压力整体呈逐渐增大的趋势。从隔舌到扩散段进口压力变化比较均匀，从扩散管进口到过渡段出口，压力增加比较明显，这跟压水室中主要是靠扩散管将80%～85%的动能转化为压力能的原理相符合。同时压力在扩散管进口处有个波动，进入扩散管后，压力出现一面高一面低的状况，即蜗壳外侧压力高于内侧压力。这是由于在扩散管进口处尺寸开始同时沿径向和轴向变化，空间扭向比较明显，液流在此处的流向改变幅度相应比较大，因此会出现压力局部降低，速度局部增大的情况，这在3所示蜗壳绝对速度矢量图中也可以反映出来。这些现象会造成一定的能量损失，因此在扩散段的设计中，可进一步考虑扩散段的长度、扩散角度，以优化扩散段的设计，提高泵的性能。 　　从压力图上还可以看出，在1断面与2断面之间明显有涡存在。涡的存在必然导致蜗壳内的能量损失。从提高水力性能的角度考虑，在泵体设计中应优化考虑蜗壳断面面积的变化规律，以便可以更好地优化设计。 　　同时由的压力图可以看出压力呈180°基本对称分布，根据径向力产生的原因，两个对称的蜗壳的径向力基本抵消。因此180°对称分布的双蜗壳结构能较好地平衡径向力。 　　3为双蜗壳绝对速度矢量分布图，由此图可以看出：绝对速度矢量呈180°对称分布，从隔舌到扩散段出口，速度矢量逐渐减小，进入扩散管后速度矢量减小更为明显，这进一步验证了压水室中80%～85%的速度能在扩散管中转变为压力能这一理论。同时在扩散管进口处速度矢量出现局部增大，并且出现速度局部增大的部位与蜗壳静压分布图中压力局部减小的部位相同，这是由双蜗壳的特殊结构造成的。因为对于单蜗壳，其扩散管多为直型锥形，尺寸只是沿着径向变化。而HD型石油化工流程双蜗壳结构的扩散段在空间是扭向的，同时沿着径向和轴向变化，因而扩散段进口处转弯比较明显，液体在此处的流向变化相应比较明显。 　　4为双蜗壳湍动能分布图，从图中可以看出，从蜗壳进口到出口，湍动能变化比较均匀，只是在蜗壳进口处和扩散管进口处出