第十周：CFD差分方法学习+ansys 后台阶流仿真案例学习

摘要

在本周，主要进行了Fluent的经典后台阶流案例的学习和对CFD中的有限差分法进行了初步的理论学习。在fluent案例中，主要学习了对两个台阶条和会发生回流的断面处进行了加密网格设置，使其更容易进行收敛的计算。对于CFD的有限差分法，在差分格式的构造，和描述差分格式精度与差分的修正方程进行了学习。

Abstract

In this week, I totally studied the classic post-step flow case of Fluent and made a preliminary theoretical study of the finite difference method in CFD. In the Fluent case, I  learned to set up the encryption mesh at the two step bars and the section where the reflow will appear, so that could make it easier to perform convergence calculations. For the finite difference method of CFD, the construction of the differential format and how to describe the modified equations the accuracy and difference of the differential format were learned.

Fluent实例部分

  该案例是关于在一个无限长的平板进行拉伸，可以对应成在一个圆管中有一个变径的形式，在流动的下游区会发生一个回流，出现回流后则会出现能量的损失，为了减少该能量的损失，则要进行网格的优化计算。在本案例中，对该模型进行了平移拉伸的方法，表述出的近二维模型。

  将模型载入fluent，将其在Z方向上拉伸1/4的台阶高度，将其变成三维模型，用三维模型代表二维。使用Fluent中将使用六面体网格绘制。对于台阶流的研究，会对台阶高度作为主要变量，总结其回流的程度，x轴为台阶高度，y轴为能量的损失。

  在Fluent中，添加边界层，使用body of influence 对两个台阶处进行加密，将网格尺寸设置成13mm，一般设置为台阶的1/10，若最后无法收敛，则继续进行加密。在对两个发生回流的平面出进行面网格加密，尺寸设定为0.65。

  

将网格划分模式设定为仅剩曲率，相当于安最大（小）值进行网格划分的设定，最长不可超过边长，否则网格将会出现不均匀的现象，将其设定为5。该案例中由于存在梯度的原因，对其进行网格加密可达成更好的收敛。

网格划分如下图：

进入求解器设置，由于该案例网格数较少，故可以选用k-omega SST模型，进行精密计算。SST模型中在y+值1附近（至少低于5）有较好的表现，倘若y+值很大，则选用k-epsilon模型可能会更好。材料设定为air，湍流强度和水力直径如下图所示： 

对于水力直径的选择，应设定为2倍的y方向上的高度。进口速度设定为41.7，其余设置默认，进行300次的迭代运算，在fluent中，会将湍流粘性比控制在小于50w的数下，结果为下图：

可观察到湍流速度并未达到稳定状态，需将连续性方程的残差调试至0.00001，再次进行运算，由于连续性方程的稳定最难达到，故只需调整其残差，其余项必定收敛稳定，调试结果如下：

可见在进行100次运算后，各项参数均达到稳定状态，模拟收敛。以下为流场压力、速度云图，在速度云图中，可看到其回流部分的速度方向和大小：

CFD部分

关于cfd学习的部分见下图：

总结

  在本周的fluent实例中，了解到后台阶流模型是一种经典的湍流分析模型，由于其存在梯度特性，且会发生回流，故在实际工程问题中如何减少其能量的损失便需要将发生回流的部分进行网格加密和调整其网格最大值的操作，且为使最终结果收敛，可将最难收敛的项即连续性方程的残差调整，使其稳定，故其余项便达到稳定状态。在CFD的学习中，只是将有限差分法的概念进行了一个开头，在下周将会对其展开讨论。