CFD计算速度的分析.doc

CFD计算速度的分析 CFD计算旋转流时旋转中心与几何中心不重合的速度变换 贾继斌 陈功 (中国石油大学(北京)化工学院，北京102249) 摘 要：基于CFD计算旋转流的三维速度计算值是以空间的几何中心为基准给定的。但旋转流由于流动的不稳定性或空间结构的非轴对称性导致旋转流的旋转中心偏离空间的几何中心。这就使参照旋转中心为基准的三维速度分布不同于基于CFD计算的三维速度分布，两者之间存在着一定的差别。在实际应用中一般分析问题是建立在以旋转中心为基准的三维速度分布上的，因此需要将参考几何中心CFD计算的速度变换为参考旋转中心的速度。本文以旋风分离器为实例分析了由于基准不同所产生的速度变化，给出了CFD计算旋流时旋转中心与几何中心不重合的速度变换的方法。结论也同样适用于其他类似旋转流的计算分析。 关键词：CFD，旋流，旋转偏心，速度变换，旋风分离器 中图分类号：TQ 051.8 文献标识码 A 文章编号 引 言 旋转流是一种工业上广泛应用的流动形式，如旋风分离器，水力旋流器，旋转射流等。一般基于CFD计算旋转流的三维速度计算值是以空间的几何中心为基准坐标给定的。但由于旋转流动过程的不稳定性或几何结构上非轴对称等原因，旋转流的旋转中心偏离空间的几何中心[1]，存在着旋转偏心的现象，即旋转流的旋转中心与几何中心不在一个轴线上，两个中心不重合。若以不同的中心为参考基准，旋转流在空间上同一点的三维速度是不同的，存在着一定的差别，互相直接对等替换必定为后续的计算或分析带来较大的误差。在实际应用中，由于分析问题的特点，在对旋转流流场进行分析中有时需要以旋转中心为基准的三维速度分布，此时应对基于CFD计算的旋转流三维速度值进行修正，将参考几何中心计算的流场变换为参考旋转中心的流场。 本文以蜗壳式旋风分离器流场为CFD计算对象，重点分析旋风分离器内旋转流的旋转中心与旋风分离器的几何中心不在一个轴线上所导致的速度分量之间的差异，给出了将基于几何中心计算的流场换算到参考旋转中心的流场的方法。结论方法也同样适用于其他有关类似旋转流的计算分析。 1 基于几何中心的CFD计算旋转流时的速度分布 1.1数学模型及方法 1.1.1几何模型 模拟旋风分离器尺寸如图1所示，进气口尺寸176mm×73mm，筒体内径300mm，升气管内径96mm。为保证充分发展条件的成立，将出口管路加长。 图1 旋风分离器结构示意图 1.1.2湍流模型 旋风分离器内流场是复杂的三维、湍流、强旋流场，具有很强的各向异性特点，为此采用Fluent软件计算流场时，选择能较好地反映湍流各向异性的RSM（Reynolds stress model）模型。压力梯度项采用PRESTO!（pressure staggering option）方法进行预处理 1.1.3 边界条件 1.1.3.1 入口边界条件 旋风分离器入口边界条件采用速度入口条件。流体为常态的空气, 入口速度Vi＝=20m/s，湍流脉动动能＝0.735 m2/s2 ，＝＝0.49 m2/s2， m2/s2，= 17.021 m2/s3， 其中，＝0.09，＝1.42， 为入口处的平均流速，I是湍流强度，DH为入口水力直径，为湍流特征尺寸。取I=0.05，，DH =0.116m。 1.1.3.2 出口边界条件 出口边界条件按充分发展管流条件处理，所有变量在出口截面处轴向梯度为零，即。这种处理方法对出口截面的选择要求比较高,一般只有当出口区域有一平直段且离开回流区已经较远时采用。为此在计算中将旋风分离器的出口管路加长，以保证充分发展条件的成立。 1.1.3.3 壁面边界条件 壁面处采用无滑移边界条件，默认壁面粗糙度为0.5。壁面效应是旋涡和湍流的主要来源,因此近壁区的处理对数值求解结果的准确性有显著影响。由于在靠近固体壁面的区域内，流底层的粘性作用增强而湍流扩散相对减弱，作用于高雷诺数下的湍流输运方程已不能严格有效。近壁网格点用标准壁面函数近似处理。 1.1.4差分格式 在Fluent中，认的差分格式是一阶迎风格式。论分析和计算过程中发现，阶迎风格式只具有一阶精度的截差，虽然具有良好的收敛稳定性，但数值耗散过大，尤其对于强旋转流，致使结果产生很大的误差；二阶格式具有较高的数值计算精度，但二阶格式本身具有色散性,结果会产生非物理震荡，不具有守恒性；QUICK差分格式比前两种格式都有所改善，保留了结果的守恒，又使数值结果具有二阶以上的精度的截差。QUICK差分格式最能准确的预报旋风分离器内的流场,二阶迎风格式只能在某些方面给出好的结果，另一些方面却远偏离实际情况，一阶迎风格式由于精度差，根本不能用于旋风分离器这种复杂流场的模拟。因此，本文将采用QUICK格式控制方程的离散。 1.1.5算法