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

Supersonic flow with rhoCentralFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 29, 2013, 05:24
Default Supersonic flow with rhoCentralFoam
  #1
New Member
 
Robin Debroux
Join Date: Oct 2012
Posts: 22
Rep Power: 4
robdeb is on a distinguished road
Hi everybody,

I'm trying to simulate a supersonic flow in a convergent-divergent nozzle with rhoCentralFoam.

I have a problem with my inlet boundary conditions:
I impose a totale pressure and use pressureInletVelocity for the velocity.

The problem is that I would like to impose (0 0 0) as value of the pressureInletVelocity condition but when I do this, I get a "floating point exception". When I impose a non-zero value, it works but it doesn't give me the expected results.

When I use sonicFoam with the same conditions, everything goes well. The aim of using rhoCentralFoam is to compare the two solvers, so I need some help.

Thank you.
robdeb is offline   Reply With Quote

Old   January 29, 2013, 12:28
Default
  #2
Member
 
Valentin Wibaut
Join Date: Oct 2012
Posts: 45
Rep Power: 4
vwibaut is on a distinguished road
Hello

I have the same problems with rhoCentralFoam and I also try to simulate the flow in a nozzle. But when I impose a different value for the velocity I have the same error (after more iterations).
Do you find the problem robdeb?

Regards
vwibaut is offline   Reply With Quote

Old   January 31, 2013, 05:20
Default
  #3
New Member
 
Robin Debroux
Join Date: Oct 2012
Posts: 22
Rep Power: 4
robdeb is on a distinguished road
Hi Valentin,

Yes, I get the same problem. I think the problem come from the fact I set a to big inlet pressure compare to my outlet pressure (inlet: 1100000, outlet : 101325) and the solver can't handle it. I try a setField to divide my domain and impose a pressure that decrease slowly to the outlet value but it didn't solve anything.

If anybody has a solution, it would be great.

Thank you.
robdeb is offline   Reply With Quote

Old   January 31, 2013, 12:09
Default
  #4
Senior Member
 
Mahdi Hosseinali
Join Date: Apr 2009
Posts: 124
Rep Power: 8
anishtain4 is on a distinguished road
Quote:
Originally Posted by robdeb View Post
Hi everybody,

I'm trying to simulate a supersonic flow in a convergent-divergent nozzle with rhoCentralFoam.

I have a problem with my inlet boundary conditions:
I impose a totale pressure and use pressureInletVelocity for the velocity.

The problem is that I would like to impose (0 0 0) as value of the pressureInletVelocity condition but when I do this, I get a "floating point exception". When I impose a non-zero value, it works but it doesn't give me the expected results.

When I use sonicFoam with the same conditions, everything goes well. The aim of using rhoCentralFoam is to compare the two solvers, so I need some help.

Thank you.
I believe rhoCentralFoam uses absolute pressure in the thermophysical model, so when you set it to zero it is just physically meaningless to it.

and by the way, in my personal experience rhoCenralFoam is a little bit fragile in comparison with say rhoPimpleFoam
anishtain4 is offline   Reply With Quote

Old   January 31, 2013, 12:23
Default
  #5
New Member
 
Robin Debroux
Join Date: Oct 2012
Posts: 22
Rep Power: 4
robdeb is on a distinguished road
Quote:
Originally Posted by anishtain4 View Post
I believe rhoCentralFoam uses absolute pressure in the thermophysical model, so when you set it to zero it is just physically meaningless to it.

and by the way, in my personal experience rhoCenralFoam is a little bit fragile in comparison with say rhoPimpleFoam
Thank you for your reply.
But when I set the velocity to a non-zero value, it makes more iterations, but I get the same error as vwibaut said it.
robdeb is offline   Reply With Quote

Old   January 31, 2013, 15:31
Default
  #6
Senior Member
 
Mahdi Hosseinali
Join Date: Apr 2009
Posts: 124
Rep Power: 8
anishtain4 is on a distinguished road
Check the courant number, it's a matter of wave propagation, you may play with some stuff to hinder it but if you set the case in a manner that leads to it, it will be inevitable.
anishtain4 is offline   Reply With Quote

Old   January 31, 2013, 19:17
Default
  #7
Member
 
Valentin Wibaut
Join Date: Oct 2012
Posts: 45
Rep Power: 4
vwibaut is on a distinguished road
Quote:
Originally Posted by anishtain4 View Post
Check the courant number, it's a matter of wave propagation, you may play with some stuff to hinder it but if you set the case in a manner that leads to it, it will be inevitable.
The courant number seems to be good. Moreover we give an upper limit in the controldict file so it is not our problem.
vwibaut is offline   Reply With Quote

Old   February 1, 2013, 05:23
Default
  #8
Senior Member
 
Mahdi Hosseinali
Join Date: Apr 2009
Posts: 124
Rep Power: 8
anishtain4 is on a distinguished road
this error means a variable is being out of range (divide by zero, or negative in sqrt, etc) and this usually happens when the physical modelling collapses. To see what is wrong set very small time steps so you may get many runs and have some written data before the core gets dump. Then view your flow and report what is getting out of order first?
anishtain4 is offline   Reply With Quote

Old   February 1, 2013, 15:47
Default
  #9
Member
 
Join Date: Jun 2012
Posts: 55
Rep Power: 5
maHein is on a distinguished road
What did you use as initial conditions for U? I would recommend to use values which are close to the final expected axial speed. Furthermore, you could try to steadily increase the total pressure at the inlet in order to avoid extreme pressure and velocity gradients within the first iterations.
maHein is offline   Reply With Quote

Old   February 2, 2013, 06:33
Default
  #10
Member
 
Valentin Wibaut
Join Date: Oct 2012
Posts: 45
Rep Power: 4
vwibaut is on a distinguished road
Quote:
Originally Posted by maHein View Post
What did you use as initial conditions for U? I would recommend to use values which are close to the final expected axial speed. Furthermore, you could try to steadily increase the total pressure at the inlet in order to avoid extreme pressure and velocity gradients within the first iterations.
Hi mahein,
First, thanks for your reply.
For U, I have
inlet
{
type pressureInletVelocity;
value uniform (0 0 0);
}
outlet
{
type zeroGradient;
}

I have used values wich are given by sonicFoam just before convergence and with that no more "floating point exception" error but there is another problem. The solver rhoCentralFoam doesn't give a (good) solution and never stabilize although the initial conditions are close to convergence. When I look at the solution at the different time step, I can see that the cells at the inlet have high pressure. And this pressure increases with time. I don't know if this is the cause of the no convergence.

error.jpg

Other test I made,
I use setFields to impose different pressure in the nozzle and decrease the pressure and velocity gradient for the first iterations. There is no more floating point exception but the solution is wrong!! I thtink I'll try to use a time varying condition.

Thank you very much for your help for the floating point exception
vwibaut is offline   Reply With Quote

Old   February 3, 2013, 08:32
Default
  #11
Member
 
Join Date: Jun 2012
Posts: 55
Rep Power: 5
maHein is on a distinguished road
I also experienced unphysical pressure distributions at the inlet when using "pressureInletVelocity" for U at the inlet. Try "pressureInletOutletVelocity" instead. It improved my results.
maHein is offline   Reply With Quote

Old   February 4, 2013, 09:58
Default
  #12
Member
 
Valentin Wibaut
Join Date: Oct 2012
Posts: 45
Rep Power: 4
vwibaut is on a distinguished road
Quote:
Originally Posted by maHein View Post
I also experienced unphysical pressure distributions at the inlet when using "pressureInletVelocity" for U at the inlet. Try "pressureInletOutletVelocity" instead. It improved my results.
I have tested "pressureInletOutletVelocity" :
inlet
{
type pressureInletOutletVelocity;
value uniform (0 0 0);
}

But this doesn't solve the problem
vwibaut is offline   Reply With Quote

Old   February 6, 2013, 05:30
Default
  #13
Member
 
Valentin Wibaut
Join Date: Oct 2012
Posts: 45
Rep Power: 4
vwibaut is on a distinguished road
When I replace the boundary condition for the velocity (pressureInletVelocity) by zeroGradient, I have the same problem. An idea??

For the other problem, I try with setFields but always the problem of wrong solution. Somebody has an idea about how to solve this problem?
vwibaut is offline   Reply With Quote

Old   February 11, 2013, 06:43
Default
  #14
Member
 
Valentin Wibaut
Join Date: Oct 2012
Posts: 45
Rep Power: 4
vwibaut is on a distinguished road
I see on other threads that there is problem with rhoCentralFoam when using totalPressure and/or totalTemperature. But why? I don't understand
vwibaut is offline   Reply With Quote

Old   February 18, 2013, 04:09
Default
  #15
New Member
 
Sami
Join Date: Sep 2012
Posts: 15
Rep Power: 4
abdul Sami is on a distinguished road
Plz tell me the terms, static pressure, total pressure, static temperature and total temperature with respect to convergent divergent nozzle I m so worried plzz help me I will be thankful
abdul Sami is offline   Reply With Quote

Old   February 19, 2013, 04:12
Default
  #16
Senior Member
 
Mahdi Hosseinali
Join Date: Apr 2009
Posts: 124
Rep Power: 8
anishtain4 is on a distinguished road
Quote:
Originally Posted by abdul Sami View Post
Plz tell me the terms, static pressure, total pressure, static temperature and total temperature with respect to convergent divergent nozzle I m so worried plzz help me I will be thankful
read a book like the one by anderson in gas dynamic
anishtain4 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
back pressure at exit for supersonic flow BigFrank FLUENT 6 April 24, 2014 15:13
Supersonic flow over an airfoil. FLUENT Lift don't match up with Shock Expansion Tech enriccasas FLUENT 1 August 21, 2013 06:43
Flow in Supersonic MHD Channel Antanas CFX 6 February 13, 2012 15:49
Modeling Bleed in Supersonic Flow Mohd Yousuf Ali FLUENT 0 March 17, 2007 04:51
Inviscid Drag at subsonic, subcritical Mach # Axel Rohde Main CFD Forum 1 November 19, 2001 13:19


All times are GMT -4. The time now is 19:39.