收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

含运动物体三维复杂流场数值模拟和尾流结构研究

邓见  
【摘要】: 论文首先对所使用的数值方法进行了论述:简要介绍了高精度分块耦合算法,此算法是针对不可压缩粘性流体的,采用任意曲线坐标系下原始变量形式的N-S方程进行描述。对于运动边界问题,本文引入了任意拉各朗日-欧拉法(ALE)进行处理。差分离散过程中,粘性项采用四阶精度的Pade紧致格式,对流项采用三阶迎风紧致格式离散,时间上采用四阶显式Runge-Kutta格式离散。作者所做的贡献是结合分块耦合方法的天然粗粒度并行特性,利用MPI函数库实现了并行,并且将基于分块耦合算法子域划分的任务分配策略与基于负载平衡的任务分配策略进行了比较,发现虽然分块耦合算法能够提供天然的任务分配方案,但却很难保证各进程上的负载平衡。 另一个计算方法是沉浸边界法(Immersed Boundary Method),该方法是基于笛卡儿网格系统的,因此在网格生成上可节省大量时间,而且可以大大提高计算效率,尤其适合含运动物体的流场求解。本文先归纳了前人在沉浸边界法上所做的努力,随后发展了自己的求解器,先后在交错网格和非交错网格上实现了沉浸边界法,并结合大涡模拟(LES)对高雷诺数下的湍流数值模拟进行了尝试。 用高精度分块耦合方法结合ALE对弹性支撑方柱绕流进行了数值模拟。模型简化为“弹簧-阻尼-质量”系统,柱体运动采用经典龙格-库塔方法求解。同时模拟了方柱的驰振和涡致振动,捕捉到了驰振到涡致振动的转变,捕获了涡致振动的“锁定”、“拍”、“相位开关”等现象。 对单圆柱绕流尾流三维特性进行了研究,找到了尾流向三维转捩的临界雷诺数Re=195,低于这个临界雷诺数时,尾流保持二维状态,此时的流场数值模拟可以近似为二维,但当超过这个临界雷诺数时,二维求解不再适用。当雷诺数超过临界值后,随着雷诺数的进一步提高,尾流先后经过两种过渡的三维模态,这两种模态分别被称作模态A(Mode A)和模态B(Mode B)。这两种转捩区的模态有各自不同的特征和物理机制,文中分别进行了研究。在此基础上,还研究了圆柱的旋转振荡对尾涡的影响。此研究是在Re=260下进行的,这个雷诺数下的圆柱尾流处在Mode B的状态。研究发现,当旋转振荡的振幅足够大时,能将尾流带回二维状态。 对串列双圆柱绕流三维特性进行了研究。从两个不同方面进行了分析:首先是固定雷诺数Re=220而改变上下游圆柱之间的间距(1.5≤L/D≤8),研究表明,当L/D≤3.5时,尾流呈现二维状态,而当L/D≥4时,尾流向三维转捩。据此可以判断临界间距存在于区间3.5≤L/D≤4;其次,将间距固定在L/D=3.5,来研究此间距下出现尾流三维特性的临界雷诺数,通过计算得到临界雷诺数为Re=250。在这部分的研究中,首次发现:在一个圆柱尾流场放置另一圆柱可以一定程度上抑制尾流向三维的转捩,从而推迟向湍流的转捩。 对置于方形管道里的垂直交叉双圆柱结构进行了研究,本算例采用高精度分块耦合方法的并行程序,首先给出了详细的分块处理程序,以及内部拟边界、角线和角点的处理。圆柱端部紧贴管道壁面,两圆柱间距为圆柱直径的5倍,计算雷诺数为200。研究发现,交叉结构的中心流场较复杂,上游圆柱尾流的规则卡门涡街被下游圆柱截断成复杂的缠绕结构。两圆柱间的相互影响,以及管壁的影响,形成了流场特殊的涡结构。 对无限长垂直交叉双圆柱绕流进行了研究,包括两圆柱都静止的情况以及下游圆柱强迫振动或涡致振动的情况。对于两圆柱都静止的情况,当L/D>3时,下游圆柱对上游圆柱尾流的影响只限定在下游圆柱的尾流所扫过的范围之内;当L/D<3时,下游圆柱对上游圆柱尾流的影响扩大,下游圆柱尾流扫过区上下出现两排三维流向二次涡结构。当下游圆柱做强迫振动或涡致振动时,尾流的涡结构和宏观量都发生了改变,尤其当下游圆柱弹性固定作涡致振动时,与独立的圆柱绕流相比,振幅响应的幅度降低,共振区域的带宽增加。 最后一章提出了用计算流体力学的方法来研究鱼群游动的流体动力学性质,对一些实验根本无法实现的内容进行了探索性的研究,例如:巡游状态下,特定队形对游动推进效率的影响。从前人观察得到的菱形编队中抽出三条鱼组成的基本单元,当一条鱼位于另外两条鱼中间的下游区域时,如果处于加速巡游状态,则上游鱼尾流产生的反卡门涡街使下游鱼体两侧流体产生了向前的运动,从而提高了下游鱼游推进效率,减少了功率消耗。本文首次对前人关于鱼群减阻和提效的假设进行了验证。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 范本隽,沈培玉;极坐标系边界权残法求解二维温度场问题[J];无锡轻工大学学报;2000年06期
2 泰庆华;Hybrid-Trefftz有限元简介和展望[J];力学与实践;1993年03期
3 张洪武,何扬,李锡夔;改进的广义Smith非反射边界法[J];大连理工大学学报;1991年02期
4 牛庠均;固体内任意元素的边界积分变分定理——任意裂纹开展时能量释放率的计算[J];应用数学和力学;1984年01期
5 骆振欧;Poisson方程的高精度有限差分格式[J];计算物理;1986年01期
6 徐次达;加权残数法在固体力学中的应用——我国近年来进展情况[J];应用数学和力学;1982年06期
7 徐次达;加权残数法用于固体力学近年来进展概况[J];力学进展;1982年02期
8 谢全敏,朱瑞赓,池秀文;危岩块体稳定性分析的蒙特卡罗边界法[J];灾害学;1998年02期
9 童朝锋,严以新,林跃华;长江口深水航道数学模型信息系统研究[J];海洋工程;2003年01期
10 张晓蓉,徐剑刚;日元汇率理性投机泡沫的实证检验分析[J];复旦学报(自然科学版);2005年02期
11 许永林 ,唐锦春;边界积分直接法解板的弯曲问题[J];浙江大学学报(工学版);1985年02期
12 P.A.Berge;用熔融玻璃球制作的砂岩体模拟超声速度-孔隙度关系[J];勘探地球物理进展;1995年05期
13 童卫华,姜节胜,廖巍;基于可靠性分析及频率约束的结构优化[J];机械科学与技术;1997年04期
14 黄河宁,李鉴初;散射波的二阶幅射边界条件及其应用[J];海洋学报(中文版);1988年02期
15 牛庠均;非线性固体的断裂分析[J];北京工业大学学报;1991年02期
16 蔡荫林,周健生;结构系统失效概率的近似分析[J];航空学报;1992年03期
17 彭凯,方锋,曹叔尤;在二维流动计算中应用“河床切削”技术处理动边界问题[J];水动力学研究与进展A辑;1992年02期
18 张永苍;计算二维正态分布的导函数拟合法[J];西北工业大学学报;1994年03期
19 李论;改进加权残值法收敛性的一种方法──逐次逼近法[J];西安石油学院学报(自然科学版);1995年04期
20 张剑锋,刘书;弹性波混合差分解法的透射边界条件[J];固体力学学报;1997年04期
中国重要会议论文全文数据库 前1条
1 涂小萍;;宁波市年降水量分析及消除小波变换边界影响的比较试验[A];推进气象科技创新加快气象事业发展——中国气象学会2004年年会论文集(下册)[C];2004年
中国博士学位论文全文数据库 前2条
1 周鲲鹏;基于在线动态等值和能量函数法的电力系统暂态稳定控制研究[D];武汉大学;2004年
2 田力;地下爆炸波冲击下地面结构动力响应及滑移隔震研究[D];天津大学;2004年
中国硕士学位论文全文数据库 前2条
1 白凤仙;面向工程的电磁场软件与电磁场逆问题研究[D];太原理工大学;2002年
2 周晓燕;永磁电机二维电磁场软件平台[D];沈阳工业大学;2003年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978