基于CIP方法振荡柱体绕流二维数值模拟研究

发布时间：2017-10-19 03:38

  本文关键词：基于CIP方法振荡柱体绕流二维数值模拟研究

  更多相关文章： CIP方法 振荡圆柱 浸入边界法 柱体绕流 锁定域 升力系数 涡脱

【摘要】：柱体绕流是一个复杂的工程问题,它涉及多种科学难题,如紊流、流动分离和涡脱等。准确计算出流体对柱体结构的作用荷载、涡脱频率以及锁定模态,具有重要的科学价值和工程意义。本文以柱体绕流为主要议题,着重研究了流向强迫振荡柱体绕流的水动力特性,其中涉及到主动流动控制和复杂的流固耦合作用,文中也包含了对静止柱体绕流及串列双柱体绕流问题的数值模拟研究。本文所采用的数值模型在笛卡尔直角坐标系里建立,以Navier-Stokes方程和连续性方程为控制方程,通过CIP-ZJU (Constrained Interpolation Profile method in Zhejiang University)差分格式离散对流方程,基于多相流的理论实现流固耦合同步求解,利用浸入边界法处理运动边界。首先对低雷诺数Re(Reynolds Number)下单圆柱和方柱静止绕流进行数值模拟,模拟结果与文献结果进行比较,二者吻合情况较好。通过引入压力阻力项和摩擦阻力项,重点分析Strouhal数、阻力系数、升力系数等参数随雷诺数的变化情况,探讨圆柱和方柱绕流水动力特性的差异。结果表明,在Re≈140时,方柱的分离点由后角点处转移到前角点处,并引起摩擦阻力方向的改变,而圆柱则不存在此现象；方柱绕流中压力阻力项的主导地位要比圆柱绕流中的显著,同时方柱绕流的主波长大于圆柱绕流的主波长。升力的振幅较阻力大一个数量级,且随着Re的增大,方柱振幅的增长速度远大于圆柱。然后,针对振荡柱体绕流建立数学模型,对流向强迫柱体绕流进行数值研究。分析不同振幅、不同振荡频率的振荡柱体绕流的水动力参数规律、涡脱模态以及锁定区间。发现由于振幅和振荡频率的不同,在振荡圆柱绕流和振荡方柱中均可以观察到2T的涡脱模态锁定。当振荡频率处在锁定域内时,升力系数随柱体的振荡运动具有明显的相位锁定,且升力系数和阻力系数均在锁定域内取得极大值。串列双柱体静止绕流的水动力特性相比单柱体静止绕流更加复杂。通过对Re=100,间距比s=L/D=1.5-7.0(L是柱体中心间的距离)的串列双圆柱和双方柱绕流的数值模拟研究,探讨间距对于流场的影响。上下游柱体的涡脱频率相同,且与串列双圆柱绕流相比,串列双方柱绕流的临界距离相对较大。当柱体间距L达到临界距离Lc时,水动力参数将发生跳跃。最后,对Re=100,间距比s=L/D=2.0,串列双圆柱振荡绕流的水动力特性进行数值研究。着重分析上下游柱体升力系数的频谱特性及涡脱模态。研究表明频谱由多个峰值组成：串列双圆柱静止绕流的涡脱频率f0,圆柱的振荡频率fd,谐波频,分数谐波频以及所有上述频率的总和。结果观察到的涡脱S模态与单圆柱振荡绕流相似,同时在锁定域内观察到了2P模态。
【关键词】：CIP方法 振荡圆柱 浸入边界法 柱体绕流 锁定域 升力系数 涡脱
【学位授予单位】：浙江大学
【学位级别】：硕士
【学位授予年份】：2016
【分类号】：TV131.2
【目录】：

• 致谢4-6
• 摘要6-8
• Abstract8-12
• 1 绪论12-23
• 1.1 研究背景及意义12-13
• 1.1.1 柱体绕流问题研究意义12-13
• 1.1.2 柱体振荡绕流问题研究意义13
• 1.2 国内外研究现状13-21
• 1.2.1 单柱体静止绕流研究现状13-15
• 1.2.2 串列双柱体静止绕流研究现状15-17
• 1.2.3 流向振荡柱体绕流研究现状17-19
• 1.2.4 串列双柱体振荡绕流研究现状19-20
• 1.2.5 CIP方法的发展20-21
• 1.3 研究内容21-23
• 2 基于CIP方法的流场数值求解23-32
• 2.1 控制方程23-24
• 2.2 数值方法24-31
• 2.2.1 CIP方法24-25
• 2.2.2 时间积分25-27
• 2.2.3 固体边界处理27
• 2.2.4 计算网格27-29
• 2.2.5 压力Possion方程迭代求解29-31
• 2.3 主要控制参数31
• 2.4 本章小结31-32
• 3 圆柱和方柱绕流的二维数值模拟及结果分析32-53
• 3.1 圆柱绕流的数值模型32-41
• 3.1.1 计算域的选取及网格划分32-33
• 3.1.2 边界条件和初始条件设置33
• 3.1.3 圆柱绕流的模型验证33-34
• 3.1.4 圆柱绕流参数分析34-38
• 3.1.5 圆柱绕流流场分析38-41
• 3.2 方柱绕流的数值模型41-49
• 3.2.1 计算域的选取及网格划分41-42
• 3.2.2 边界条件和初始条件设置42
• 3.2.3 方柱绕流参数分析42-46
• 3.2.4 方柱绕流流场分析46-49
• 3.3 圆柱与方柱绕流参数对比49-51
• 3.4 本章小结51-53
• 4 圆柱和方柱流向振荡绕流数值模拟53-77
• 4.1 静止流体振荡圆柱53-56
• 4.1.1 计算域的选取及网格划分53-54
• 4.1.2 边界条件和初始条件设置54
• 4.1.3 数值算例及结果分析54-56
• 4.2 振荡圆柱绕流56-62
• 4.2.1 初始参数及计算模型56-57
• 4.2.2 数值算例及结果分析57-62
• 4.3 振荡方柱绕流62-73
• 4.3.1 计算域的选取及网格划分62
• 4.3.2 边界条件和初始条件设置62-63
• 4.3.3 数值算例及结果分析63-73
• 4.4 振荡圆柱和振荡方柱绕流主要水动力参数对比73-74
• 4.5 本章小结74-77
• 5 串列双圆柱静止绕流数值模拟77-87
• 5.1 串列双圆柱绕流数值模型77-78
• 5.1.1 计算域的选取及网格划分77-78
• 5.1.2 边界条件和初始条件设置78
• 5.2 数值结果及分析78-86
• 5.3 本章小结86-87
• 6 串列双方柱静止绕流数值模拟87-97
• 6.1 串列双方柱绕流数值模型87-88
• 6.1.1 计算域的选取及网格划分87-88
• 6.1.2 边界条件和初始条件设置88
• 6.2 数值结果及分析88-96
• 6.2.1 Re=100的数值结果及分析88-95
• 6.2.2 不同雷诺数结果对比95-96
• 6.3 本章小结96-97
• 7 串列双圆柱振荡绕流数值模拟97-105
• 7.1 串列双圆柱振荡绕流的数值模型97-98
• 7.1.1 计算域的选取及网格划分97-98
• 7.1.2 边界条件和初始条件设置98
• 7.2 数值结果及分析98-103
• 7.3 本章小结103-105
• 8 结论与展望105-107
• 8.1 本文工作的总结105-106
• 8.2 进一步的工作展望106-107
• 参考文献107-117
• 附件117-118
• 个人简历118-119
• 攻读硕士期间完成的学术论文119

【相似文献】

中国期刊全文数据库 前9条

1 刘月兰;杨小雨;陈禹;刘玉喜;;CIP清洗系统在乳品行业的应用[J];哈尔滨师范大学自然科学学报;2012年03期

2 朱廷儒;彭晓东;王刚;赵小明;唐昌建;;二维高阶CIP数值方法及其应用[J];核聚变与等离子体物理;2012年02期

3 孙雪梅;王舰;;从CIP引进马铃薯资源的倍性鉴定[J];青海大学学报(自然科学版);2007年01期

4 朱琳;;立方插值拟质点(CIP)方法的改进及其应用[J];佳木斯大学学报(自然科学版);2011年01期

5 张捍卫,许厚泽,王爱生;弹性地球CIP轴的极移和岁差章动[J];北京大学学报(自然科学版);2005年05期

6 游娟;黄建林;曹莉;韩日畴;;重组昆虫病原细菌Cip蛋白基因工程菌发酵条件的初步研究[J];生物学杂志;2012年01期

7 张捍卫,许厚泽,王爱生;液核地球CIP轴的极移和岁差章动[J];山东理工大学学报(自然科学版);2005年01期

8 成善汉;谢从华;柳俊;吴艳阁;;马铃薯低温启动子CIP的克隆及其功能鉴定[J];贵州科学;2008年01期

9 ;[J];;年期

中国重要会议论文全文数据库 前2条

1 杨万林;隋启君;;CIP马铃薯资源在云南的评价、利用及前景[A];2005年全国马铃薯产业学术年会论文集[C];2005年

2 吴承金;程群;李大春;戴清堂;;CIP马铃薯资源的引进与初步评价[A];马铃薯产业与粮食安全（2009）[C];2009年

中国重要报纸全文数据库 前1条

1 闻圆圆;土耳其航空推出全新CIP休息室[N];中国民航报;2011年

中国硕士学位论文全文数据库 前2条

1 李岩;出口企业（卖方）在CIP合同下的义务和风险防范研究[D];对外经济贸易大学;2015年

2 付英男;基于CIP方法振荡柱体绕流二维数值模拟研究[D];浙江大学;2016年

，

本文编号：1058838