搜索

x

留言板

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

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

孤子内波环境下三维声传播建模

张泽众 骆文于 庞哲 周益清

引用本文:
Citation:

孤子内波环境下三维声传播建模

张泽众, 骆文于, 庞哲, 周益清

Modeling of three-dimensional sound propagation through solitary internal waves

Zhang Ze-Zhong, Luo Wen-Yu, Pang Zhe, Zhou Yi-Qing
PDF
HTML
导出引用
  • 提出了一种适用于存在孤子内波水平变化波导的高效三维水下声场计算模型. 该模型忽略反向散射, 一般情况下由于孤子内波的反向散射非常弱, 所以该模型能够提供精确的三维声场结果. 同时, 相对于双向三维耦合简正波模型, 该模型在计算效率上能够至少提高一个数量级. 除了孤子内波环境之外, 本模型还适用于存在小尺度海脊等反向散射比较弱的一般水平变化波导环境. 本文用该模型计算由KdV方程得到的孤子内波问题, 并用双向三维耦合简正波模型作为标准模型来验证本模型的计算精度. 计算结果表明本模型在反向散射比较弱的波导环境中具有非常高的计算精度.
    An accurate and numerically efficient numerical model is very important for studying the effect of internal wave on underwawter sound propagation. A full-wave, three-dimensional (3D) coupled-mode model is able to deal with the internal wave problem with satisfactory accuracy, but such a model is in general numerically inefficient. A numerically efficient 3D model is presented for sound propagation in a range-dependent waveguide in the presence of solitary internal waves in this work. The present model is a forward-marching model that neglects backscattering. In this 3D model, an efficient two-dimensional (2D) coupled-mode model, C-SNAP, is adopted to compute 2D acoustic field solutions excited by a line source. The C-SNAP is a 2D forward-marching model, which uses an energy-conserving matching condition to preserve accuracy. An appealing aspect of C-SNAP is that its efficiency is competitive with that of the existing parabolic equation model. The integral transform technique is used to extend C-SNAP to a 3D model, where a complex integration contour is used for evaluating the wavenumber integral. A brief review of C-SNAP and formulation of the present 3D model are given. The forward-marching models are primarily suitable for treating the range-dependent problems with weak backscattering, such as with a slowly varying bathymetry. Since in general the backscattering from internal wave is weak, which is also validated numerically in this work, the present model is able to address the problem of sound propagation through internal wave with satisfactory accuracy. At the same time, it achieves an efficiency gain of at least an order of magnitude over that of full two-way, 3D model. In addition to the internal wave, the present model is also suitable for solving the general range-dependent problems where backscattering is weak, such as in the presence of a bottom ridge of a small height. Numerical simulations are also provided to validate the present model, where a two-way, 3D model serves as the benchmark. The numerical results show that the effect of the internal wave on the acoustic field is negligible for the region between the source and the internal wave. However, the effect is significant on the other side of the internal wave. A more interesting observation is the angular dependence of the interference pattern induced by the internal wave.
      通信作者: 骆文于, lwy@mail.ioa.ac.cn
    • 基金项目: 国家自然科学基金(批准号: 11774374)资助的课题
      Corresponding author: Luo Wen-Yu, lwy@mail.ioa.ac.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 11774374)
    [1]

    王宁, 张海青, 王好忠, 高大治 2010 声学学报 35 38

    Wang N, Zhang H Q, Wang H Z, Gao D Z 2010 Acta Acust. 35 38

    [2]

    马树青, 杨士莪, 朴胜春, 李婷婷 2009 振动与冲击 28 73Google Scholar

    Ma S Q, Yang S E, Piao S C, Li T T 2009 J. Vib. Shock 28 73Google Scholar

    [3]

    Badiey M, Katsnelson B G, Lynch J F, Pereselkov S, Siegmann W L 2005 J. Acoust. Soc. Am. 117 613Google Scholar

    [4]

    Luo J, Badiey M, Karjadi E A, Katsnelson B, Tskhoidze A, Lynch J F, Moum J N 2008 J. Acoust. Soc. Am. 124 66Google Scholar

    [5]

    Lin Y T, Duda T F, Lynch J F 2009 J. Acoust. Soc. Am. 126 1752Google Scholar

    [6]

    Colosi J A 2008 J. Acoust. Soc. Am. 124 1452Google Scholar

    [7]

    Yang T C 2014 J. Acoust. Soc. Am. 135 610Google Scholar

    [8]

    Siegmann W L, Kriegsmann G A, Lee D 1985 J. Acoust. Soc. Am. 78 659Google Scholar

    [9]

    Lee D, Botseas G, Siegmann W L 1992 J. Acoust. Soc. Am. 91 3192Google Scholar

    [10]

    Collins M D, Chin-Bing S A 1990 J. Acoust. Soc. Am. 87 1104Google Scholar

    [11]

    Ferla C M, Porter M B, Jensen F B 1993 C-SNAP: Coupled SACLANTCEN Normal Mode Propagation Loss Model (La Spezia: SACLANT Undersea Research Center) pp1–46

    [12]

    Luo W Y, Yang C M, Qin J X, Zhang R H 2012 Sci. Chin. Phys. Mech. Astron. 55 572Google Scholar

    [13]

    杨春梅, 骆文于, 张仁和, 秦继兴 2014 声学学报 39 295

    Yang C M, Luo W Y, Zhang R H, Qin J X 2014 Acta Acust. 39 295

    [14]

    Luo W Y, Zhang R H 2015 Sci. China-G: Phys. Mech. Astron. 58 1

    [15]

    Evans R B 1983 J. Acoust. Soc. Am. 74 188Google Scholar

    [16]

    Collins M D, Westwood E K 1991 J. Acoust. Soc. Am. 89 1068Google Scholar

    [17]

    Jensen F B, Kuperman W A, Porter M B, Schmidt H 2011 Computational Ocean Acoustics (2nd ed.) (New York: Springer) pp337–445

    [18]

    Buckingham M J 1989 J. Acoust. Soc. Am. 86 2273Google Scholar

    [19]

    Luo W Y, Yu X L, Yang X F, Zhang Z Z, Zhang R H 2016 Chin. Phys. B 25 124309Google Scholar

    [20]

    Jensen F B 1998 J. Acoust. Soc. Am. 104 1310Google Scholar

  • 图 1  水平变化波导(红色虚线)的阶梯近似

    Fig. 1.  Stair step approximation of a sloping bottom (red dashed line)

    图 2  反傅里叶变换的积分围线(圆圈代表简正波的本征值, 即水平波数)

    Fig. 2.  Complex integration contour for evaluation of the wavenumber integral. The dots indicate horizontal wavenumbers of the normal modes.

    图 3  二维声速剖面示意图 (a)无孤子内波的声速剖面; (b)孤子内波位于4 km处的声速剖面; (c)孤子内波位于4 km处, 海脊位于6 km处的声速剖面

    Fig. 3.  Sound speed fields considered in this paper: (a) The background sound speed field; (b) in the presence of an internal wave soliton, centered at range 4 km from the source; (c) in the presence of both an internal wave soliton centered at 4 km and a cosine-bell shaped bottom ridge centered at 6 km from the source.

    图 4  孤子内波环境下DGMCM2D和C-SNAP在70 m深度传播损失曲线 蓝色实线和红色虚线分别表示无孤子内波时DGMCM2D和C-SNAP的计算结果, 绿色和枚红色虚线分别表示孤子内波波包在4 km时DGMCM2D和C-SNAP的计算结果, 黑色虚线表示孤子内波波包位置

    Fig. 4.  Transmission loss results for the internal solitary wave problem computed by DGMCM2D and C-SNAP. The blue solid line and red dashed line are the results by DGMCM2D and C-SNAP for the case without internal waves, respectively, and the green and magenta dashed lines are the results by DGMCM2D and C-SNAP for the case with a soliton located at range 4 km, respectively. The black dashed line indicates the center of the soliton.

    图 5  70 m深度水平面上的三维传播损失结果 (a) DGMCM3D的结果; (b)本文提出的三维模型的结果, 黑色虚线代表孤子内波波包的位置

    Fig. 5.  Three-dimensional transmission loss results in the horizontal plane at depth 70 m computed by (a) DGMCM3D and (b) the present 3D model. The center location of the internal wave is indicated by dashed black lines.

    图 6  纵向距离$ y = 0 \;{\rm{km}}$, 深度z = 70 m, 随距离传播的损失曲线, 蓝色实线和红色虚线分别为DGMCM3D和本文提出的三维模型的结果

    Fig. 6.  Transmission loss lines versus range along the cross-range $ y = 0 \;{\rm{km}}$ at depth z = 70 m computed by DGMCM3D (the blue, solid curve) and the present 3D model (the red, dashed curve).

    图 7  海底存在海脊的水平变化波导示意图

    Fig. 7.  Geometry of a range-dependent waveguide with a bottom ridge.

    图 8  孤子内波和海脊同时存在时深度70 m处的二维传播损失结果, 其中蓝色实线为DGMCM2D结果, 红色虚线为C-SNAP结果

    Fig. 8.  Two-dimensional transmission loss results at a depth of 70 m for the problem involving a solitary internal wave as well as a bottom ridge by DGMCM2D (the blue, solid curve) and C-SNAP (the red, dashed curve).

    图 9  同时存在孤子内波和海脊时深度70 m水平平面上的三维传播损失 (a) DGMCM3D结果; (b)本文提出的三维模型的结果. 黑色和红色虚线分别表示孤子内波波包和海脊中心的位置

    Fig. 9.  Three-dimensional transmission loss results in the horizontal plane at depth of 70 m in the presence of a solitary internal wave as well as a bottom ridge computed by (a) DGMCM3D and (b) the present 3D model. The center locations of the internal wave and the ridge are indicated by dashed black and red lines, respectively.

    Baidu
  • [1]

    王宁, 张海青, 王好忠, 高大治 2010 声学学报 35 38

    Wang N, Zhang H Q, Wang H Z, Gao D Z 2010 Acta Acust. 35 38

    [2]

    马树青, 杨士莪, 朴胜春, 李婷婷 2009 振动与冲击 28 73Google Scholar

    Ma S Q, Yang S E, Piao S C, Li T T 2009 J. Vib. Shock 28 73Google Scholar

    [3]

    Badiey M, Katsnelson B G, Lynch J F, Pereselkov S, Siegmann W L 2005 J. Acoust. Soc. Am. 117 613Google Scholar

    [4]

    Luo J, Badiey M, Karjadi E A, Katsnelson B, Tskhoidze A, Lynch J F, Moum J N 2008 J. Acoust. Soc. Am. 124 66Google Scholar

    [5]

    Lin Y T, Duda T F, Lynch J F 2009 J. Acoust. Soc. Am. 126 1752Google Scholar

    [6]

    Colosi J A 2008 J. Acoust. Soc. Am. 124 1452Google Scholar

    [7]

    Yang T C 2014 J. Acoust. Soc. Am. 135 610Google Scholar

    [8]

    Siegmann W L, Kriegsmann G A, Lee D 1985 J. Acoust. Soc. Am. 78 659Google Scholar

    [9]

    Lee D, Botseas G, Siegmann W L 1992 J. Acoust. Soc. Am. 91 3192Google Scholar

    [10]

    Collins M D, Chin-Bing S A 1990 J. Acoust. Soc. Am. 87 1104Google Scholar

    [11]

    Ferla C M, Porter M B, Jensen F B 1993 C-SNAP: Coupled SACLANTCEN Normal Mode Propagation Loss Model (La Spezia: SACLANT Undersea Research Center) pp1–46

    [12]

    Luo W Y, Yang C M, Qin J X, Zhang R H 2012 Sci. Chin. Phys. Mech. Astron. 55 572Google Scholar

    [13]

    杨春梅, 骆文于, 张仁和, 秦继兴 2014 声学学报 39 295

    Yang C M, Luo W Y, Zhang R H, Qin J X 2014 Acta Acust. 39 295

    [14]

    Luo W Y, Zhang R H 2015 Sci. China-G: Phys. Mech. Astron. 58 1

    [15]

    Evans R B 1983 J. Acoust. Soc. Am. 74 188Google Scholar

    [16]

    Collins M D, Westwood E K 1991 J. Acoust. Soc. Am. 89 1068Google Scholar

    [17]

    Jensen F B, Kuperman W A, Porter M B, Schmidt H 2011 Computational Ocean Acoustics (2nd ed.) (New York: Springer) pp337–445

    [18]

    Buckingham M J 1989 J. Acoust. Soc. Am. 86 2273Google Scholar

    [19]

    Luo W Y, Yu X L, Yang X F, Zhang Z Z, Zhang R H 2016 Chin. Phys. B 25 124309Google Scholar

    [20]

    Jensen F B 1998 J. Acoust. Soc. Am. 104 1310Google Scholar

  • [1] 张沐安, 王进卿, 吴睿, 冯致, 詹明秀, 徐旭, 池作和. 多孔介质内气泡Ostwald熟化特性三维孔网数值模拟.  , 2023, 72(16): 164701. doi: 10.7498/aps.72.20230695
    [2] 罗勇, 杨党国, 武从海, 李虎, 张树海, 吴军强. 三维空腔流动波系建模及模态演化.  , 2022, 71(19): 194301. doi: 10.7498/aps.71.20220922
    [3] 李森清, 张肖, 林机. 熔融耦合器中耦合模式与新型孤子结构.  , 2022, 71(23): 234207. doi: 10.7498/aps.71.20221273
    [4] 曹明鹏, 吴晓鹏, 管宏山, 单光宝, 周斌, 杨力宏, 杨银堂. 基于对偶单元法的三维集成微系统电热耦合分析.  , 2021, 70(7): 074401. doi: 10.7498/aps.70.20201628
    [5] 原青云, 孙永卫, 张希军. 基于电荷守恒定律的航天器内带电三维仿真简化模型.  , 2019, 68(19): 195201. doi: 10.7498/aps.68.20190631
    [6] 杨雄, 程谋森, 王墨戈, 李小康. 螺旋波等离子体放电三维直接数值模拟.  , 2017, 66(2): 025201. doi: 10.7498/aps.66.025201
    [7] 陆昌根, 朱晓清, 沈露予. 三维边界层内诱导横流失稳模态的感受性机理.  , 2017, 66(20): 204702. doi: 10.7498/aps.66.204702
    [8] 秦继兴, Katsnelson Boris, 彭朝晖, 李整林, 张仁和, 骆文于. 三维绝热简正波-抛物方程理论及应用.  , 2016, 65(3): 034301. doi: 10.7498/aps.65.034301
    [9] 陈茂林, 夏广庆, 毛根旺. 多模式离子推力器栅极系统三维粒子模拟仿真.  , 2014, 63(18): 182901. doi: 10.7498/aps.63.182901
    [10] 宋天明, 杨家敏. 三维柱腔内辐射输运的一维模拟.  , 2013, 62(1): 015210. doi: 10.7498/aps.62.015210
    [11] 张兰强, 顾乃庭, 饶长辉. 大气湍流三维波前探测模式层析算法分析.  , 2013, 62(16): 169501. doi: 10.7498/aps.62.169501
    [12] 王新鑫, 樊丁, 黄健康, 黄勇. 双钨极耦合电弧数值模拟.  , 2013, 62(22): 228101. doi: 10.7498/aps.62.228101
    [13] 刘磊, 费建芳, 黄小刚, 程小平. 大气-海浪-海流耦合模式的建立和一次台风过程的初步试验.  , 2012, 61(14): 149201. doi: 10.7498/aps.61.149201
    [14] 张霞萍. 强非局域空间三维光孤子短程相互作用.  , 2011, 60(3): 034211. doi: 10.7498/aps.60.034211
    [15] 李鹏, 赵建林, 张晓娟, 侯建平. 三角结构三芯光子晶体光纤中的模式耦合特性分析.  , 2010, 59(12): 8625-8631. doi: 10.7498/aps.59.8625
    [16] 郝保良, 肖刘, 刘濮鲲, 李国超, 姜勇, 易红霞, 周伟. 螺旋线行波管三维频域非线性注波互作用的计算.  , 2009, 58(5): 3118-3124. doi: 10.7498/aps.58.3118
    [17] 蔡鑫伦, 黄德修, 张新亮. 基于全矢量模式匹配法的三维弯曲波导本征模式计算.  , 2007, 56(4): 2268-2274. doi: 10.7498/aps.56.2268
    [18] 董慧媛, 刘 楣, 吴宗汉, 汪 静, 王振林. 由介质球构成的三维光子晶体能带结构的平面波研究.  , 2005, 54(7): 3194-3199. doi: 10.7498/aps.54.3194
    [19] 刘敬伟, 陈少武, 余金中. 一种分析三维楔脊形光波导与光纤耦合的方法.  , 2005, 54(1): 6-11. doi: 10.7498/aps.54.6
    [20] 朱学光, 匡光力, 谢纪康. 三维洞穴环天线的耦合计算.  , 2000, 49(10): 1978-1981. doi: 10.7498/aps.49.1978
计量
  • 文章访问数:  8486
  • PDF下载量:  93
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-04-02
  • 修回日期:  2019-08-20
  • 上网日期:  2019-10-01
  • 刊出日期:  2019-10-20

/

返回文章
返回
Baidu
map