留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于体积力模型的潜艇应急上浮运动数值模拟分析

魏可可 高霄鹏

魏可可, 高霄鹏. 基于体积力模型的潜艇应急上浮运动数值模拟分析[J]. 中国舰船研究, 2021, 16(2): 1–8 doi:  10.19693/j.issn.1673-3185.01879
引用本文: 魏可可, 高霄鹏. 基于体积力模型的潜艇应急上浮运动数值模拟分析[J]. 中国舰船研究, 2021, 16(2): 1–8 doi:  10.19693/j.issn.1673-3185.01879
WEI K K, GAO X P. Numerical simulation of emergency surfacing motion of submarine based on volumetric force model[J]. Chinese Journal of Ship Research, 2021, 16(2): 1–8 doi:  10.19693/j.issn.1673-3185.01879
Citation: WEI K K, GAO X P. Numerical simulation of emergency surfacing motion of submarine based on volumetric force model[J]. Chinese Journal of Ship Research, 2021, 16(2): 1–8 doi:  10.19693/j.issn.1673-3185.01879

基于体积力模型的潜艇应急上浮运动数值模拟分析

doi: 10.19693/j.issn.1673-3185.01879
基金项目: 国家自然科学基金资助项目(51479207)
详细信息
    作者简介:

    魏可可,男,1992年生,博士生。研究方向:船舶水动力学。E-mail:1229726215@qq.com

    高霄鹏,女,1971年生,博士,副教授。研究方向:船舶水动力学。E-mail:18627093611@qq.com

    通讯作者:

    魏可可

  • 中图分类号: U661.3

Numerical simulation of emergency surfacing motion of submarine based on volumetric force model

  • 摘要:   目的  为了探索潜艇应急上浮运动的规律,基于体积力的螺旋桨简化模型,开展不同航速下潜艇的应急上浮运动数值模拟。  方法  对艇体水下定深直航、上浮以及上浮出水进行数值模拟,对比分析不同航速下潜艇上浮过程中艇的螺旋桨转速与航速的匹配,以及上浮运动时间和艇体姿态的变化。  结果  结果显示,随着转速的增大,潜艇获得稳定航速的时间越短,在水下定深直航运动的时间也越短;无论是纵倾角还是横倾角,其第1次发生变化的时刻、产生峰值的时刻等随转速的增大呈减小趋势,同时随着转速的增大,纵倾角的来回震荡波动时间会越来越大,而横倾角的来回震荡波动时间则越来越小。  结论  相关计算方法和研究结果可为潜艇应急上浮运动研究提供一定的实际参考价值。
  • 图  1  艇体三维模型

    Figure  1.  Three dimensional model of hull

    图  2  计算域及边界条件

    Figure  2.  Computational domain and boundary conditions

    图  3  艇体整体网格及域剖面网格

    Figure  3.  Overall grid and domain profile grid of hull

    图  4  不同漂角下艇体纵向水动力的试验值与数值计算值的对比

    Figure  4.  Comparison between experimental and numerical values of longitudinal hydrodynamics of the hull at different drift angles

    图  5  不同转速下艇体阻力和螺旋桨推力时历曲线

    Figure  5.  Time history curves of hull resistance and propeller thrust in different rotate speed

    图  6  不同转速下艇体航速的时历曲线

    Figure  6.  Time history curves of hull velocity in different rotate speed

    图  7  不同转速下艇体的垂向位移运动时历曲线

    Figure  7.  Time history curves of vertical displacement motion in different rotate speed

    图  8  不同转速下艇体纵倾角时历曲线

    Figure  8.  Time history curves of trim angle in different rotate speed

    图  9  不同转速下艇体横倾角时历曲线

    Figure  9.  Time history curves of heeling angle in different rotate speed

    表  1  模型主参数

    Table  1.   Main parameters of model

    参数数值
    艇长L/m4.356
    艇宽B/m0.508
    水下排水体积V/m30.708
    水下排水质量M/kg704.045
    重心位置/m(2.009,0,0.002)
    惯量/(kg·m−2)(22.6,693.562,693.562)
    下载: 导出CSV

    表  2  DTMB 4382螺旋桨主参数

    Table  2.   Main parameters of DTMB 4382 propeller

    参数数值
    剖面翼x型NACA
    直径D/mm270
    毂径比d/D0.2
    盘面比Ae/Ao0.725
    侧斜角/(°)36
    有否纵倾
    设计进速系数0.889
    下载: 导出CSV

    表  3  计算工况

    Table  3.   Calculation conditions

    工况上浮力与艇体排水
    质量的比值/%
    艇体重/kg上浮力/kg螺旋桨转速
    /(r·min−1)
    15668.835.2300
    25668.835.2400
    35668.835.2500
    45668.835.2600
    下载: 导出CSV

    表  4  潜艇定深直航运动计算结果

    Table  4.   Calculational results of fixed depth direct sailing motion of sumbmarine

    转速/(r·min−1)t1/st2/sU/(m·s−1)
    3006.5499.980.97
    4006.1590.051.35
    5005.7484.111.72
    6005.6170.022.10
    下载: 导出CSV

    表  5  不同转速下纵倾角的分析

    Table  5.   Analysis of trim angle in different rotate speed

    转速
    /(r·min−1)
    纵倾角第1次
    产生变化时刻/s
    纵倾角第1峰值
    产生时刻/s
    纵倾角第1次
    峰值/(°)
    来回震荡波动
    时长/s
    300100.2109−2229.8
    4009299−2347
    5009196−2247
    6007175−2754
    下载: 导出CSV

    表  6  不同转速下横倾角的分析

    Table  6.   Analysis of heeling angle in different rotate speed

    转速
    /(r·min−1)
    横倾角第1次
    产生变化时刻/s
    横倾角第1峰值
    产生时刻/s
    横倾角第1次
    峰值/(°)
    来回震荡波动
    时长/s
    30010110828.587
    400929922.379
    500889818.173
    600707822.269
    下载: 导出CSV
  • SMALLWOOD D A, WHITCOMB L L. Adaptive identification of dynamically positioned underwater robotic vehicles[J]. IEEE Transactions on Control Systems Technology, 2003, 11(4): 505–515. doi:  10.1109/TCST.2003.813377
    WATT G D. Estimating the hydrodynamic characteristics of the dorado remote minehunting vehicle[R]. [S. l.]: Defence Research and Development Canada, 2002.
    KIM J, KIM K, CHOI H S, et al. Estimation of hydrodynamic coefficients for an AUV using nonlinear observers[J]. IEEE Journal of Oceanic Engineering, 2002, 27(4): 830–840. doi:  10.1109/JOE.2002.805098
    WATSON K P, WEBSTER J S, CRANE J W, et al. Prediction of submersible maneuvering performance at high incidence angles[C]//Proceedings of OCEANS'93. Victoria, BC, Canada: IEEE, 1993: 289–294.
    贾欣乐, 杨盐生. 船舶运动数学模型——机理建模与辨识建模[M]. 大连: 大连海事大学出版社, 1999.

    JIA X L, YANG Y S. Mathematical model of ship motion[M]. Dalian: Dalian Maritime University Press, 1999 (in Chinese)..
    戴余良, 俞科云. 潜艇应急上浮运动稳定性与分叉分析[J]. 上海交通大学学报, 2010, 44(10): 1400–1404, 1413.

    DAI Y L, YU K Y. Stability and bifurcations analysis of submarine during emergency ascent[J]. Journal of Shanghai Jiaotong University, 2010, 44(10): 1400–1404, 1413 (in Chinese).
    陈志法, 戴余良, 林雄伟, 等. 潜艇水下高速运动稳定性分析[J]. 舰船科学技术, 2015, 37(3): 15–20. doi:  10.3404/j.issn.1672-7649.2015.03.003

    CHEN Z F, DAI Y L, LIN X W, et al. Motion stability analysis of high-speed submarines under water[J]. Ship Science and Technology, 2015, 37(3): 15–20 (in Chinese). doi:  10.3404/j.issn.1672-7649.2015.03.003
    王晓玢, 孙尧, 莫宏伟. 潜艇操纵运动分叉突变特性[J]. 工程力学, 2009, 26(10): 252–256.

    WANG X F, SUN Y, MO H W. Bifurcation and catastrophe characteristics in submarine motion[J]. Engineering Mechanics, 2009, 26(10): 252–256 (in Chinese).
    CARRICA P M, ISMAIL F, HYMAN M, et al. Turn and zigzag maneuvers of a surface combatant using a URANS approach with dynamic overset grids[J]. Journal of Marine Science and Technology, 2013, 18(2): 166–181. doi:  10.1007/s00773-012-0196-8
    CARRICA P M, SADAT-HOSSEINI H, STERN F. CFD analysis of broaching for a model surface combatant with explicit simulation of moving rudders and rotating propellers[J]. Computers & Fluids, 2012, 53: 117–132.
    CHASE N, MICHAEL T, CARRICA P M. Overset simulation of a submarine and propeller in towed, self-propelled and maneuvering conditions[J]. International Shipbuilding Progress, 2013, 60(1/2/3/4): 171–205.
    周广礼, 欧勇鹏, 高霄鹏. 潜艇上浮出水运动及粘性流场直接数值模拟方法[J]. 海军工程大学学报, 2017, 29(4): 86–92.

    ZHOU G L, OU Y P, GAO X P. Numerical simulation method of surfacing maneuvers and flow field of submarine[J]. Journal of Naval University of Engineering, 2017, 29(4): 86–92 (in Chinese).
    周广礼, 董文才, 欧勇鹏. 潜艇应急上浮六自由度运动及黏性流场数值模拟[J]. 国防科技大学学报, 2017, 39(2): 199–206. doi:  10.11887/j.cn.201702030

    ZHOU G L, DONG W C, OU Y P. Numerical simulation of six degree of freedom motion and viscous flow for submarine's emergency ascent[J]. Journal of National University of Defense Technology, 2017, 39(2): 199–206 (in Chinese). doi:  10.11887/j.cn.201702030
    廖欢欢, 庞永杰, 李宏伟, 等. 基于重叠网格技术的潜艇应急上浮空间运动的数值模拟[C]//第二十七届全国水动力学研讨会文集. 南京, 2015.

    LIAO H H, PANG Y J, LI H W, et al. Numerical simulation based on the overset grid technique about submarine emergency ascent space movement[C]//The 27th National Hydrodynamics Symposium. Nanjing, 2015 (in Chinese).
    WILCOX D C. Formulation of the k-w turbulence model revised[J]. AIAA Journal, 2008, 46(11): 2823–2838. doi:  10.2514/1.36541
    MENTER F R, KUNTZ M, LANGTRY R. Ten years of industrial experience with the SST turbulence model[C]//Turbulence, Heat and Mass Transfer. Antalya: Begell House, 2003.
    FELDMAN J. DTNSRDC Revised standard submarine equations of motion[R]. Taylor Naval Ship Research and Development Center, Netherlands, 1979.
  • [1] 张海洋, 谷海涛, 林扬, 孙原, 高浩, 冯萌萌2020: 无动力运载器倾斜爬升式上浮特性分析 , 中国舰船研究, 15, 38-47. doi: 10.19693/j.issn.1673-3185.01729
    [2] 孙斌, 黄祥兵2020: 潜艇大型集体逃生舱上浮速度控制仿真计算 , 中国舰船研究. doi: 10.19693/j.issn.1673-3185.01847
    [3] 谭骏怡, 胡俊华, 陈国明, 杨健, 葛阳2019: 水空跨介质航行器斜出水过程数值仿真 , 中国舰船研究, 14, 104-121. doi: 10.19693/j.issn.1673-3185.01519
    [4] 章奇, 司梦婷, 王东旭, 王安振, 李灿, 王振, 罗自学, 娄春, 程强2019: 多光谱反演火焰温度及烟黑体积分数研究 , 中国舰船研究, 14, 42-48. doi: 10.19693/j.issn.1673-3185.01413
    [5] 陈康, 姚志崇, 周恩东, 夏彬2019: 基于体积力的海水自流循环系统流速确定方法 , 中国舰船研究, 14, 70-76. doi: 10.19693/j.issn.1673-3185.01480
    [6] 陈柱, 徐国华, 王冠学, 刘炎, 翟云峰, 郑煜2018: AUV主动应急自救机制与策略 , 中国舰船研究, 13, 120-127. doi: 10.19693/j.issn.1673-3185.01027
    [7] 叶礼裕, 王超, 郭春雨, 常欣2018: 潜艇破冰上浮近场动力学模型 , 中国舰船研究, 13, 51-59. doi: 10.3969/j.issn.1673-3185.2018.02.007
    [8] 李日杰, 陈庆龙, 常文田2018: 潜艇应急上浮运动的CFD预报 , 中国舰船研究, 13, 36-42. doi: 10.19693/j.issn.1673-3185.01070
    [9] 魏成柱, 易宏, 李英辉2017: 新概念高速穿梭艇系列船型及其直航性能 , 中国舰船研究, 12, 12-21. doi: 10.3969/j.issn.1673-3185.2017.02.002
    [10] 陈林, 熊海峰, 魏青2016: 深V船型大型化应用分析 , 中国舰船研究, 11, 9-13. doi: 10.3969/j.issn.1673-3185.2016.05.002
    [11] 舒敏骅, 陈科, 尤云祥, 胡天群, 刘恒2016: 动网格方法在螺旋桨非定常轴承力数值计算中的应用 , 中国舰船研究, 11, 25-31. doi: 10.3969/j.issn.1673-3185.2016.03.005
    [12] 王慧婷, 毕毅2016: 基于体积力法的全附体KCS型船模PMM运动数值模拟 , 中国舰船研究, 11, 29-37,66. doi: 10.3969/j.issn.1673-3185.2016.04.005
    [13] 董润鹏, 董文才, 姚朝帮2015: 深V型船零速横摇改善措施 , 中国舰船研究, 10, 15-20. doi: 10.3969/j.issn.1673-3185.2015.06.003
    [14] 傅江妍, 陈炜, 葛慧晓, 陈韬颖2015: 加装舭龙骨对千吨级深V船型静水阻力和横摇运动的影响 , 中国舰船研究, 10, 22-26. doi: 10.3969/j.issn.1673-3185.2015.05.004
    [15] 吴召华, 陈作钢, 代燚, 邓德衡2013: 基于体积力法的船舶回转运动水动力数值研究 , 中国舰船研究, 8, 12-19. doi: 10.3969/j.issn.1673-3185.2013.04.003
    [16] 唐友刚, 陶海成, 王榕, 曹菡2012: 系泊油轮与海上平台的碰撞力分析 , 中国舰船研究, 7, 36-40. doi: 10.3969/j.issn.1673-3185.2012.04.007
    [17] 姚熊亮, 张成, 孙士丽2010: 考虑可压缩性及惯性力的油膜力研究 , 中国舰船研究, 5, 33-40. doi: 10.3969/j.issn.1673-3185.2010.06.007
    [18] 冯奇, 迟卫2010: 舰船尾迹气泡上浮运动的一种表征模型 , 中国舰船研究, 5, 27-29. doi: 10.3969/j.issn.1673-3185.2010.03.006
    [19] 钱家昌, 彭旭, 陈明2008: 机械设备对基座激励力估算方法研究 , 中国舰船研究, 3, 48-50. doi: 10.3969/j.issn.1673-3185.2008.01.011
    [20] 夏翔, 李水才, 李慧敏2006: 表面桨产生垂向力的原因初析 , 中国舰船研究, 1, 68-70. doi: 10.3969/j.issn.1673-3185.2006.06.016
  • 加载中
图(10) / 表 (6)
计量
  • 文章访问数:  51
  • HTML全文浏览量:  13
  • PDF下载量:  1122
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-01-07
  • 修回日期:  2020-08-18
  • 网络出版日期:  2020-12-10

基于体积力模型的潜艇应急上浮运动数值模拟分析

doi: 10.19693/j.issn.1673-3185.01879
    基金项目:  国家自然科学基金资助项目(51479207)
    作者简介:

    魏可可,男,1992年生,博士生。研究方向:船舶水动力学。E-mail:1229726215@qq.com

    高霄鹏,女,1971年生,博士,副教授。研究方向:船舶水动力学。E-mail:18627093611@qq.com

    通讯作者: 魏可可
  • 中图分类号: U661.3

摘要:   目的  为了探索潜艇应急上浮运动的规律,基于体积力的螺旋桨简化模型,开展不同航速下潜艇的应急上浮运动数值模拟。  方法  对艇体水下定深直航、上浮以及上浮出水进行数值模拟,对比分析不同航速下潜艇上浮过程中艇的螺旋桨转速与航速的匹配,以及上浮运动时间和艇体姿态的变化。  结果  结果显示,随着转速的增大,潜艇获得稳定航速的时间越短,在水下定深直航运动的时间也越短;无论是纵倾角还是横倾角,其第1次发生变化的时刻、产生峰值的时刻等随转速的增大呈减小趋势,同时随着转速的增大,纵倾角的来回震荡波动时间会越来越大,而横倾角的来回震荡波动时间则越来越小。  结论  相关计算方法和研究结果可为潜艇应急上浮运动研究提供一定的实际参考价值。

English Abstract

魏可可, 高霄鹏. 基于体积力模型的潜艇应急上浮运动数值模拟分析[J]. 中国舰船研究, 2021, 16(2): 1–8 doi:  10.19693/j.issn.1673-3185.01879
引用本文: 魏可可, 高霄鹏. 基于体积力模型的潜艇应急上浮运动数值模拟分析[J]. 中国舰船研究, 2021, 16(2): 1–8 doi:  10.19693/j.issn.1673-3185.01879
WEI K K, GAO X P. Numerical simulation of emergency surfacing motion of submarine based on volumetric force model[J]. Chinese Journal of Ship Research, 2021, 16(2): 1–8 doi:  10.19693/j.issn.1673-3185.01879
Citation: WEI K K, GAO X P. Numerical simulation of emergency surfacing motion of submarine based on volumetric force model[J]. Chinese Journal of Ship Research, 2021, 16(2): 1–8 doi:  10.19693/j.issn.1673-3185.01879
    • 潜艇的应急上浮是指潜艇当发生舱室破损进水、舵卡等事故时,采取排除部分或者全部主压载水舱内的水,从而使潜艇迅速浮至水面的一种紧急操纵措施。为保证潜艇应急上浮的安全性和稳定性,准确预报潜艇在应急上浮过程中的运动特性十分有必要。

      和正常的上浮相比,潜艇应急上浮中艇的排水方式不同且排水速率快,艇体姿态变化幅度大。常见的潜艇应急上浮运动预报方法是通过拘束摸试验[1-4]或参数辨识[5]等方法求解水动力系数,然后再基于六自由度运动模型进行预报。戴余良等[6-7]应用非线性系统运动稳定性与分叉理论,以及同伦延拓数值计算方法,分析了潜艇应急上浮运动的稳态响应及运动稳定性。王晓玢等[8]通过对潜艇垂直面操纵运动进行非线性建模,利用分叉与突变理论方法对潜艇上浮运动中由静态分叉和动态分叉引发的状态突变进行了分析。潜艇水动力系数的求解较复杂且有一定的难度,而上述预报方法主要是针对潜艇某自由度上的运动,并不能真实反映潜艇在六自由度完全释放情况下的应急上浮运动特性,且预报精度还有待提高。

      随着计算机技术的快速发展,计算流体力学在船舶运动研究中的应用越来越广泛。Carrica等[9-10]基于自主研发的流体软件CFDShip-Iowa V4,实现了不同工况下潜艇六自由度运动的数值模拟计算,并对其操纵性进行了较高精度的预报。Chase等[11]采用重叠网格技术,对Suboff潜艇模型的直航拖曳和自由自航等运动进行了数值模拟计算。周广礼等[12-13]基于RANS方程及流体体积模型,采用整体动网格技术,对Suboff模型的应急上浮六自由度运动及其出水运动予以了研究。廖欢欢等[14]采用流体体积法(volume of fluid,VOF)六自由度求解器和重叠网格技术,研究了潜艇在一定潜深下的应急上浮六自由度运动,并对潜艇上浮过程中的姿态和水动力进行了分析。

      以上大多学者对潜艇应急上浮运动的研究都没有考虑螺旋桨以及初始航速等问题的嵌入,为此,本文将基于STAR-CCM+软件平台,采用基于体积力模型的虚拟桨盘,针对潜艇在不同航速下的应急上浮运动进行数值模拟,并对比分析潜艇在应急上浮运动中的姿态和上浮运动时间等,用以为潜艇在实际中的上浮运动操控提供一定的借鉴。

    • 本文将以美国泰勒研究中心提供的潜艇试验模型Suboff为研究对象。国内外针对该模型开展了大量的拖曳及流场测量试验,其三维模型图如图1所示,模型主尺度参数如表1所示。

      图  1  艇体三维模型

      Figure 1.  Three dimensional model of hull

      表 1  模型主参数

      Table 1.  Main parameters of model

      参数数值
      艇长L/m4.356
      艇宽B/m0.508
      水下排水体积V/m30.708
      水下排水质量M/kg704.045
      重心位置/m(2.009,0,0.002)
      惯量/(kg·m−2)(22.6,693.562,693.562)
    • 应用粘流CFD方法开展数值计算的基础是纳维-斯托克斯(Navia-Stokes,N-S)方程,然而受限于计算能力,当前还不能通过直接求解N-S方程(direct numerical simulation,DNS)来获得潜艇周围复杂的粘性流场。将N-S方程中的湍流脉动项进行时均化处理可以极大地减小计算量,所得雷诺平均方程(RANS)的具体表达如式(1)所示。然而,N-S方程在时均化的过程中增加了含未知量时均的雷诺应力项,通过引入湍流模型,可使方程封闭。目前,该方法已在舰艇流体力学理论研究及工程实践中得到广泛应用。

      $$\left\{ \begin{aligned} & \frac{{\partial {{\bar u}_i}}}{{\partial t}} \!+\! {{\bar u}_j}\frac{{\partial {{\bar u}_i}}}{{\partial {x_j}}}\! = \! -\! \frac{1}{\rho }\frac{{\partial \bar p}}{{\partial {x_i}}} \!+ \!\frac{1}{\rho }\frac{\partial }{{\partial {x_j}}}\left(\mu \frac{{\partial {{\bar u}_i}}}{{\partial {x_j}}} \!-\! \rho \overline {{{u'}_i}{{u'}_j}}\right) + {f_i} \\& \frac{{\partial {{\bar u}_i}}}{{\partial {x_i}}} = 0 \end{aligned} \right.$$ (1)

      式中:$\rho $为流体密度;$\mu $为流体粘性;$p$为静压;${f_i}$为单位质量的质量力;${u_i}$${u_j}$分别为$i,j$点上$x$方向的速度分量;t为时间;${x_i}$${x_j}$分别为$i,j$点上$x$方向的位移分量;上加线“―”表示参量时均化;$ - \rho \overline {{{u'}_i}{{u'}_j}} $为雷诺应力。

      开展计算时,控制方程离散均采用二阶迎风格式,并应用分离式解法(segregated flow)对离散后的方程组进行求解。压力速度耦合迭代采用SIMPLE算法,涉及非稳态计算时,设定计算过程中的时间步长为$\Delta t$且固定不变,时间迭代采用二阶形式,即

      $$\frac{{{\rm{d}}\chi }}{{{\rm{d}}t}} = \frac{{3 \times \left[ {\chi (n + 1) - \chi (n)} \right] - \left[ {\chi (n) - \chi (n - 1)} \right]}}{{2 \times \Delta t}}$$ (2)

      式中:$\;\chi = \int\limits_{\Delta v} {\rho u\phi {\rm{d}}V} $,为变量$\rho u\phi $在控制体积上的积分,其中$V$为合速度,$u$$x$方向的速度分量,$\phi $是广义变量;$\;\chi (n)$为第$t = n \times \Delta t$时刻$\;\chi $的量值;$n$为时间步数。

      采用RANS方法开展潜艇水动力计算时,湍流模型的选取是确保计算精度的关键要素,众多学者已针对该问题开展了大量研究,但所得结论存在较大的不同,且湍流模型的适用性也与网格布局及操作者的数值计算经验相关。综合各类文献资料,本文选取Wilcox $k - \omega $与SST $k - \omega $湍流模型进行优选分析。

      Wilcox $k - \omega $湍流模型最早由Wilcox于1988年提出,后经过不断的修正,模拟精度及其适用范围均得到了提升,STAR-CCM+软件平台中就提供有最新的修正模型,具体表达形式如下:

      $$\left\{ \begin{aligned} & \frac{\partial }{{\partial t}}(\rho k) + \frac{\partial }{{\partial {x_j}}}(\rho {u_j}k) = \rho {\tau _{ij}}\frac{1}{\rho }\frac{{\partial {u_i}}}{{\partial {x_j}}} - {\beta ^ * }\rho k\omega +\\&\qquad \frac{\partial }{{\partial {x_j}}}\left[\left(\mu + {\sigma ^ * }\frac{{\rho k}}{\omega }\right)\frac{{\partial k}}{{\partial {x_j}}}\right] \\& \frac{\partial }{{\partial t}}(\rho \omega ) + \frac{\partial }{{\partial {x_j}}}(\rho {u_j}\omega ) = \alpha \frac{\omega }{k}\rho {\tau _{ij}}\frac{1}{\rho }\frac{{\partial {u_i}}}{{\partial {x_j}}} - \beta \rho {\omega ^2} + \\&\qquad{\sigma _{\rm{d}}}\frac{\rho }{\omega }\frac{{\partial k}}{{\partial {x_j}}}\frac{{\partial \omega }}{{\partial {x_j}}} + \frac{\partial }{{\partial {x_j}}}\left[\left(\mu + \sigma \frac{{\rho k}}{\omega }\right)\frac{{\partial \omega }}{{\partial {x_j}}}\right] \end{aligned} \right.$$ (3)

      式中:$k$为湍动能;$\omega $为特殊湍动能耗散;其余参数的具体定义及参数值参见文献[15]。

      Menter等[16]结合标准$k - \varepsilon $$k - \omega $湍流模型,提出了剪切应力模型SST $k - \omega $,式(4)给出了该模型的最新修正形式,具体推导过程及参数定义详见文献[15]。

      $$\left\{ \begin{aligned} & \frac{\partial }{{\partial t}}(\rho k) + \frac{\partial }{{\partial {x_j}}}(\rho {u_j}k) = \widetilde {{P_k}} - {\beta ^ * }\rho k\omega + \frac{\partial }{{\partial {x_j}}}\left[\left(\mu + {\sigma _k}{\mu _{\rm{t}}}\right)\frac{{\partial k}}{{\partial {x_j}}}\right] \\& \frac{\partial }{{\partial t}}(\rho \omega ) + \frac{\partial }{{\partial {x_j}}}(\rho {u_j}\omega ) = \alpha \rho {S^2} - \beta \rho {\omega ^2} + \\&\qquad 2(1 - {F_1})\rho {\sigma _{\omega 2}}\frac{1}{\omega }\frac{{\partial k}}{{\partial {x_j}}}\frac{{\partial \omega }}{{\partial {x_j}}} + \frac{\partial }{{\partial {x_j}}}\left[(\mu + {\sigma _\omega }{\mu _{\rm{t}}})\frac{{\partial \omega }}{{\partial {x_j}}}\right] \end{aligned} \right.$$ (4)

      式中,$\widetilde {{P_k}} = \min \left[{\mu _{\rm{t}}}\dfrac{{\partial {u_i}}}{{\partial {x_j}}}\left(\dfrac{{\partial {u_i}}}{{\partial {x_j}}} + \dfrac{{\partial {u_j}}}{{\partial {x_i}}}\right),10{\beta ^ * }\rho k\omega \right]$

    • 计算域选取为方形域,除域尾端边界设置为压力出口外,其余边界条件均设置为速度入口,域大小及计算域与自由液面的初始相对位置如图2所示。

      图  2  计算域及边界条件

      Figure 2.  Computational domain and boundary conditions

      进行网格划分时,为能更好地捕捉两相流自由液面,对自由液面可能通过的区域进行了网格加密。此外,为准确模拟潜艇出水后艇体周围的流场,对模型首部及指挥室围壳周围网格进行了细化处理。同时,为满足网格贴体的要求,对近壁面区域进行了局部加密,以确保艇体壁面y+在30~300之间[12]。当y+在30~200之间时,表征第1层网格布置于对数层,可通过假定的壁面函数对近壁面流动进行模拟。该方法不直接求解粘性底层与过渡层内的流动,适用性较强,可极大地减少近壁面网格数量。艇体表面、域剖面及近壁面域的网格如图3所示,总网格数为700万。

      图  3  艇体整体网格及域剖面网格

      Figure 3.  Overall grid and domain profile grid of hull

    • 通过对螺旋桨入流面流速进行加权积分,获得桨的进流速度,然后代入到已知螺旋桨水动力曲线中,预估得到推力及扭矩,接着,通过添加源项的方式,对设定的体积力域内的流体施加等效的力及力矩,该方法即为体积力模型。体积力模型不仅能够将力及力矩施加于潜艇上,而且能更为准确地模拟由螺旋桨引起的潜艇尾流场变化,相较于直接力模型,其计算量大,但计算效率远优于开展的真实螺旋桨计算。

      为了重点分析潜艇的运动特征,提高计算效率,螺旋桨用基于体积力模型的虚拟桨代替。该虚拟桨的参数采用DTMB 4383标准桨的,模型主参数如表2所示。

      表 2  DTMB 4382螺旋桨主参数

      Table 2.  Main parameters of DTMB 4382 propeller

      参数数值
      剖面翼x型NACA
      直径D/mm270
      毂径比d/D0.2
      盘面比Ae/Ao0.725
      侧斜角/(°)36
      有否纵倾
      设计进速系数0.889
    • 针对漂角$\beta = 2^\circ \sim 14^\circ $、间隔$2^\circ $的Suboff模型的无因次纵向水动力X'进行数值计算,计算结果与试验值[17]的对比如图4所示。

      图  4  不同漂角下艇体纵向水动力的试验值与数值计算值的对比

      Figure 4.  Comparison between experimental and numerical values of longitudinal hydrodynamics of the hull at different drift angles

      图4可看出,不同漂角下艇体纵向水动力数值计算结果与试验值吻合良好,证明该方法可行。因此,选取该方法进行后续的数值仿真计算。

    • 为探讨不同航速对潜艇上浮运动的影响规律,选取上浮力为总排水量的5%、螺旋桨转速分别为300,400,500,600 r/min的工况进行数值模拟计算,具体计算工况如表3所示。

      表 3  计算工况

      Table 3.  Calculation conditions

      工况上浮力与艇体排水
      质量的比值/%
      艇体重/kg上浮力/kg螺旋桨转速
      /(r·min−1)
      15668.835.2300
      25668.835.2400
      35668.835.2500
      45668.835.2600
    • 假定潜艇上浮之前在水下做定深直航运动,为获取模型水下稳定直航状态,设定初始直航航速U=0,在螺旋桨推力的作用下,潜艇沿x方向作直航运动。为避免在初始计算阶段出现计算结果不稳定的现象,在t=0~4 s时间段内将艇体作静态处理,待螺旋桨推力及流场计算结果稳定后,即在t1时刻,将x方向运动完全释放,潜艇开始直航运动。

      不同转速下螺旋桨推力及艇体阻力的时历曲线如图5所示。由图可知,艇体阻力及螺旋桨推力随时间逐渐向自航值逼近,且二者随时间的绝对变化率较为相近,航速随时间逐步增大并在t2时刻后趋于稳定;随着转速的增大,螺旋桨获得的推力以及艇体的阻力值越大,其运动稳定的时间就越小,最终稳定的航速值也越大,这从图6即可看出。

      图  5  不同转速下艇体阻力和螺旋桨推力时历曲线

      Figure 5.  Time history curves of hull resistance and propeller thrust in different rotate speed

      图  6  不同转速下艇体航速的时历曲线

      Figure 6.  Time history curves of hull velocity in different rotate speed

      表4给出了不同转速下潜艇做定深直航运动时的计算结果。由表4可知,4个转速下t1值间差别较小,在t2时刻是随转速的增大而减小,转速越大,需要稳定的时间就越短;转速为300~600 r/min时,最终稳定的定深下的航速分别为0.97,1.35,1.72和2.1 m/s。

      表 4  潜艇定深直航运动计算结果

      Table 4.  Calculational results of fixed depth direct sailing motion of sumbmarine

      转速/(r·min−1)t1/st2/sU/(m·s−1)
      3006.5499.980.97
      4006.1590.051.35
      5005.7484.111.72
      6005.6170.022.10
    • 图7给出了不同转速下艇体垂向位移运动的时历曲线。由图可知,艇体做应急上浮至出水一般会经历水下定深直航运动—上浮运动—水面航行这3种状态的运动。当螺旋桨转速为300 r/min时,其水下定深直航时间为101 s,上浮时间8 s,上浮出水的时刻为108 s;当转速为400 r/min时,其水下定深直航时间为92 s,上浮时间9 s,上浮出水的时刻为101 s;当转速为500 r/min时,其水下定深直航时间为91 s,上浮时间9 s,上浮出水的时刻为103 s;当转速为600 r/min时,其水下定深直航时间为71 s,上浮时间12 s,上浮出水的时刻为83 s。由此可知,随着转速的增大,潜艇获得稳定航速的时间越短,艇体在水下做定深直航运动的时间也就越短;不同转速下艇体上浮的时间虽然有所差别,但其差值量级较小,上浮运动的时间主要还是受上浮力的影响;同时,潜艇的上浮运动是一个非线性的复杂的六自由度运动,当上浮力一致时,螺旋桨转速越大,其上浮运动速度不一定越大。

      图  7  不同转速下艇体的垂向位移运动时历曲线

      Figure 7.  Time history curves of vertical displacement motion in different rotate speed

    • 图8所示为不同转速下艇体纵倾角的时历曲线,表5给出了纵倾角第1次产生变化的时刻、发生峰值时刻、峰值大小以及来回震荡波动的时间。由图8表5可知,当转速为300 r/min时,纵倾角第1次产生变化的时刻是100.2 s,第1次产生峰值的时刻是109 s,其峰值为−22°,纵倾角来回震荡波动的时间是29.8 s;当转速为400 r/min时,纵倾角第1次产生变化的时刻是92 s,第1次产生峰值的时刻是99 s,其峰值为−23°,纵倾角来回震荡波动的时间是47 s;转速为500 r/min时,纵倾角第1次产生变化的时刻是91 s,第1次产生峰值的时刻是96 s,其峰值为−22°,纵倾角来回震荡波动的时间是47 s;转速为600 r/min时,纵倾角第1次产生变化的时刻是71 s,第1次产生峰值的时刻是75 s,其峰值为−27°,纵倾角来回震荡波动的时间是54 s。结果显示,随着转速的增大,纵倾角第1次产生变化时刻和发生峰值的时刻呈减小趋势,而纵倾角来回震荡波动的时间则越来越大,也即转速越大,艇体出水时受自由液面扰动越小。纵倾角的峰值变化主要发生在潜艇上浮运动过程中,来回震荡波动主要发生在潜艇出水时刻。对比3.2节有关上浮的时刻可知,无论转速多大,潜艇第1次倾角的峰值都是发生在潜艇上浮运动过程中,由对比分析可知,随着转速的增大,纵倾角的峰值并未呈现出明显的变化规律。

      表 5  不同转速下纵倾角的分析

      Table 5.  Analysis of trim angle in different rotate speed

      转速
      /(r·min−1)
      纵倾角第1次
      产生变化时刻/s
      纵倾角第1峰值
      产生时刻/s
      纵倾角第1次
      峰值/(°)
      来回震荡波动
      时长/s
      300100.2109−2229.8
      4009299−2347
      5009196−2247
      6007175−2754

      图  8  不同转速下艇体纵倾角时历曲线

      Figure 8.  Time history curves of trim angle in different rotate speed

    • 图9所示为不同转速下艇体横倾角的时历曲线,表6给出了横倾角第1次产生变化的时刻、发生峰值时刻、峰值大小以及来回震荡波动的时间。由图9表6可知,当转速为300 r/min时,横倾角第1次产生变化的时刻是101 s,第1次产生峰值的时刻是108 s,其峰值为28.5°,横倾角来回震荡波动的时间是87 s;当转速为400 r/min时,横倾角第1次产生变化的时刻是92 s,第1次产生峰值的时刻是99 s,其峰值为22.3°,横倾角来回震荡波动的时间是79 s;当转速为500 r/min时,横倾角第1次产生变化的时刻是88 s,第1次产生峰值的时刻是98 s,其峰值为18.1°,横倾角来回震荡波动的时间是73 s;当转速为600 r/min时,横倾角第1次产生变化的时刻是70 s,第1次产生峰值的时刻是78 s,其峰值为22.2°,横倾角来回震荡波动的时间是69 s。结果显示,随着转速的增大,横倾角第1次产生变化时刻和发生峰值的时刻以及来回震荡波动的时间越来越小。横倾角的峰值变化主要发生在潜艇上浮运动过程中,来回震荡波动主要发生在艇出水时刻。对比3.2节有关上浮的时刻可知,无论转速是多少,潜艇第1次横倾角的峰值均发生在潜艇上浮运动过程中,随着转速的增大,横倾角的峰值并未呈现出明显的变化规律。

      图  9  不同转速下艇体横倾角时历曲线

      Figure 9.  Time history curves of heeling angle in different rotate speed

      表 6  不同转速下横倾角的分析

      Table 6.  Analysis of heeling angle in different rotate speed

      转速
      /(r·min−1)
      横倾角第1次
      产生变化时刻/s
      横倾角第1峰值
      产生时刻/s
      横倾角第1次
      峰值/(°)
      来回震荡波动
      时长/s
      30010110828.587
      400929922.379
      500889818.173
      600707822.269
    • 本文基于体积力模型的虚拟桨盘模型,对不同转速下潜艇的水下定深直航、六自由度上浮、出水以及水面航行进行了数值模拟,并对比分析了不同转速下潜艇螺旋桨与转速的匹配以及其姿态等,该方法可为潜艇的应急上浮运动提供一定的参考。

      通过对比分析不同转速下艇体上浮运动过程中的姿态,可得出:随着转速的增大,潜艇获得稳定航速的时间越短,艇体在水下做定深直航运动的时间也越短;不同转速下其上浮运动的时间相差不大;转速越大,无论是纵倾角还是横倾角,其倾角第1次产生变化时刻和发生峰值的时刻会越来越小;同时随着转速的增大,纵倾角的来回震荡波动时间会越来越大,而横倾角的来回震荡波动时间则越来越小,且倾角的峰值随转速的变化并未呈现出明显的变化规律。

参考文献 (17)

目录

    /

    返回文章
    返回