基于重叠网格的航行体水动力特性计算_韩守根.pdf

第 卷 第 期兵 器 装 备 工 程 学 报 年 月 收稿日期：；修回日期：基金项目：国家自然科学基金面上项目（）作者简介：韩守根（），男，硕士研究生，主要从事计算流体力学研究，：。通信作者：杨岩（），男，工学硕士，助教，主要从事兵器发射理论与技术、力学、结构强度设计等研究。：基于重叠网格的航行体水动力特性计算韩守根，杨 岩，刘丙鑫，陶 钢（南京理工大学 能源与动力工程学院，南京；南京航空航天大学 直升机传动技术重点实验室，南京）摘要：利用重叠网格技术结合实验方法，研究了航行体跨界质过程中的阻力来源。将试验结果与仿真结果进行对比，得出该过程中的水动力及水动力系数变化趋势以及原因，结果表明：航行体在水下速度由 下降到 ，水动力以及系数在航行体即将出水时刻达到最大值，分别为 和 。出水过程水动力及力系数主要受水冢现象、浮力变化以及附连水质量变化影响而呈现波动并逐渐趋向于零，粘滞力是航行体完全出水后主要的流体阻力来源。该研究为水航行体减阻设计提供了参考依据。关键词：流体动力；多相流；重叠网格；阻力计算；数值模拟本文引用格式：韩守根，杨岩，刘丙鑫，等 基于重叠网格的航行体水动力特性计算 兵器装备工程学报，（）：，：，（）：，中图分类号：文献标识码：文章编号：（），（，；，）：，：；引言潜射航行体的发射是一个涉及多因素干扰，目前为应对各种复杂环境需求，航行体朝着变深度发射、简化结构、提高发射可靠性等趋势发展，其中发射可靠性一直是其重点研究方向之一。以往的研究往往只关注水下过程，忽略了出水过程的阻力动态演变。因此出水过程的阻力来源对于水下航行体的减阻设计是非常必要的。利用实验方法与重叠网格技术结合，揭示了航行体出水过程中水冢现象、浮力变化以及附连水质量变化对于水动力及水动力系数影响的内在规律，为工程设计提供了指导。数学模型 模型 模型中动量方程为求解整个区域内统一的动量方程，作为结果的速度场由各相共享，即：（）（）（）（）相间界面跟踪通过求解一项或多项容积比率的连续方程来完成。对第 相方程为：（）（）（?）（）式（）中：?为第 相对第 相的质量传输；?为 相对第 相的质量传输；为源项；为容积比率。水动力系数水动力系数是水下航行体结构设计的重要参数，一般而言，水动力系数大多通过水动力测量试验和理论方法计算相结合获得，由于航行体水动力系数很多，想要全部辨识是不切实际的，本文选取最常用的无量纲化水动力系数来建立数学模型，有：（）（）式（）中：为水的密度；为水深 处航行体最大纵向速度，即为质心速度在纵向方向的分量；为航行体特征横截面积；为航行体轴向长度。航行体质心的流体力合力为：（）（）式（）中：（）为表面压力积分项；为黏性力积分项。实验方案实验前航行体与发射平台在室内进行吊装，吊装过程如图 所示，吊装完成后，将航行体和发射平台移至湖内水下，竖直固定在湖内水下。为获得航行体表面压力，实验时在航行体表面布置了多个压力监测点，测点 为航行体头部驻点，测点 分布于头部母线段，测点 分布在柱段，测点 位于尾段。航行体几何尺寸及测点具体位置如图 所示。图 航行体发射平台室内吊装工作示意图 图 航行体几何尺寸及测点位置 试验当天无风，湖面基本无浪，航行体从水下 处垂直起射直至射出水面。其水下弹射出筒过程如图 所示，射出水面过程如图 所示。图 航行体水下垂直弹射过程 图 航行体出水过程 韩守根，等：基于重叠网格的航行体水动力特性计算 实验中航行体从初始体到完全出水等过程时序如图 所示。航行体尾部离开发射筒时，高压发射气体在冲出筒口的过程中将会对流场产生了剧烈扰动，形成筒口效应，这将会使得此时的数据误差较大、可信度低，数据分析时要避免这一段的筒口效应影响，因此本文数值模拟重点为模拟筒口效应消失到航行体完全出水这一段过程，并对这一过程的流场以及水动力系数进行分析讨论。图 出水过程时序图 数值模拟采用重叠网格技术，利用（）湍流模型、模型计算当前时刻瞬态流场，通过对航行体表面压力积分，获得当前时刻航行体所受水动力，在此基础上计算航行体的水动力系数等；利用迭代获得下一时刻的航行体边界和网格，进而求解下一时刻的瞬态流场，接着计算下一时刻航行体的加速度、速度和位移等得到航行体运动特性。上述过程不断循环，直至航行体完全出水为止。为简化模型和计算，做如下假设：）流体常物性，不可压缩；）流体在航行体运动之前是静止的；）航行器是刚体，其外形关于、平面对称；）不考虑其偏航运动；）不考虑空泡的影响。流体控制方程流体力学的控制方程为 个基本的守恒方程，即连续性方程（质量守恒方程）、动量守恒方程（运动方程）和能量守恒方程，其中能量守恒方程是包含有热交换的流动问题必须满足的基本方程，因为本文所做数值计算中涉及到的流体（海水）是不可压缩流体，故将其忽略，数值计算时只需要考虑质量守恒方程和动量守恒方程即可。考虑到本文中的流体处于湍流状态，所以还需遵守附加的湍流输运方程。和 等人分别将 和 模型应用于近壁面和远流场区。本文应用的湍流模型为 湍流输运方程。连续性方程：（）非定常 方程：（）（）()（）湍流动能方程：（）（）（）（）湍流动能耗散率方程：（）（）（）（）式（）（）中：为流体粘度；为湍流动能普朗特数；为湍流耗散普朗特数；和 分别为 和 的耗散；为湍流耗散率的源项；为交叉扩散项；为湍流黏度。几何模型及边界和初始条件将整个计算域划分为 个区域：一个是数值波浪水池区域（背景区域）；另一个为航行体区域（子区域）。航行体区域设置相对比较简单，对于包裹体（航行体的重叠区域）的外部边界设置为重叠网格边界。由于没有波浪和海流，因此不需要考虑波浪和海流的进出口，边界设置相对简单，波浪水池四周设置为对称面，底部设置为滑移壁面，顶部设置为压力出口。边界条件设定如图 所示。对于模拟输入的初始条件，我们采用了实验测定的 时的速度值，即：。图 边界条件设定 重叠网格模型重叠网格的划分主要分为 个区域，即背景区域和航行体区域（子区域）。背景区域在自由表面附近划分很细的网格，用来捕捉航行体穿越自由液面时自由液面的细节变化，在竖直方向划分出与子区域尺寸一致的切割体网格，用来更好地与子网格装配。背景计算域网格划分完成后，对子区域进行网格划分，并对其近壁面网格细化。子网格与背景网格重叠，子网格通过插值向背景网格传递计算结果，加快计算效率和提高计算精度。网格划分如图 所示。兵 器 装 备 工 程 学 报：图 背景区域和航行体区域网格划分 结果与分析将实验数据与仿真计算的结果整理，绘制速度曲线如图 所示，航行体表面压力曲线图如图 所示。从速度曲线可知，航行体在水下速度下降显著，水下 内速度实验值由 下降至 ，计算值由 下降至 。从压力曲线分布可知，头部驻点压强随着航行体速度减小而降低；头锥段压力，因为绕流作用和环境压力的下降而呈现出较稳定的下降趋势，后头部出水，驻点压强呈现明显的下降趋势，出水前达到最大压力为 。仿真计算结果与实验吻合较好，验证了用重叠网格方法处理航行体运动的准确性和可靠性。图 航行体运动速度曲线 图 航行体表面压力曲线 取航行体在水下 截面压力云图如图 所示，由图 可以看出，航行体运动对整个流场压力的分布影响很小，其主要影响集中在航行体周围，可以看出明显的高压和低压区。其中最明显的高压区包括航行体顶部高压区，达到 ，是阻力的主要来源，低压区主要表现在航行体肩部附近的低压区和尾端，分别达到 和 。图 时航行体表面压力云图 航行体出水过程数值仿真如图 所示，其中 指的是航行体头部出水时刻；是肩部出水时刻；是因为肩部到主段之间承接处出水时刻；指的是航行体尾部刚离开自由液面时刻；时刻拐点是航行体出水后完全摆脱附连水时刻。由图 可知，航行体出水产生水冢，还会穿越水面而形成水冢破碎，之后随着航行体出水的距离增加，水冢还会增长一段时间，在达到最大水冢之后，水面开始沿着航行体表面回落，完全出水后仅有少量附连水附着于航行体表面，附连水引起的粘滞力是航行体完全出水后主要的流体阻力来源。图 航行体出水过程数值仿真 同时，绘制水动力曲线如图 所示，水动力系数曲线如图 所示。其中纵轴代表力的大小和方向。由图、图 可以看出，从航行体头部到完全出水过程中，航行体水动力及系数在 时刻达到最大值，分别为 和 。在整个过程中，水动力及力系数主要是沿 轴方向变化明显，方向波动很小，基本可以忽略。在航行体出水过程中，水动力和力系数波动变化拐点主要出现在整个航行体出水过程，水下段水动力和力系数变化相对稳定。但从筒口效应到出韩守根，等：基于重叠网格的航行体水动力特性计算水过程中，水动力及系数在 最大，这主要是因为水的静压所导致的。图 水动力曲线 图 水动力系数曲线 结合图，分析航行体出水过程水力和力系数变化曲线上各个拐点出现的原因。时刻出现的拐点是因为航行体头部出水导致的，这是航行体出水过程的开始时刻，水动力以及水动力系数达到最大值（从头部接触水面到完全出水）；时刻出现的拐点是由于肩部出水导致受力发生变化；时刻拐点是因为肩部到尾部之间有段承接处几何形状改变导致的；时刻拐点是由于航行体尾部出水，这是航行体出水过程的结束时刻。此时航行体已经出水，所受浮力变为零，出水时刻速度较低，其空气阻力又特别小，可以忽略不计，但此时航行体水动力不为零，这是由于航行体出水时其表面所附带的水沾湿产生的粘滞力导致的，粘滞力是航行体刚出水后主要的流体阻力来源。时刻拐点是因为航行体出水后完全摆脱水的影响，而空气阻力又特别小，可以忽略不计，因此所受到的流体力及力系数趋于零。结论）航行体在水下时顶部压力是航行体阻力的主要来源力。水动力及系数呈现较为稳定的上升趋势，而速度下降显著，计算值由 下降 。）出水过程水动力及力系数主要受水冢现象、浮力变化以及附连水质量变化影响而呈现波动并逐渐趋向于零。）粘滞力是航行体完全出水后主要的流体阻力来源。参考文献：倪火才 潜载导弹水下发射技术的发展趋势分析 舰载武器，（）：，（）：，（）：，：，（）：，（）：，（）：，（）：崔可 航行体出筒以及水下运动的水动力研究 哈尔滨：哈尔滨工业大学，：，（）：，（）：（下转第 页）兵 器 装 备 工 程 学 报：，（）：，：杜飞，姚舜才，安坤，等 新型直爪混合励磁爪极电机设计与研究 机电工程，（）：，（）：倪有源，何强，陈俊华 爪极发电机负载漏磁计算 微特电机，（）：，（）：，张继鹏，陈鹏，苏锦智，等 一种三维磁路永磁电机的集中参数磁路模型 电机与控制学报，（）：，（）：，（）：，张卓然，王东，花为 混合励磁电机结构原理、设计与运行控制技术综述及展望 中国电机工程学报，（）：，（）：，科学编辑 徐永能 博士（南京理工大学副教授、硕导）责任编辑 胡君德（上接第 页），（）：，周超杰，洪亮，张周康，等 重叠网格在多浮体结构 中的应用 兵器装备工程学报，（）：，（）：孙士丽，孙义龙，胡竞中，等 圆柱体出水运动的自由面效应及水冢现象分析 计算物理，（）：，（）：，（）：杨晓光，陈焕龙，刘华坪，等 导弹水下运动及出水过程的三维流场仿真 弹道学报，（）：，（）：科学编辑沈超 博士（中国船舶科学研究中心工程师）责任编辑 徐佳忆杨宇哲，等：车用爪极发电机新型三维分段等效磁路法研究