六十六、Fluent离心泵旋转流场模拟

本文模拟离心泵泵水的过程，泵叶轮有5个叶片，旋转速度为1450转/分。蜗壳出口的质量流量为90 kg/s。入口表压为0 pa。模拟将确定泵产生的压头，表示流体的总压力增加

3.1 导入网格：

使用Fluent软件打开Chapter66.msh网格文件，文件在本文末尾链接资源内。注意，本案例为3D模型，打开Fluent软件时需勾选3D

导入mesh文件后，点击Display，将显示模型

3.2 修改模型尺寸

本案例模型尺寸不需要修改，打开Scale mesh查看模型尺寸是否正确

3.3 求解器设置

4.1 湍流模型设置

不打开能量方程

湍流方程选择k-omega-SST模型，与其他双方程模型相比，如k-e模型，k-omega-SST湍流模型对于叶轮机械中的流动分离模拟更加准确

4.2 设置water材料属性

在材料设置界面，双击air，打开材料设置界面。点击Fluent Database

在Fluent Fluid Material栏下滑到最下面，找到water-liquid，点击下面的copy，即将液体water复制到材料中

本案例有两个计算域，分布为impeller叶轮和volute蜗壳

5.1 impeller设置

双击impeller，弹出如下界面

Material Name选择water-liquid

勾选Frame Motion

Rotational Velocity下的Speed (rpm)设置为1450

Rotation-Axis Direction设置z轴为1

单击apply

注：Frame Motion和Mesh Motion的适用性

Frame Motion用于在相对于固定坐标系定义流体运动时使用。当流体的运动简单且可以由几个参数定义时。例如，如果想模拟绕旋转圆柱的流动，可以使用Frame Motion。

Mesh Motion用于在相对于网格定义流体运动时使用。当流体的运动复杂且不能由几个参数定义时。例如，如果想模拟发动机中活塞周围的流动，可以使用Mesh Motion。

简单来说，Frame Motion是坐标系动，Mesh Motion是网格动，复杂性Frame Motion<Mesh Motion<yynamic mesh

5.2 volute设置

双击volute，弹出如下界面

Material Name 选择water-liquid

其他的保持默认即可

本案例的边界很多，但是对于internal和interface这里边界基本不需要设置，只需要对inlet和wall边界进行设置即可

6.1 wall边界设置

impeller-hub为叶轮轮毂，因此会随着叶轮一起运动，其相对于叶轮间的流体可认为是静止的。

inblok-shroud为入口的壁面，相对于绝对坐标系是静止的。

对impeller-hub边界进行设置，双击impeller-hub

Wall Motion栏选择Moving Wall，Motion栏选择Relative to Adjacent Cell Zone，下面勾选Rotational。其他的保持默认。

解释一下：

Wall Motion栏选择Moving Wall表示impeller-hub是移动壁面，不是固定不动的；

Motion栏是对这个移动壁面的运动进行具体的设置。

选择Relative to Adjacent Cell Zone，表示下面的运动设置都是相对于相邻计算域的.

举个例子，比如下面设置选择Rotational，设置速度speed为100，表示的并不是impeller-hub旋转速度是100r/min，而是相对于相邻计算域（impeller是1450r/min），旋转速度是100r/min，即1550r/min

勾选Rotational，表示运动是旋转运动，因此需要设置旋转轴的坐标和方向。

由于impeller-hub相对于impeller是静止的，因此Speed设置速度是0。

对inblock-shroud边界进行设置，双击inblock-shroud

Wall Motion栏选择Moving Wall，Motion栏选择Absolute，下面勾选Rotational。其他的保持默认。

由于inblock-shroud相对于绝对参考系是静止的，因此Speed设置速度是0。

这里Wall Motion设置为Stationary Wall即固定壁面也是可以的，上述设置只是为演示需要。

6.2 inlet设置

inlet边界条件设置为压力入口

将Supersonic/Initial Gauge Pressure设置为-100pa。这个压力只是初始化压力，对计算结果无影响。参考文章五十二、Fluent瞬态可压缩流动

6.3 outlet设置

设置mass-flow-outlet-11为质量流量出口。Fluent默认

mass-flow-outlet-11为质量流量入口，需要进行更改

双击mass-flow-outlet-11，设置Mass Flow Rate为90

打开Turbo Model涡轮模型，

Domain-Turbo Model-Enable

创建turbo interfaces.

Domain → Turbo Model → Turbo Create...

Mesh Interface：输入nps

Interface Zones Side 1：选择 interface-impeller-outlet

Interface Zones Side 2：选择interface-volute-outlet

Interface Options：勾选General Turbo Interface

General Turbo Interface：选择No Pitch-Scale

8.1 求解方法

Spatial Discretization-Gradient 选择Green-Gauss Node Based

勾选High Order Term Relaxation

勾选Pseudo Transient，参考文章四十三、Fluent增强收敛性-伪瞬态计算

8.2 求解控制

Solution-Solution Controls

保持默认即可

9.1 监测蜗壳出口的总压力

Solution → Reports → Definitions → New → Surface Report → Mass-Weighted Average...

双击Reports Definitions，弹出下图窗口

依次点击New，Surface Report，Mass-Weighted Average可弹出下图窗口。

name：输入p-out

Field Variable：选择Pressure...和Total Pressure

Surfaces：选择Outlet下的mass-flow-outlet-11

勾选Report-Plot，在Fluent会显示曲线图

检测mass-flow-outlet-11面的总压力，以质量平均方式计算出口平均总压力。

9.2 检测入口总压

Solution → Reports → Definitions → New → Surface Report → Mass-Weighted Average...

重复8.1的步骤，直到弹出下图

name：输入p-in

Field Variable：选择Pressure...和Total Pressure

Surfaces：选择inlet下的inlet

勾选Report-Plot，在Fluent会显示曲线图

检测inlet面的总压力，以质量平均方式计算入口平均总压力。

9.3 叶片blade检测

Solution → Reports → Definitions → New → Surface Report → Integral...

双击Reports Definitions，弹出下图窗口

依次点击New，Surface Report，Integral可弹出下图窗口

name：输入p-blade

Field Variable：选择Pressure...和Static Pressure

Surfaces：选择wall下的blade

勾选Report-Plot，在Fluent会显示曲线图

检测blade面的静压力，对blade进行静压力积分计算。

9.4 创建泵头表达式

Solution → Reports → Definitions → New → Expression...

双击Reports Definitions，弹出下图窗口

单击Expression，输入公式

name：输入head

公式栏输入：({p-out} - {p-in}) / (998.2[kg/m^3] * 9.81[m/s^2])

其中的998.2为水的密度，9.81为重力加速度，这里是计算压力差产生的泵压头。

勾选Report Plot，表示在图形界面显示曲线图

勾选Print to Console，表示在控制台输出值

混合初始化

单击More Settings...，对混合初始化进行高级设置。

勾选Absolute

用于设置初始速度是绝对速度（absolute）还是相对速度。只有当问题涉及移动参照系moving reference frames或者sliding meshes网格时，此选择才有意义。如果没有，则这两个选项是等效的。

勾选Use Specified Initial Pressure on Inlets

正常混合初始化是无法设置初始物理场的，勾选此选项，表示使用inlet压力作为初始化压力。

回到混合初始化界面，单击initialize

Number of Iterations设置为500，

Timescale Factor设置为10，这是伪瞬态设置，可参考文章四十三、Fluent增强收敛性-伪瞬态计算

单击 Calculate

12.1泵产生的水头

12.2 出口总压和压头曲线

12.3 速度云图

12.4 压力云图

下载链接见附件