# Transient axial rotor/stator convergence issue?

 Register Blogs Members List Search Today's Posts Mark Forums Read

March 16, 2010, 12:12
Transient axial rotor/stator convergence issue?
#1
New Member

Nicola Viscanti
Join Date: Apr 2009
Posts: 5
Rep Power: 8
Hi all!

I made a transient analysis of an axial rotor/stator. First of all I made a steady analysis with frozen rotor interface and I reached a good convergence after 250 iterations (between 10^-5 and 10^-6 RMS).
Then I launched the transient analysis using the steady results like initial values.
I have the same rotor/stator blades number so I used 1:1 and I have a pitch ratio that's exactly 1.
I set the total time = (2pi/rotor blades)/omega where omega is rotational speed in rad/s and I used 50 time steps.
1) max coeff. loop =10 (default)
2) max coeff. loop =30

B.C. are:
Inlet -> Total pressure and Total temperature, turbulence intensity 5%
Outlet-> Average static pressure
k-epsilon

I enclose 2 plots:
1) massflow goes through the stage
2) axial force on rotor blade

I expected a constant massflow. What are oscillationts in the last time steps? Did my transient analysis reach convergence?
I expected a periodic behaviour for axial force (the same value or very similar at 250 and 300 timestep) because I have pitch ratio exactly 1. I'm right? Plots are the same using max coeff. loop 10 or 30

Thank you very much!
Attached Images
 Massflow.jpg (51.8 KB, 46 views) Axial_force.jpg (62.5 KB, 32 views)

 March 16, 2010, 13:41 #2 New Member   Join Date: Apr 2009 Posts: 14 Rep Power: 8 Your calculation probably hasn't been stabilized. It takes several cycles to develop a repeated pattern. By the way, specifying pressure at both inlet and outlet is not a robust setup.

 March 16, 2010, 18:14 #3 New Member   Join Date: Jun 2009 Posts: 8 Rep Power: 8 Hi Nicola, I am having a bit of trouble understading the set-up that you have discribed. 1. You say that you have the same number of blades on the rotor and the stator how many is that? 2. Are you modeling a full sector or a full-wheel? 3. The calcuation: 2*Pi [rad/revolution]/Omega[rad/s]/blades[blade/revolution] = [s/blade] = or blade period.....but you called this the total time? 4. What is the 50 steps for? 1 revolution or 1 blade pass? 5. Your figure show iteration...should this be time step? i usually set one of these problems up with 10-20 time steps per blade pass. 3-5 iterations per time step. Things will begin to develop into a regular period after a few revolutions. If you do not get good convergence (or accuracy of the transient features) then half the time step. The time step requirement depends on the macine and the strength of the coupling between the rotor and stator. Very sharp features such as shocks, separations, etc. from an upstream component need to be resolved with smaller time steps than "smoother interactions". Why would you expect a constant mass flow rate with Ptotal in and Pstatic out on the stage? Think of it this way: as you solve the transient interaction you are clocking and then solving the full transient equations...then clocking again right? What would you get if you just clocked and solved to steady state and then clocked and solved again......slightly different mass flows because it is a different configurations....right? For boundary conditions it depends what operating point you are trying to solve but try to set boundary conditions for what is "most constant". On the flat part of the curve mass flow inlet and statc pressure outlet works well.....on the steep part and near choking (of it is compressible) then Ptotal inlet and Pstatic out works well. If you are having problems solving at a difficult point try solving an easier one for a few rotations and then "moving the solution over" to the point you want. All of the flow dynamics associated with rotation will be developed and you only need to solve for what is additionally challenging due to the rotor-stator interaction at the desired point. I hope this helps, let us know how it works out..........Bak_Flow

 March 17, 2010, 05:15 #4 New Member   Nicola Viscanti Join Date: Apr 2009 Posts: 5 Rep Power: 8 Hi bak_flow, I have 73 stator blades and 73 rotor blades so I used a stationary domain made by 1 stator blade and periodic surfaces and 1 rotor blade with periodic surfaces. Pitch ratio is exactly 1! Rotational speed=314.16 rad/s Rotor pitch modelled = 1*(2pi/73)= 0.086 rad Time to pass through 1 pitch= 0.086 rad/ 314.16 rad/s=2.7397e-04 s time steps 50 -> 2.7397e-04 /50=5.4794e-06 s It's a HP steam turbine. Mach numbers are low! Why is massflow different at 250 and 300 time steps? The rotor blade is in the same position! I'm not sure about massflow has to be constant for all timesteps but I guess I must have the same values at 250 and 300! Why is axial force different at 250 and 300 time steps? Thanks

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post siw CFX 5 October 30, 2010 05:45 icesniffer CFX 1 August 8, 2009 07:25 Sonja CFX 0 December 8, 2008 09:34 Adam CFX 1 April 12, 2007 11:34 mahesh prakash Main CFD Forum 1 January 21, 1999 14:45

All times are GMT -4. The time now is 05:49.