《Journal of Fluid Mechanics》:Energy dissipation mechanisms in an acoustically driven slit
编辑推荐:
研究人员量化了平面波通过大长径比狭缝几何时,入射声能转化为涡旋运动和粘性耗散的过程。研究人员在入射声压级(ISPL)、斯特劳哈尔数(St)和雷诺数(Re)构成的宽广施加参数空间内进行了直接数值模拟(DNS)。谱本征正交分解(SPOD)在每个频率下产生能量排序的
研究人员量化了平面波通过大长径比狭缝几何时,入射声能转化为涡旋运动和粘性耗散的过程。研究人员在入射声压级(ISPL)、斯特劳哈尔数(St)和雷诺数(Re)构成的宽广施加参数空间内进行了直接数值模拟(DNS)。谱本征正交分解(SPOD)在每个频率下产生能量排序的相干结构,研究人员据此构造了谱动能(KE)和粘性损失(VL)的逐模态场,以研究声吸收机制。在ISPL=150 dB时,当St≤4St0时,声-水动力能量转换最高,对应于有效Keulegan–Carpenter数(Kc)大于40。在该区域,三维模拟显示主导流动响应是二维的,狭缝角附近的振荡剪切层卷起形成脱落涡。VL占KE贡献的20%–60%。对于更大的St,Stokes层的约束产生X形近缝模态,使能量输入减少约50%。Re的影响取决于振幅。在150 dB时,较大的Re对应抑制宽带波动并锐化谐波峰。在ISPL=120 dB时,边界层保持附着,涡脱落较弱,吸收随粘度单调变化,且Re和St依赖性变得相当。在所有条件下,超过99%的VL集中在狭缝口附近的紧凑区域。KE–VL谱识别了狭缝几何中增强或抑制声阻尼的区域,为基于声学的设计提供了物理上可解释的基础。
声波与流体在受限几何中的相互作用是流体力学的基础问题,对声阻尼、噪声控制和声衬设计具有重要意义。狭缝谐振器因其简单且蕴含丰富物理而成为典型模型。然而,以往研究多将声波通过狭缝时的阻尼与背腔共振响应耦合,很少将其作为独立机制考察。为此,研究人员在无背腔、末端消声的二维大长径比狭缝中开展直接数值模拟(DNS),聚焦狭缝口局部的声能耗散机制。研究基于入射声压级(ISPL)、斯特劳哈尔数(St)和雷诺数(Re)三个无量纲参数,覆盖从线性粘性吸收到强非线性涡脱落的宽广工况,旨在阐明声能转化为涡动能(KE)并最终经粘性损失(VL)耗散的过程。该论文发表在《Journal of Fluid Mechanics》。
研究方法方面,研究人员使用GPU加速的可压缩流求解器MFC,求解可压缩Navier–Stokes方程,采用五阶加权本质无振荡(WENO)格式进行空间离散,以总变差递减三阶Runge–Kutta格式推进时间,并用浸没边界法处理壁面。数据库包含42个二维算例,涵盖2个ISPL(120和150 dB)、7个频率(对应St/St
0=1至12)和3个粘度(对应Re=Re
0、Re
0/2、Re
0/3);另在特定工况进行展向长度为11d和22d的三维模拟,验证流动的二维特性。分析采用谱本征正交分解(SPOD),在每个频率提取能量排序的相干结构,进而构造逐模态的谱KE和VL场。
研究结果如下。在3.1 设置(Set-up)中,研究人员构建的参数空间覆盖低马赫数(Ma≤0.05)下的层流到弱非线性流动,所有算例均保持声波传播所需的压缩性。在3.2 数据(Data)中,瞬时涡量场显示边界层在狭缝角卷起、周期性脱落,高频时涡结构被限制在狭缝附近;VL场在狭缝口最大,且120 dB下VL比150 dB低约三个量级。在3.3 施加控制参数与诱导狭缝响应之间的映射(Mapping between imposed control parameters and induced in-slit response)中,研究人员定义了基于狭缝峰值速度的狭缝St和Re,发现150 dB下三条Re曲线可塌缩为单一函数,表明高振幅吸收由狭缝St决定;120 dB下则保持Re依赖性,粘性效应主导。在3.4 声吸收(Acoustic absorption)中,功率吸收系数显示高ISPL下St为主控参数,低ISPL下Re与St作用相当;反射系数随St单调增加,透射系数减小,狭缝阻抗呈现非线性增长。
在4.1 谱本征正交分解(Spectral proper orthogonal decomposition)中,SPOD将流场分解为能量排序模态,为后续KE–VL分析提供框架。在4.2 代表性SPOD特征值谱和模态(Representative SPOD eigenvalue spectrum and modes)中,150 dB、St=4St
0、Re=Re
0工况的特征值谱在基频及谐波处出现峰值,领先模态占95%能量,空间结构集中于狭缝四角,反映声能向相干涡运动转换。在4.3 参数空间中的谱特征(Spectral features across the parameter space)中,较低Re抑制宽带波动并锐化谐波峰;高频时出现非整数倍频峰,表明存在非线性三波相互作用;高St下模态变为X形,边界层变薄。在4.4 三维效应(Three-dimensional effects)中,2D与3D模拟的吸收系数相差约11%,展向波数分解显示能量几乎全在k
z=0分量,证明主导流动实质为二维。在4.5 耗散机制(Dissipation mechanism)中,150 dB下VL占KE的20%–60%;速度最大梯度与VL场空间对齐,确认边界层动力学是粘性耗散主因;所有工况超过99%的VL集中在狭缝口附近区域。在4.6 ISPL依赖性(The ISPL dependence)中,从120 dB升至150 dB时,狭缝Re从O(10
1–10
2)增至O(10
3),主导模态从局部受限变为大尺度涡脱落,KE成为主要耗散通道,宽带分量随振幅增大而增强。
讨论部分指出,150 dB下总VL在St=4St
0附近最大,随后在St≥6St
0时下降约50%,该行为可通过Keulegan–Carpenter数(K
c)解释:K
c>40对应强涡卷起,促进剪切层耗散;高St下出现非谐波峰,提示三波相互作用转移能量。研究结论可翻译为:本研究建立了严格的数值和谱分析框架,量化了薄型大长径比狭缝中入射声能向相干动能和粘性损失的转换。二维DNS覆盖宽参数空间,三维模拟验证了主要物理可由二维模型捕获。SPOD揭示了声–KE–VL交换由St、Re和ISPL耦合控制;所有条件下超过99%的VL来自近狭缝区域,领先SPOD模态占基频能量的95%以上。在150 dB、St≤4St
0、K
c≥40时,振荡边界层卷起大尺度涡,声能转换最强;St>8St
0时模态受限,耗散降低约50%。这些结果为声衬和噪声控制设计提供了物理可解释的依据。