海绵层模型

海绵层模型提供了在直接噪声模拟中抑制外部边界附近的虚假反射干扰的方法。

直接噪声模拟求解可压缩流体流中的时间相关流变化。这些流变化将与流进行对流(对于涡干扰),或者以声速向外传播到流域的外边界(对于声干扰)。对于高质量气动声学模拟,数值算法必须允许这些干扰平滑地穿过外部边界,而不在内部反射。

某些边界(如压力出口边界)提供非反射条件。但是,非反射边界条件的效率取决于干扰到达边界的倾角 - 倾角切向分量越大,反射抑制效率越低。

为了克服这种与倾角的相关性,Simcenter STAR-CCM+ 提供了海绵层模型,可以通过抑制指定厚度层(即海绵层)内的非稳态流变化实现边界处的非反射条件。该体积方法可防止声波和涡干扰到达边界,从而消除虚假反射。

海绵层模型将以下阻尼项添加到连续性、动量和能量方程的右侧:

1. EQUATION_DISPLAY
S s l = - σ s l ( ϕ - ϕ r e f )
(4744)

其中:

  • ϕ { ρ , ρ v , ρ E } 表示已传输的量。 ρ 为密度, v 为速度矢量, E 为单位质量的总能量。
  • σ s l 海绵层阻尼系数
  • ϕ r e f 参考温度

海绵层阻尼系数

海绵层阻尼系数 σ s l 用于海绵区域中阻尼量的平滑过渡。该系数计算为边界距离 d 的函数:

2. EQUATION_DISPLAY
σ s l ( d ) = σ max { 1 cos [ 2 π ( 1 d w ) ] }
(4745)

其中 w 为海绵层厚度,为模型系数

σ max 为计算的最大阻尼系数:

3. EQUATION_DISPLAY
σ max = c w log ( 1 α 100 )
(4746)

其中 α 是相对抑制量,为模型系数

参考求解

参考求解 ϕ r e f 是模型阻尼海绵区域中波动所参考的求解。参考求解基于移动时间平均,以显式形式给出,如下所示:

4. EQUATION_DISPLAY
ϕ r e f n + 1 = ϕ r e f n + Δ t Δ T ( ϕ ϕ r e f n )
(4747)

其中

  • n + 1 为当前时间级别, n 为上一个时间级别。
  • Δ t 为时间步长。

Δ T 为局部时间过滤器宽度,计算如下:

5. EQUATION_DISPLAY
Δ T = w c
(4748)

其中 c 为声速。

模型系数

α w
90.0 1.0 m