《Ocean Engineering》：A numerical study of waterborne noise propagation in air across a thin solid layer
编辑推荐：
研究人员建立了一个基于卷积完美匹配层（Convolutional Perfectly Matched Layer, CPML）吸收边界与二阶精度有限体积离散（Finite Volume Method, FVM）的二维/三维声波方程数值求解框架，用以模拟由流体动
研究人员建立了一个基于卷积完美匹配层（Convolutional Perfectly Matched Layer, CPML）吸收边界与二阶精度有限体积离散（Finite Volume Method, FVM）的二维/三维声波方程数值求解框架，用以模拟由流体动力湍流场（通过Lighthill类比给出源项）激发的声压波在具有不同密度和声速比的多层介质（重流体–固体–轻流体，如水–钢/玻璃钢–空气）中的传播过程。文中引入非维数化分析识别出控制问题的一组无量纲数（马赫数Ma、阻抗比ρc、几何比L/δ等），并在界面处采用二阶后向差分格式的连续性条件求解交界面压力与通量。通过与单层无界域中点源解析解及平面波在两层界面的反射/透射系数理论对比验证了算法的二阶时空精度与CPML衰减效果。随后将该方法分别用于单频点声源（monopole）与壁分辨大涡模拟（Wall-Resolved Large Eddy Simulation, WR-LES）平面槽道流提供的宽带Lighthill源项，考察不同固体层厚度及材料（钢、玻璃钢）下声能在空气侧的透射特性，并以近似弹性波散射矩阵估算等效流体模型引起的S波缺失误差并予以修正。
论文解读：《Ocean Engineering》刊载 "Numerical simulation of hydroacoustic noise propagation across layered media with fluid-dynamic source coupling via CPML-absorbing boundaries and finite volume discretization"
一、研究背景与意义
水下结构（如船体、海洋平台）内部或外部湍流边界层、机械振动等流激噪声（flow-induced noise / hydroacoustic noise）可经固体围护结构传递至空气中，涉及水–固–气三层介质间巨大的声阻抗失配（acoustic impedance mismatch, Z=ρc）。现有解析解仅适用于半无限空间或单一界面，对含薄板中间层的多层系统尤其是由宽带湍流源激励的情形缺乏普适解，且全弹性波动（Navier–Lamé方程）计算代价高昂。因此有必要发展一种兼顾精度与效率的分层声传播数值模型，量化材料属性与几何参数对跨介质声透射的影响，为舰船噪声评估与隔声设计提供依据。
本文研究人员提出将固体层简化为只传纵波（P-wave）的等效流体（equivalent fluid approximation），结合卷积完美匹配层（CPML）处理外场辐射边界，用有限体积法（FVM）在结构化网格上求解含Lighthill源项的波动方程，并通过无量纲分析提炼控制参数，系统考察点声源与湍流槽道流两种典型激励下钢板与玻璃钢（fiberglass）夹层对声能向空气侧传输的滤波作用，并以弹性波散射矩阵近似修正S波缺省带来的偏差。
二、主要技术方法概述
研究人员采用以下关键技术路线开展研究：(1) 在二阶线性化欧拉方程基础上导出常密度介质中含CPML辅助变量(ψi, ζi)与阻尼因子(σi, αi)的修正声波方程（Acoustic Wave Equation with CPML），时间推进用二阶中心差分显式格式，空间离散用有限体积法（FVM）在六面体单元上积分，交界面压力连续性(p1=p2)与法向速度连续(1/ρ1?np=1/ρ2?np)用二阶后向/向前差分求解界面压力pint；(2) CPML参数按Collino & Tsogka (2001)、Komatitsch & Martin (2007)设置，其中σi∝(xi/Li)2ln(1/R)，R=10?7，αi=πf0(1?xi/Li)；(3) 无维数化引入特征长度δ（源尺度或λ=c/f）、流速尺度uf，得控制无量纲数组：Ma=uf/c0、ρ2/ρ1、c2/c1、L1,2,3/δ等；(4) 验证阶段对比无界域单极子解析解pm(r,t)=?ikρcQmei(ωt?kr)/(4πr)及水–气平面波R/T系数；(5) 分层算例设水(c=1500 m/s, ρ=1000 kg/m3)–钢(c=5900 m/s, ρ=7800 kg/m3)或玻璃钢(c=3000 m/s, ρ=2500 kg/m3)–空气(c=300 m/s, ρ=1 kg/m3)，点源频率25–2000 Hz，湍流源来自Reτ=395平面槽道流WR-LES（WALE亚格子模型，OpenFOAM PISO算法），Lighthill张Tij=ρuiuj?σij+[(p?p0)?c02(ρ?ρ0)]δij在低Ma下仅保留ρuiuj项；(6) 对等效流体模型忽略横波(S-wave)的误差用入射/透射能量流角积分与散射矩阵估计整体误差εsteel≈33%，εfiberglass≈20%，修正prms,c=√((1?εi)prms2)。
三、研究结果
3.1 Unbounded and homogeneous domain（无界均匀域验证）
采用不同CPML厚度(LCPML/λ=0.10–0.20)对水中50 Hz单极子进行模拟，对比解析解得到相对误差<1%，加倍CPML厚度仅降低误差约1%，最终选用LCPML/λ=0.133为折衷。非均匀网格拉伸测试误差亦低于均匀网格情形，证实CPML与FVM联合实施有效，算法具二阶时空收敛精度（通过MMS法验证）。
3.2 Two-layer domain（两层介质界面验证）
在水–气平面波垂直入射情形下，数值解与理论反射/透射系数吻合良好（误差<1%），证明交界面二阶差分格式能正确复现压力连续与法向质点速度连续条件，为三层模拟奠定基础。
5.1 Point-like source in layered media（分层介质中点声源传播）
对水中单频点源穿越钢/玻璃钢薄板进入空气的系列模拟表明：(a) 固体–空气界面因巨大阻抗失配引起声压骤跌落；(b) 随频率升高（即L2/δ增大）透射声压下降，钢因更大Z失配较玻璃钢衰减更显著；(c) 能量流透射系数TI在δ/L2≈93.75附近急剧下降，存在临界板厚使长波"透明"区终结；(d) 经S波修正后的透射强度钢较玻璃钢低约一个数量级，说明材料选择显著影响跨介质传声。
5.2 Flow-induced noise in layered media（流激噪声在分层介质中传播）
5.2.1 Acoustic similarity and the prototype case（声学相似律与原型缩放）
通过无量纲分析建立模型尺度（h=0.17 m, um=1 m/s, νm=2.45×10?5m2/s, cm=150 m/s, ρm=100 kg/m3 "metafluid"）复现实船Ma≈6.7×10?3与阻抗比，允许Re偏低但声学无量纲数组一致。
5.2.2 Flow-dynamic source（流体动力声源）
Reτ=395槽道流WR-LES给出Lighthill项时间谱，其呈低频平台–中频隆起(f*≈40–70, ～ω2放大)–高频衰减形态，源于?i?j(uiuj)在波数域对应kikj乘速度谱卷积。均方根Lighthill源在缓冲层(y+≈10–30)最强，近壁是主要发声区。
5.2.3 Acoustic analysis（声学分析结果）
声压谱在水侧(a1)呈?7/3（惯性子区）与?11/3（剪切相互作用区）两段幂律衰减；在空气侧(a2,a3)仅保留?7/3段，?11/3段消失且高频(f*>100–110)被大幅抑制，证实固体层具低通滤波效应。板增厚使水侧反射增强（水侧谱幅上升）而空气侧透射压降低。同厚度下钢较玻璃钢透射压更低，与单频点源结论一致，说明阻抗失配主导跨介质传声强弱，材料属性调制滤波幅度。
四、讨论与结论翻译
研究表明所建CPML–FVM分层声传播模型经解析解校验具二阶精度，可有效模拟流激噪声跨水–固–气介质的传递。固体层起低通滤波器作用，透射强依赖于板厚与材料阻抗，钢因更大ρc失配较玻璃钢衰减更强，临界板厚对应δ/L2≈93.75处透射强度陡降。湍流边界层噪声主要源自近壁区，频谱保留?7/3惯性子区衰减而在透过固体层后丢失?11/3剪切相互作用区成分。提出的a posteriori S波修正因子使等效流体模型结果趋近全弹性解。该框架为水下结构辐射噪声评估及隔声层设计提供了计算可行、物理合理的数值工具。
Concluding Remarks（结论原文意译）： 本研究探讨了噪声在多介质层中的产生与传播机制，常见于船舶水下线源经结构传入大气的工程问题。研究人员开发并验证了基于多维有限体积法与CPML吸收边界的声波方程数值方法，建立了识别跨介质传声控制参数的无量纲框架，并将其用于水体中单频点源与湍流流场宽带源经薄固体层向空气传播的模拟，比较了钢与玻璃钢两种固体介质，引入对固体中横波未模拟的a posteriori修正。点源情形下固体层呈频率选择性滤波器，透射依赖频率与材料属性，钢因更大阻抗失配在高频段衰减更强，玻璃钢透射略高且频率选择性弱；强度透射系数存在临界厚度使透射骤降，两材料临界值相近。槽道流情形确认近壁为主要发声区，固体层具低通滤波性，高频成分被抑制，钢与玻璃钢透射差异与点源一致。所提数值方法可稳健分析强阻抗失配下分层介质水声噪声问题，后续可拓展至变Ma数、更多介质层及不同噪声源情形，或与Navier–Lamé方程耦合用于振–声（vibro-acoustic）问题。
