ansys电机边界条件（Fluent电子风扇效能及噪音仿真攻略）

ansys电机边界条件（Fluent电子风扇效能及噪音仿真攻略）通过Fluent Mesh制作Interface第三阶段瞬时计算(Transient Flow)：ANSYS Fluent做噪声仿真属于瞬态仿真(Transient)流程(细节后面将描述)，需注意下列工况及流程：第一阶段：稳态计算(Steady Flow)：第二阶段瞬时计算(Transient Flow)：從MRF改成Mesh Motion

导读：大家好，我是阿朱Simulation，刚刚成为了一名仿真秀专栏作者，这也是我第一篇仿真科普文章。本文介绍如何利用ANSYS Fluent进行电子风扇效能及噪音的仿真，关键技术已显示在流程中，按照本案例流程，可顺利且有效进行风扇流体仿真，如有不当欢迎大家批评指正。

电子产品中风扇占据一个重要的散热角色，尤其电子风扇之可靠度要求倍增，消费者无法接受经常性地更换风扇，且电子风扇供货商也无法承受此成本的提升，更无法承受市场上面临淘汰的困境，所以对电子风扇的设计往往将通过仿真进行前期开发。

本案例模型藉由产品的要求及自身未来的发展，进行了仿真分析与验证。主要着重在单体风扇设计，从叶形设计到CFD的性能特征预测，及后段的噪音评估等，都是产品设计的指标。本文内容包含了风扇性能曲线及噪音的分析流程。

一、电子风扇性能曲线分析

本例模型采用如上Blower风扇。

ANSYS Fluent做噪声仿真属于瞬态仿真(Transient)流程(细节后面将描述)，需注意下列工况及流程：

• 转速：3600 rpm
• 网格数：1 662 569个

第一阶段：稳态计算(Steady Flow)：

• k-omega SST紊流模式
• 计算1000 iterations
• Broadband Noise Model

第二阶段瞬时计算(Transient Flow)：從MRF改成Mesh Motion

• LES紊流模式(Large Eddy Simulation)
• 计算500 time step (60 iterations/ time step )到动态稳定

第三阶段瞬时计算(Transient Flow)：

• LES紊流模式(Large Eddy Simulation)
• 计算500 time step (60 iterations/ time step )撷取声压随时间之变化数据
• FW-H Noise Model

通过Fluent Mesh制作Interface

若用户想采流程化Fluent Mesh来建立网格，因New Fluent Mesh还不直接支持Interface的边界条件，因此学习者可把MRF旋转流体区及周围其它计算域分两次制作网格(耗费时间差不多)，之后再于Fluent中进行Domain/Append/Append Case File(可读入*.msh)。

要建立一个可与风洞相比拟的环境，我们首先就要将模型处理成风洞的实验状态，如下图。我们把风扇出口的风道管给建立出，再将周围适当的空气域包覆整个主要器件。

接着将MRF 旋转区及周围空气区做Mesh Interfaces，将共面位置给连接起。

第一阶段计算稳态：计算噪声不应直接从瞬态进行，因在瞬态抽取噪声数据前，要达稳定需要一段时间；因此我们要采先进行稳态计算，可快速先达到稳态平衡。默认动量计算采k-omega SST。

空气密度部份，若气动速度非常快，例如转速8000RPM或甚至超过上万RPM等，我们要考虑采Incompressible-idea-gas，其余工况下可采Constant。

转速及转轴设置部份，在旋转流体区勾选Frame Motion，设置Rotational Velocity=3600 rev/min(RPM)；Rotation-Axis Direction设置的是风扇转动方向，采右手定则，姆指为轴心，其余四指为叶片转动方向。

边界条件设置，模型有两个主要边界，一个是风道出口，一个是周围环境。实务上，在风道出口是接入风洞腔体中，由风洞来控制出口压力；因此在Fluent软件也是类似操作，入口是自由一大气压边界，而关键是控制出口的压力变化。出口压力上升，通过风扇的流量就会降低，便可形成风扇性能区线(P-Q Curve)。

第二阶段动态稳定平衡：此阶段涉及到噪声，因此风扇必须转动，我们由此切换设置为Mesh Motion。因要计算噪声涡损，改以瞬态模型做计算，并切换为LES紊流模式(Large Eddy Simulation)，当达到动态压力持续稳定后，才进入第三阶段。

二、电子风扇噪音的分析流程

第三阶段FWH计算噪声：计算噪声前，我们要知道所要计算噪声的频率到什么范围。人耳一般能听到的频率在20 000Hz，而高频率如蝙蝠所能听到的频率落在100 000Hz；一般采用无响室测试规范，频率范围约采在20~50 000Hz。

CFD中仿真移动件的压力与速度波动通常是近场(near field)噪音，而近场位置通常会最大的噪音发生之处，大多数CFD应用中只要将计算域延伸将可进行近场分析与预测，远场距离对象一段距离，通常会在CFD计算域以外，而近场的压力与速度波动为决定远场噪音的重要因素，这些波动将成为远场噪音预测的声源。

ANSYS Fluent要计算完整噪声采用FWH Ffowcs Williams & Hawkings model

通过FWH模型计算噪声需要特别注意网格密度及质量，我们可以知道一个声波的波长尺度，从声波波长大小要去控制网格密度。从理论来探讨，一般最少在一个波长中至少要有5个网格点，计算越高频噪声所需要的网格点通常一个波长范围会高达20个，这也是造成噪声计算量特别大的一个原因；同样地，若不精准控制波长区间网格数，准确性就会明显下降。

启动FWH噪声模型：

定义噪声源为转动壁面(Define Sources)：

接着定义要放Speaker的位置，也就是要监控噪声的位置，这里建立两个噪声监控位置(Define Receivers)，在风扇本体上下个0.1m处。

定义时间步长。我们将要计算频率达20 000Hz的噪声问题，而我们知道频率为时间的倒数，因此一个声波行进距离为1/20 000Hz = 0.00005s，若要精准计算每个声波波长上网格点的数据，则至少要除以5，也就是DT = 0.00001s。实务上更复杂的问题，时间步长可能会更小。

由于前面两个阶段我们已达成动态平衡，接下来风扇每转一圈对比于每转更多圈噪成的声压结果，理论上会一致，因此学习者只要确保你所要监控的位置在风扇转动一圈后，声压能够传递到即可。若监控的位置过远，你就需要让风扇转动更多圈，让波形能够传递到该位置。

3600RPM → 60 rev./s → 1/60 s/rev. = 0.016666s (转动一圈需要的秒数)

Number of Time Steps = 0.016666s/0.00001s = 1667

运行完后的残差如下：

通过FFT转换

• 上方0.1m处噪音为 41.22 dB
• 下方0.1m处噪音为 39.87 dB

以上内容就是本次分享的全部内容。