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

HELP! Small Turbomachinery Diffuser Turbulence Modeling

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 27, 2012, 14:37
Default HELP! Small Turbomachinery Diffuser Turbulence Modeling
  #1
New Member
 
James Berlin
Join Date: Apr 2012
Posts: 3
Rep Power: 4
jkbg38@mail.missouri.edu is on a distinguished road
Hello, I am trying to model a small jet engine (70mm centrifugal compressor) in an attempt to evaluate diffuser performance to compare one design to another. The fluid volume is from the outlet of the compressor to the inlet of the combustor. I am having trouble getting the model to converge which I think is attributable to using the default turbulence parameters. Can anybody suggest a range of suitable parameters for turbulence and or any suggestions for the rest of the model?

My enabled models are as follows:
All y + Wall Treatment
SST (Menter) K Omega
K Omega Turbulence
Reynolds Averaged Navier Stokes
Turbulent
Coupled Energy
Ideal Gas
Coupled Flow
Gas
Steady
Three Dimensional


It is a rotational periodic interface about 24 degrees

Velocity Inlet (from compressor)
Flow Direction: (121.3,337.6,0) Cylindrical Coordinate frame
Static Temperature: 372.2239
Turbulence Intensity: 0.1
Turbulent Length Scale: 0.01m
Velocity Magnitude: 358.8m/s


Pressure Outlet (To Combustor)
Pressure: 0 (Reference Pressure 2.0942BAR)
Static Temperature: 409.861K
Turbulent Intensity 0.01
Turbulent Viscosity Ratio: 10
I tried Target Mass Flow Rate Option but it caused sudden divergence after 200 iterations

jkbg38@mail.missouri.edu is offline   Reply With Quote

Old   April 28, 2012, 01:45
Default
  #2
Senior Member
 
Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 442
Rep Power: 10
LuckyTran is on a distinguished road
Quote:
Originally Posted by jkbg38@mail.missouri.edu View Post
I tried Target Mass Flow Rate Option but it caused sudden divergence after 200 iterations
What version of Star-CCM are you using?

I did not find any glaring errors in your setup yet. However, your modelling is fairly complex with a coupled solver and ideal gas. How are you initializing the simulation? You will need a fairly good initialization for the simulation to not diverge.

I'm assuming you've already tried lowering the under relaxation factors.

Try running without all the complicated setups right off the bat. For starters I would turn off the target mass flow rate option initially and let the solution settle a bit. Next I would start with cutting back on ideal gas and use constant density. Help the solver as much as you can!
LuckyTran is offline   Reply With Quote

Old   April 28, 2012, 08:04
Default
  #3
Senior Member
 
Join Date: Oct 2009
Location: Germany
Posts: 637
Rep Power: 11
abdul099 is on a distinguished road
Again the most important question is: HOW DID YOU JUDGE CONVERGENCE? How long did you run it and how big is the model? Are your sure, your mesh is suitable to capture all important flow phenomena?

A lot of people complain about a not converging simulation but don't judge convergence the right way.
abdul099 is offline   Reply With Quote

Old   April 28, 2012, 15:50
Default
  #4
New Member
 
James Berlin
Join Date: Apr 2012
Posts: 3
Rep Power: 4
jkbg38@mail.missouri.edu is on a distinguished road
Lucky,
Im using Version 6.04. By initializing the solution do you mean initial conditions? I have:
pressure = 0
Static Temp = 300K
Turbulence Intensity = 0.01
Turbulent Velocity Scale = 1m/s
Turbulent Viscosity Ratio= 10
Velocity=0

I have not tried any under relaxation factors.

Abdul, the volume mesh is 2420377 cells, 15413936 faces and the residual plot is below. Also, a pressure scaler plot will continue to change as the iterations continue. As for if my mesh is suitable for capturing all important flow phenomena, unfortunately I have no idea, but I would appreciate some insight.
jkbg38@mail.missouri.edu is offline   Reply With Quote

Old   April 28, 2012, 16:56
Default
  #5
Senior Member
 
Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 442
Rep Power: 10
LuckyTran is on a distinguished road
Quote:
Originally Posted by jkbg38@mail.missouri.edu View Post
Lucky,
Im using Version 6.04. By initializing the solution do you mean initial conditions? I have:
pressure = 0
Static Temp = 300K
Turbulence Intensity = 0.01
Turbulent Velocity Scale = 1m/s
Turbulent Viscosity Ratio= 10
Velocity=0
Your initialization is terrible. They are nothing like what the flow is even at the inlet. Your flow is incoming at 358 m/s but you initialize it with 0 m/s? Your inlet static temperature is 370+ K and you initialized it with 300K. It is no surprise that it diverged, you tried to use target mass flow rate yet you initialized it with no flow!

The geometry is a bit complex, so giving a good initialization will be hard without writing your own field functions. But as I said before, you need to help the solver anyway you can.

You need to monitor the solution periodically and check that they even make sense.

And you should run it longer

Last edited by LuckyTran; April 29, 2012 at 10:38.
LuckyTran is offline   Reply With Quote

Old   April 28, 2012, 19:23
Default
  #6
Super Moderator
 
Ryne Whitehill
Join Date: Aug 2009
Posts: 313
Rep Power: 8
rwryne is on a distinguished road
Quote:
Originally Posted by LuckyTran View Post
Your initialization is terrible. They are nothing like what the flow is even at the inlet. Your flow is incoming at 358 m/s but you initialize it with 0 m/s? Your inlet static temperature is 370+ K and you initialized it with 300K. It is no surprise that it diverged, you tried to use target mass flow rate yet you initialized it with no flow!

The geometry is a bit complex, so giving a good initialization will be hard without writing your own field functions. But as I said before, you need to help the solver anyway you can.

You need to monitor the solution periodically and check that they even make sense.
This.

If it works, grid sequencing initialization is excellent. That said, it doesnt always work for me depending on my geometry.

Make sure you try ramps on your courant number & relaxation factors
rwryne is offline   Reply With Quote

Old   April 29, 2012, 07:38
Default
  #7
Senior Member
 
Join Date: Oct 2009
Location: Germany
Posts: 637
Rep Power: 11
abdul099 is on a distinguished road
Well, it doesn't look too bad. Of course, better initialisation would help, but it also should work without doing anything.

Maybe it just needs more than 1500 iterations? At least the residual plot doesn't look too bad since the residuals are not exploding but still slightly decreasing.
Maybe it's an unsteady flow while your're running steady? (There are not many really steady flows in the real world). That would cause the pressure plots to change forever, although important engineering values would be okay when averaging them over some iterations.

The target mass flow option is useless when using a velocity inlet. You fix the incoming mass to a specific value, calculated from inlet velocity, area of your inlet boundary and density. Now imagine, that doesn't meet the desired outlet mass flow (might be due to numerical issues since most of the values are just single precision). The solver might ramp the outlet pressure up (or down) to reach it's target mass flow, but that could never be achieved without violating the mass conservation. So it will continue to rise (or lower) the outlet pressure, giving you unphysical values or sudden divergence.

Therefore this option should only be used when using a pressure outlet in combination with a stagnation inlet or something similar, giving an additional degree of freedom to the solver. And even when using the target mass flow option, it shouldn't be switched on unless the flow pattern more or less meets the expected flow pattern. Otherwise it could get unstable.

For the initialisation, I wouldn't bother too much with this, since your geometry is too complicated to reliable predict the flow behaviour. It's not your job to know the solution to initialise the simulation, it's the solver's job to get the solution from some initial values. Of course, a better initialisation will help to solver to get the solution quicker and more stable, but for a steady flow, it should even work with a poor initialisation.
I would just ramp up the inlet velocity with a field function over some hundred iterations. That should help enough to stabilise the solution, but of course it will take longer to get a converged state.
Grid sequencing mentioned by rwryne might also be a good idea, but it's not guaranteed to work well for internal flows. It was desinged to be used mainly with external flows. The more complex and winding your domain the more propably it will not work.

Maybe you can set the static initialisation temperature to meet some expected average temperature in your domain. But don't bother too much about this. Ramping up the velocity and / or Courant number should help enough.

And please monitor some engineering values like inlet pressure, volume averaged turbulent kinetic energy etc. and wait until they are leveled out or start oscillating around a constant value while the residuals don't ddrop anymore. This will be the point to judge the simulation to be converged.
abdul099 is offline   Reply With Quote

Old   May 1, 2012, 10:52
Default
  #8
New Member
 
James Berlin
Join Date: Apr 2012
Posts: 3
Rep Power: 4
jkbg38@mail.missouri.edu is on a distinguished road
Should I keep going or call this good?
Attached Images
File Type: jpg cfdforum.jpg (83.0 KB, 38 views)
File Type: jpg cfdforum1.JPG (72.2 KB, 37 views)
File Type: jpg cfdforum2.JPG (56.6 KB, 34 views)
jkbg38@mail.missouri.edu is offline   Reply With Quote

Old   May 1, 2012, 13:17
Default
  #9
Senior Member
 
Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 442
Rep Power: 10
LuckyTran is on a distinguished road
Quote:
Originally Posted by jkbg38@mail.missouri.edu View Post
Should I keep going or call this good?
Look at your actual pressure and kinetic energy plots. They are not converged by any means, there are large and rampant changes in their values.
Verdict: Not converged
LuckyTran 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
turbulence modeling in turbofan Saima Main CFD Forum 5 November 21, 2010 17:47
Turbulence Modeling meshkati FLUENT 0 January 23, 2010 16:04
Turbulence Modeling in CD-ROM? huygen julyend FLUENT 0 December 1, 2008 07:11
turbulence modeling error at a stagnation point erdem CD-adapco 4 August 9, 2006 10:44
CFD Modeling of Two-phase Flow in Small Dia.Tubes Eric Poindexter Main CFD Forum 2 September 22, 2000 09:21


All times are GMT -4. The time now is 04:50.