ANSYS教学算例集FL 离心泵叶轮气蚀模拟课件

1 ©

2018

ANSYS,

Inc. August

30,2023 ANSYS

Confidential离心泵叶轮气蚀模拟摘要在分析泵系统的性能时，气蚀出现会显著影响泵的性能。气蚀是由于局部压力低于当前温度下液体的饱和蒸汽压力，产生的汽化现象。产生气蚀时，管路容易充满气泡，导致泵的性能曲线下降、能耗上升，还会产生振动和噪音，不仅影响泵的使用寿命和可靠性，而且对周边产生噪音污染。因此仿真计算可以帮助设计工程师在产品设计时模拟水泵流量压力变化时的气蚀性能，从而指导工程师通过结构调整提高水泵的抗气蚀性能。本算例使用ANSYS

Fluent

19.0软件，对离心泵内部流动进行仿真，文档内包含详细的网格导入、模型选择、材料物性、边界条件、求解参数、后处理的设置。通过仿真计算获得泵内流场的形态及气蚀状况。ANSYS

Confidential2©

2018

ANSYS,

Inc.August

30,2023案例描述ANSYS

Confidential3©

2018

ANSYS,

Inc.August

30,2023本算例仿真对象为一离心泵，包含五个叶片，以2160rpm运作。流体通过进口管路进入到泵中，通过叶轮的做功，使流体最终沿着径向方向流出。仿真中假设流动介质为水，流动为稳态。流体参数如下：液体密度：1000kg/m3液体粘度：0.001Pa

s蒸汽密度：0.01927kg/m3蒸汽粘度：8.8E-6Pa

s液体饱和蒸汽压力：3540Pa本次计算假设液体以总流量210m3/hr进入进口，计算中将折算成速度入口边界。液体进入进口假设为充分发展的湍流，进口湍流强度为5%。出口边界条件设置

为均一的压力出口条件，进口静压取决于出口绝对压力。在当前case中，出口压力设置为500000Pa。操作步骤ANSYS

Confidential4©

2018

ANSYS,

Inc.August

30,20233.1.创建工作目录并启动Fluent在硬盘上创建英文名称的文件夹（例如Cavitation），将网格文件centrif-pump.msh.gz拷贝至该目录下。启动Fluent

19.0，在Fluent

Launcher中，

Dimension选择3D，Display

Options中勾选Display

Mesh

After

Reading和Workbench

Color

Scheme，Processing

Options选择Serial，设置使用单核计算（用户可以根据现有的硬件资源和License授权酌情选择合适的并行数），更改WorkingDirectory路径至该网格文件目录下，点击OK启动Fluent

19.0。操作步骤3.2.导入网格文件菜单中点击【File】>【Read】>【Mesh…】，选取网格文件centrif-pump.msh.gz，点击OK导入网格。此时，图形界面中可以查看导入的网格。ANSYS

Confidential5©

2018

ANSYS,

Inc.August

30,2023操作步骤General一般设置在最左侧的树中，鼠标左键双击【Setup】>【General】，进行网格相关的操作以及选择求解器。缩放网格树中点击【Setup】>【General】>【Mesh】>【Scale…】，View

Length

Unit

In选择

mm，以毫米为单位查看当前模型的尺寸范围是否正确。ANSYS

Confidential6©

2018

ANSYS,

Inc.August

30,2023操作步骤3.3.2.调整转速单位树中点击【Setup】>【General】>【Units】，在【Set

Units】界面下，选择

Quantities为angular-velocity，选定右侧Units单位为rpm（即以“转/分”单位显示和设置转速）。ANSYS

Confidential7©

2018

ANSYS,

Inc.August

30,2023操作步骤3.4.Models模型设置在最左侧的树中，鼠标左键双击【Setup】>【Models】，进行物理模型设置。3.4.1.选择湍流模型树中【Setup】>【Models】>【Viscous】，鼠标左键双击Viscous，Model中选择k-epsilon，k-epsilonModel中选择Standard，Near-WallTreatment中选择Standard

Wall

Functions，点击OK完成设置。ANSYS

Confidential8©

2018

ANSYS,

Inc.August

30,2023操作步骤3.4.2.混合多相流模型设置树中【Setup】>【Models】>【Multiphase】，鼠标左键双击Multiphase，Model中选择Mixture，去选Slip

Velocity

[1]。保持Number

of

Eulerian

Phases中数值为2，点击

OK完成设置。[1]当第二相相比起第一相密度小很多时，可以不考虑滑移速度。ANSYS

Confidential9©

2018

ANSYS,

Inc.August

30,2023操作步骤3.5.Materials材料设置在最左侧的树中，鼠标左键双击【Setup】>【Materials】，进行材料物性设置。3.5.1.定义液态水属性在Materials的Task

Page中选中点击Fluent

Database…选项按钮，在弹出的FluentDatabaseMaterials对话框中选择water-liquid，点击Copy按钮，将液态水参数拷贝至Materials下。ANSYS

Confidential10©

2018

ANSYS,

Inc.August

30,2023操作步骤回到Create/Edit

Materials，在Properties下设置Density(kg/m3)值为1000，设置

Viscosity(kg/m-s)值为0.001。点击Change/Create选项。3.5.2.定义水蒸气属性在Materials的Task

Page中选中点击Fluent

Database…选项按钮，在弹出的FluentDatabaseMaterials对话框中选择water-vapor，点击Copy按钮，将水蒸气参数拷贝至Materials下。ANSYS

Confidential11©

2018

ANSYS,

Inc.August

30,2023操作步骤回到Create/Edit

Materials，在Properties下设置Density(kg/m3)值为0.01927，设置

Viscosity(kg/m-s)值为8.8e-6。点击Change/Create选项。ANSYS

Confidential12©

2018

ANSYS,

Inc.August

30,2023操作步骤3.6.混合相定义在最左侧的树中，鼠标左键双击【Setup】>【Models】>【Phases】，进行混合项定义。选择phase-1–PrimaryPhase，点击Edit…，在PhaseMaterial选择water-liquid，点击OK确认设置。ANSYS

Confidential13©

2018

ANSYS,

Inc.August

30,2023操作步骤选择phase-2–PrimaryPhase，点击Edit…，在PhaseMaterial选择water-vapor，点击OK确认设置。回到Phases面板，点击Interaction…，在弹出的Phase

Interaction面板中，设置

Number

of

Mass

Transfer

Mechanisms为1，选择Mass

Transfer栏下的To

Phase为phase-2，设置Mechanism为cavitation。ANSYS

Confidential14©

2018

ANSYS,

Inc.August

30,2023操作步骤设置空化模型后，弹出空化模型设置窗口。选择Model为Schnerr-Sauer，CavitationProperties下Vaporization

Pressure设置为3540

pascal，点击OK确认设置。ANSYS

Confidential15©

2018

ANSYS,

Inc.August

30,2023操作步骤3.7.CellZone

Conditions设置在最左侧的树中，鼠标左键双击【Setup】>【CellZoneConditions】，进行体网格区域条件的设置。在Cell

Zone

Conditions的Task

Page中点击Operating

Conditions…，选择

Operating

Pressure为0

pascal，点击OK确认。ANSYS

Confidential16©

2018

ANSYS,

Inc.August

30,2023操作步骤在CellZoneConditions的TaskPage中选中fluid点击【Edit…】，在弹出的Fluid窗口中，勾选FrameMotion，在ReferenceFrame标签栏下选择Rotation-AxisDirection为0,0,1。（Fluent中转动方向按照右手定则定义，Z=1代表绕z轴正方向旋转。）在Rotational

Velocity中输入Speed值为2160

rpm。点击OK退出Fluid设置对话框。ANSYS

Confidential17©

2018

ANSYS,

Inc.August

30,2023操作步骤3.8.Boundary

Conditions设置在最左侧的树中，鼠标左键双击【Setup】>【BoundaryConditions】，进行边界条件的设置。ANSYS

Confidential18©

2018

ANSYS,

Inc.August

30,2023操作步骤3.8.1.定义速度进口边界条件在BoundaryConditions的TaskPage中，鼠标左键单击选中inlet，点击【Edit…】，弹出Velocity

Inlet窗口，Momentum选项卡中设置Direction

Specification

Method为Magnitude，NormaltoBoundary，设置VelocityMagnitude为7.0455m/s，设置Supersonic/Initial

Gauge

Pressure为450000。设置Turbulence下SpecificationMethod为Intensity

and

Hydraulic

Diameter，在TurbulentIntensity为5%，Hydraulic

Diameter为103mm。点击OK完成设置。ANSYS

Confidential19©

2018

ANSYS,

Inc.August

30,2023操作步骤回到BoundaryConditions的TaskPage中，保持选中inlet，切换Phase中的phase-2，点击【Edit…】。在打开的对话框中Multiphase标签栏下设置Volume

Fraction为0。点击OK确认设置。ANSYS

Confidential20©

2018

ANSYS,

Inc.August

30,2023操作步骤3.8.2.定义出口边界条件在BoundaryConditions的TaskPage中，鼠标左键单击选中pressure-outlet-stator，确认其Type是pressure-outlet，如果不是，从Type的下拉菜单中切换选择。确认将

Phase项设置调整为mixture，点击【Edit…】，弹出Pressure

Outlet窗口，在

Momentum选项卡中，Gauge

Pressure设置为450000。点击OK完成设置。回到Boundary

Conditions的Task

Page中，保持选中outlet，切换Phase中的phase-2，点击【Edit…】。ANSYS

Confidential21©

2018

ANSYS,

Inc.August

30,2023操作步骤在打开的对话框中Multiphase标签栏下设置BackflowVolumeFraction为0。点击OK确认设置。ANSYS

Confidential22©

2018

ANSYS,

Inc.August

30,2023操作步骤3.8.3.周期边界条件设置在BoundaryConditions的TaskPage中，鼠标左键单击选中periodic.10，确认Phase项为mixture，点击【Edit…】，弹出Periodic窗口，设置PeriodicType为Rotational。点击OK确认设置。同理，将periodic.11边界设置成同样的旋转周期类型。ANSYS

Confidential23©

2018

ANSYS,

Inc.August

30,2023操作步骤ANSYS

Confidential24©

2018

ANSYS,

Inc.August

30,20233.8.1.定义壁面边界条件1.定义出口轮毂、机匣及进口机匣边界（hub-outer、shroud-outer、inlet-shroud），这三个边界由于在物理上不旋转，因此需要将其定义为静止运动边界。在BoundaryConditions的TaskPage中，鼠标左键单击选中hub-outer，确认其Type是wall，如果不是，从Type的下拉菜单中切换选择。点击【Edit…】，弹出Wall窗口，在Momentum选项卡中，选择Wall

Motion为Moving

Wall，选择Motion为Absolute，Rotational，Speed设置为0rpm，点击OK完成设定。操作步骤回到BoundaryConditions的TaskPage中，点击Copy…打开CopyConditions设置窗口，在From

Boundary

Zone中选择hub-outer，在To

Boundary

Zones中选择inlet-shroud、shroud-outer，点击Copy完成复制边界设置。2.定义转子轮毂、机匣及叶片边界（hub、shroud、blade），这三个边界由于是随着转动部件，因此需要设置为旋转边界。在Boundary

Conditions的Task

Page中，鼠标左键单击选中blade，确认其Type是wall，如果不是，从Type的下拉菜单中切换选择。点击【Edit…】，弹出Wall窗口，在

Momentum选项卡中，选择Wall

Motion为Moving

Wall，选择Motion为Relative

toAdjacent

Cell

Zone，Rotational，Speed设置为0rpm，点击OK完成设定。ANSYS

Confidential25©

2018

ANSYS,

Inc.August

30,2023操作步骤回到BoundaryConditions的TaskPage中，点击Copy…打开CopyConditions设置窗口，在From

Boundary

Zone中选择blade，在To

Boundary

Zones中选择hub、shroud，点击Copy完成复制边界设置。ANSYS

Confidential26©

2018

ANSYS,

Inc.August

30,2023操作步骤ANSYS

Confidential27©

2018

ANSYS,

Inc.August

30,20233.9.Solution

Methods求解方法设置在最左侧的树中，鼠标左键双击【Solution】>【Methods】，进行求解方法的设置。设置的Scheme为Coupled，设置Spatial

Discretization中的Gradient为Least

SquaresCellBased，Pressure为PRESTO!，Momentum为SecondOrderUpwind，其余为FirstOrderUpwind。操作步骤在最左侧的树中，鼠标左键双击【Solution】>【Controls】，进行松弛因子的设置。选择Under-RelaxationFactors下TurbulentKineticEnergy为0.5，TurbulentDissipation

Rate为0.5，Turbulent

Viscosity为0.1。ANSYS

Confidential28©

2018

ANSYS,

Inc.August

30,2023操作步骤3.9.1.流场域内出口流量监视设置树中【Solution】>【Monitors】>【ReportPlots】，鼠标左键点击New…，进行自定义物理量的监视设置。在弹出的New

Report

Plot窗口中，点击New…，选择Surface

Monitor>Area-Weighted

Average，并同时在Surface选项栏中选择进口面inlet，勾选Report

Plot和Print

to

Console，点击OK。回到New

Report

Plot窗口，点击OK。ANSYS

Confidential29©

2018

ANSYS,

Inc.August

30,2023操作步骤3.10.Initialization初始化树中【Solution】>【Initialization】，鼠标左键双击Initialization，在TaskPage中确认InitializationMethods选择为StandardInitialization，点击Computefrom下拉菜单，选择inlet，点击Initialize按钮进行初始化。ANSYS

Confidential30©

2018

ANSYS,

Inc.August

30,2023操作步骤3.11.Run

Calculation运行计算树中【Solution】>【Run

Calculation】，鼠标左键双击Run

Calculation，在TaskPage中的Number

of

Iterations输入1500。点击Calculate开始运行计算。运行计算开始后，可以在视图窗口中看到残差曲线和进口平均压力曲线随迭代步数的变化。ANSYS

Confidential31©

2018

ANSYS,

Inc.August

30,2023操作步骤ANSYS

Confidential32©

2018

ANSYS,

Inc.August

30,2023操作步骤计算进行了364步迭代后弹出Information窗口显示计算已经完成。点击OK完成计算。计算所需要进行的迭代步数由以下三个因素共同决定：1. Convergence

Conditions中定义的平均值的监视条件；2. Residual

Monitors中定义的残差曲线的收敛标准；3. Run

Calculation中定义的Number

of

Iterations。以上三个条件任意一个达到满足，计算即停止。当前的计算是因为达到最小残差量级而停止。ANSYS

Confidential33©

2018

ANSYS,

Inc.August

30,2023操作步骤3.12.Post-processing后处理在最左侧的树中，使用【Results】中的相关工具，对计算结果进行后处理。3.12.1.检查进出口流量平衡树中【Results】>【Reports】>【Fluxes】，鼠标左键双击Fluxes，弹出FluxReports窗口。保持默认选项Options中Mass

Flow

Rate，选择Boundary中inlet、outlet，点击Compute按钮显示这些边界上的流量。在NetResults(kg/s)中计算得到进出口净流量为-0.01997364

kg/s，通常计算域内流入流出差值与总流量的比值小于

0.5%可认为质量守恒，此处远小于0.5%可认为本次计算质量流量平衡。ANSYS

Confidential34©

2018

ANSYS,

Inc.August

30,2023操作步骤3.12.2.显示叶轮、轮毂面上的压力树中【Results】>【Graphics】>【Contours】，鼠标左键双击Contours，弹出

Contours窗口，在Options中勾选Filled，保持Contours

of中Pressure…>StaticPressure选项，在Surfaces对话栏中选择blade、hub，点击Save/Display显示压力分布。ANSYS

Confidential35©

2018