第十一周：CFD有限差分法+排气歧管内流场案例仿真

摘要

在本周中，我主要对CFD中的有限差分法部分、拉格朗日插值以及牛顿插值进行了初步的学习。了解到什么是迎风型差分格式和如何对复杂网格进行差分处理，何为守恒型的差分格式。以及拉格朗日的插值多项式和牛顿插值多项式。关于Fluent的实例运用中，本周进行了关于排气歧管内流场的前、后处理，以及求解器设置。

Abstract

In this week, I mainly studied the finite difference method part of CFD, Lagrange interpolation, and Newtonian interpolation，and learned what the windward differential format is, how to deal with the  complex meshes in differentially process, and what was the conserved differential format. and Lagrangian interpolated polynomials and Newtonian interpolated polynomials. In the case of Fluent, this week I worked on the pre- and post-processing of the flow field inside the exhaust manifold, as well as the solver settings.

Ansys 部分

前处理

关于排气歧管的前处理，主要关注于其面网格、边界层、体网格的设置，下图为本次进行内流场模拟的模型：

面网格绘制将最大尺寸设定为6mm，最小尺寸设置为0.6mm，初步进行面网格的绘制，绘制结果如下：

由于本次案例的网格总数较少，且无需特别加密网格的部分，故不添加局部尺寸，对内流场进行边界封堵，确定内流场的位置，将进出口封堵起来，使软件自动识别其内部流场区域进行描述展示，并检查是否有未被划出的部分：

内流场网格绘制情况：

流体区域：

固体区域：

划分好流场后，进行边界层设置，设定为3层，增加率为1.2，后按软件给与的默认值进行体网格的绘制，完成绘制后发现网格质量为0.07<0.1，网格质量不合格，排除其体网格绘制中精度问题后，返回面网格部分进行重新绘制：

将面网格最小精度调整为0.4mm，最大尺寸不变，其余设置相同与上次操作，进行体网格设置，发现网格质量为0.16>0.1，增加了约四万个网格，网格质量合格。

体网格剖面图如下，前处理部分完成：

求解器设置

本次案例中，以空气为流体，采用k-epsilon、realizable、SWF模型，该选择稳定性好，发散概率小，易于收敛，适于新手使用。只关注其速度与压力的关系，故不勾选能量模型。且重力的影响较少，故不选则重力选项。

三个进口速度都设置为10m/s，垂直入射，出口表压由于舍入误差和该案例压降存在几十pa，相对于大气压强过于小，故设置为0，如若需要调整工作条件，可在物理模型-工作条件中对模拟环境进行调整。

对运算进行初始化，基于每个网格以迭代的方式满足设定方程所需的值，将残差调整为0.0001，进行300次迭代运算，求解器设置完成。

后处理

进行186次运算后，判定其达到收敛，以下为其加权平均面积分布图：

对于判断收敛，有三个条件：

（1）残差足够小

（2）通量到达守恒

（3）所观测值曲线图达到稳定

查看报告，查看进出口的通量是否守恒：

发现最终结果基本无差距，故本次运算到达收敛。流场壁面压力分布如下：

速度分布如下：

速度矢量分布：

迹线图：

CFD 部分

见下图：

总结

本周Fluent部分使我受益良多，明白并不是按照案例进行一步步的运算即可出现稳定的结果，要时刻关注其网格绘制情况，若不注意即会酿成大错，且对于模型后处理的部分也有更深刻的理解。在CFD基础理论中，明白了有限差分法在遇到复杂网格处理时的劣势，需将复杂网格通过雅可比变换等方式化成均匀网格才可进行计算，以及在存在有间断处的情况，如何使用守恒型的差分格式对其进行描述。