汽车空气动力学

高速行驶时,气动噪声主导发动机、排气和轮胎发出的噪声。

气动噪声来自经过风窗玻璃立柱、后视镜和车底的气流以及穿过管道的流体(包括风扇噪声)。 本节介绍以下推荐的汽车空气动力学建模方法:

推荐的过程

  1. 对于基于 RANS 的稳态空气动力学情况,应用网格分辨率准则。 确保平均流正确至关重要。 仅当车尾气动声学无关紧要时,才会折衷考虑尾流区网格加密。
  2. 绘制网格截止频率,来判断声源区域中的加密是否求解 3-5 kHz 区域中的频率。 实际上,需要尺寸小于或等于 2 mm 的网格单元。 适当地加密并将求解映射到加密的网格。
    可以从声学相关性较大的区域中标识声源区域。 将 Proudman 用于体积源:

    将 Curle 用于表面源:

  3. 确保近场中的网格足以局部传播声音。 例如,在 A 柱、侧视镜和侧窗玻璃区域中使用小于或等于 3.5 mm 的网格单元传播 5 kHz 信号。
  4. 选择适用于以下约束中要求最苛刻(即,最小)的约束的时间步长:频率求解、显式 CFL 或隐式 CFL(具体视情况而定)和库朗数。
  5. 从用于判断先前要求的初始稳态模拟重启瞬态工况。 将湍流模型从 RANS 切换到 DES 或 LES。 使用二阶时间离散和基于中心差分的空间离散。
  6. 选择并微调数值参数,确保收敛在每个时间步内。 满足以下所有(而非部分)条件时,可以推断求解已收敛在时间步内:
    1. 在时间步结束时单调减少或最小化残差。
    2. 监视器保持不变。 监视器可以为点压力、升力或曳力。

推荐的网格设置

根据条件,选择网格单元尺寸 Δ

条件

注释/描述

公式

典型值

平均流

使用外部稳态空气动力学准则。 加密近壁区域以捕捉体分离。

Δmin= 2 mm

y+ = 1(对于低雷诺数壁面处理)

湍流积分尺度,Lintegral

度量耗散涡的“典型”尺寸。 从稳态近似,Δ0.5Lintegral

Lintegral=Cμ0.75k1.5ϵ

Lintegral > 10 mm

Δ2 mm

与湍流结构关联的频率

网格截止频率,基于相对准则中推荐的网格尺寸的局部湍流波动。 从稳态近似。

FMCF=23k2Δ

Δ= 4 mm 可捕捉 1000-1500 Hz 范围内的频率。

压力波传播

基于二阶空间和时间离散,每声学波长 λ 20 个网格单元(高阶空间格式更少)。 a 为声速。

λ=aF
Δλ20

Δ3.5 mm。基于在空气中传播 F = 5000 Hz 压力波的要求。

推荐的时间步

根据条件,选择时间步长 Δt

条件

注释/描述

公式

典型值

F 频信号

奈奎斯特采样定理表明,频率为 F 的信号可通过以大于 2F 的速率采样来分辨。 对于二阶空间格式,测试结果表明求解波幅值需要约 15 个点。

ΔtF115F,

推荐

ΔtF = 2.22E-5 s。

基于捕捉要求

F = 3000 Hz。

对流和声学

基于对流速度 U 和声速(声压波以局部声速 a 传播)。

库朗 (CFL) 数基于 (U+a)

显式:

Δtexp=Δ(U+a)

隐式:

Δtimp=CFLΔtexp

Δtexp = 5.0E-6 s。

基于

Δ = 2.0 mm,

U = 50 m/s,

a = 350 m/s。

Δtimp = 5.0E-5 s。

基于 CFL = 10。

仅对流

基于对流速度 U

对流库朗数仅基于 U

ΔtCourant=ΔU

ΔtCourant = 4.0E-5 s。

基于

Δ = 2.0 mm,

U = 50 m/s。

扩散时间尺度

对于局部(涡流)粘度 Veff 最大的边界层最为重要。

对声学现象无关紧要

Δtdiff=Δnormal2Veff

Δtdiff = 1.0E-6 s。

基于

Δnormal= 5.0E-6 m,

Veff = 2.0E-5 m2/s。

推荐的求解器设置

可以为慢速(“常规”)网格和时间步设置、快速(“主动”)设置或默认设置选择求解器设置。

速度 URF

压力 URF

每个时间步的迭代次数

默认

0.8

0.2

5

常规

Δt = 1.0E-4 到 1.0E-5 s。

Δmin = 1-4 mm。

0.7

0.7

4-10

主动

Δt = 2.0E-5 s。

Δmin = 2 mm。

1.0

0.9

3