1、ANSYS,Inc.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitating Centrifugal PumpWS5-1April 28,2009Inventory#002599IntroductionThe Purpose of the tutorial is to model cavitation in a centrifugal pump,which involves the use of a rotation domain and the cavitation model.The p
2、roblem consists of a five blade centrifugal pump operating at 2160 rpm.The working fluid is water and flow is assumed to be steady and incompressible.Due to rotational periodicity a single blade passage will be modeled.The initial flow-field will be solved without cavitation.It will be turned on lat
3、er.ANSYS,Inc.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitating Centrifugal PumpWS5-2April 28,2009Inventory#0025991.Start Workbench and save the project as centrifugalpump.wbpj2.Drag CFX into the Project Schematic from the Component Systems toolbox3.Start CFX-Pre by dou
4、ble clicking Setup4.When CFX-Pre opens,import the mesh by right-clicking on Mesh and selecting Import Mesh ICEM CFD5.Browse to pump.cfx56.Keep Mesh units in m7.Click OpenWorkbench ANSYS,Inc.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitating Centrifugal PumpWS5-3April 28
5、,2009Inventory#002599Modifying the material properties:1.Expand Materials in the Outline tree2.Double-click Water 3.On the Material Properties tab change Density to 1000 kg/m3 4.Change Dynamic Viscosity to 0.001 kg m-1 s-1 under Transport Properties5.Click OKCreating Working FluidsANSYS,Inc.Propriet
6、ary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitating Centrifugal PumpWS5-4April 28,2009Inventory#002599Setting up the Fluid Domain1.Double-click on Default Domain2.Under Fluid and Particle Definitions,delete Fluid 1 and then create a new Fluid named Water Liquid3.Set Material to
7、Water4.Create another new Fluid named Water Vapour5.Next to the Material drop-down list,click the“”icon,then the Import Library Data icon(on the right of the form),and select Water Vapour at 25 C under the Water Data objectClick OK6.Back in the Material panel,select Water Vapour at 25 CClick OKANSYS
8、,Inc.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitating Centrifugal PumpWS5-5April 28,2009Inventory#002599Setting up the Fluid Domain7.Set the Reference Pressure to 0 Pa8.Set Domain Motion to Rotating9.Set Angular Velocity to 2160 rev min-110.Switch on Alternate Rotatio
9、n Model11.Make sure Rotation Axis under Axis Definition is set to Global Z11.Switch to the Fluid Models tab,and set the following:12.Turn on Homogeneous Model in the Multiphase section13.Under Heat Transfer set the Option to Isothermal,with a Temperature of 25 C14.Set Turbulence Option to Shear Stre
10、ss Transport15.Click OKANSYS,Inc.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitating Centrifugal PumpWS5-6April 28,2009Inventory#002599Inlet Boundary Condition1.Insert a boundary condition named Inlet2.On the Basic Settings tab,set Boundary Type to Inlet3.Set Location to
11、 INLET4.Set Frame Type to Stationary5.Switch to the Boundary Details tab6.Specify Mass and Momentum with a Normal Speed of 7.0455 m/s7.Switch to the Fluid Values tab8.For Water Liquid,set the Volume Fraction to a Value of 19.For Water Vapour,set the Volume Fraction to a Value of 010.Click OKANSYS,In
12、c.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitating Centrifugal PumpWS5-7April 28,2009Inventory#002599Outlet Boundary Condition1.Inset a boundary condition named Outlet2.On the Basic Settings tab,set Boundary Type to Opening3.Set Location to OUT4.Set Frame Type to Stat
13、ionary5.Switch to the Boundary Details tab6.Specify Mass and Momentum using Entrainment,and enter a Relative Pressure of 600,000 Pa7.Enable the Pressure Option and set it to Opening Pressure8.Set Turbulence Option to Zero Gradient9.Switch to the Fluid Values tab10.For Water Liquid,set the Volume Fra
14、ction to a Value of 111.For Water Vapour,set the Volume Fraction to a Value of 012.Click OKANSYS,Inc.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitating Centrifugal PumpWS5-8April 28,2009Inventory#002599Periodic Interface1.Click to create an Interface,and name it Periodi
15、c2.Set the Interface Type to Fluid Fluid3.For Interface Side 1,set the Region List to DOMAIN INTERFACE 1 SIDE 1 and DOMAIN INTERFACE 2 SIDE 1(use the“”icon and the Ctrl key)4.For Interface Side 2,set the Region List to DOMAIN INTERFACE 1 SIDE 2 and DOMAIN INTERFACE 2 SIDE 25.Set the Interface Models
16、 option to Rotational Periodicity6.Under Axis Definition,select Global Z7.Set Mesh Connection Option to 1:18.Click OKANSYS,Inc.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitating Centrifugal PumpWS5-9April 28,2009Inventory#002599Wall Boundary Conditions1.Insert a boundar
17、y condition named Stationary2.Set it to be a Wall,using the STATIONARY location3.On the Boundary Details tab,enable a Wall Velocity and set it to Counter Rotating Wall4.Click OK5.In the Outline Tree,right-click on the Default Domain Default boundary and rename it to MovingThe default behavior for th
18、e Moving boundary condition is to move with the rotating domain,so there is nothing that needs to be setANSYS,Inc.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitating Centrifugal PumpWS5-10April 28,2009Inventory#002599Initialization1.Click to initialize the solution2.On t
19、he Fluid Settings form,set Water Liquid Volume Fraction to Automatic with Value,and set the Volume Fraction to 13.Set Water Vapour Volume Fraction to Automatic with Value,and set the Volume Fraction to 04.Click OKANSYS,Inc.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitat
20、ing Centrifugal PumpWS5-11April 28,2009Inventory#002599Solver Control1.Double click Solver Control in the Outline tree2.Set Timescale Control to Physical timescaleA commonly used timescale in turbomachinery is 1/omega,where omega is the rotation rate in radians per second.You can use an expression t
21、o determine a timestep from this.In this case,2/omega will be used to achieve faster convergence.3.Enter the following expression in the Physical Timescale box:1/(pi*2160 min-1)4.Set Residual Target to 1e-55.On the Advanced Options tab,turn on Multiphase Control,then turn on Volume Fraction Coupling
22、 and set the Option to Coupled6.Click OKANSYS,Inc.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitating Centrifugal PumpWS5-12April 28,2009Inventory#002599Output Control1.Double Click on Output Control in the Outline tree2.On the Monitor tab,turn on Monitor Options3.Under
23、Monitor Points and Expressions,create a new object and call it InletPTotalAbs4.Set Option to Expression5.Specify the following expression:massFlowAve(Total Pressure in Stn Frame)Inlet6.Create a new object called InletPStatic,and set Option to Expression 7.Specify the following expression:areaAve(Pre
24、ssure)Inlet8.Click OKANSYS,Inc.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitating Centrifugal PumpWS5-13April 28,2009Inventory#002599Solver1.Close CFX-Pre and switch to the Workbench Project window2.Save the project3.Now double click on Solution in the Project Schematic
25、 to start the Solver Manager4.When the Solver Manager opens,click Start Run5.When the solution has completed,close the Solver Manager and return to the Project window6.Save the projectANSYS,Inc.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitating Centrifugal PumpWS5-14Apr
26、il 28,2009Inventory#002599Post-processing1.View the results in CFD-Post by double clicking Results in the Project Schematic2.Insert a Contour by clicking3.For the Location,click ,expand Regions and then select BLADE4.Set Variable to Absolute Pressure from the extended list5.Set Range to Global6.On t
27、he Render tab switch off Lighting and Show contour Lines7.Click ApplyANSYS,Inc.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitating Centrifugal PumpWS5-15April 28,2009Inventory#002599Post-processing9.Insert another Contour on the HUB location,using the variable Absolute P
28、ressure coloured by Local Range.Turn off Lighting and Show Contour Lines.10.Insert another Contour on the SHROUD location,using the variable Absolute Pressure coloured by Local Range.Turn off Lighting and Show Contour Lines.The minimum pressure is above the Saturation Pressure of 2650 Pa for Water h
29、ere.In the next step,the outlet pressure will be reduced enough to initiate Cavitation.ANSYS,Inc.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitating Centrifugal PumpWS5-16April 28,2009Inventory#002599Adding another Analysis1.Close CFD-Post and return to the Project Schem
30、atic2.Click the arrow next to the A cell and select DuplicateA new CFX project is created as a copy of the first3.Change the name of the new Simulation to Cavitation4.Use the arrow next to the A cell to Rename it to No Cavitation5.Save the Project6.Double-click Setup for the Cavitation simulation to
31、 open CFX-PreANSYS,Inc.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitating Centrifugal PumpWS5-17April 28,2009Inventory#002599Physics Modifications1.Edit the Default Domain2.On the Fluid Pair Models tab set Mass Transfer to Cavitation 3.Set Option to Rayleigh Plesset4.Tu
32、rn on Saturation Pressure5.Set a Saturation Pressure of 2650 Pa6.Click OK7.Edit the Outlet Boundary Condition8.On the Boundary Details tab,set the Relative Pressure to 300,000 Pa9.Click OKMost cavitation solutions should be performed by turning cavitation on and then successively lowering the system
33、 pressure over several runs to more gradually induce cavitation.To speed up this workshop,a sudden change in pressure is introduced.Note that this approach may not be suitable for modelling some industrial cases.ANSYS,Inc.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitati
34、ng Centrifugal PumpWS5-18April 28,2009Inventory#002599Physics Modifications1.Edit Solver Control2.Set the Max.Iterations to 1503.Set the Residual Target to 1e-44.Click OK5.Close CFX-Pre and save the project6.In the Project Schematic,drag cell A3 onto cell B3The non-cavitating solution will be used a
35、s the initial guess for the cavitating solution7.Double-click Solution for the Cavitation systemIn the Solver Manager note that the initial conditions have been provided from the project schematic8.Click Start RunANSYS,Inc.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitat
36、ing Centrifugal PumpWS5-19April 28,2009Inventory#002599Cavitation SolutionThere is a significant spike in residuals,in part due to the outlet pressure difference,but also due to the fact that the absolute pressure is low enough to induce cavitation.1.When the run completes,close the Solver Manager a
37、nd return to the Project Schematic2.Save the project3.Double-click Results for the Cavitation project to openCFD-PostANSYS,Inc.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitating Centrifugal PumpWS5-20April 28,2009Inventory#002599Post-processing1.If it is not enabled,tur
38、n on visibility for the Wireframe and turn off visibility for any User Locations and Plots2.Create an XY Plane at Z=0.01 m3.Colour it by Absolute Pressure(the variable is available in the Extended List by clicking ).Use a Global RangeThe minimum absolute pressure is equivalent to the Saturation Pres
39、sure specified earlier,which is a strong hint that some cavitation has occurred4.Change the Colour Variable to Water Vapour.Volume Fraction5.Change the Colour Map to Blue to WhiteANSYS,Inc.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitating Centrifugal PumpWS5-21April 28
40、,2009Inventory#002599Post-processing1.Turn off visibility for Plane 12.Create a Volume using the Isovolume method3.Set the Variable to Water Vapour.Volume Fraction4.Set Mode to Above Value,and enter a value of 0.55.To view 360 degrees of the model,double-click Default Transform6.Uncheck Instancing I
41、nfo from Domain7.Set#of copies to 58.Set#of Passages to 59.Click OKANSYS,Inc.Proprietary 2009 ANSYS,Inc.All rights reserved.Workshop SupplementWS5:Cavitating Centrifugal PumpWS5-22April 28,2009Inventory#002599Post-processingThe main area of cavitation exists between the suction side of the blade and
42、 the shroud in this geometry.A secondary area of cavitation is just behind the leading edge of the blade on the pressure sideFurther steps to try:1.Calculate torque on the BLADE using the function calculator(hint,use the extended region list to find the BLADE,and use Global Z axis)2.Plot velocity Ve
43、ctors on Plane 1,using the variableWater Liquid.Velocity in Stn.Frame3.Calculate the mass flow through the pump(hint:use the function calculator to evaluate massFlow at the Outlet region)4.Using a similar method to step 2,calculate the drop in Total Pressure from Inlet to Outlet5.Plot Streamlines,starting from the Inlet location
