变加速运动圆球附加质量与阻力仿真分析.docVIP

变加速运动圆球附加质量与阻力仿真分析 　　作者简介： 姚保太（1973―），男，河南获嘉人，高级工程师，研究方向为机械制造及仿真计算，(E??mail)yaobaotai@163.com 　　通信作者： 康宁（1963―），女，辽宁金州人，副教授，博士，研究领域方向为流体仿真计算，(E??mail)kangning@0引言 　　附加质量和阻力计算是水下物体运动分析的关键问题.［1］以往对附加质量的计算主要是根据势流理论［1??4］，如在无界无黏流体中球体等简单形状物体的附加质量存在解析解,复杂形状物体的附加质量通过估算法或三维面元法求得.近年来，随着计算流体力学的发展，开始用数值模拟的方法计算物体的附加质量，发展基于动网格技术计算附加质量的新方法.［5］数值方法有更广泛的应用空间，不仅可以计算理想流体中的附加质量［5??6］，还可以计算黏性流体中的附加质量［7??10］. 　　人们已经熟知无界无黏流体中物体加速运动时的附加质量只与物体形状及流体密度有关，与运动加速度无关，目前还缺少系统地研究在有界区域中流体黏性和加速度对做变加速运动物体附加质量的影响.由于工程中实际问题的复杂性，比如水下导弹发射问题等，因此有必要研究在有界区域黏性流体中物体做变加速运动的附加质量问题. 　　本文采用数值模拟方法和动网格技术，计算管体域中黏性不可压流体中圆球做变加速运动时的附加质量、阻力和阻力因数，深入分析和研究在管体域中圆球变加速运动对附加质量、阻力和阻力因数的影响. 　　1计算模型和计算方法 　　1.1控制方程组 　　水中运动物体满足非定常不可压N??S方程， 　　连续方程： ??vi??xi=0(1)动量方程：??vi??t+vj??vi??xj=Fi－1ρ??p??xi+μρ?氮?xj??vi??xj(2)式中：vi和vj为速度；F为质量力；p为压力；ρ为流体密度；μ为流体动力黏性系数. 　　1.2计算域、边界条件和网格划分 　　1.2.1管体域 　　管体域由管体、喇叭口和水域组成，圆球直径为0.254 m，圆球在管体域中运动前的计算域见图1.通海水口边界条件为压强进口，出口为压强出口，压强均为57 m水深处压强574 826.8 Pa；根据圆球平移运动特点，动网格模型选择动态网格层变方法，圆球表面为无滑移运动壁面条件；管体和喇叭口为无滑移静止壁面条件，水域边界为滑移静止壁面条件.圆球在管体域中运动前、后网格分布见图2和3.采用六面体结构网格，网格数约为58万个. 　　 　　图 1管体域中圆球运动前的计算域 　　Fig.1Computation domain in tube domain before ball 　　movement 　　 　　图 2管体域中圆球运动前的网格分布 　　Fig.2Mesh distribution of tube domain before ball movement 　　 　　图 3管体域中圆球运动后的网格分布 　　Fig.3Mesh distribution of tube domain after ball movement 　　 　　1.2.2小计算域 　　作为对比，对圆球在小计算域中的运动情况进行计算.小计算域为圆柱形，取较大尺度，直径为3.6 m，长度为14 m.纵向对称面网格分布见图4.六面体结构网格，网格数约为139万个.计算域左侧为压强入口边界条件，右侧为压强出口条件，相对压强为0；其他边界和动网格模型同管体域. 　　 　　图 4小计算域纵向对称面的网格分布 　　Fig.4Mesh distribution in longitudinal symmetry plane of 　　small computation domain 　　 　　1.2.3大计算域 　　为验证附加质量仿真方法的可靠性，对圆球在大计算域中的运动情况进行计算.大计算域仍为圆柱形，直径加大到14.4 m，长度加大到20.86 m.采用六面体结构网格，网格数约为264万个.边界条件和动网格模型同管体域.1.3圆球变速运动参数 　　1.3.1先加速后减速 　　为模拟发射体在管体内直至出管的运动情况，对圆球先加速后减速运动情况进行计算.加速阶段模拟发射时在主动力作用下的运动情况，减速阶段模拟在惯性力作用下直至出管的运动情况，具体参数见表1. 　　表 1圆球先加速后减速的运动参数 　　Tab.1Movement parameters when ball is accelerated motion first and then decelerated初速度/ 　　(m/s)加速度/ 　　(m/s2)时间/s末速度/ 　　(m/s)位移/ 　　m加速阶段