设计仿真 | 新型风扇气动噪声组合分析方法

首先进行非稳态CFD仿真进行。计算域由静止场和包围风扇的旋转域组成，如下图1所示。该模型捕捉了所有细节，例如风扇孔口、风扇护罩、线圈和百叶窗。在本研究中，线圈被建模为多孔区域，并应用滑动网格方法来计算Actran气动声学模拟所需的非定常CFD结果。旋转域（风扇）的旋转频率为1118RPM。例子中的时间步长为0.0005s。此次模拟，0.8秒的总时间确保所求最小频率远小于37.2Hz（叶片通过频率）。

图1：旋转域包围风扇叶片，静止域包含流动障碍物和多孔线圈

Lighthill声类比有两种源的计算方式，分别是在域的整个体积上和在源的表面上执行源的计算。如果采用前一种方法，源计算需要整个域上的CFD信息，但在后一种方法中，我们只需要在单个表面而不是体积上读取速度信息（以及不可压缩模拟情况下的密度），从文件管理的角度来看，这是一个很大的优势。本研究中风扇作为唯一的主要声源，为了加速CFD模拟仅导出包围风扇的表面，即转子-定子界面静态侧的CFD数据。CFD求解器采用480个核心的仿真时间接近40小时。最耗时的部分是在每个时间步下写入数据这个步骤，这会减慢求解速度。

本例基于Actran2021.1，采用Lighthill面源方法，将声源映射到声学网格上，并完成时域气动声源转成频域的计算。这项研究的新颖之处在于这一步，其中采用两个互补的离散傅里叶变换（DFT）设置来精确计算线谱音调和宽频噪声，同时避免了由于采样时间有限而在高频下出现不切实际的声压级波动：

① 对于叶片通过频率(BPF)及其谐波引起的线谱音调噪声，使用最小二乘法在整个采样时间内定义并完成第一个DFT。该方法强制提取用户设置的频率。在本研究中，设置为从BPF开始并提取BPF高达1000Hz的每个谐波；

② 对于宽频带噪声，时域数据样本分解为多个较小的样本（多重离散傅里叶变换），这些样本彼此重叠50%，并对每个样本应用Hanning窗以平滑每个子样本之间的过渡。在本研究中，原始时域数据被分为78个重叠样本，每个样本的持续时间为0.02秒，正好对应50个CFD时间步长。子采样持续时间为0.02s，最小频率和频率步长固定为50Hz。选择这些参数是为了获得不包含BPF或其谐波的频率列表，因此仅提取宽频带噪声的结果。此方法的缺点是没有对CFD的全部结果进行利用。

通过上面的两个步骤，便获得了两组气动噪声：

①一个DFT得到BPF及其谐波的噪声。

②多重DFT方法得到的宽频带噪声。使用两组互补的结果，并使用脚本进行合并，就可以获得组合气动声学仿真的总体频率响应，如图2所示。

图2：音调和宽频带结果相结合以获得模拟的频率响应

图3：接收器相对于装置的放置示意图

图4：使用Lighthill类比的传统混合方法（蓝色）、新组合方法（红色）与测试方法（绿色）的SPL比较图

查看完整原文请搜索：Ali Nikparto, Alexis Talbot, Asad Sardar, Clement Roy. Outdoor unit split fan noise prediction using a hybrid lighthill analogy method. NOISE-CON. 2022: 727-736.