超空泡航行体流体动力延迟特性数值模拟研究.pdf
《超空泡航行体流体动力延迟特性数值模拟研究.pdf》由会员分享,可在线阅读,更多相关《超空泡航行体流体动力延迟特性数值模拟研究.pdf(10页珍藏版)》请在咨信网上搜索。
1、 年 月第 卷第 期西 北 工 业 大 学 学 报 :收稿日期:基金项目:国家重点实验室基金()资助作者简介:刘喜燕(),西北工业大学博士研究生,主要从事航行器水动力学研究。通信作者:袁绪龙(),西北工业大学副教授,主要从事跨介质航行力学研究。:超空泡航行体流体动力延迟特性数值模拟研究刘喜燕,袁绪龙,罗凯,杜兆星(西北工业大学 航海学院,陕西 西安)摘 要:由于超空泡形态的延迟效应,超空泡航行体在机动航行过程中的流体动力会展现出极强的非定常特性,给航行体流体动力的获取、运动性能的分析以及控制方案的设计带来了巨大的挑战,因此开展航行体超空泡状态下的流体动力特性和形成机理的研究具有重要的意义。采用
2、动态网格技术,建立了超空泡航行体连续变攻角运动的三维数值模型,并通过水洞试验结果验证了该模型的准确性,仿真研究了不同速度、摆动频率和空化器预置舵角下超空泡航行体非定常空泡发展过程和流体动力特性变化规律。研究表明:超空泡航行体非定常运动状态下的流体动力具有延迟特性,采用常规航行体流体动力计算方法预报结果存在较大误差;空泡延迟效应和附体空泡的动态发展是导致超空泡航行体非定常流体动力延迟特性形成的根本原因;超空泡航行体流体动力延迟特性随着速度的增大而显著增强、摆动频率增大而减弱,空化器预置舵角导致航行体在正负攻角摆动时呈现不同的延迟特性。关 键 词:超空泡航行体;非定常流体动力;延迟特性;附体空泡;
3、数值模拟中图分类号:文献标志码:文章编号:()由于空泡延迟效应的存在,当超空泡航行体改变攻角时,新生成的空泡需经过一段时间才能到达航行体尾部。空泡状态的更新相对于航行体姿态的变化存在一定的时间滞后,沾湿面积改变,从而引起流体动力的变化。空泡与航行体之间相对位置关系的实时变化,导致超空泡航行体机动航行过程中的流体动力展现出极强的复杂性和非定常特性,从而使超空泡航行体弹道预报和稳定控制充满挑战。因此,研究超空泡航行体流体动力的延迟特性,具有十分重要的理论与工程应用价值。目前国内外学者围绕超空泡航行体机动航行过程中的空泡形态和非定常流体动力特性进行了大量研究。李魁彬等基于空泡截面独立扩张原理、空泡镜
4、像模型和航行体动力学方程,推导了恒定推力作用和匀变速情况下超空泡形态的解析解,探索了记忆效应对空泡形态的影响,获得了空化数、空化器尺寸和加速度对记忆效应的影响规律。周景军等基于 方法研究了超空泡航行体加速和操舵过程中的流体动力特性,获得了航行体空化器、尾舵操舵、通气量、重力效应以及航行体攻角等参数变化过程中的空泡形态及航行体瞬态流体动力特性变化规律。李杰等通过建立超空泡航行体自由运动的非定常流场求解数值模拟方法,对尾拍流体动力的形成和演化规律进行了研究,发现尾拍力与空泡壁面形状、航行速度及尾部沾湿面积等密切相关。陈诚通过建立考虑空泡延迟效应的超空泡航行器动力学模型,对航行器的尾部位置力特性进行
5、了仿真研究,得到了空化器舵角导致流体动力出现附加攻角的规律。王威、何广华等采用数值方法分别研究了垂直方向和水平方向的周期性来流扰动对通气超空泡航行体流体动力特性的影响,分析了周期性来流的波长和波幅对空泡形态和流体动力的影响规律。罗凯等通过建立定量求解通气超空化的三维数值模型,研究了通气量变化对空泡形态的影响,获得了通气量突变会使得空泡特征呈现不同程度时间迟滞特性的规律。西 北 工 业 大 学 学 报第 卷综上所述,公开的文献主要集中在考虑空泡延迟效应对动力学建模和流体动力特性的研究,缺乏对非定常流体动力延迟特性形成机理的研究。目前,课题组已通过搭建连续变攻角测力试验平台,定量分析了航行速度、摆
6、动频率及空化器预置舵角等因素对超空泡航行体连续变攻角过程中延迟时间特性的影响。在此研究基础上,本文通过建立连续变攻角运动数值模型,结合对空泡与流体动力变化的分析,进一步揭示不同影响因素下非定常流体动力延迟特性的形成机理,为超空泡航行体动力学建模和机动航行控制技术提供参考。计算模型与方法 连续变攻角运动模型建立图 给出了连续变攻角运动模型坐标系示意图。设置地面坐标系为,初始时刻航行体质心位置为地面坐标系原点,且 平面与航行体的对称面重合,航行体中轴线与 轴重合,空化器指向 轴负方向。设置航行体的转轴中心 位于航行体质心 后 处(水洞试验中尾支撑转动位置),旋转轴正方向垂直于对称面向外,建模过程中
7、不考虑重力的影响。航行体的摆动规律为 ()()()()式中:为摆动角度;为摆动的最大角度,本文中;为摆动频率;为角速度。图 坐标系定义 控制方程本文通过数值模拟和水洞试验验证相结合的方法开展研究。受制于试验条件,水洞中的水速最多只能达到 ,形成自然空化能力较弱,因此采用主动通气实现超空泡流动。借鉴目前通气超空泡数值模拟研究经验,采用流体体积函数()多相流模型进行水、汽两相流动问题模拟研究。)连续性方程()()式中:,为密度;为体积分数;下表,分别代表混合相、水相和空气相;为混合密度的平均速度,其表达式为 ;为第 相体积分数;为第 相的密度;为第 相的速度,。)动量方程()()()式中:为压力;
8、为源项;为剪切应力,为混合物黏性,为第 相动力学黏度。)湍流方程本文选取 湍流模型,它在更广泛的流动领域中具有更好的鲁棒性和求解精度。()()()()()()式中:为湍流动能;为湍流耗散率;和分别为 和 的有效扩散系数;和 分别为湍流导致的 和 的耗散;为交叉扩散项;,为湍流产生项。计算模型与边界条件设置航行体模型如图 所示,主要由圆盘空化器、圆锥段、圆柱段以及尾喷管组成。图 航行体模型示意图第 期刘喜燕,等:超空泡航行体流体动力延迟特性数值模拟研究航行体空化器直径 ,航行体圆柱段直径 ,尾喷管直径 ,全长 ,锥段长 ,航行体质心位置距离空化器前端 。计算域及边界条件设置如图 所示,计算域为长
9、方体,考虑流域边界效应对空泡发展的影响,结合作者对壁面效应的研究经验,计算域尺寸长、宽、高分别为,航行体放置于对称面上,空化器头部距离计算域左侧为(为空泡最大长度),尾喷管端面距离计算域右侧为。计算域左、上、下面设置为速度入口,前、后面以及航行体设置为壁面,右侧面设置为压力出口。图 计算域及边界条件本文采用 软件进行结构化网格划分,重点对空化器、航行体肩部及尾流空泡区的网格进行加密处理。通过研究 套不同单元数量的网格(万,万,万)对计算精度和效率的影响,发现网格数量由 万增加到 万,最大空泡长度增加,最大空泡直径增大,而当网格数量由 万增加到 万,最大空泡长度增加,最大空泡直径增加,综合分析,
10、最终确定网格数量为 万,网格划分细节如图 所示。图 网格划分细节 计算方法本文中超空泡航行体做受迫连续变攻角运动,模型不仅具有相对于水速的向前速度,还同时有转动速度。采用动计算域法实现航行体的位置和姿态更新,用 控制超空泡航行体按照给定的正弦规律进行摆动,同时将计算获得的航行体受力结果输出至文件。超空泡航行体流场与运动耦合仿真基于 软件实现,压力速度耦合算法选用 算法,压力项离散格式选用 格式,动量、体积分数、湍流动能和湍流耗散率的离散选用一阶迎风格式,计算时间步长保持在。计算流程:内,航行体保持静止,使空泡充分发展;内航行体连续往复摆动若干周期,然后航行体停止运动;最后 空泡恢复原状态。计算
11、结果与分析 连续变攻角测力水洞试验系统依托西北工业大学高速水洞完成连续变攻角测力试验平台的搭建。试验系统如图 所示,由攻角调节与测量系统、通气流量控制与测量系统、流体动力测量系统、压力测量系统、图像采集系统及模型主体组成。模型以尾支撑方式安装于水洞工作段内的电动转台上,试验开始前,使模型保持在 攻角位置上。通过 编写综合控制软件,由单片机控制步进电机,实现模型按设定规律连续转动。通气管路、测压管路和天平导线从导流罩内孔引出至洞体外,压缩空气经流量控制器接至模型通气管。图 试验系统构成 数值模型验证)试验结果重复性验证为减小试验中通气流量和水洞水速脉动变化带来的影响,确保结果准确性,试验工况均做
12、 组西 北 工 业 大 学 学 报第 卷重复试验。图 给出水速 ,频率,预置舵角 工况时,模型摆动方向受力的数据处理前后对比情况。由图 可以看出,组数据的重合度较好,试验的重复性良好,可以说明试验结果具有较高的可靠性。图 试验重复性结果)数值结果与试验结果对比验证图 给出了通气流量 为 ,空化数为 时的试验与仿真计算结果对比。由图 可以看出,参照试验空泡结果,仿真空泡最大直径相对误差为,空泡最大长度相对误差为,说明建立的数值模型基本合理。图 水洞试验与数值仿真结果对比图 给出了连续变攻角运动过程中摆动方向所受流体动力相对于运动相位之间的延迟效应示意图。计算各个运动周期 内攻角曲线峰值与流体动力
13、峰值之间的时间差,取多个周期内时间差值求均值作为该工况下的流体动力延迟时间。为了方便对比,对延迟时间 进行无量纲化处理,定义为延迟系数,其表达式为()式中:为无量纲参数;,其中,为航行体总长;为参考速度。图 延迟时间示意图图 给出了来流速度 分别为,、摆动频率 均为 、保持相同通气量 、预置舵角 等不同工况水洞试验获得的延迟系数 与数值仿真的对比结果。由图 可以看出,仿真结果与试验结果变化规律一致,延迟系数随着速度增加而增大。数值计算得到的延迟系数与试验结果的最大相对误差为,进一步证明建立的数值模型具有较好的准确性。图 延迟系数 随速度 的变化规律 速度影响分析图 给出了水速 ,空化数,改变攻
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 空泡 航行 流体 动力 延迟 特性 数值 模拟 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。