留言板

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

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

不同外形声呐的水动力及流噪声数值模拟

邓云 傅何琪 广超越 陈威 林永水

邓云, 傅何琪, 广超越, 等. 不同外形声呐的水动力及流噪声数值模拟[J]. 中国舰船研究, 2021, 16(4): 1–7 doi: 10.19693/j.issn.1673-3185.02051
引用本文: 邓云, 傅何琪, 广超越, 等. 不同外形声呐的水动力及流噪声数值模拟[J]. 中国舰船研究, 2021, 16(4): 1–7 doi: 10.19693/j.issn.1673-3185.02051
DENG Y, FU H Q, GUANG C Y, et al. Numerical simulation on hydrodynamic and flow noise of different shape sonars[J]. Chinese Journal of Ship Research, 2021, 16(4): 1–7 doi: 10.19693/j.issn.1673-3185.02051
Citation: DENG Y, FU H Q, GUANG C Y, et al. Numerical simulation on hydrodynamic and flow noise of different shape sonars[J]. Chinese Journal of Ship Research, 2021, 16(4): 1–7 doi: 10.19693/j.issn.1673-3185.02051

不同外形声呐的水动力及流噪声数值模拟

doi: 10.19693/j.issn.1673-3185.02051
基金项目: 中央高校基础研究基金资助项目(WUT:2019IVA059);博士科研启动基金资助项目(WUT:40120268)
详细信息
    作者简介:

    邓云,男,1998年生,硕士生。研究方向:水动力及流噪声。E-mail:1569168520@qq.com

    陈威,男,1988年生,博士,研究方向:海洋工程水动力学。E-mail:whutcw01@126.com

    通讯作者:

    陈威

  • 中图分类号: U661.44; U666.7

Numerical simulation on hydrodynamic and flow noise of different shape sonars

  • 摘要:   目的  为降低传统圆形声呐外表面流噪声,设计了3种不同的声呐外形(圆形、椭圆形以及方形),并对其外表面水动力和流噪声进行了研究。  方法  基于Fluent软件中的标准$k - \varepsilon$湍流计算模型和Lighthill声类比方法对3种不同外形声呐外表面的流场和声场进行分析。  结果  结果显示:方形的声呐升力系数幅值和阻力系数均值最大,其次是圆形,最小的是椭圆形;各方案声呐水动力差异是由于边界层分离点及尾涡的不同引起的。对比流噪声结果发现,圆形的总声压级最大,椭圆形的最小。圆形和方形呈现了“正8字”的偶极子声源特性,声辐射最大值在垂直于来流方向;而椭圆呈现“倒8字”的偶极子声源特性,最大值在水平来流方向。  结论  可用椭圆形形声呐替代传统的圆形声呐,结果作为声呐外形设计的参考。
  • 图  1  (a) 不同外形的声呐计算域图

    Figure  1.  Calculation domain of sonar with different shapes

    图  2  不同外形声呐升阻力系数

    Figure  2.  Lift and drag coefficient of sonar with different shapes

    图  3  不同外形声呐频谱

    Figure  3.  Spectrum diagram of sonar with different shapes

    图  4  不同外形声呐的涡量图

    Figure  4.  Vorticity picture of sonar with different shapes

    图  5  圆形声呐监测点的定义图

    Figure  5.  Definition diagram of monitoring points for circular sonar

    图  6  不同外形声呐的声压频谱曲线

    Figure  6.  Sound pressure picture of sonar with different shapes

    图  7  不同外形声呐声场外辐射指向特性

    Figure  7.  Directivity of radiation field of sonar with different shapes

    表  1  $Re = 3\;900$的圆柱绕流计算结果和参数

    Table  1.   Calculation results and parameters of flow around cylinder at $Re = 3\;900$

    算例$\Delta x/$mm$\Delta t/$s$\overline {{C_{\rm{D}}}} $$St$
    A1 0.06 0.000 1 1.10 0.229
    A2 0.000 2 1.08 0.227
    A3 0.000 5 1.08 0.225
    A4 0.1 0.000 1 1.09 0.234
    A5 0.000 2 1.08 0.237
    A6 0.000 5 1.07 0.237
    A7 0.3 0.000 1 1.09 0.232
    A8 0.000 2 1.08 0.239
    A9 0.000 5 1.05 0.236
    端木玉[19] 1.53 0.233
    詹昊[20] 1.26 0.220
    下载: 导出CSV

    表  2  各测点对应的总声压

    Table  2.   The total sound pressure corresponding to each measuring point

    测点总声压级/dB
    圆形方形椭圆
    1156.36136.94122.74
    2134.39129.52123.51
    3133.49124.15121.02
    4145.17146.98112.55
    5145.18146.99112.57
    695.7090.9944.75
    下载: 导出CSV
  • [1] 刘明星, 许欣然, 冯亮, 等. 高航速下艇艏声呐平台水动力噪声特性研究[C]//2019年全国声学大会论文集. 深圳: 中国声学学会, 2019: 117–118.

    LIU M X, XU X R, FENG L, et al. The characteristics research of hydrodynamic noise in bow sonar area at high speed[C]//Proceedings of the 2019 National Acoustics Conference. Shenzhen: Acoustic Society of China, 2019: 117–118 (in Chinese).
    [2] 张娟, 李天匀, 刘敬喜, 等. 空间不均匀湍流边界层激励下声呐腔自噪声统计能量分析[J]. 中国舰船研究, 2007, 2(6): 24–29.

    ZHANG J, LI T Y, LIU J X, et al. SEA analytical model of sonar's self-noise induced by turbulent boundary layer[J]. Chinese Journal of Ship Research, 2007, 2(6): 24–29 (in Chinese).
    [3] 张小锋, 刘国庆, 赵成, 等. 声类比水下圆柱绕流声学特性研究[J]. 舰船科学技术, 2019, 41(6): 120–124, 128.

    ZHANG X F, LIU G Q, ZHAO C, et al. Research on acoustic properties of flow around underwater cylinders based on acoustic analogue[J]. Ship Science and Technology, 2019, 41(6): 120–124, 128 (in Chinese).
    [4] 王建春, 吴乘胜, 王星, 等. 中等雷诺数方柱绕流的直接数值模拟及涡系分析[J]. 船舶力学, 2019, 23(8): 893–905.

    WANG J C, WU C S, WANG X, et al. Direct numerical simulation of medium Reynolds number flow past square cylinder and analysis of vortex structures[J]. Journal of Ship Mechanics, 2019, 23(8): 893–905 (in Chinese).
    [5] 彭华康, 李天匀, 朱翔. 湍流脉动压力下椭球声呐腔水动力自噪声分析[J]. 中国舰船研究, 2011, 6(4): 37–42. doi: 10.3969/j.issn.1673-3185.2011.04.008

    PENG H K, LI T Y, ZHU X. Self noise of hydrodynamic component for ellipsoid-shaped sonar cavity induced by turbulent fluctuating pressure[J]. Chinese Journal of Ship Research, 2011, 6(4): 37–42 (in Chinese). doi: 10.3969/j.issn.1673-3185.2011.04.008
    [6] KO S H, NUTTALL A H. Analytical evaluation of flush-mounted hydrophone array response to the corcos turbulent wall pressure spectrum[J]. The Journal of the Acoustical Society of America, 1991, 90(1): 579–588. doi: 10.1121/1.401232
    [7] KO S H. Performance of various shapes of hydrophones in the reduction of turbulent flow noise[J]. The Journal of the Acoustical Society of America, 1993, 93(3): 1293–1299. doi: 10.1121/1.405414
    [8] KO S H. Performance of shaded, rectangular hydrophones in the reduction of turbulent flow noisea)[J]. The Journal of the Acoustical Society of America, 1993, 94(3): 1681–1687. doi: 10.1121/1.408141
    [9] 王斌, 汤渭霖, 范军. 水听器非轴线布放时的拖线阵流噪声响应[J]. 声学学报, 2008(5): 402–408.

    WANG B, TANG W L, FAN J. Flow induced noise calculations for non-axially distributed hydrophones in towed arrays[J]. Acta Acustica, 2008(5): 402–408 (in Chinese).
    [10] 王晓林, 上官经邦, 王世闯, 等. 拖曳线列阵拖曳噪声抑制的试验研究[J]. 声学技术, 2017, 36(4): 123–126.

    WANG X L, SHANGGUAN J B, WANG S C, et al. Experimental research on towing noise suppression of towed line array[J]. Technical Acoustics, 2017, 36(4): 123–126 (in Chinese).
    [11] 刘明星, 许欣然, 夏铁坚. 潜艇水动力噪声对声呐声基阵影响分析[J]. 声学与电子工程, 2020(1): 18–22.

    LIU M X, XU X R, XIA T J. Analysis of the influence of submarine hydrodynamic noise on sonar acoustic array[J]. Acoustics and Electronics Engineering, 2020(1): 18–22 (in Chinese).
    [12] 龙双丽, 聂宏, 许鑫. 不同雷诺数下圆柱绕流气动噪声数值模拟[J]. 声学技术, 2011, 30(2): 111–116.

    LONG S L, NIE H, XU X. Numerical simulation of noise induced by flow around a cylinder at different Reynolds number[J]. Technical Acoustics, 2011, 30(2): 111–116 (in Chinese).
    [13] LYRINTZIS A S. Surface integral methods in computational aeroacoustics: from the (CFD) near-field to the (Acoustic) far-field[J]. International Journal of Aeroacoustics, 2003, 2(2): 95–128. doi: 10.1260/147547203322775498
    [14] 唐科范, FRANKE J. 用解NS方程和FW-H积分的混合方法计算圆柱绕流噪声[J]. 水动力学研究与进展, 2009, 24(2): 190–199.

    TANG K F, FRANKE J. Numerical simulation of noise induced by flow around cylinder using the hybrid method with the solutions of NS equation and FW-H integration[J]. Chinese Journal of Hydrodynamics, 2009, 24(2): 190–199 (in Chinese).
    [15] LUO K H, LAI H. A hybrid LES-acoustic analogy method for computational aeroacoustics[M]//LAMBALLAIS E, FRIEDRICH R, GEURTS B J, et al. Direct and Large-Eddy Simulation VI. Dordrecht: Springer Netherlands, 2006: 537-544.
    [16] 何鸿涛. 圆柱绕流及其控制的数值模拟研究[D]. 北京: 北京交通大学, 2009.

    HE H T. Numerical simulation to characteristics and control of flow around a circular cylinder[D]. Beijing: Beijing Jiaotong University, 2009 (in Chinese).
    [17] 郑力铭. ANSYS Fluent 15.0流体计算从入门到精通[M]. 北京: 电子工业出版社, 2015.

    ZHENG L M. ANSYS Fluent 15.0 fluid calculation from entry to proficiency[M]. Beijing: Publishing House of Electronics Industry, 2015 (in Chinese).
    [18] AI Y H, FENG D K, YE H K, et al. Unsteady numerical simulation of flow around 2-D circular cylinder for high Reynolds numbers[J]. Journal of Marine Science and Application, 2013, 12(2): 180–184. doi: 10.1007/s11804-013-1183-0
    [19] 端木玉, 万德成. 雷诺数为3900时三维圆柱绕流的大涡模拟[J]. 海洋工程, 2016, 34(6): 11–20.

    DUAN M Y, WAN D C. Large-eddy simulation of the flow past a cylinder with Re=3 900[J]. The Ocean Engineering, 2016, 34(6): 11–20 (in Chinese).
    [20] 詹昊, 李万平, 方秦汉, 等. 不同雷诺数下圆柱绕流仿真计算[J]. 武汉理工大学学报, 2008, 30(12): 129–132.

    ZHAN H, LI W P, FANG Q H, et al. Simulation calculation of flow around a cylinder under different Reynolds numbers[J]. Journal of Wuhan University of Technology, 2008, 30(12): 129–132 (in Chinese).
    [21] 林平. 钝体绕流数值模拟[D]. 杭州: 杭州电子科技大学, 2015.

    LIN P. The numerical simulation of turbulent flow over a blunt body[D]. Hangzhou: Hangzhou Dianzi University, 2015 (in Chinese).
    [22] 李玲, 刘沛清, 邢宇, 等. 亚临界雷诺数圆柱绕流远场气动噪声实验研究[J]. 北京航空航天大学学报, 2016, 42(5): 977–983.

    LI L, LIU P Q, XING Y, et al. Far-field aeroacoustic experimental study of flow around a circular cylinder at subcritical Reynolds number[J]. Journal of Beijing University of Aeronautics and Astronautics, 2016, 42(5): 977–983 (in Chinese).
  • [1] 柏铁朝, 许建, 王国栋, 余琨, 胡笑寒.  近冰面航行潜艇阻力及绕流场分析 . 中国舰船研究, doi: 10.19693/j.issn.1673-3185.01975
    [2] 邓超, 常宗瑜, 冯展霞, 张嘉坤.  波浪动力船水翼的推进性能分析 . 中国舰船研究, doi: 10.19693/j.issn.1673-3185.01672
    [3] 解学参, 常晟铭, 王超, 郭春雨.  开设内槽翼型的水动力性能数值模拟 . 中国舰船研究, doi: 10.19693/j.issn.1673-3185.01789
    [4] 张磊, 周志勇, 张雨新, 黄少锋.  基于数值模拟的螺旋桨吸气及水动力性能研究 . 中国舰船研究, doi: 10.19693/j.issn.1673-3185.01719
    [5] 李士强, 肖昌润, 曹植珺.  基于STAR-CCM+的潜艇尾流场及水动力数值分析 . 中国舰船研究, doi: 10.19693/j.issn.1673-3185.01216
    [6] 王晓聪, 桂洪斌, 刘洋.  三维有限长圆柱绕流数值模拟 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2018.02.004
    [7] 鲁利, 熊鹰, 王睿.  RANS, DES和LES对螺旋桨流噪声预报的适用性分析 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2017.06.007
    [8] 伍宏亮, 周其斗, 吕晓军, 孟庆昌.  基于湍流脉动压力的波数-频率谱预报流噪声 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2017.06.006
    [9] 于汉, 李清, 杨德庆.  水面舰船粘性流场和流噪声的数值计算 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2017.06.004
    [10] 王春旭, 吴崇建, 陈乐佳, 邱昌林, 熊济时.  流致噪声机理及预报方法研究综述 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2016.01.008
    [11] 李环, 刘聪尉, 吴方良, 陈灿.  水动力噪声计算方法综述 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2016.02.011
    [12] 王康, 刘国庆, 王强, 张咏鸥, 李旭.  水下航行器水动力噪声分离预报 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2015.04.005
    [13] 王冶, 徐筱欣.  船用换热器三维流场数值模拟 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2013.04.014
    [14] 方昭昭, 赵丙乾, 金武雷, 朱仁传.  顶浪中船舶水动力计算与运动模拟研究 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2013.05.006
    [15] 高菊, 陈美霞.  偶极子源与力激励作用下截顶锥形壳振动与声辐射研究 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2011.01.009
    [16] 卢锦国, 梁中刚, 吴方良, 周轶美.  水下航行体回转水动力数值计算研究 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2011.06.002
    [17] 詹金林, 卢晓平, 崔焰, 王家林.  基于势流理论的船舶操纵水动力分析 . 中国舰船研究, doi: 10.3969/j.issn.16733185.2011.05.009
    [18] 柏铁朝, 梁中刚, 周轶美, 敬军.  潜艇操纵性水动力数值计算中湍流模式的比较与运用 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2010.02.005
    [19] 刘明静, 马运义, 吴军.  潜艇艏舵绕流场的数值模拟 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2010.05.008
    [20] 孙荣, 吴晓光, 姜治芳, 吴启锐, 卢晓晖.  带自由面船体绕流场数值模拟 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2008.02.001
  • 加载中
图(7) / 表 (2)
计量
  • 文章访问数:  15
  • HTML全文浏览量:  5
  • PDF下载量:  2
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-07-30
  • 修回日期:  2020-12-07
  • 网络出版日期:  2021-03-23

不同外形声呐的水动力及流噪声数值模拟

doi: 10.19693/j.issn.1673-3185.02051
    基金项目:  中央高校基础研究基金资助项目(WUT:2019IVA059);博士科研启动基金资助项目(WUT:40120268)
    作者简介:

    邓云,男,1998年生,硕士生。研究方向:水动力及流噪声。E-mail:1569168520@qq.com

    陈威,男,1988年生,博士,研究方向:海洋工程水动力学。E-mail:whutcw01@126.com

    通讯作者: 陈威
  • 中图分类号: U661.44; U666.7

摘要:   目的  为降低传统圆形声呐外表面流噪声,设计了3种不同的声呐外形(圆形、椭圆形以及方形),并对其外表面水动力和流噪声进行了研究。  方法  基于Fluent软件中的标准$k - \varepsilon$湍流计算模型和Lighthill声类比方法对3种不同外形声呐外表面的流场和声场进行分析。  结果  结果显示:方形的声呐升力系数幅值和阻力系数均值最大,其次是圆形,最小的是椭圆形;各方案声呐水动力差异是由于边界层分离点及尾涡的不同引起的。对比流噪声结果发现,圆形的总声压级最大,椭圆形的最小。圆形和方形呈现了“正8字”的偶极子声源特性,声辐射最大值在垂直于来流方向;而椭圆呈现“倒8字”的偶极子声源特性,最大值在水平来流方向。  结论  可用椭圆形形声呐替代传统的圆形声呐,结果作为声呐外形设计的参考。

English Abstract

邓云, 傅何琪, 广超越, 等. 不同外形声呐的水动力及流噪声数值模拟[J]. 中国舰船研究, 2021, 16(4): 1–7 doi: 10.19693/j.issn.1673-3185.02051
引用本文: 邓云, 傅何琪, 广超越, 等. 不同外形声呐的水动力及流噪声数值模拟[J]. 中国舰船研究, 2021, 16(4): 1–7 doi: 10.19693/j.issn.1673-3185.02051
DENG Y, FU H Q, GUANG C Y, et al. Numerical simulation on hydrodynamic and flow noise of different shape sonars[J]. Chinese Journal of Ship Research, 2021, 16(4): 1–7 doi: 10.19693/j.issn.1673-3185.02051
Citation: DENG Y, FU H Q, GUANG C Y, et al. Numerical simulation on hydrodynamic and flow noise of different shape sonars[J]. Chinese Journal of Ship Research, 2021, 16(4): 1–7 doi: 10.19693/j.issn.1673-3185.02051
    • 拖曳式线列阵声呐的主要作用是利用声波在海洋中传播能量衰减的能力弱且能远距离传送等特点,对目标进行定位、警戒等。拖曳式声呐由水听器和护套2个部分组成,当声呐工作时,声呐与海水间的相对运动会产生湍流边界层,引起压力起伏,产生流噪声。近年来,拖曳式列阵声呐的流噪声问题受到了研究者们的广泛关注[1-2]

      研究发现,声呐外形很大程度上会影响水动力及流噪声性能,故学者们从2个方面展开了一些研究。在水动力数值研究方面,张晓峰等[3]以三维刚性圆柱为研究对象,探讨了不同雷诺数、间距比和排列方式下单、双圆柱的水动力特性,研究结果为圆形声呐流场特性的精确预报提供了一定的指导。王建春等[4]分析方柱绕流的复杂涡系结构,为研究方形声呐外辐射奠定了理论基础。彭华康等[5]运用统计能量法建立了水下椭球壳外形的声呐腔模型,考虑在椭球声呐湍流脉动压力空间分布不均匀的条件下计算得到不同航速时声呐罩表面湍流边界层的分离点,以为水下椭圆形声呐自噪声预报和控制提供参考。

      在流噪声数值研究方面,Sung等[6-8]分析了声呐内部不同外形的水听器对降噪的影响,发现矩形和圆形的降噪效果优于菱形和三角形。汤渭霖等[9]研究了流噪声产生的机理,指出流噪声的强度大小与水听器结构设计存在着很大的关系。王晓林等[10]从扩大水听器和护套表面距离以及结构减振优化的角度提出了拖曳声呐流噪声的抑制方法。刘明星等[11]对不同声呐基阵类型的水动力噪声进行仿真分析,得到了球阵水动力噪声要小于圆柱阵水动力噪声的结果。由于目前的研究一般是通过设计声呐内部的水听器外形来降低流噪声在声呐内部的耦合效应,所以在降低声呐护套外部的湍流边界层压力起伏引起的流噪声研究方面则显得不足。因此,有必要从抑制流噪声的角度对声呐护套外形进行合理的声学设计。

      然而,开展声呐流噪声的相关研究需从不同声呐外形的角度对其流场和声场进行分析,以得到不同外形声呐的声场和流场分布规律,这不仅在声呐结构设计方面具有重要的研究价值,在提高声呐工作性能方面也具有重要的实用价值。

      鉴此,为降低传统圆形声呐外表面流噪声,本文设计了3种不同外形(圆形、方形、椭圆形)的声呐,拟通过对比声呐外表面的流场和声场分布规律,降低传统的圆形声呐外表面流噪声,合理优化声呐的设计。

    • 研究过程中,描述流场的无量纲数主要包括雷诺数$Re$、阻力系数${C_{\rm{D}}}$、升力系数${C_{\rm{L}}}$、斯特劳哈尔数St。斯特劳哈尔数等于当地惯性力与迁移惯性力之比,而雷诺数等于惯性力与黏性力之比。具体定义式如下:

      $$ St = \frac{{f_{\rm{v}}D}}{U} $$ (1)
      $$ Re = \frac{{UD}}{\upsilon } $$ (2)

      式中,$f_{\rm{v}}$为涡脱落频率;$D$为圆柱直径;$U$为来流速度;$\upsilon $为流体动力黏度。

      升力系数CL和阻力系数CD定义式如下:

      $$ {C_{\rm{D}}} = \frac{{{F_{\rm{D}}}}}{{0.5\rho \mathop U\nolimits^2 A}} $$ (3)
      $$ {C_{\rm{L}}} = \frac{{{F_{\rm{L}}}}}{{0.5\rho \mathop U\nolimits^2 A}} $$ (4)

      式中,${F_{\rm{D}}}$为圆柱所受到的阻力;${F_{\rm{L}}}$为圆柱受到的升力;$\rho $为流体密度;$A$为圆柱的迎流面积,在二维情况下为圆柱的直径。

    • 在进行近场流场计算后,取声呐周围为声源,得到近场强度大小和分布规律,并将声场相关的量传至远场。在整个传递过程中忽略流动和波动对声场的影响,龙双丽等[12]使用基于FW-H方程[13-14]的积分外推方法计算出了外流场的声学信息,如式(5)所示。FW-H方程的右边3项代表声辐射源,其中,第1项代表流体本身的湍流应力,具有四极子特性;第2项代表施加在某些界面上非稳定力的散度,具有偶极子特性;第3项代表进入到流体中的非稳定质量流,具有单极子特性[15]。由于声场计算过程中时间具有滞后性,会造成声场结果的不可靠性,所以,为得到准确声场信息,应在稳定的流场中进行声场计算。

      $$ \begin{split} & \frac{1}{{{a_0}^2 }}\frac{{\partial^2 p' }}{{\partial^2 t}} - \nabla^2 p' = \frac{{\partial^2 }}{{\partial x_i\partial x_j}}\{ T_{ij}H(f)\} - \\[-3pt]&\;\;\;\; \frac{\partial }{{\partial x_i}}\{ [P_{ij}{{\mathit{\boldsymbol{n}}}}_j + \rho u_i(u_{\rm{n}} - v_{\rm{n}})\omega (f)\} + \\[-3pt]&\qquad \frac{\partial }{{\partial t}}\{ [\rho_ 0v_{\rm{n}} + \rho u_i(u_{\rm{n}} - v_{\rm{n}})\omega (f)\} \end{split} $$ (5)

      式中,${{\mathit{\boldsymbol{n}}}}_j $为控制面上的单位外法向矢量;$u_i$$x_i$方向的速度分量;$u_{\rm{n}}$为流体在$f = 0$面上的法向速度;$v_{\rm{n}}$物面速度方向的法向分量;$\rho_ 0 $为流体密度参考值;$\omega (f)$为Dirichlet 函数;$H(f)$为Heaviside广义函数;$\nabla^2 $为拉普拉斯算子;$\mathop p\nolimits' $为远场声压;$a_{0}$为远场的声速;$T_{ij}$为Lighthill应力张量;$P_{ij}$为可压缩流体应力张量。

    • 本文主要研究二维情况下不同外形的声呐表面水动力及流噪声的分布规律。建模过程从几何建模、设置物理参数、边界条件以及求解方式等几个方面来进行。圆形声呐建模过程是以坐标原点为圆心,圆心距离上、下边界分别为$10D$,距离左、右边界分别为$10D$$25D$,建模过程取圆柱直径$D = 10$ mm,方形声呐建模将其中的圆形变为方形即可。椭圆的流场特性与椭圆的长、短轴的比值有关,为了研究方便,取长、短轴比值为2∶1,即$a = 2D$$b = D$图1所示为计算域。

      图  1  (a) 不同外形的声呐计算域图

      Figure 1.  Calculation domain of sonar with different shapes

      模拟中,流体介质为水,其密度为998.2 $\mathop {{\rm{kg}}\cdot{\rm{m}}}\nolimits^{ - 3} $,运动黏度为$\mathop {1.004 \times 10}\nolimits^{ - 6}\; \mathop {\rm{m}}\nolimits^2\cdot\mathop {\rm{s}}\nolimits^{ - 2} $,本文雷诺数取为3 900,该雷诺数下流动状态为湍流[16],湍流模型选择标准$k - \varepsilon $模型。相较于其他湍流模型,湍流$k - \varepsilon $模型具有稳定性、经济性和比较高的计算精度等优点[17]。左侧边界为速度入口,右侧边界为出口边界,上、下边界为对称边界,内部圆柱为无滑移静止壁面。为加快求解的收敛,本文选择SIMPLEC压力速度耦合方式,该求解方式相比SIMPLE,在四边形网格上具有更好的收敛性。

    • 为验证计算结果的合理性,以圆形声呐流场模型为例,对不同网格密度和时间步长的组合进行验证[16,18],计算参数和结果如表1所示。表中,$\Delta x$为圆柱周围最小的网格尺寸,$\Delta t$为时间步长,$\overline {{C_{\rm{D}}}}$为平均阻力系数。

      表 1  $Re = 3\;900$的圆柱绕流计算结果和参数

      Table 1.  Calculation results and parameters of flow around cylinder at $Re = 3\;900$

      算例$\Delta x/$mm$\Delta t/$s$\overline {{C_{\rm{D}}}} $$St$
      A1 0.06 0.000 1 1.10 0.229
      A2 0.000 2 1.08 0.227
      A3 0.000 5 1.08 0.225
      A4 0.1 0.000 1 1.09 0.234
      A5 0.000 2 1.08 0.237
      A6 0.000 5 1.07 0.237
      A7 0.3 0.000 1 1.09 0.232
      A8 0.000 2 1.08 0.239
      A9 0.000 5 1.05 0.236
      端木玉[19] 1.53 0.233
      詹昊[20] 1.26 0.220

      表1可见,3种网格计算出的$\overline {{C_{\rm{D}}}} $均能够较好地满足数值计算结果的要求。从计算效率的角度而言,网格尺寸过小会造成计算效率偏低,故计算时采用了相对较大的0.3 mm网格尺寸。在确定网格尺寸之后,根据$Courant$数选择合适的时间步长。$Courant$数定义如下:

      $$ Courant = \frac{{U.\Delta t}}{{\Delta x}} $$ (6)

      $Courant$数表示的是在一个时间步长内一个流体质点可以穿过的网格数。若选择的时间步长合理,则$Conrant$数需小于1,如此才能确保计算的精度。若选择了中等长度的时间步长0.000 2 s,$Courant$数为0.26。此时,结果与文献[19-20]比较吻合,这表明结果及建立的模型是有效的。因此,所有模型均选择A8算例进行水动力及流噪声的分析。

    • 分析3种不同外形的声呐阻力、升力系数,结果如图2所示。随着时间的推移,圆形声呐${C_{\rm{L}}}$振幅逐渐增加,${C_{\rm{D}}}$均值先增后减,然后缓慢增加到稳定状态,最后呈现出${C_{\rm{L}}}$幅值在0.85、${C_{\rm{D}}}$均值在1.08的周期性单一频率振动。而方形的水动力系数较大,${C_{\rm{L}}}$幅值为1.31,${C_{\rm{D}}}$均值为1.80。与圆形和方形外形的声呐相比,椭圆形的声呐水动力系数都较小,${C_{\rm{L}}}$幅值为0.05,${C_{\rm{D}}}$均值为0.43。

      图  2  不同外形声呐升阻力系数

      Figure 2.  Lift and drag coefficient of sonar with different shapes

      图3所示为3种外形声呐的升力频谱特性曲线。其中,圆形的升力频率(涡脱落频率)为9.26 Hz,其斯特劳哈尔数为0.236;方形的升力频率(涡脱落频率)为5.82 Hz,其斯特劳哈尔数为0.148;椭圆形的升力频率(涡脱落频率)为10.80 Hz,其斯特劳哈尔数为0.275。通过比较发现,方形声呐的斯特劳哈尔数最小,椭圆的斯特劳哈尔数最大。说明方形涡旋脱落周期较长,椭圆的涡旋脱落周期较短。

      图  3  不同外形声呐频谱

      Figure 3.  Spectrum diagram of sonar with different shapes

      图4所示为3种外形声呐的涡量图。声呐周围表面的涡旋脱落会引起表面压强的周期性变化,使升力和阻力周期性变化,且因正、负涡交替脱落的涡强度相等,造成升力系数均值为0。水平端前、后两侧受到的压力方向一致,造成阻力系数均值不为0[21],由图4可见,在相同计算域内,圆形与椭圆外形的声呐存在4对脱落涡旋,方形存在2对脱落涡旋,这说明在圆形与椭圆外形下涡脱落频率相差不大,但与方形涡脱落频率相比,相差较大,约为方形涡脱落频率的2倍。可见,涡脱的差异造成了升力频率不一。 另外,方形边界层分离位置在方形的前端上、下两个角点位置,而圆形边界层分离的位置则在中间且稍偏后的位置,椭圆形的分离位置相比圆形,位置更靠后一些,且边界层拉长造成的尾涡长度大于圆形。可见,分离点的差异造成了压差阻力的不同。因此,相比较而言,方形声呐的阻力系数最大,椭圆形的最小。

      图  4  不同外形声呐的涡量图

      Figure 4.  Vorticity picture of sonar with different shapes

    • 在圆形外形的声呐周围定义6个测点,位置分别为1(0,0),2(D,0),3(−D,0),4(0,D),5(0,−D),6(10D, 0)。其中,2和3号测点的位置是关于垂直轴对称,4和5号测点是关于水平轴对称。测点的具体位置如图5所示。方形和椭圆外形的声呐周围测点的位置与圆形的相同。

      图  5  圆形声呐监测点的定义图

      Figure 5.  Definition diagram of monitoring points for circular sonar

      图6所示为3种不同外形声呐的声压频谱曲线。各测点声压最大值所对应的频率为涡脱落频率,且各声压峰值处的频率均为涡脱落频率的整数倍。测点6距离声呐距离最远,其声压值最小。如表2所示,3种声呐的总声压级中内部的总声压级最大,距离最远的测点的最小,距离中心等距的测点(例如,测点2,3,4和5)总声压级相差不大。对比3种不同外形声呐,椭圆的总声压级较小。

      图  6  不同外形声呐的声压频谱曲线

      Figure 6.  Sound pressure picture of sonar with different shapes

      表 2  各测点对应的总声压

      Table 2.  The total sound pressure corresponding to each measuring point

      测点总声压级/dB
      圆形方形椭圆
      1156.36136.94122.74
      2134.39129.52123.51
      3133.49124.15121.02
      4145.17146.98112.55
      5145.18146.99112.57
      695.7090.9944.75

      图7所示为3种外形的声呐外辐射指向特性图。由图可见,各外形声呐外辐射指向都呈现出一个“8”字形,说明流噪声的主要来源为表面交替涡脱落产生的非定常脉动力形成的偶极子声源。其中,圆形和方形声呐的外辐射特性规律类似,都呈现出“正8字”形的偶极子声源特性,与李玲等[22]给出的结果相似,即声辐射横向最大、水平来流方向最小。这是因为横向声辐射主要源于涡脱落产生的非定常脉动力,而水平来流方向主要源于涡脱落中湍流应力,从图3中也可以看出,圆形和方形声呐的横向的升力幅值较大。此外,从表2也可以发现,等距的横向测点4和5的总声压级要大于水平测点2和3的总声压级。与圆形和方形声呐不同,椭圆形声呐外辐射特性呈现出“倒8字”形,声场外辐射的最大值发生在水平来流方向,而与来流方向成70°,外辐射值最小,对应的升力幅值也较小(见图2),且等距的水平测点2和3的总声压级要大于横向测点4和5的总声压级(见表2)。

      图  7  不同外形声呐声场外辐射指向特性

      Figure 7.  Directivity of radiation field of sonar with different shapes

    • 本文采用$k - \varepsilon $湍流模型与声类比方法对3种不同外形的声呐水动力及流噪声进行了研究,得到如下结论:

      1)水动力方面,方形边界层分离位置靠前,圆形的为中部偏后,椭圆的比圆形的更靠后一些。边界层分离位置以及尾涡的差异造成方形阻力系数均值和升力系数幅值最大,椭圆最小。

      2)流噪声方面,柱体的声场具有一定的对称性;距离声源越近,总声压级越大;相较于圆形和方形,椭圆总声压级最小;不同外形声呐外辐射都呈现“8”字形,但指向方位有所差异。圆形和方形声呐的声辐射最大值在横向,椭圆最大值在水平来流方向。

      本文的研究成果在一定程度上有助于声呐外形的声学设计。

参考文献 (22)

目录

    /

    返回文章
    返回