水下稳态涡流场声传播的数值模拟和实验研究.pdf
《水下稳态涡流场声传播的数值模拟和实验研究.pdf》由会员分享,可在线阅读,更多相关《水下稳态涡流场声传播的数值模拟和实验研究.pdf(7页珍藏版)》请在咨信网上搜索。
1、声波穿过涡流场时声传播特性会发生改变,可由此判断涡流场的存在及其特征参数。文章首先推导了二维稳态涡流场的射线微分方程组,实现了声波通过涡流场的声线轨迹模拟。然后开展水池超声传感器实验,采用现代计算流体力学(Computational Fluid Dynamics,CFD)仿真方法获取涡流场参数。最后将通过涡流场声信号的时延结果,与实验测量值进行对比分析。结果表明,对于涡核半径为30 mm,最大切向速度为1.2 ms-1的涡流场,数值模拟误差小于10%,验证了基于射线声学的水下稳态涡流场声传播数值模拟的有效性。关键词:射线声学;稳态涡流场;声传播;时延中图分类号:TB566 文献标志码:A 文章
2、编号:1000-3630(2023)-04-0419-07Numerical simulation and experimental study of sound propagation in underwater steady vortex fieldYU Min1,FAN Dingfan1,ZHANG Ye2,LIU Hang1(1.School of Naval Architecture,Ocean and Energy Engineering,Wuhan University of Technology,Wuhan 430074,Hubei,China;2.Systems Engine
3、ering Research Institute,CSSC,Beijing 100094,China)Abstract:When sound waves pass through vortex field,the propagation characteristics of sound waves are changed,and hereby,the existence and characteristic parameters of vortex field can be judged.In this paper,the ray differential equations of two-d
4、imensional steady vortex field are derived,and the trajectory of sound wave passing through vortex field is simulated.Then,the ultrasonic sensor experiment in pool is conducted and the parameters of vortex field are obtained by modern CFD simulation method.Finally,the numerical simulation results of
5、 acoustic signal delay in vortex field are compared with the experimental measurement results,and it is shown that the numerical simulation error is less than 10%for the vortex field with a core radius of 30 mm and a maximum tangential velocity of 1.2 ms-1,which verifies the validity of numerical si
6、mulation of sound propagation in steady underwater vortex field based on ray acoustics.Key words:ray acoustics;steady vortex field;acoustic propagation;time delay0引 言声波探测作为一种直接的、非侵入式的和全局式的测量方法,可用于水下旋涡特性的测量1-2。声波在穿过涡流场时,其传播路径发生变化3-4,由此产生声波相位的变化,相位变化量和涡流场的形态特征密切相关。因此,发展基于涡流场声传播特性发展水下涡流场探测识别技术成为可能,在军事和
7、工程中具有应用前景。使用声波探测方法进行涡特性的计算首先在超声领域得到了应用。1992年,文献5-6提出了在非均匀介质(涡流场)中,引入光学中相位共轭镜的方法对目标的反射声进行聚焦,以达到定位的目的。1997年,Roux等7通过实验从接收信号的幅值与相位差两个方面研究声信号穿过涡流场后产生的变化。从1998年开始,文献8-10先后使用多种实验装置对水下涡流场进行深入研究,指出涡流同时改变了声波的传播方向和传播速度。2001年,Pagneux等11基于射线声学研究了二维空间中的理想涡对高频声波的散射。2005年,Rosny等12针对低马赫数涡流场引起的声信号相位跳变量较小、不易测量的问题,采用了
8、TRM技术来增大相位跳变量、降低干扰,取得了较好的效果。2023年,陈雨晨等13针对介质运动对声传播的影响,建立一种利用高斯波束法求解亚音速运动介质中的声传播模型,认为在海流会明显影响声的传播。上述文献研究中,多采用理想涡流速分布公式,甚少通过现代计算流体力学(Computational 引用格式:喻敏,樊丁繁,张烨,等.水下稳态涡流场声传播的数值模拟和实验研究J.声学技术,2023,42(4):419-425.YU Min,FAN Ding-fan,ZHANG Ye,et al.Numerical simulation and experimental study of sound prop
9、agation in underwater steady vortex fieldJ.Technical Acoustics,2023,42(4):419-425.DOI:10.16300/ki.1000-3630.2023.04.002收稿日期:2023-03-02;修回日期:2023-06-27基金项目:国家自然科学基金(面上基金,52271349)。作者简介:喻敏(1978),女,湖北武汉人,博士研究生导师,副教授,研究方向为水声定位技术,水下目标识别。通信作者:喻敏,E-mail:2023 年声学技术Fluid Dynamics,CFD)方法获取涡流场流速参数。文献13虽然采用射线声学
10、方法对运动介质中的声传播进行研究,但没有涉及复杂流场中的声传播问题。本文针对声波通过水下稳态涡流场的数值模拟方法及实验展开研究,首先根据稳态涡流场中的射线微分方程数值模拟声线传播轨迹,然后开展水池超声传感器实验,最后将数值模拟和实验测量的通过涡流场的声信号时延量进行对比验证,证明了本文通过射线声学方法数值模拟水下稳态涡流场声传播的有效性。1涡流场声线轨迹数值模拟涡流场声传播特性的改变,主要表现为通过涡流场的声波振幅与相位的变化14。涡流场声传播效应示意图如图1所示,顺流区声速加速,逆流区声速减速,最终导致在接收端声波的相位超前或滞后。射线声学是几何声学的近似理论,适用于解决高频情况下的声传播问
11、题。基于射线理论可以有效模拟声波通过涡流场的声线轨迹及信号变化,具有直观、计算效率高的特点。为模拟声波通过涡流场的声线轨迹,根据稳定移动介质中的程函方程15,推导出可进行声线轨迹模拟计算的稳态涡流场射线声学微分方程组:ddx=(r)ucosdydx=tan(1)式中:(r)为涡度;表示声线掠射角;u表示涡流场速度。声线轨迹数值仿真示意图如图2所示。在理论计算过程中,对于某一条声线s,固定每步迭代中x方向的增量Dx,通过式(1)获取y方向的增量Dy,进而确定该声线的轨迹。声线通过涡流场后的数值模拟轨迹如图3所示。由图3可知,通过涡流场的声线掠射角发生偏移,越靠近涡核中心其掠射角变化越大,并在涡核
12、后方形成声线焦散区。仿真中,涡流场切向速度(简称:流速)r分布采用二维稳态兰金涡流(Rankine vortex):v(r)=2ara ra21r ra(2)式中:r表示涡流场某一测点到涡核中心点的距离,a为涡流场流速最大处涡核特征半径;=2acMa,定义为涡流环量;c=1 500 ms-1为水中声传播速度;Ma为马赫数,定义为涡流最大切向速度与声速之比。本文数值模拟结果与文献16实验结果进行对比,在相同条件设置下,声信号相位差计算结果如图4所示。由图4可见,相位差的数值模拟结果与图1 涡流场声传播效应示意图Fig.1 Schematic diagram of sound propagatio
13、n effect in vor-tex field图2 声线轨迹数值仿真示意图Fig.2 Schematic diagram of sound ray numerical simulation-40-2002040图3 a=10 mm,Ma=4.610-2时的声线传播轨迹Fig.3 Sound ray propagation trajectory when a=10 mm,Ma=4.610-2图4 声信号相位差本文数值模拟结果与文献16实验结果的对比Fig.4 Comparison between the numerical simulation result in this paper an
14、d the experimental result from reference 16 regarding the phase difference of acoustic signal420第 4 期喻敏等:水下稳态涡流场声传播的数值模拟和实验研究文献实验结果整体上吻合较好。为了进一步验证数值模拟方法,本文设计并搭建了超声测量涡流场实验平台,开展相关实验验证。2实验装置及涡流场参数2.1实验仪器设备实验设备实物图如图5(a)所示。实验通过螺旋桨产生稳态涡流场,螺旋桨直径d=60 mm,固定在尺寸为2 m0.6 m0.8 m的亚克力长方体水池中,实验布置如图5(b)所示。螺旋桨位于整涡管顶部,
15、由电机带动旋转产生涡流。涡流流出整涡管后,使用超声传感器发射和接收通过涡流的声信号。为保证两个超声传感器中轴在同一条直线上,超声传感器刚性固定在同一移动滑台上。水池装有导轨和滑尺,便于精确标定传感器位置。针对螺旋桨高速旋转产生的振动,在实验中将造涡装置和测量装置分别固定在两组滑台上,在螺旋桨与其固定滑台联接处使用螺丝螺母固定并加垫片以减少振动。仪器设备型号及参数如表1所示。2.2涡流场CFD仿真由于实验中的涡流场速度分布不易获取,本文基于现代流体仿真软件Fluent,根据实验布置对涡流场分布进行数值仿真。本文使用Design Modeler建模软件建立三维模型。计算域的划分与边界条件的设置如图
16、6所示。采用ICEM前处理软件划分网格,使用六面体结构化网格对整体区域进行划分,并在旋转域表面进行局部加密,以更好地模拟近壁流动。螺旋桨所在的旋转域采用稳态(Moving Reference Frame,MRF)方法,并绕z轴以1 400 rmin-1的转速旋转。基于有限体积法求解RANS方程,湍流模型选用标准k-模型,调用标准壁面函数进行计算。采用三维双精度基于压力隐式求解器,耦合方式采用 SIMPLE 算法,离散格式采用二阶迎风格式。2.3涡流场参数定义整涡管出口所在平面为z=0平面,分别在整涡管出口下端z=10、30和50 mm处取截面,得到涡流场速度分布矢量图如图7所示,其中z=10
17、mm图5 实验设备实物图及实验布置示意图Fig.5 Photo of the experimental equipment and schematic di-agram of experiment layout表1实验仪器设备型号及参数Table 1Model and parameters of experimental equipment实验设备螺旋桨整涡管开发板发射/接收传感器示波器型号定做定做DYSTM32型DYW-2M-01T型TDS2014C型参数4桨叶,直径60 mm,可调电源电压624 V亚克力材质,长240 mm,内径65 mm,外径70 mm总线宽度32位,最高频率168 M
18、Hz直径35 mm,波束宽度:41,可用量程0.032 m4通道,带宽100 MHz,最高采样率2.0 GHz图6 数值模拟计算域Fig.6 Numerical simulation computing domain4212023 年声学技术截面上的流速分布云图如图8所示。由图7、8可以看到涡流场在流出整涡管后向外消散,符合自由场涡流运动规律。涡流流速最大位置对应涡核半径尺寸,由图8可知,z=10 mm平面上的最大速度umax=1.258 ms-1,判断涡核半径a=30.6 mm。根据z=10、30和50 mm处的流场速度数据可得到不同截面的涡流场主要参数,如表2所示,为基于射线声学的涡流场声
- 配套讲稿:
如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。