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

Compressible flow simulation at mach .7

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 15, 2009, 14:46
Default Compressible flow simulation at mach .7
  #1
New Member
 
Matthew
Join Date: Apr 2009
Posts: 12
Rep Power: 8
mbeals is on a distinguished road
This has probably been asked before, but a search didn't seem to help me all that much.

I have a model that I'm attempting to run in compressible mode, but it's diverging horribly.

My initial conditions are based off of a incompressible run which is converged to under 1x10^-4 continuity with the second-order solver.

This solution was run with these conditions:

operating pressure: 218 hPa
air density: 0.3507 Kg/m^3
velocity inlet @ 225 m/s on front face
outflow on rear face
symmetry for surrounding faces
wall boundary for the body of the object

For the compressible runs, I set all of the domain boundaries to pressure-far-field with:
gauge pressure: 0
temp: 216 K
mach #: 0.7
momentum along +x vector

I then patch the domain, setting domain temp to 216K

I'm running the simulation pressure based, cell based, and inviscid

I've tried density based as well and it does the same thing.

It begins warning that it's limiting pressure in certain nodes to 1, then it detects divergence and starts dropping the Courant number and timestep for some cells. Finally the temperature in some cells starts to spike and has to be limited. The whole time residuals are spiking. If I let it continue, it finally dies with a NaN error in the thermal solution.

Any ideas of what could be the culprit or what to change/modify to make it happy?

Last edited by mbeals; April 15, 2009 at 15:22.
mbeals is offline   Reply With Quote

Old   April 16, 2009, 02:17
Default
  #2
Member
 
Join Date: Apr 2009
Posts: 38
Rep Power: 8
doki is on a distinguished road
there might be several reasons for the problem.
firstly you probably have to use density based/implicit solver.
second, the outer boundary should be set to pressure far field
check the operating conditions and set the operating pressure to ambient pressure.
you do not need to patch the temperature for the domain.
the mesh for the case should be somehow a fine mesh.
the total pressure is the sum of static and DYNAMIC pressure. so consider it!
doki is offline   Reply With Quote

Old   April 16, 2009, 12:26
Default
  #3
New Member
 
Matthew
Join Date: Apr 2009
Posts: 12
Rep Power: 8
mbeals is on a distinguished road
All of my outer boundaries are set to pressure far field and a gauge pressure of 0 hPa. I've also set the thermal component to the ambient air temp (216 K) and the momentum vector to point in the +x direction.

The operating conditions have the ambient pressure set to the proper value.

I did attempt another run using K-omega turbulence for viscosity and the density based solver and it was MUCH more stable. There was a rather large initial spike in residuals, especially the energy, but it ran steadily for 1000 iterations last night.

Now the problem is that it ran for 15 hours and residuals dropped less then an order of magnitude.

Any thoughts on how to get it to converge faster, or should I just let is churn all weekend?
mbeals is offline   Reply With Quote

Old   April 17, 2009, 07:44
Default
  #4
Member
 
Join Date: Apr 2009
Posts: 38
Rep Power: 8
doki is on a distinguished road
that is sth. related to the mesh size and also the system RAM and CPU.
if your system has a kind of multi-core CPU, you can use the parallel computing ability of the software and get smaller times per iteration.
also, generally you can change the courant No. (control/solution panel) and set it to greater values. of course that depends to the case and greater than 5 values for the number may cause the case to diverge!
doki is offline   Reply With Quote

Old   April 17, 2009, 10:29
Default
  #5
New Member
 
Matthew
Join Date: Apr 2009
Posts: 12
Rep Power: 8
mbeals is on a distinguished road
Right now I'm running the solution on 4 cores (of a 16 core cluster) with 62 gig of ram. I don't have enough licenses to use more cores. The mesh I'm working with is 6.5 million elements.

Besides courant number, is there anything else I can tweak? Any of the under relaxation factors?

edit: I may have found something. I used a incompressible solution converged with a second order flow solver and I didn't change it back to first order before I started iterating. I switched it over to 1st order flow and it seems to be converging fast. I'll have to let it run a few hours to be sure though.

Last edited by mbeals; April 17, 2009 at 12:29.
mbeals is offline   Reply With Quote

Old   April 18, 2009, 05:35
Default
  #6
Member
 
Pratik Mehta
Join Date: Mar 2009
Posts: 32
Rep Power: 8
pratikmehta is on a distinguished road
hi,

I am sorry if i have missed something but from the discussion above it is slightly unclear what is your geometry . You said you intilisalised the compressible run with the inputs of the incompressible flow ( which is ok )
here is something i really did not get it

"My initial conditions are based off of a incompressible run which is converged to under 1x10^-4 continuity with the second-order solver.

This solution was run with these conditions:

operating pressure: 218 hPa
air density: 0.3507 Kg/m^3
velocity inlet @ 225 m/s on front face ??? ( did u input this for incompresible case set up ??, if yes then why
outflow on rear face
symmetry for surrounding faces
wall boundary for the body of the object
"

What turbulence model did u use ??


Thanks for your reply

ciao
pratik
pratikmehta is offline   Reply With Quote

Old   April 18, 2009, 10:49
Default
  #7
New Member
 
Matthew
Join Date: Apr 2009
Posts: 12
Rep Power: 8
mbeals is on a distinguished road
When I switched to compressible, I changed all of my boundaries to pressure-far-field. I input the ambient air temp, kept the gauge pressure at 0, set the mach number to 0.7 and made sure the momentum vector pointed in the +x direction. So the initial solution for the compressible run consisted of the incompressible solution for the inner cells of the domain and new boundary conditions.

I am using the K-omega turbulence model with the default values.
mbeals 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
Compressible Fluid Flow in COMSOL Multiphysics BBG Main CFD Forum 1 November 19, 2008 15:05
compressible flow computation amv Main CFD Forum 5 June 27, 2003 07:27
Compressible Flow Modelling? yeo FLUENT 4 March 7, 2003 08:08
simulation of compressible flow svlesh Phoenics 0 March 7, 2002 06:53
compressible channel flow.. R.D.Prabhu Main CFD Forum 0 July 17, 1998 17:23


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