Fluent仿真实例-大涡模拟大风吹过圆柱体的噪声
流体流过圆柱体产生的噪声案例描述：空气以69.2m/s的速度吹向直径为1.9cm的圆柱体，用Fluent仿真此时产生的噪声。基于圆柱体直径的Reynolds数大概是90000。其他尺寸参数见下图。对于声学仿真，推荐使用LES湍流模型，因为LES模型求解所有涡旋尺度比网格尺度大的涡旋，能较好预测到噪声。1、启动软件并导入网格1.1启动Fluent软件，选择2D双精度版本，单核求解。1.2导入网格文件“cylinder2d.msh.gz”，网格下载在文章底部。为了改善求解速度，将网格重新读取编录，操作：Mesh->Reorder->Domain在文本窗口中显示Fluent采用了ReverseCuthill-McKee方法进行。2、求解器设置3、模型设置3.1湍流模型-大涡LES模型在2D求解器中，LES模型是隐藏的，就是你打开湍流模型面板是找不到的。在文本窗口中输入下面命名“(rpsetvar'les-2d?'#t)”，键盘回车键。命令输入要英文状态，括号也要输入，还有一点就是不能复制黏贴输入，只能手动敲键盘输入才有效，本人亲测过了，Fluent版本是15.0。再次打开湍流模型，就发现LES已经出现可选了。此时会弹出一个warning提示框，点击OK即可。4、边界条件4.1inlet边界，边界类型为velocity-inlet。4.2outlet边界，边界类型为pressure-outlet。保留默认设置。5、求解设置5.1离散方案设置。5.2松弛因子设置，将pressure松弛因子调到0.7。5.3残差设置。5.4初始化5.5此时，保持case和data文件，命名为“cylinder2dt0.00”。设置计算步数，先计算20步。5.6迭代完20步后，设置面监控升阻力。先设置阻力监控具体参数设置如下。设置升力监控。设置计算参考值。保持case和data文件，文件名“cylinder2dt0.00”，覆盖之前保持的。5.7输入迭代步数4000.点击calculate计算。计算完后流动的时间是t=0.02s。此时，升力系数如下。此时，阻力系数如下。保持case和data文件，命名为“cylinder2dt0.02”。6、设置声学模型操作：Models->Acoustics->Edit...将压力松弛因子设置为1.点击计算，将会继续迭代4000步，计算完后，保持case和data文件，命名为“cylinder2dt0.04”。修改声学模型的参数，操作：Models->Acoustics->Edit...将SourceCorrelationLength设置为0.095m，此值等于圆柱体直径0.019m的5倍。如果是3D模型计算，此值不需要设置。设置声学信号。操作：Solution->RunCalculations->AcousticSignals...设置两个信号接收器receivers。其中一个在Y方向的35D（D为圆柱体直径）就是-0.665m处，另外一个在128D就是-2.432m处。计算和读取声学信号进两个接收器中。7、声学后处理7.1显示声压随时间变化图。选择两个接收器。点击plot按钮，显示声压随时间变化曲线。可以看到，红色线（接收器2）的时间相对黑色线（接收器1）平移到大约5e-3s，而且也弱很多。这是因为接收器2比接收器1离圆柱体更加远。7.2傅里叶转换。操作：Results–>Plots–>FFT。参数设置如下接收器1的声压谱显示如下。来源：CFD饭圈
