莱特希尔声类比

分类:耦合分析·声学耦合 | 2026-10-01 改写

理论:把流动方程改写成声方程

概述

🙋

喷气式飞机为什么那么吵?

🎓

因为高速射流与周围空气混合时的湍流在发声。莱特希尔于1952年把流动方程变形,左边改写为静止空气中声波的方程,右边为湍流产生的声源项。右边的声源是称为四极子的形式,做量纲分析可知发出的声功率与速度的八次方成正比。速度加倍声音变为256倍,即大24 dB。所以射流速度只要稍微降低,噪声就大幅下降。现代喷气发动机用大风扇让大量空气缓慢流动产生推力,降低射流速度以变得安静。

莱特希尔方程

$$ \frac{\partial^2\rho'}{\partial t^2} - c_0^2\nabla^2\rho' = \frac{\partial^2 T_{ij}}{\partial x_i\,\partial x_j},\qquad T_{ij} = \rho u_i u_j + (p' - c_0^2\rho')\delta_{ij} - \tau_{ij} $$

ρ′为密度脉动,c₀为静止空气声速,T_ij为莱特希尔应力张量(由雷诺应力、熵脉动、粘性应力组成)。右边的二阶导数表示四极子声源。

八次方定律与声效率

$$ W = K\,\frac{\rho_0 U^8 D^2}{c_0^5},\quad K \approx 3\times10^{-5},\qquad L_W = 10\log_{10}\frac{W}{10^{-12}\,\mathrm{W}},\qquad \eta = \frac{W}{\tfrac12\rho_0 U^3 A} \propto M^5 $$

K为实验确定的系数,D为射流直径,A为射流截面积,M为射流马赫数。适用于亚声速冷射流的近似。

Coffee Break 随笔闲谈

喷气客机时代之初

莱特希尔发表声类比的1952年,也是世界第一架喷气客机彗星号投入运营的年份。喷气发动机噪声在机场周边成为大问题,八次方定律表明降低射流速度是最有效的对策。这后来发展为用风扇输送大量空气的高涵道比涡扇发动机,同样推力下噪声大幅下降。声类比至今仍是用CFD预测噪声的方法基础。

算例

例1:直径0.5 m射流的声功率

射流速度马赫数声功率声功率级射流动能声效率100 m处声压级
100 m/s0.290.02 W102.8 dB0.12 MW1.6×10⁻⁷51.8 dB
200 m/s0.584.9 W126.9 dB0.94 MW5.1×10⁻⁶75.9 dB
300 m/s0.87124 W140.9 dB3.18 MW3.9×10⁻⁵90.0 dB

100 m处声压级是假设向所有方向均匀扩散的值(实际射流噪声在后方斜向最强)。速度降低一成为−3.7 dB,加倍为+24 dB。

例2:计算声音的网格(300 m/s)

斯特劳哈尔数频率波长网格(每波长10格)
0.2120 Hz2.86 m286 mm
1.0600 Hz0.57 m57 mm
🙋

流动能量只有0.004%变成声音,却这么响。

🎓

射流携带的能量高达3 MW,所以0.004%也有124 W。人耳能听到弱到10⁻¹² W/m²的声音,124 W是非常大的声音。这种微小比例也使声音预测变难。声音的脉动远小于流动脉动,用同一计算直接求流动和声音,声音会淹没在数值误差里。因此常用分工的混合方法:只对湍流部分用LES等精细计算求声源,向远处的传播用莱特希尔方程或FW-H法等积分计算。流动计算也需要能分辨关心的最高频率声源的网格和时间步长。

噪声预测流程

  1. 用八次方定律等简单公式把握噪声大小和速度效应。
  2. 根据关心的频率确定网格大小和时间步长。
  3. 用LES等非定常计算求解湍流区域,求声源。
  4. 用FW-H法或莱特希尔积分计算远处观测点的声音。
  5. 与实测噪声按频率、方向比较。
Coffee Break 随笔闲谈

“把排气速度降低一成”

某工厂排气管道出口气流噪声在厂界处略微超标。也有加装消声器的方案,但先用八次方定律估算,出口速度降低一成噪声约降3.7 dB,就能达标。把出口截面扩大约一成降低速度后,实测噪声几乎与计算一样下降了。这是在采取大规模对策之前估算速度效应发挥作用的例子。

常见错误

错误与对策

错误后果对策
在同一网格上直接求解流动和声音声音淹没在误差中混合方法
用RANS求声源得不到脉动用LES等非定常
网格粗缺少高频按频率定网格
忽略指向性观测点声音出错按方向评估
超声速用八次方定律高估确认适用范围
🙋

我想了解相关内容。

关联模拟器

在本领域的交互式模拟器中体验理论

模拟器列表

相关领域

结构分析电磁界分析热分析
本文评价
谢谢回答!
有帮助
需要更
详细
错误
报告
有帮助
0
需要更详细
0
错误报告
0