新型通海阀流噪声的数值模拟对比研究.pdf

第45卷第1 6 期2023年8 月舰船科学技术SHIP SCIENCEANDTECHNOLOGYVol.45,No.16Aug.,2023新型通海阀流噪声的数值模拟对比研究贾文尖，蒋鹏，史文祥，娄燕鹏(中国船舶集团有限公司第七二五研究所，河南洛阳47 1 0 0 0)摘要：角式双球阀作为一种新型通海阀，针对其内部流道结构在不同流速下产生的噪声问题，基于计算流体动力学（CFD）与直接边界元法（BEM）对其进行流噪声数值模拟研究，分析2 种不同结构形状角式双球阀的流噪声声压级以及其不同流速下的流噪声变化规律。结果表明，2 种结构声压级均随频率的升高而降低，低流速时，双球阀低频特性更明显；开放式流道结构流噪声声压级明显高于封闭流道结构，由于封闭流道结构内部流体较好的流通性，其噪声总声压较开放式结构下降约2 6.9%，降幅为1 7 dB(A)，封闭结构设计可以有效降低通海阀内流道流体脉动压力与流噪声，从而为新型通海阀的降噪优化提供一种新思路。关键词：角式双球阀；通海阀；流噪声；声压级；流速中图分类号：TH134文章编号：1 6 7 2-7 6 49(2 0 2 3)1 6-0 0 48-0 7Comparative study on numerical simulation of flow noise of a new type of sea valveAbstract:Angle double ball valve is a new type of sea valve,aiming at the noise problem caused by its internal flowchannel structure at different flow rates,the flow noise numerical simulation of sea double ball valves is carried out based oncomputational fluid dynamics(CFD)and direct boundary element method(BEM).The flow noise sound pressure level oftwo kinds of angle double ball valves with different structure shapes and the variation law of flow noise at different flowrates are analyzed.The results show that the sound pressure level of the two structures decreases with the increase of fre-quency,and the low frequency characteristic of double ball valve is more obvious at low flow rate.The sound pressure levelof the flow noise of the open channel structure is obviously higher than that of the closed channel structure.Due to the goodfluid fluidity of the closed channel structure,the total noise sound pressure level of the closed channel structure is reduced byabout 26.9%compared with the open structure,with a decrease of 17dB(A).The closed structure design can effectively re-duce the fluid pulsating pressure and flow noise in the inner channel of the sea valve,which provides a new idea for the noisereduction optimization of the new type sea valve.Key words:angle double ball valve;sea valve;flow noise;sound pressure level;flow velocity0引言通海阀作为船舶管路系统重要的控制单元，是海水管路系统中的第一道阀门，但同时也是管路系统中主要的噪声源。由于通海阀的节流作用，阀内紊流程度高，流体对流道的冲刷严重，流体介质产生的压力脉动会直接辐射产生流噪声，流噪声的频率具有典型的宽频特性2 ，传播距离远，范围广，是影响水下装备隐蔽性的关键因素之一。随着噪声探测技术的飞速发展，有效降低阀门的噪声，提升装备隐蔽性已成收稿日期：2 0 2 2-0 8-2 5作者简介：贾文尖(1 9 9 7-），男，硕士研究生，研究方向为振动与噪声。文献标识码：AJIA Wen-jian,JIANG Peng,SHI Wen-xiang,LOU Yan-peng(The 725 Research Institute of CSSC,Luoyang 471000,China)doi:10.3404/j.issn.1672-7649.2023.16.010为研究的主要课题3。目前，针对海水管路阀件噪声的研究大多数基于对截止类阀、调节阀及管路进行声学数值模拟分析。方超等4 基于单向流固耦合，采用声学边界元法对某通海阀流噪声及流激振动噪声进行对比分析研究，结果表明流激振动噪声完全潼没在流噪声中，对该系统通海阀噪声进行治理时应该优先考虑流噪声。徐国栋等5 以船舶海水冷却系统通海阀出口管路为对象，采用FEM/AML技术，进行水下辐射噪声分析，对比不第45卷同通海阀水下深度、出口管径对流噪声的影响，结果表明海水冷却系统通海阀水下流噪声以低频为主。桂瞬丰等6 采用CFX软件对某型通海阀的噪声进行数值模拟分析，结果表明海阀低频噪声声能较大，高频部分所占比重较小。聂欣等7 根据LES和Lighthill声类比理论，采用流体力学软件Fluent和声学仿真软件ACTRAN，研究计算了直简笼式阀门在相同流量不同开度下阀门内外部的流噪声特性。结果表明，当阀门开度降低时，阀门内的湍流强度、速度和压力波动均相应增强，导致阀门内外的整体声压分贝变强，阀门内部噪声能量也增强，阀门噪声具有偶极子特性。WeiLin等8 基于CFD，应用FW-H模型来研究高压减压阀中的流噪声特点，发现进口压力对声压级有很大影响，最大马赫数接近1 时，低频的声压级会非常高。LiuKai等9 提出一种基于计算流体力学和Lighthill声学类比的理论的混合方法来模拟管路流噪声，结果显示该方法模拟结果与实验数据吻合较好。近几年，采用角式双球阀作为通海阀的方案已被提出并实现初步应用，但针对角式双球阀噪声方面的分析研究与数据不完善，因此研究其噪声特性对其后续的结构优化及应用具有重要的科学意义和实际价值。本文以2 种新型角式双球阀结构为研究对象，利用CFD技术及直接边界元法（BEM），通过对双球阀全开时内部流场进行数值模拟，联立Fluent软件和LMSVirtual.Lab软件对双球阀的流噪声进行计算研究，分析不同流速与流噪声频谱特性间的关系，获得噪声声压级的数值大小，提供数据支撑，为进一步提高通海阀的降噪优化水平优化提供数据支撑和理论指导。1角式双球阀的结构设计1.1结构组成和设计参数角式双球阀采用一体化双球体结构，主要由阀体、阀盖、球体、阀座、等零部件组成。开放流道角式双球阀结构如图1 所示，封闭流道角式双球阀三维模型方案图如图2 所示。1.2流道结构设计特点1）开放式流道角式双球阀。为角球阀与直通球阀连通的常规设计，一道阀与二道阀不成封闭流通结构，部分介质流出一道阀后进入球与阀体的空腔内；同时角球阀采用9 0 转角球结构，转角半径较小。2）封闭流道角式双球阀。通过采用角球球体中心与出口流道具有偏心距设计，增加弯头使一道阀与二道阀内流道实现封闭连通；同时也增加了角球流道的贾文尖，等：新型通海阀流噪声的数值模拟对比研究一道阀9 0 阀盖角球二道阀闽体直球阀座阀座档环图1 开放式流道结构图Fig.1 Structure diagram of unclosed flow channel一道阀加长弯头角球阀盖二道阀阀座直球阀体图2 封闭流道结构图Fig.2 Structure diagram of closed flow channel转弯半径，约为开放流道角球转弯半径的1.5倍。2流噪声数学理论基础2.1流场计算模型2.1.1k-8模型各流速下雷诺数均大于30 0 0，由GB/T30832-2014阀门流量系数和流阻系数试验方法A.3.1测试管道内流体的雷诺数高于30 0 0，确认该流动为端流，因此稳态分析时选择标准k-s模型。1 9 7 2 年，Launder 等1 0 提出了标准k-8流模型，在一方程模型的基础上引人了关于端流耗散率的方程，k方程为精确方程，8 方程是由经验公式推导而来。标准k-模型假定流场全是端流，分子间的粘性可以忽略，该模型计算简单，结果精度高，在工程实践中得到广泛应用。标准k-模型表达式分别为：k方程a(kp)(kul)ot0 x1Gk+Gb-Px-YM+S2,：方程49阀座档环ok从+OXF+oxfk(1)500(08)十otC18(G+C.CG)-Cap%+1.k式中：SE、S k 为湍能耗散率及端动能的源项；Ci、C2、C3为设定的经验系数这里取Ci=1.44、C2=1.9 2、C3=0.09；YM为可压缩流体中脉动扩张量；G,为由于浮力引起的湍动能k的产生项；Gk为由于速度梯度所引起的k的产生项。端动能k和耗散率。的端流Prandtl取值分别为1.0、1.3。2.1.2大涡模拟瞬态分析采用大涡模拟（LargeEddySimulation，LES），即在一定的简化条件下直接求解N-S方程，可将湍流流场中的涡流分为不同尺度的的涡流，小涡采用时均化处理。通过近似的亚格子模型进行考虑，大涡保留与时间相关项而可以通过数值方法直接求解，大涡模拟所得端流解包含丰富的频率信息，更有利于分析其频谱特性1-1 2 。2.2声学理论1952年，Lighthill在研究喷气噪声时基于流体力学中的N-S方程提出了声比拟方法，得到了著名的莱特希尔方程1 3-1 4。后来Curlel15 在Lighthill 方程进行了扩展，得到了考虑固体壁面的气动声学声类比方程，该方程表达了由固体壁面由压力脉动引起的偶极子噪声。运用Curle的积分解公式，有壁面边界存在时的流噪声问题得到很好的求解。针对此问题，Ffowcs-Wil-liams与Hawkings16 进一步考虑了旋转壁面的情况，推导出了FW-H方程，为流噪声预测提供了理论依据。FW-H方程为：828t2axiXj综上，可看出噪声声源为单极子、偶极子和四极子，在声学模拟中，单极子存在于流体速度低时的不稳定状态，一般不作考虑。流体遇到壁面由于脉动压力产生的偶极子源与动量能量交换产生的四极子源辐射声功率与流速分别成6 次方或8 次方关系。因此，当流体马赫数较低时，不需要考虑四极子声源，研究重点应关注在壁面脉动压力构成的偶极子声源上1 7 。3数值模拟及分析3.1流场数值模拟运用Fluent软件进行稳态及瞬态流场模拟，流场舰船科学技术分析过程中为了让前后的流体充分发展，利用Ansys/从+xixf第45卷+Qz.xr-,pound(f)。Workbench前处理模块建立流体介质的三维流体域模型，为了更准确模拟流经通海阀2 种不同流道时流体(2)的流动状态，根据GB/T30832-2014阀门流量系数和流阻系数试验方法，通海阀的前端人口选取5D的流场通道，后端出口选取1 0 D的流场通道。为保证计算精度，在模型处理环节上尽量保留几何细节，网格划分采用四面体，检验网格质量并进行网格的无关性验证。流场分析过程中，流体状态为三维不可压缩，属流流动，本次模拟介质为水，根据通海双球阀实际工况，入口为速度边界条件，分别取1 m/s、2 m/s、3m/s、4m/s；出口为压力边界条件，为固定压力5MPa，收敛精度为0.0 0 1，壁面为无滑移面。模拟计算步骤为：首先应用标准k-模型对流场进行稳态求解，流动达到稳定状态后，再改用LES大涡模拟进行瞬态运算1 8-1 9 。这种方法不需要舍去非定常过程中初始部分结果，更节省时间，而且已被证明可靠。本文计算中时间步长取0.0 0 0 2 5s，采样总时长为0.1 S。3.2流噪声数值模拟将前述流场计算数据保存，提取流体域与双球阀耦合区域壁面压力脉动信息，将模拟计算的不同流速下流固耦合面流体压力脉动信息以CGNS格式导人LMSVirtual.Lab软件中，基于单向流固耦合，将流场外壁面的压力脉动传递至通海阀流道内壁面，经快速傅里叶变换，对传递后的壳体内壁面的压力脉动时域信息进行数据转换，得到频域信息，声学计算流程如图3所示。本文研究流体马赫数较低，因此将此压力脉动作为壁面偶极子声源激励，采用直接边界元法计算流体流经通海阀产生的流噪声，频率为1 0 2 0 0 0 Hz。(3)CFD计算CGNS偶极子声源快速傅里叶变换(FFT)时域频域图3声学计算流程Fig.3 Acoustic calculation process为了后续实验测量流经通海阀的流噪声声压值，根据标准GJB4058-2000舰船设备噪声、振动测量方法的规定：测量表面与基准体的测量距离应不小于1m，一般取1 m，故设定监测点的位置如图4所示，通过SI、S 2、S 3、S 4四个监测点分析双球阀的流噪声特性。边界元法(BEM)声学计算S第45卷Fig.43D model of double ball valve and pipeline3.3流噪声影响分析3.3.1不同流速下声压频响特性曲线分析根据上述模拟实验条件，对2 种角式双球阀进行声学响应分析，得到监测点流噪声声场信息，将不同7560(V)aP/里453015贾文尖，等：新型通海阀流噪声的数值模拟对比研究SSS,图4双球阀及管道三维模型S4(V)aP/单453015051速度下各监测点的声压级数据进行整理，定量对比通海双球阀2 种结构在不同流速下的噪声特性。不同流速下通海双球阀的声压频率特性曲线如图5所示。分析可知，通海阀在不同速度下的频谱规律具有明显的宽频特性。4个监测点不同频率下声压级的极大值区域与极小值区域呈不规则交错分布，但封闭、开放式2 种结构的频响曲线总体变化一致，声压级均随着频率的升高呈下降趋势，4个监测点的最大峰值声压级频率集中在2 0 0 40 0 Hz，即低频声能量所占比重大，中高频声能量所占比重较小。进一步验证在实际使用环境中，阀门流噪声多呈低频性。通过上述分析可知，2 种结构各监测点的声压随频率的总体变化趋势一致，为探究声压频响特性曲线随速度的变化规律，提取4个流速下在监测点S2处的声压频响特性曲线，如图6 所示。可知，在S2处，频率为1 0 2 0 0 0 Hz范围内,速度不同时各个频率下的声757560S.S460(V)aP/军4530150S2S3S4(a)封闭结构速度1 m/s7560(v)aP/王4530150(b)开放结构速度1 m/s7560S3(V)AP/王单S4530150(c)封闭结构速度2 m/sSI75SSS4SS60S(V)P/王里4530150(d)开放结构速度2 m/s7560(V)AP/单4530150200(g)封闭结构速度4 m/s图5不同监测点的声压频响曲线对比Fig.5 Comparison of sound pressure frequency response at different monitoring points(e)封闭结构速度3m/s7560S4()P/王单4530150600100014001800频率/Hz()开放结构速度3m/sSS2S3S4200600100014001800频率/Hz(h)开放结构速度4m/s52(v)ap/402002008060(V)aP/单40200200图6 不同速度声压频响曲线对比图Fig.6Comparison of sound pressure frequency response atdifferent speeds压级相差较大，变化明显。表明2 种结构声压频响曲线的波动情况在不同流速下会有很大的区别。为了分析2 种结构频响特性曲线的差别，对比双球阀2 种流道结构4个流速下在监测点S2的频响特性曲线，如图7 所示。可知，监测点S2处，开放式流道结构声压级整体处于封闭结构上方，二者的峰值声压级对应的频率不同，极大值极小值区域波动的幅度也不同。速度为1 m/s时，封闭流道结构较大的峰值频率集中在40 0 Hz附近，而开放式流道峰值极大值频率则集中在2 0 0 Hz、40 0 H z、1 50 0 H z 附近，且其声压值均大于封闭结构，封闭、开放式结构声压级对应的最高峰值为35.6 6 dB(A)、53.36 d B(A)；速度从1 m/s变化到4m/s，封闭结构的峰值极大值点增多，但是在频率低于6 0 0 Hz时，2 种结构的波动趋势较为一致，频率越大，开放式结构波动越剧烈。在监测点S2，2 种结构在低频段，声压级下降幅度大；8 0 0 1 2 0 0 Hz，声压级较为稳定；超过1 50 0 Hz后声压级迅速降低；速度为1 m/s、2 m/s 时，随着频率增大，2 种结构声压值急剧降低，而3m/s、4m/s 时压降较缓，声压级降低一定值后即在小范围内波动；表明低速时，双球阀流噪声的低频特性更为明显。3.3.2总声压级对比分析对不同频率下各场点声压值进行加权运算得到噪舰船科学技术8060一封闭1 m/s-封闭2 m/s60封闭3m/s封闭4m/s600100014001800频率/Hz开放1 m/s开放2 m/s开放3m/s-开放4m/s600100014001800频率/Hz第45卷一V=1m/s封闭-V=1 m/s开放45(V)aP/单301506045(V)AP/里301507560(v)aP/王单45301507560(v)aP/王单4530150200图7 2 种结构声压频响曲线对比图Fig.7 Comparison of sound pressure frequency response attwo structures声总体强度的声压级2 0 ，总声压级计算式为：LpAs=101g(10*0+0%.+10%)式中：ni为频率i处的声压值；N为总采样频率数。通过对不同监测点的声压进行叠加，双球阀在102000Hz频率范围内的总声压级变化如图8 所示，得到不同流速下2 种结构总声压级如表1 所示。由图8(a)可知，封闭流道结构下，速度1 m/s、2m/s、3m/s 时总声压级均低于50 dB(A)，4m/s 时总声压级略高于50 dB(A)，接近于开放式结构速度为1m时的总声压级；当流场速度增加时，各监测点处流噪声总声压级的数值会相应地增加，但同一流速距一V=2m/s封闭-V=2m/s开放V=3m/s封闭V=3m/s开放V=4m/s封闭V=4m/s开放6001000 14001800频率/Hz(4)第45卷(V)P/田单6050-401图8 总声压级随位置、速度的变化曲线Fig.8 Variation curve of total sound pressure levels with positionand velocity表1 总声压级对比dB(A)Tab.1 Comparison of overall sound pressure level between thetwo structures速度/ms-结构SiS2封闭36.36 5539.8 8 7 9 37.52 3540.0 0 7 51开放54.7 6 0 6 56.1 9 2 0 54.9 40 7 50.8 0 9 1封闭41.1 1 59 42.48 6 541.6 357 38.9 8 562开放57.0 9 6 455.5549 57.1 52 1 48.9 2 59封闭46.1 7 7 2 46.432 8 47.1 2 8 2 47.57 1 73开放6 5.6 2 436 9.1 8 8 7 6 6.36 6 36 3.1 8 8 4封闭54.8 6 1 556.42 6 354.8 8 2 357.1 6 8 84开放7 5.6 2 2 1 7 5.9 6 6 7 7 6.59 7 7 7 0.2 7 53阀1 m处各监测点的总声压级差别不大，在小范围内波动，趋于常值；开放式流道结构在监测点S4会总声压值出现下降，表明由于内部结构不贯通，流噪声声压会在结构改变位置出现明显变化。由图8（b）可知，速度从1 m/s增大到4m/s过程中，监测点声压级总体呈增大趋势，表明随着速度的增加，雷诺数随之增加，通海双球阀壁面湍流边界层的不规则流动也变得更加剧烈，脉动压力也增强，因贾文尖，等：新型通海阀流噪声的数值模拟对比研究90+封团1 m/s开放1 m/s封闭4m/s封团2 m/s并效2 m/s并放4m/s+封闭3m/s并放3m/s75604530180封闭S1开放S1封闭S2并放S270封闭S3开放S3封闭S4开放S42速度/ms-1(b)随速度变化曲线S3S4总声压级平均降幅26.953此声场中的由脉动压力引起的流噪声声压级也增大。图中封闭结构各速度的总声压值远小于开放式结构的值，且封闭结构4m/s的总声压值也小于开放结构速度1 m/s的声压值。从表1 可看出，封闭结构相对开放式结构的总声压级平均降幅比例约为2 6.9%，降幅为1 7 dB(A)，压降明显，能够有效抑制阀门噪声。分析可知，封闭23监测点(a)随位置变化曲线344流道相较于开放式流道结构，流体介质始终处在阀门各部件所构成的自带流道内，不会进入内部的空腔，避免介质冲刷和内腔涡流，提高了阀门的流通性，有效避免泥沙及介质颗粒的沉积，不会造成阀门的卡涩、堵死。因此封闭结构流体介质流动阻力小，而开放式流道流体介质进人空腔，形成涡流，流体脉动压力大，阻滞严重，因此直接辐射流噪声声压级较大。4结语本文基于流噪声理论，对2 种不同流道结构通海双球阀在实际应用工况下进行流噪声数值模拟，得到结论如下：1）封闭流道增大转弯半径，增加弯角，虽增加了加工难度，但提高了阀门的流通能力。2）通过2 种结构不同速度频响曲线对比，可看出封闭与开放式2 种结构下的频响曲线总体变化一致，声压级均随着频率的升高呈下降趋势，最大峰值声压频率集中在2 0 0 40 0 Hz。3）通过对同一监测点不同速度的频响曲线的对比分析，表明各流速下开放式流道结构流噪声声压级明显高于封闭流道结构，且速度越大，压降越不明显，低流速时，双球阀低频特性更明显。4）通过对比分析总声压级的计算结果，初步揭示了总声压级的变化规律，即流道结构改变会导致监测点的声压级出现明显波动，且流场速度增加，监测点总声压级也呈增大趋势。5）与开放式流道结构对比，封闭流道结构降低约为2 6.9%的噪声声压，降幅为1 7 dB(A)，能够降低通海阀流道内流体脉动压力与流噪声声压，从而为新型通海阀的结构设计优化提供一种新的思路。参考文献：【1 】曾杏.钛合金通海阀结构及性能优化研究D.哈尔滨：哈尔滨工程大学,2 0 1 5.2郭涛.管路的流致振动及噪声研究D.武汉：华中科技大学,2012.543王开春，马洪林，赵凡，等.围壳形状对潜艇流致噪声的影响计算.空气动力学学报，2 0 1 8,36(5)：7 7 4-7 7 9.WANG Kaichun,MA Honglin,ZHAO Fan,et al.Numericalsimulation for effects of submarine fairwater shape on flow-in-duced noise.Acta Aerodynamica Sinica,2018,36(5):774779.【4方超，马士虎，蔡标华,等.基于BEM的通海阀流噪声与流激振动噪声数值模拟对比研究1 .舰船科学技术，2 0 1 8,40(3):3034.FANG Chao,MA Shi-hu,CAI Biao-hua.Numerical simula-tion contrastive study on flow noise and flow induced vibra-tion noise of sea valve based on BEMJ.Ship Science andTechnology,2018,40(3):30-34.5徐国栋，孙启，郑荣等.基于FEM/AML的船舶海水冷却系统出海管路流噪声预报.舰船科学技术，2 0 2 1,43(7):102-106XU Guo-dong,SUN Qi,ZHENG Rong flow noise predictionof sea water cooling system discharge pipe based onFEM/AML JJ.Ship Science and Technology,2021,43(7):102-1066桂瞬丰,幸福堂，李群燕，等.基于大涡模拟的通海阀噪声分析J.武汉科技大学学报(自然科学版),2 0 1 5(2)：1 2 9-1 33.GUI Shunfeng,XING Futang,LI Qunyan,et al.Noise analy-sis of a sea suction valve based on LESJ.Journal of WuhanUniversity of Science and Technology,2015(2):129-133.【7】聂欣，张童伟，朱羊羊，等.直筒笼式阀门流噪声数值模拟.热力发电,2 0 1 7,46(6):8 0-8 7.NIE X,ZHANG T W,ZHU Y Y,et al.Numerical simulationon flow noise of enterprises valveJ.Thermal Power Genera-tion,2017,46(6):8087.8 WEI L,ZHU G,QIAN J,et al.Numerical simulation of flow-induced noise in high pressure reducing valveJ.PLoS ONE,2015,10(6):e0129050.9 LIU K,QI S,et al.Flow-induced noise simulation using de-tached eddy simulation and the finite element acoustic analogymethodJ.Advances in Mechanical Engineering,2016,8(7):1-8.10 LAUNDER B E,SPALDING D B.Lectures in mathematical model of turbulenceJ.Academic Press,1972.11】徐，王德忠，王志敏，等.核电站主蒸汽隔离阀气流诱发振动与噪声的数值分析.原子能科学技术，2 0 1 0,44(1):48-53.XU Zheng,WANG De-zhong,WANG Zhi-min.et al Numeri-cal analysis for vibration and noise due to airflow of mainsteam isolation valve in nuclear power stationJ.Atomic Ener-舰船科学技术gy Science and Technology,2010,44(1):48-53.12】李树勋，赵子琴,张云龙，等.高温高压过热蒸汽疏水阀消声减振研究.振动与冲击,2 0 1 1,30(1 0):1 1 6-1 2 1.LI Shuxun,ZHAO Ziqin,ZHANG Yunlong.Noise elimina-tion and vibration reduction for a superheat steam trap withhigh temperature and high pressureJ.Journal of Vibration andShock,2011,30(10):116-121.13 LIGHTHILL M J.On sound generated aerodynamically.I.General TheoryJ.Proceedings of the Royal Society of Lon-don.Series A.Mathematical and Physical Science,1952,211(1107):564-587.14 LIGHTHILL.M J.On sound generated aerodynamically.II.turbulence as a source of soundJ.Proceedings of the RoyalSociety of London.Series A.Mathematical and Physical Sci-ence.1954,222(1148):1 3215 CURLE N.The influence of solid boundaries upon aerodynam-ic soundJ.Proceedings of the Royal Society A,1955,231(1187):505-514.16 FFOWCS WILLIAMS J.E,HAWKINGS D.L.Sound genera-tion by turbulence and surface in arbitrary motionJ J.Philo-sophical Transactions of the Royal Society of London,SeriesA,1969,264:321 342.17】王春旭,吴崇健，陈乐佳，等.流致噪声机理及预报方法研究综述.中国舰船研究,2 0 1 6,1 1(1):57-6 7.WANG Chun-xu,WU Chong-jian,CHEN Le-jia,et al.A com-prehensive review on the mechanism of flow-induced noiseand related predictions methodsJ.Chiese Journal of Ship Re-search,2016,11(1):57-67.18吕世金，俞孟萨，李东升.水下航行体水动力辐射噪声预报方法研究.水动力学研究与进展A辑，2 0 0 7,2 2(4):475-482.LYU S J,YU M S,LI D S.Prediction of hydrodynamic radia-tion noise of underwater vehicleJ.Journal of Hydrodynamics(Ser.A),2007,22(4):475482.19湛含辉，朱辉,陈津端，等.9 0 弯管内二次流(迪恩涡)的数值模拟.锅炉技术,2 0 1 0,41(4):1-5.ZHAN Hanhui,ZHU Hui,CHEN Jinduan,et al.Numericalsimulation of secondary flow(Dean vortices)in 90curvedtube.Boiler Technology,2010,41(4):1-5.20李树勋，康云星，孟令旗等，多级降压疏水调节阀流致噪声数值模拟研究J.振动与冲击，2 0 2 0,39(1 4)：1 1 6-1 2 1.LI Shuxun,KANG Yunxing,MENG Lingqi,et al.Flow-in-duced noise of a high pressure drop control valveJ.Journal ofVibration and Shock,2020,39(14):116-121.第45卷