Fluent教程案例4-非牛顿流体三维环空内流动.pdfVIP

实验四、三维环空流动的数值模拟

在石油工程，环空内的流动是最常见的一种流动，本实验模拟环空内的牛顿流体和非牛顿流体的流动。

本文旨在学习非牛顿流体模拟的设置，辅助线法构建网格和移动（旋转）壁面条件的应用。

1物理模型

三维环空管长5米，外圆半径0.5m，内小圆半径0.1m，小圆偏心距为0.1米。

流体介质：非牛顿流体。

Inlet：流速入口2m/s

Outlet：流出outflow

2数值模拟原理

方程求解：采用双精度求解器，定常流动，层流，SIMPLEC算法。

3建立模型

3.1首先建立三维水平放置环空的几何模型Geometry。

如果不利用辅助线而是直接对偏心环空进行网格构造则产生不好的网格。如下图是对偏心圆无辅

助线直接绘制的网格，网格质量差。故本文采用添加辅助线构建合理化网格。

1）利用geometry/face/createrealcircularface生成同心大小圆。

将小圆x方向移动0.1m，形成偏心圆。

2）为了改善环空网格，利用move/copyvertices生成新节点，利用节点添加过两圆心的辅助线。

连接两节点，生成辅助线。

将辅助线扫略（sweep），向z轴正方向sweep5个单位大小生成辅助面。其中SweepEdges面板中Vector默

认的Magnitude是1m，需要调整到5米。

3）利用面的布尔运算，将小圆从大圆中减去，Face/SubtractRealFaces，生成Face1，得到偏心圆面。

4）为改善网格将得到的偏心面用辅助面分割，再Sweep形成计算域的三维环空。

将分割后的两个面选中做扫略（Sweep）成三维体。Geometry/volume/sweepFcae，其中SweepEdges面板

中Vector默认的Magnitude是1m，需要调整到5米。

3.2生成网格，由边到面网格到体网格。.

1）设置大圆和小圆的边节点数intervalcount为25，辅助线部分节点数intervalcount10，完成边网格

设置，选中2个面，利用Elements默认Quad，Type：Submap点击应用完成面网格生成。

再将z方向的四条边设置100网格节点，即可去生成体网格。选中2个半环空体，按默认参数生成体网格。

3.3设置边界条件：左侧两面为速度进口，右侧两面为出口outflow；设置内部表面的固体壁面为innerwall，

其余的默认为固体壁面wall。

3.4保存文件生成huankong.dbs，输出huankong.msh等相关文件。

4在FLUENT中求解非牛顿流体流动

4.1FLUENT设置:3d双精度

除以下设置为非牛顿层流所必须设置的外，其余选项和三维圆管相同，不再详述。

1）Viscous设置

环空内流动为层流，Viscous设置为Laminar。

2）材料设置为非牛顿流体，需要自定义材料，点击Materials面板右侧的User-DefinedDatabase…,

输入材料名称”Newnonnewton“，点击OK，在弹出的“User-DefinedDatabaseMaterials”对话框

中单击“New”按钮，在name中输入材料名称，Type中选择fluid,如下图，

在可用属性（AvailableProperties）中将密度和粘性分别选中，单击Edit按钮，将密度设为

1200kg/m3，

再选中Viscous，可用属性选择非牛顿幂律模型，单击Edit按钮，将ConsistencyIndex设为1.8，

PowerLaw为1.5，单击ok，再单击Apply按钮，回到“User-DefinedDatabaseMaterials”对话框，

选中Non-n-powerlaw,单击Copy按钮，添加Non-n-powerlaw材料到当前工程，即将流动中的流体

定义为Non-n-powerlaw。

3）Boundary边界条件设置

设置速度入口，为2m/s，

Outlet设置为自由出口Outflow。