CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > CFX

Transient axial rotor/stator convergence issue?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 16, 2010, 12:12
Default Transient axial rotor/stator convergence issue?
  #1
New Member
 
Nicola Viscanti
Join Date: Apr 2009
Posts: 5
Rep Power: 7
Nicola Viscanti is on a distinguished road
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.
I made 2 analysis:
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

Please any suggest is welcome!

Thank you very much!
Attached Images
File Type: jpg Massflow.jpg (51.8 KB, 46 views)
File Type: jpg Axial_force.jpg (62.5 KB, 32 views)
Nicola Viscanti is offline   Reply With Quote

Old   March 16, 2010, 13:41
Default
  #2
New Member
 
Join Date: Apr 2009
Posts: 14
Rep Power: 7
longbow is on a distinguished road
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.
longbow is offline   Reply With Quote

Old   March 16, 2010, 18:14
Default
  #3
New Member
 
Join Date: Jun 2009
Posts: 8
Rep Power: 7
Bak_Flow is on a distinguished road
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
Bak_Flow is offline   Reply With Quote

Old   March 17, 2010, 05:15
Default
  #4
New Member
 
Nicola Viscanti
Join Date: Apr 2009
Posts: 5
Rep Power: 7
Nicola Viscanti is on a distinguished road
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
Nicola Viscanti is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Best practice for transient simulations? siw CFX 5 October 30, 2010 05:45
transient simulation of a rotating rectangle icesniffer CFX 1 August 8, 2009 07:25
transient statistics for turbo components Sonja CFX 0 December 8, 2008 09:34
Transient vs Steady State Adam CFX 1 April 12, 2007 11:34
difference between false and true transient mahesh prakash Main CFD Forum 1 January 21, 1999 14:45


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