Tutorial 11. Using Sliding Meshes Introduction The analysis of turbomachinery often involves the examination of the unsteady effects due to flow interaction between the stationary components and the rotating blades. In this tutorial, the sliding mesh capability of FLUENT is used to to analyze the unsteady flow in an axial compressor stage. The rotor-stator interaction is modeled by allowing the mesh associated with the rotor blade row to rotate relative to the stationary mesh associated with the stator blade row. This tutorial demonstrates how to do the following: • Create periodic zones. • Set up the unsteady solver and boundary conditions for a sliding mesh simulation. • Set up the grid interfaces for a periodic sliding mesh model. • Sample the time-dependent data and view the mean value. Prerequisites This tutorial assumes that you are familiar with the menu structure in FLUENT and that you have completed Tutorial 1. Some steps in the setup and solution procedure will not be shown explicitly. Problem Description The model represents a single-stage axial compressor comprised of two blade rows. The first row is the rotor with 16 blades, which is operating at a rotational speed of 37,500 rpm. The second row is the stator with 32 blades. The blade counts are such that the domain is rotationally periodic, with a periodic angle of 22.5 degrees. This allows you to model only a portion of the geometry, namely, one rotor blade and two stator blades. Due to the high Reynolds number of the flow and the relative coarseness of the mesh (both blade rows are comprised of only 13,856 cells total), the analysis will employ the inviscid model, so that FLUENT is solving the Euler equations. c Fluent Inc. January 31, 2007 11-1
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Tutorial 11. Using Sliding Meshes
Introduction
The analysis of turbomachinery often involves the examination of the unsteady effectsdue to flow interaction between the stationary components and the rotating blades. Inthis tutorial, the sliding mesh capability of FLUENT is used to to analyze the unsteadyflow in an axial compressor stage. The rotor-stator interaction is modeled by allowingthe mesh associated with the rotor blade row to rotate relative to the stationary meshassociated with the stator blade row.
This tutorial demonstrates how to do the following:
• Create periodic zones.
• Set up the unsteady solver and boundary conditions for a sliding mesh simulation.
• Set up the grid interfaces for a periodic sliding mesh model.
• Sample the time-dependent data and view the mean value.
Prerequisites
This tutorial assumes that you are familiar with the menu structure in FLUENT and thatyou have completed Tutorial 1. Some steps in the setup and solution procedure will notbe shown explicitly.
Problem Description
The model represents a single-stage axial compressor comprised of two blade rows. Thefirst row is the rotor with 16 blades, which is operating at a rotational speed of 37,500rpm. The second row is the stator with 32 blades. The blade counts are such that thedomain is rotationally periodic, with a periodic angle of 22.5 degrees. This allows youto model only a portion of the geometry, namely, one rotor blade and two stator blades.Due to the high Reynolds number of the flow and the relative coarseness of the mesh(both blade rows are comprised of only 13,856 cells total), the analysis will employ theinviscid model, so that FLUENT is solving the Euler equations.
1. Download sliding_mesh.zip from the Fluent Inc. User Services Center or copyit from the FLUENT documentation CD to your working folder (as described inTutorial 1).
2. Unzip sliding_mesh.zip.
axial comp.msh can be found in the sliding mesh folder created after unzippingthe file.
3. Start the 3D (3d) version of FLUENT.
Step 1: Grid
1. Read in the mesh file axial comp.msh.
File −→ Read −→Case...
2. Check the grid.
Grid −→Check
FLUENT will perform various checks on the mesh and will report the progress inthe console. Pay particular attention to the minimum volume, and make sure thisis a positive number.
Warnings will be displayed regarding unassigned interface zones, resulting in thefailure of the grid check. You do not need to take any action at this point, as thisissue will be rectified when you define the grid interfaces in a later step.
(a) Retain the default selections from the Surfaces selection list.
(b) Click Display and close the Grid Display panel.
(c) Rotate the view to get the display shown in Figure 11.2.
ZY
X
GridFLUENT 6.3 (3d, pbns, lam)
Figure 11.2: Rotor-Stator Outline Display
The inlet to the rotor mesh is colored blue, the interface between the rotor and statormeshes is colored yellow, and the outlet of the stator mesh is colored red.
5. Use the text user interface to change zones rotor-per-1 and rotor-per-3 fromwall zones to periodic zones.
(a) Press <Enter> in the console to get the command prompt (>).
1. Specify air (the default material) as the fluid material, using the ideal gas law tocompute density.
Define −→Materials...
(a) Retain the default entry of air in the Name text entry field.
(b) Select ideal-gas from the Density drop-down list in the Properties group box.
(c) Retain the default values for all other properties.
(d) Click Change/Create and close the Materials panel.
As reported in the console, FLUENT will automatically enable the energy equation,since this is required when using the ideal gas law to compute the density of amaterial.
(b) Click OK to close the Operating Conditions panel.
Since you have set the operating pressure to zero, you will specify the boundarycondition inputs for pressure in terms of absolute pressures when you define themin a later step. Boundary condition inputs for pressure should always be relative tothe value used for operating pressure.
4. Set the boundary conditions for the outlet (stator-outlet).
(a) Enter 1.08 atm for the Gauge Pressure.
(b) Enable the Radial Equilibrium Pressure Distribution option.
(c) Click the Thermal tab and enter 288 K for Backflow Total Temperature.
(d) Click OK to close the Pressure Outlet panel.
Note: The momentum settings and temperature you input at the pressure outletwill be used only if flow enters the domain through this boundary. It is impor-tant to set reasonable values for these downstream scalar values, in case flowreversal occurs at some point during the calculation.
5. Retain the default boundary conditions for all wall zones.
Note: For wall zones, FLUENT always imposes zero velocity for the normal velocitycomponent, which is required whether or not the fluid zone is moving. Thiscondition is all that is required for an inviscid flow, as the tangential velocityis computed as part of the solution.
1. Create a periodic grid interface between the rotor and stator mesh regions.
Define −→Grid Interfaces...
(a) Enter int for Grid Interface.
(b) Enable Periodic in the Interface Type group box.
Enabling this option, allows FLUENT to treat the interface between the slidingand non-sliding zones as periodic where the two zones do not overlap.
(c) Select rotor-interface from the Interface Zone 1 selection list.
Note: In general, when one interface zone is smaller than the other, it isrecommended that you choose the smaller zone as Interface Zone 1. Inthis case, since both zones are approximately the same size, the order isnot significant.
(d) Select stator-interface from the Interface Zone 2 selection list.
(e) Click Create and close the Grid Interfaces panel.
2. Check the grid again to verify that the warnings displayed earlier have been re-solved.
6. Run the calculation for one revolution of the rotor.
Solve −→Iterate...
(a) Enter 6.6666e-6 s for the Time Step Size.
This time step represents the length of time during which the rotor will rotate1.5 degrees. Since the periodic angle of the rotor is 22.5 degrees, the passingperiod of the rotor blade will equal 15 time steps, and a complete revolution ofthe rotor will take 240 time steps.
(b) Enter 240 for the Number of Time Steps.
(c) Retain the default setting of 20 for the Max Iterations per Time Step in theIteration group box.
(d) Click Iterate.
The solution will converge in approximately 3100 iterations.
The residuals jump at the beginning of each time step and then fall at least two tothree orders of magnitude. Also, the relative convergence criteria is achieved beforereaching the maximum iteration limit (20) for each time step, indicating the limitdoes not need to be increased.
7. Examine the monitor histories for the first revolution of the rotor (Figures 11.4,11.5, and 11.6).
Figure 11.6: Static Pressure at the Interface During the First Revolution
The monitor histories show that the large variations in flow rate and interface pressurethat occur early in the calculation are greatly reduced as time-periodicity is approached.
8. Save the case and data files (axial comp-0240.cas and axial comp-0240.dat).
File −→ Write −→Case & Data...
! When the sliding mesh model is used, you must save a case file whenever adata file is saved. This is because the case file contains the grid information,which is changing with time.
Note: For unsteady-state calculations, you can add the character string %t to thefile name so that the iteration number is automatically appended to the name(e.g., by entering axial comp-%t for the File Name in the Select File dia-log box, FLUENT will save files with the names axial comp-0240.cas andaxial comp-0240.dat).
9. Rename the monitor file names in preparation for further iterations.
By saving the monitor histories under a new file name, the range of the axes willautomatically be set to show only the data generated during the next set of iterations.This will scale the plots so that the fluctuations are more visible.
(a) Click the Define... button for monitor-1 to open the Define Surface Monitorpanel.
i. Enter monitor-1b.out for the File Name.
ii. Click OK to close the Define Surface Monitor panel.
(b) Similarly, change the File Name for monitor-2 and monitor-3 to be monitor-2b.outand monitor-3b.out respectively.
(c) Click OK to close the Surface Monitors panel.
Extra: Instead of creating a new file for each monitor, you could have adjustedthe ranges of the axes to make the fluctuations visible, and then allowed thedata from the next set of iterations to be appended onto the original monitorfiles. To do this, click the Axes... button in the Define Surface Monitor panelof each monitor, disable the Auto Range option, enter in appropriate values inthe Range group box, and click Apply.
10. Continue the calculation for 720 more time steps to simulate three more revolutionsof the rotor.
Solve −→Iterate...
! Calculating three more revolutions will require significant CPU re-sources. Instead of calculating the solution, you can read a data file(axial comp-0960.dat.gz) with the precalculated solution for this tu-torial. This data file can be found in the folder where you found the meshfile.
The calculation will run for approximately 9,100 more iterations.
11. Examine the monitor histories for the next three revolutions of the rotor to verifythat the solution is time-periodic (Figures 11.7, 11.8, and 11.9).
Note: If you read the provided data file instead of iterating the solution for threerevolutions, the monitor histories can be displayed by using the Plot/File...menu option. Simply click the Add button in the File XY Plot panel, select oneof the monitor histories in the Select File dialog box, click OK, and then clickPlot.
Figure 11.9: Static Pressure at the Interface During the Next 3 Revolutions
Extra: Note that the y-axis for Figure 11.7 does not show enough significant figuresto fully display the values of the mass flow rate. You can display the full valuesby using the Plot/File... menu option. Use the Add... button in the File XY Plotpanel to select monitor-1b.out, and click the Axes... button to open the Axes -File XY Plot panel. There, select Y from the Axis list, enter 6 for Precision,and click Apply. Finally, click Plot in the File XY Plot panel to redisplay themonitor history.
12. Save the case and data files (axial comp-0960.cas and axial comp-0960.dat).
File −→ Write −→Case & Data...
13. Change the File Name for monitor-1, monitor-2, and monitor-3 to be monitor-1c.out,monitor-2c.out, and monitor-3c.out, respectively (as described in a previousstep), in preparation for further iterations.
In the next two steps you will examine the time-averaged values for the mass flow ratesat the inlet and the outlet during the final revolution of the rotor. By comparing thesevalues, you will verify the conservation of mass on a time-averaged basis for the systemover the course of one revolution.
1. Examine the time-averaged mass flow rate at the inlet during the final revolutionof the rotor (as calculated from monitor-1c.out).
(b) Click the Plot/Modify Input Signal... button to open the Plot/Modify InputSignal panel.
i. Examine the values for Min, Max, Mean, and Variance in the Signal Statis-tics group box.
ii. Close the Plot/Modify Input Signal panel.
(c) Select the folder path ending in monitor-1c.out from the Files selection list.
(d) Click the Free File Data button.
2. Examine the time-averaged mass flow rate at the outlet during the final revolutionof the rotor (as calculated from monitor-2c.out), and plot the data.
Plot −→FFT...
(a) Click the Load Input File... button to open the Select File dialog box.
i. Select All Files from the Files of type drop-down list.
ii. Select monitor-2c.out from the list of files.
iii. Click OK to close the Select File dialog box.
(b) Click the Plot/Modify Input Signal... button to open the Plot/Modify InputSignal panel.
i. Examine the values for Min, Max, Mean, and Variance in the Signal Statis-tics group box.
Note that the outlet mass flow rate values correspond very closely withthose from the inlet, with the mean having approximately the same ab-solute value but with opposite signs. Thus, you can conclude that massis conserved on a time-averaged basis during the final revolution of therotor.
ii. Click the Set Defaults button.
iii. Click Apply/Plot to display the mass flow rate at the outlet (Figure 11.10).
(c) Select wall from the Surface Types selection list.
Scroll down the Surface Types selection list to find wall.
(d) Click Display and close the Contours panel.
(e) Rotate the view to get the display shown in Figure 11.11.
Contours of Mean Static Pressure (atm) (Time=7.9999e-03)FLUENT 6.3 (3d, dbns imp, unsteady)
1.34e+00
1.29e+00
1.23e+00
1.18e+00
1.12e+00
1.07e+00
1.01e+00
9.60e-01
9.06e-01
8.52e-01
7.98e-01
7.43e-01
6.89e-01
6.35e-01
5.81e-01
5.27e-01
4.72e-01
4.18e-01
3.64e-01
3.10e-01
2.55e-01 ZY
X
Figure 11.11: Mean Static Pressure on the Outer Shroud of the Axial Compressor
Shock waves are clearly visible in the flow near the outlets of the rotor and stator,as seen in the areas of rapid pressure change on the outer shroud of the axialcompressor.
This tutorial has demonstrated the use of the sliding mesh model for analyzing unsteadyrotor-stator interaction in an axial compressor stage. The model utilized the density-based solver in conjunction with the unsteady, dual-time stepping algorithm to computethe inviscid flow through the compressor stage. The solution was calculated over timeuntil the monitored variables displayed time-periodicity (which required several revolu-tions of the rotor), after which time-averaged data was collected while running the casefor the equivalent of one additional rotor revolution (240 time steps). The Fast FourierTransform (FFT) utility in FLUENT was employed to determine the time averages fromstored monitor data. Although not described in this tutorial, you can further use theFFT utility to examine the frequency content of the unsteady monitor data (in this case,you would observe peaks corresponding to the passing frequency and higher harmonicsof the passing frequency).
Further Improvements
This tutorial guides you through the steps to reach a second-order solution. You maybe able to obtain a more accurate solution by adapting the grid. Adapting the grid canalso ensure that your solution is independent of the grid. These steps are demonstratedin Tutorial 1.