关于非线性水声波的近似色散关系分析.pdf

International Journal of Fluid Dynamics 流体动力学流体动力学,2023,11(3),94-102 Published Online September 2023 in Hans.https:/www.hanspub.org/journal/ijfd https:/doi.org/10.12677/ijfd.2023.113009 文章引用文章引用:付鲁丹,王苹.关于非线性水声波的近似色散关系分析J.流体动力学,2023,11(3):94-102.DOI:10.12677/ijfd.2023.113009 关于非线性水声波的近似色散关系分析关于非线性水声波的近似色散关系分析 付鲁丹付鲁丹1,2，王，王 苹苹1,2*1青岛科技大学数理学院，山东 青岛 2青岛科技大学数学与跨学科科学研究所，山东 青岛 收稿日期：2023年8月21日；录用日期：2023年8月28日；发布日期：2023年9月28日 摘摘 要要 目前在非线性水声波的研究中，为了简化计算通常假设色散关系是线性的。然而线性色散关系式对大振目前在非线性水声波的研究中，为了简化计算通常假设色散关系是线性的。然而线性色散关系式对大振幅波是不恰当的。为此，幅波是不恰当的。为此，本文研究了弹性板覆盖下的有限深可压单层流体中非线性水声波的色散关系。本文研究了弹性板覆盖下的有限深可压单层流体中非线性水声波的色散关系。假设流体是无粘、可压缩的，假设流体是无粘、可压缩的，且且流动是无旋的流动是无旋的。我们构建了控制方程和表示水动力、弹性力和惯性力之我们构建了控制方程和表示水动力、弹性力和惯性力之间关系的边界条件，并间关系的边界条件，并得到得到水声波水声波的的一种近似一种近似非线性色散关系非线性色散关系，对水声波模态的特性进行分析对水声波模态的特性进行分析，讨论了讨论了弹性板厚度等弹性板厚度等重要物理因素重要物理因素对波动传播特性的影响对波动传播特性的影响。该研究该研究为为极地海洋资源探测、水下目标探测、海底极地海洋资源探测、水下目标探测、海底地震和海啸预警等工程实际问题提供了地震和海啸预警等工程实际问题提供了理论参考理论参考。关键词关键词 弹性板弹性板，非线性水声波非线性水声波，色散关系色散关系 Study on Approximate Dispersion Relation of Nonlinear Hydroacoustic Wave Ludan Fu1,2,Ping Wang1,2*1College of Mathematics and Physics,Qingdao University of Science and Technology,Qingdao Shandong 2Institute of Mathematics and Interdisciplinary Science,Qingdao University of Science and Technology,Qingdao Shandong Received:Aug.21st,2023;accepted:Aug.28th,2023;published:Sep.28th,2023 Abstract Currently,in the research of nonlinear hydroacoustic waves,it is commonly assumed that the dis-persion relation is linear in order to simplify calculations.However,the linear dispersion relation *通讯作者。付鲁丹，王苹 DOI:10.12677/ijfd.2023.113009 95 流体动力学 is not suitable for large-amplitude waves.So,this study investigates the dispersion relation of nonlinear hydroacoustic waves in a finite depth compressible single-layer fluid covered by an elastic plate.The fluid is assumed to be inviscid,compressible,and irrotational.We construct the governing equations and boundary conditions that represent the relationships between hydrody-namic,elastic,and inertial forces,and obtain an approximate nonlinear dispersion relation for hydroacoustic waves.The characteristics of hydroacoustic wave modes are analyzed,and the ef-fects of important physical factors such as the thickness of the elastic plate on wave propagation properties are discussed.This research provides theoretical references for engineering practical problems such as polar marine resource detection,underwater target detection,seafloor seismic and tsunami warning.Keywords Elastic Plate,Nonlinear Hydroacoustic Wave,Dispersion Relation Copyright 2023 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 水声波的传播特性对于极地海洋环境中的水下目标探测、地震引发海啸的机理研究和预警机制具有重要影响。目前，关于线性水声波的传播特性研究已有了大量的重要成果。然而，在极地特殊海洋环境中存在着大面积冰层覆盖、垂向密度分层和底部粘性等复杂因素，非线性水声波在这些环境中的传播特性、衰减、反射和吸收机制等问题仍然是未解决的关键科学问题，而且对于极地环境中的海洋资源探测、水下目标探测、海底地震和海啸预警具有重要的参考价值。Mikhin 和 Morozov(2008)1分析了线性水声波在冰盖覆盖下的传播特性。S.Michele(2020)2等人利用线性薛定谔方程来描述大的时空间尺度上水声波的传播特性。Zhang 等(2015)3对冰洋环境中水声波传播进行了数值模拟和理论分析。刘亚东(2022)4探究表面弹性冰层与受扰动影响的可压缩流体之间的声弹性问题，利用 Green 函数得到波形的积分形式解。线性模型将海洋中的冰层视为弹性板，并运用水弹性理论来分析波浪作用下冰层的响应。这种线性模型在描述水声波与冰层的相互作用行为方面提供了有价值的结果。它能够有效预测波浪在冰层上的传播特性，并为水声工程设计和应用提供重要的参考。然而，线性模型对实际海况不再适用。首先，线性模型仅考虑了波浪和冰层之间的线性相互作用。在复杂的海洋环境中，大振幅波浪和孤立波等非线性效应可能对冰层的响应产生显著影响。其次，线性模型在考虑冰层的复杂情况方面有所限制。实际的冰层可能具有不均匀的厚度、刚度和密度分布，也可能存在断裂、缝隙以及部分融化等复杂情形。现有的线性模型很难描述这些复杂的冰层特征。因此，研究非线性水声波问题对于更准确地预测冰层的行为和保证工程的可靠性具有重要意义。张世功等人(2016)5针对一维非线性水声波传播问题进行了有限元仿真和实验研究，推导了一维非线性声波方程的有限元形式；此外，Yang(2018)6首次利用同伦分析方法(HAM)对声重力波的稳态共振三元相互作用进行了非线性解析近似方法的研究。研究发现，共振水声波引起的海底微震较非共振情况下所引起的海底微震更大。Kadri(2016)7对非线性声重力波三元共振的时间反转类比进行了研究；最近，Kadri和 Wang(2021)8在柱坐标下研究了声重力波非线性三元相互作用的近似解，将 Yang 的结论扩展到了三Open AccessOpen Access付鲁丹，王苹 DOI:10.12677/ijfd.2023.113009 96 流体动力学 维范围。除了共振现象，水声波也可以由海底地震引发。胡晨彤(2020)9推导出海底振动与海水运动之间关系的理论解，并通过使用 ADINA 计算结果和测量的地震数据验证了这个理论解。Abdolali、Kadri和 Parsons(2018)10等学者研究了水声波在深度变化的水中传播时底部压力和色散关系的变化。本文通过考虑有限深单层无粘无旋弱可压流体和在水面上覆盖无限长弹性板的情况，建立波动方程以及控制方程和边界条件。在已有的线性色散关系的基础上，通过引入一个参数来修正色散关系，从而得到一种近似的非线性色散方程，并进一步详细分析了实际极地海洋环境中水声波的传播特性，包括频率差异和非线性增强的影响，从而对极地海洋资源的探测、地震引发海啸机理的研究和预警机制的改进具有重要意义。2.数学模型数学模型 建立如图 1 直角坐标系 oxz，其中 oz 轴以垂直向上为正，ox 轴与静止时水面重合，水平向右为正。厚度为 d 的无限长弹性板漂浮在水面上，吃水为 0。假设有限深流体是无粘的、可压的且流动是无旋的，则速度势(),x z t=满足波动方程()222210,0,hzct=(1)其中2为 Laplace 算子，c 为流体中的声速，t 为时间，h 为流体深度。Figure 1.Coordinate diagram of nonlinear interaction between underwater sound and infinite elastic plate in a single layer fluid 图图 1.单层流体中水声与无限长弹性板的非线性相互作用坐标示意图 在弹性板和流体交界的表面(),zx t=，运动学边界条件和动力学边界条件分别为 0,txxz+=(2)210,2pgt+=(3)其中 g 是重力加速度，(),p x t是板下表面压力，为流体密度。假设弹性板是均匀的薄弹性材料，且流体和板不分离。在流体底部，根据固壁的非渗透条件有()0,.zhz=(4)3.近似的非线性色散关系近似的非线性色散关系 给定2222,2,ennDTddfDTkggc=，Abdolali、Kadri 和 Parsons 10等人得到线性色散关系如下 付鲁丹，王苹 DOI:10.12677/ijfd.2023.113009 97 流体动力学 ()()()242tan1,1tannnnnnnghDkTkdh=+(5)其中为水声波的频率，()1,2,3,nkn=为水声波的波数，T 为侧向应力，e为弹性板的密度。在复杂的海洋环境中，大振幅波浪的非线性效应对水声波的传播产生显著影响，真实的波频率与线性波频率存在偏差，为较准确地描述非线性水声波在复杂海洋环境中的传播特性，本文引入,=(6)其中是一个略大于 1 的常数。增大时，水声波的非线性增强。将(6)代入(5)得()()()2242tan1,1tannnnnnnghDkTkdh=+(7)给定常数,eg hc E v T d 的值，公式(7)是一个 f 关于nk的复杂超越方程，直接求得 f 与nk的色散关系比较困难，我们需要对公式(7)做一些变换。根据 Abdolali、Kadri 和 Parsons 10等人的工作，我们引入 2,aagKKd (8)它解释了重力相对于弹性板振动的影响。当1aK 时，有()1nhn，令()1,nhn=+(9)是一个无限接近于 0 的量。色散方程变为()()()2422tan11,1tan1+=+nnnngnDkTkdn (10)整理得()2422tan1,1tan=+nnnngDkTkd (11)对方程(11)进行()tan 麦克劳林级数展开并保留主导项，整理后色散关系变为：()24221,1nnnngDDkTkd=+(12)又由(9)得()()()()()242211,11=+nnnnnnghnDkTkdhn (13)由于色散方程(13)的右端是一个关于n的偶函数，所以本文仅考虑正轴的情况。由2222nnkc=知，当，则()1,2,3,nknN=是实数，对应 N 个水声波模态的波数。付鲁丹，王苹 DOI:10.12677/ijfd.2023.113009 98 流体动力学 令2339.8 m/s,1500 m,1024 kg*m,917 kg*m,1500 m/s,5 GPa,0.3,0=eghcEvT，代入色散方程(11)可以得到水声波的波数频率关系。首先，考虑对水声波传播特性的影响。分别取1=与1.01=，代入色散方程(11)，得到水声波的波数频率关系(见图 25)，可见同频率下水声波波数随着的增大而逐渐增大；频率较低时对水声波波数的影响较大，频率较高时对水声波波数敏感度较小。而且，频率较低时波数随频率变化较慢，频率较高时波数随频率变化较快。而且，从这些图中可以看出，我们的结果与 Abdolali、Kadri 和 Parsons 10等人得到的线性色散关系结果大体一致，近似非线性色散关系就是对该线性色散关系线性修正后获得，从而说明我们的结论是准确的和科学的。Figure 2.Wave number-frequency relation of the second mode hydroacoustic wave when the thickness of the elastic plate is 40 meters 图图 2.弹性板厚度为 40 m 时第二模态水声波的波数频率关系 Figure 3.Wave number-frequency relation of the second mode hydroacoustic wave when the thickness of the elastic plate is 60 meters 图图 3.弹性板厚度为 60 m 时第二模态水声波的波数频率关系 付鲁丹，王苹 DOI:10.12677/ijfd.2023.113009 99 流体动力学 Figure 4.Wave number-frequency relation of the third mode hydroacoustic wave when the thickness of the elastic plate is 40 meters 图图 4.弹性板厚度为 40 m 时第三模态水声波的波数频率关系 Figure 5.Wave number-frequency relation of the third mode hydroacoustic wave when the thickness of the elastic plate is 60 meters 图图 5.弹性板厚度为 60 m 时第三模态水声波的波数频率关系 其次，考虑弹性板 d 对水声波传播特性的影响。基于上述色散方程表示的色散关系，分别取弹性板厚度为20,40,60,80 md=，取1=，得到水声波的波数频率关系(见图 6、图 7)，对于水声波前三阶模态，在频率较低时，随着板厚度的增加，同一频率下的波数会增加，这意味着波长随着板厚度的增加而缩短，这时波的传播受到弹性板厚度的影响较大；在频率较高时，弹性板厚度对波的传播的影响逐渐减小，这意味着随着频率的增加，波数的变化对弹性板厚度的敏感性减弱。最后，考虑水声波模态对传播特性的影响，分别取水声波模态为2,3,4n=，20d=，1.01=和1=。得到水声波模态对传播特性的影响(见图 8、图 9)。对于水声波前三阶模态，随着模态增大，同一频率下的波数会减小；频率较低时，不同模态对应的波数间差异较大，随着频率增大，这种差异逐渐减小，可见在频率较低时，同频率下波数受模态的影响较大，而在频率较高时，同频率下波数受模态的影响较小。付鲁丹，王苹 DOI:10.12677/ijfd.2023.113009 100 流体动力学 Figure 6.Wave number-frequency relation of water acoustic wave when 1=图图 6.1=时水声波的波数频率关系 Figure 7.Wave numbe-frequency relation of water acoustic wave when 1.01=图图 7.1.01=时水声波的波数频率关系 Figure 8.Influence of water acoustic wave modes on propagation characteristics for 1=图图 8.1=时水声波模态对传播特性的影响 付鲁丹，王苹 DOI:10.12677/ijfd.2023.113009 101 流体动力学 Figure 9.Influence of water acoustic wave modes on propagation characteristics for 1.01=图图 9.1.01=时水声波模态对传播特性的影响 4.结论结论 在势流理论的框架下，本文研究了无限长弹性板覆盖下的有限深可压单层流体中非线性水声波的色散关系。主要结论如下：1)非线性影响：引入一个参数来修正线性色散关系，以考虑大振幅波浪对水声波的非线性影响。频率较低时对水声波波数的影响较大，频率较高时对水声波波数影响较小。2)弹性板厚度影响：通过分析弹性板厚度对水声波传播特性的影响，发现弹性板的厚度增加会导致同一频率下的波数增加。在高频率下，波的传播受到冰层厚度的影响相对较小，其他因素可能更加显著地影响波的传播行为。3)水声波模态影响：研究了不同水声波模态对传播的影响。发现随着模态增加，同频率下的波数减小，频率较低时模态变化对波数影响较大，而高频时影响较小。这些研究结果对于理解非线性水声波在复杂海洋环境中的传播特性具有重要意义。这不仅对极地海洋资源的探测、地震引发海啸的机理研究和预警机制的改进提供了参考，还为工程应用中的水声传播提供了理论指导。未来的研究可以进一步探讨更多复杂环境中的非线性水声波传播特性，以更全面地揭示其行为和影响。参考文献参考文献 1 Mikhin,M.and Morozov,A.(2008)Acoustic Wave Propagation in an Ice-Covered Ocean.The Journal of the Acous-tical Society of America,123,951-963.2 Michele,S.and Renzi,E.(2019)Effects of the Sound Speed Vertical Profile on the Evolution of Hydroacoustic Waves.Journal of Fluid Mechanics,883,A28.https:/doi.org/10.1017/jfm.2019.907 3 Zhang,J.,Li,Z.and Jiang,W.(2015)Theoretical and Numerical Investigations of Acoustic Propagation in Ice-Covered Oceans.The Journal of the Acoustical Society of America,137,2520-2530.4 刘亚东.海洋表面冰盖、密度分层、底部粘性对水声波传播特性的影响D:硕士学位论文.上海:上海大学,2022.5 张世功,吴先梅,张碧星,安志武.一维非线性声波传播特性J.物理学报,2016,65(10):8.6 Yang,X.,Dias,F.and Liao,S.(2018)On the Steady-State Resonant Acoustic-Gravity Waves.Journal of Fluid Me-chanics,849,111-135.https:/doi.org/10.1017/jfm.2018.422 7 Kadri,U.(2016)Generation of Hydroacoustic Waves by an Oscillating Ice Block in Arctic Zones.Advances in Acous-tics and Vibration,2016,Article ID:8076108,1-7.https:/doi.org/10.1155/2016/8076108 付鲁丹，王苹 DOI:10.12677/ijfd.2023.113009 102 流体动力学 8 Kadri,U.and Wang,Z.(2021)Approximate Solution of Nonlinear Triad Interactions of Acoustic-Gravity Waves in Cylindrical Coordinates.Communications in Nonlinear Science and Nu-Merical Simulation,93,Article 105514.https:/doi.org/10.1016/sns.2020.105514 9 胡晨彤.考虑海水可压性的海底震动与海水运动关系及其在海啸生成中的应用D:硕士学位论文.上海:上海交通大学,2020.10 Abdolali,A.,Kadri,U.,Parsons,W.and Kirby,J.T.(2018)On the Propagation of Acoustic-Gravity Waves under Elastic Ice Sheets.Journal of Fluid Mechanics,837,640-656.https:/doi.org/10.1017/jfm.2017.808