第十九周：FLUENT学习周报

目录

摘要

Abstract

Fluent仿真实例部分

一、关于同心圆的切向速度与解析解的对比

1.基本描述

1.1几何模型

1.2解析解描述

2.网格划分

3.求解器设置

二、同心圆热转换

1.基本描述

2.求解器设置

Fluent基础知识部分

一、边界层

二、Y+值

三、网格质量的判断

总结

---

摘要

在本周中，我主要进行了两个案例的仿真模拟练习和Fluent中的基础知识补充，两个案例使用同样的几何建模，仿真模拟的主要内容为:对两个同心圆柱体间的定常层流进行模拟,观察其切向的速度与解析解的对比和其热转换的观察。关于Fluent基础知识的补充部分，包括了：边界层的概念和分布；Y+值的概念和选择方式以及网格质量的判断。

Abstract

In this week, I have mainly maked the two case simulation exercises and the basic knowledge supplement in Fluent, both cases use the same geometric modeling, the main content of the simulation is to simulate the stationary laminar flow between two concentric cylinders, observe the tangential velocity compared to the analytical solution, and observe the thermal transformation. A supplementary section on the basics of Fluent includes: the concept and distribution of boundary layers; The concept and selection of Y+ value and the judgment of mesh quality.

Fluent仿真实例部分

一、关于同心圆的切向速度与解析解的对比

1.基本描述

  对两个同心圆柱体间的定常层流进行模拟。内部圆筒以恒定的角速度旋转引起流动，外部圆筒保持静止。流动是稳态的。利用层流的解析方程可以计算出不同截面的切向速度。这些数值可用于与仿真结果进行比。

1.1几何模型

外部圆柱的半径为46.28mm，内部圆柱为17.8mm，流体为密度1kg/m，粘度为0.0002kg/m-s的空气，外部圆柱保持不动，内部圆柱以1rad/s的角速度进行旋转。

1.2解析解描述

  由于问题中几何体为对称形式，故只存在其切向的速度u0不为0，在纵向的速度为0，其切向速度为半径的函数，不考虑其温度，且外壁为静止，故角速度为0。

  根据一个或两个圆柱体的稳定角速度，在两个问心圆柱体之间保持稳定的流动。设内柱半径r0，角速度 w0，温度T0，外柱半径r1，角速度 w1，温度T1。这种几何结构使得只有切问速度u0不为零，变量、T和p为半径r的函数。在极坐标下的运动方程简化如下：

可得出其切向动量方程解的形式为：

从边界条件中得到C1和C2的值，得出其解析解的形式如下：

2.网格划分

  在Design Modeler中完成几何建模后，转入Ansys mesh对网格进行绘制，选择内外壁调用尺寸调整（Sizing）选项，选择分区数量，将其分区数量改为300，在高级中选择硬（hard），使网格的划分只随着分区数量进行，划分后网格如下图所示：

后适用网格划分方法命令（Method）进行调整，将网格全调整为四边形网格，调整后网格如下：

发现还存在三角形区域，故调用面网格划分，使其四边形网格在面上均匀分布，划分完成后网格划分如下：

网格节点数为12900，单元数为12600，完成网格划分。其相关网格质量指标如下：

3.求解器设置

进入Fluent求解器设置，检查网格质量后，选择计算模型为层流模型，将材料的密度按按照上述描述进行更改，具体设置如下：

设置边界条件，内壁圆柱体角速度为1rads，具体设置如下：

使用SIMPLE方法进行耦合，采用混合初始化，进行400步迭代，结果如下所示：

生成一个等值面在半径的位置，观察内壁半径切向上的速度变化，导出其在y方向上的速度变化，导入EXCEL表格中，进行线性拟合的操作，具体设置如下：

切向速度变化：

在Excel中的拟合结果：

二、同心圆热转换

1.基本描述

再对圆环热转换进行仿真模拟，其基本条件为：有两个同心的圆柱体，它们以相同的20rps的角速度向相反的方向旋转。内筒是在305K。而外筒的温度为295K。这两个气缸之间的空间充满了甘油。网格划分同上。

2.求解器设置

求解器设置中，打开能量方程，并选择甘油为流体材料，其具体数值如下：

对内、外壁圆柱体进行运动设置，将角速度设置为-20rads，温度改为305k，壁厚改为0.001m，外壁角速度设定为-20rads，为逆时针旋转，温度为295k，完成边界条件设置。

(inner cylinder的速度应为20rads，为正数，表示为顺时针旋转）

求解方式采用耦合求解（coupled），具体设置如下：

设置残差为10-6，进行混合初始化后开始计算，迭代次数为1000，计算结果如下：

速度云图：

速度矢量图如下，可以看到在两圆柱体相交处，顺逆时间旋转的速度方向变化：

温度云图：

Fluent基础知识部分

一、边界层

二、Y+值

三、网格质量的判断

总结

  在本周里找到了同样运用同心圆讨论的案例，对于Fluent的操作更加熟悉，且加强了原本不理解的FLuent中的术语，网格划分作为仿真模拟中最重要的方面，关于如何判断网格质量和如何绘制一幅好的网格，这些基础知识必不可少。在下周，我会加强后处理即结果展示部分和流体力学的学习，加深在实际操作中运用知识的能力。