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

Surge simulation

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 25, 2016, 16:55
Default Surge simulation
  #1
New Member
 
Piotr
Join Date: Jan 2016
Posts: 6
Rep Power: 3
piotr1987 is on a distinguished road
Hi !
I am going to perform surge simulation of centrifugal impeller with vaneless diffuser. I managed to calculate steady case and noticed big turbulance inside impeller. I also tried harmonics unsteady model, however it seems to be just for rotor/stator interaction (correct me if i am wrong) - results don't change with time. Is it possible to model surge in NUMECA turbo ? Which unsteady model should I use ? Does anybody have tutorial of such case ?
Best regards,
Piotr
piotr1987 is offline   Reply With Quote

Old   March 28, 2016, 05:19
Default
  #2
Member
 
DarylMusashi's Avatar
 
Holger Dietrich
Join Date: Apr 2011
Location: Germany
Posts: 92
Rep Power: 7
DarylMusashi is on a distinguished road
Hi Piotr,
you are right, non-linear harmonics unsteady method does not work in your case, because you have no rotor-stator or rotor-rotor combination(s) in your machine.
Besides the NLH-method you can use the standard unsteady method.
DarylMusashi is offline   Reply With Quote

Old   March 29, 2016, 02:15
Default
  #3
New Member
 
Jessica Wei
Join Date: Mar 2016
Posts: 6
Rep Power: 2
wkjshon is on a distinguished road
Hi,

Thank you for your sharing.
Do you by any chance know if the NLH-method could work in a compressor with a volute, or it may only work in blade rows?



Quote:
Originally Posted by DarylMusashi View Post
Hi Piotr,
you are right, non-linear harmonics unsteady method does not work in your case, because you have no rotor-stator or rotor-rotor combination(s) in your machine.
Besides the NLH-method you can use the standard unsteady method.
wkjshon is offline   Reply With Quote

Old   March 29, 2016, 16:52
Default
  #4
Member
 
DarylMusashi's Avatar
 
Holger Dietrich
Join Date: Apr 2011
Location: Germany
Posts: 92
Rep Power: 7
DarylMusashi is on a distinguished road
It works in a volute, too (compressor and turbine). But make sure to mesh just one blade passage of each row. The blade passage of each row needs periodic boundaries in circumferential direction. If you are interested in the background of this please ask, I could try to draw a sketch to explain the reason.
DarylMusashi is offline   Reply With Quote

Old   April 6, 2016, 15:45
Default
  #5
New Member
 
Piotr
Join Date: Jan 2016
Posts: 6
Rep Power: 3
piotr1987 is on a distinguished road
Quote:
Originally Posted by DarylMusashi View Post
Hi Piotr,
you are right, non-linear harmonics unsteady method does not work in your case, because you have no rotor-stator or rotor-rotor combination(s) in your machine.
Besides the NLH-method you can use the standard unsteady method.
Daryl,
Should I use domain scalling method or phased-lagged one ?
piotr1987 is offline   Reply With Quote

Old   April 6, 2016, 16:59
Default
  #6
Member
 
DarylMusashi's Avatar
 
Holger Dietrich
Join Date: Apr 2011
Location: Germany
Posts: 92
Rep Power: 7
DarylMusashi is on a distinguished road
Hi Piotr,
the advantage of the phase lagged method is a reduced CPU time needed for unsteady computations. But it is only allowed for one-stage machines (stator-rotor) or 1.5-stage machines with the same number of stators.

Domain scaling is the only available unsteady method for you then.
DarylMusashi is offline   Reply With Quote

Old   April 7, 2016, 03:10
Default
  #7
New Member
 
Piotr
Join Date: Jan 2016
Posts: 6
Rep Power: 3
piotr1987 is on a distinguished road
Quote:
Originally Posted by DarylMusashi View Post
Hi Piotr,
the advantage of the phase lagged method is a reduced CPU time needed for unsteady computations. But it is only allowed for one-stage machines (stator-rotor) or 1.5-stage machines with the same number of stators.

Domain scaling is the only available unsteady method for you then.

Daryl,
I am really grateful for your advices. I haven't checked Rotating Boundary Condition in Boundary Condition option. I used steady result for initial solution. In control variables i have to input physical time step. In NUMECA manual i found following formula:

As i understand omega is rotational speed (in my case 9000RPM), Nperiods (number of blades=21).
Could you please tell me how to calculate NOFROT ?
piotr1987 is offline   Reply With Quote

Old   April 7, 2016, 16:26
Default
  #8
Member
 
DarylMusashi's Avatar
 
Holger Dietrich
Join Date: Apr 2011
Location: Germany
Posts: 92
Rep Power: 7
DarylMusashi is on a distinguished road
Hi Piotr,
you are right, because you have no rotor-stator interface you have to calculate the physical timestep size manually. The following questions lead to the desired duration of the physical time step.

1. What is your rotational speed in [rad/s]?
2. How far in [rad] must one blade rotate for one blade passage?
3. How long does this take?
4. With how many time steps do you want to resolve one blade passage?

If you are lazy have a look at the attached file.
Attached Images
File Type: jpg physical_timestep.JPG (50.2 KB, 16 views)

Last edited by DarylMusashi; April 9, 2016 at 16:51.
DarylMusashi 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
Simulation of high pressure diesel injector - all phases compressible with cavitation fivos CFX 4 July 30, 2015 06:48
Huge file sizes when Running VOF simulation aarratia FLUENT 0 May 8, 2014 12:27
Simulation of a complex wing in solidworks flow simulation niels1900 FloEFD, FloWorks & FloTHERM 6 April 20, 2011 10:44
GUI crash and simulation engine still running RPJones FLOW-3D 2 November 9, 2010 09:18
Continuous vs interrupted simulation sega OpenFOAM Running, Solving & CFD 4 November 3, 2008 15:29


All times are GMT -4. The time now is 10:02.