高雷诺数下静止管道绕流的数值分析.pdfVIP

大 庆 石 油 学 院 学 报 第31卷 第6期 2007年12月 OF PETROLEUMINSTITUTE V01．31No．6 Dec． 2007 JOURNALDAQING 高雷诺数下静止管道绕流的数值分析 龙飞飞，戴 光，付 冰 (大庆石油学院机械科学与工程学院，黑龙江大庆 163318) 摘要：根据流体力学原理，分析静止管道的旋涡脱落}根据管道旋涡脱落原理，建立高雷诺数下静止管道的物理模 型及数学模型I通过数值分析，得出在绕流问题中雷诺数是决定流型的主要因素，在不同的雷诺数下，旋涡脱落也随之不 同；升力系数和阻力系数曲线在涡稳定脱落后呈周期性变化，幅值不变，且阻力的脉动频率是升力的2倍．得到涡脱落的 位置与管壁切应力的关系以及绕流时管道的压力分布情况． 关键词：管道；旋涡脱落}雷诺数；升力系数}阻力系数；压力分布 中圈分类号：0353．4 文献标识码：A 文章编号：1000—1891(2007)06—0090一05 0 引言 长输管道穿越江河时，如果河床受到冲刷，易使管道裸露并悬空，形成梁式结构．裸露管道在水流作 用下，水动力作用强，水流在管道后面形成尾流区，在一定流速时，旋涡发放频率可能接近管道的自振频 率，从而出现垂直于水流方向，类似共振的“频率锁定”现象．具体表现为管道以接近旋涡自然脱落频率作 横向振动，或以接近2倍旋涡自然脱落频率作纵向振动时，管道与尾流之间发生强烈的非线性相互作用， 旋涡将以管道的自振频率脱落．这种“频率锁定”将导致沿水流方向的振动和拖曳力增大，产生很大的交 变应力．其中，管道悬空段横向涡激振动的振幅远大于顺流向涡激振动的振幅，使管道产生疲劳破坏[1]． 软件，模拟均匀流中的静止管道绕流的旋涡脱落． 1 圆柱绕流的基本特征 流体绕过圆柱体时，物体表面边界层流动在逆压梯度下分离，或在角点处分离，形成向下游拖曳的剪 切层．圆柱体两侧的剪切层相互作用，在一定的雷诺数范围产生规则的旋涡结构，诱发顺流体流动方向和 垂直于流体流动方向的作用力．旋涡的形成和脱落呈现周期性，表征这一特性的相似准则数为斯特鲁哈 (Strouhal)数Sr．斯特鲁哈数随雷诺数变化，对于不同的柱体，其变化规律不同嘲． 2 基本方程 不可压缩黏性流在层流状态下，采用二维不可压黏性流体的原始变量的N二-S方程、连续方程和流函 数涡量方程‘6|．二维不可压黏性流体的原始变量的N—S方程和连续方程矢量形式为 2H， (1) p+口v 挈+(H·v)=F一』v o‘ D V·H一0。 式中：H为速度矢量；口为运动黏性系数；F为单位质量流体的质量力；P为压强；JD为水的密度． 收稿日期：2007一Ol一18}审稿人：邹龙庆；编辑：任志平 作者筒介：龙飞飞(1978一)，男，博士生，讲师，主要从事声发射检测与结构完整性评价方面的研究 ·90· 第6期 龙飞飞等：高雷诺数下静止管道绕流的数值分析 3 计算模型 模拟水流冲击作用下，Re为4×105的管道周围的流 场．此时，管道表面上的边界层转变为湍流，分离点后移， 阻力显著下降，旋涡的发放不规则．对复杂绕流，采用非稳 态RNGk一￡模型求解方程，获得旋涡模拟结果．迭代中压 力速度耦合采用SIMPLE算法，湍动能和湍流耗散率采用 二阶精度，以获得较为准确的结果[7-9]．均匀水流中静止 管道绕流的网格划分见图1． 图1管道绕流的网格划分 3．1非稳态计算时间步长 长： T 1