CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

rhoSimpleFoam for compressible fluid through an orifice

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

LinkBack Thread Tools Display Modes
Old   February 2, 2010, 09:45
Default rhoSimpleFoam for compressible fluid through an orifice
New Member
Willem van der Meer
Join Date: Jan 2010
Posts: 3
Rep Power: 9
WillemvdM is on a distinguished road

Currently I am trying to solve a problem consisting of a compressible gas flow through an orifice. I am using rhoSimpleFoam for this simulation and has been struggling with this solver for weeks now.

My problem definitely does not lie with the relaxation factors as discussed in quite a few threads in this forum because adjusting this only varies when the solution blows up.

I have also tested the simulation with different solvers like rhoPisoFoam etc with no success.

After a thorough scan a my most recent blow-up a realised that it is temperature related. As soon as it reaches 0C it blows up because the thermophysical values can not be computed. The weird thing for me is that the temperature is not suppose to vary to much, but in my simulation is varies consideribly where the velocities are the highest. For all my simulations I have used the following thermoType:

hPsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>

I used this thermophysical model because it was the only one that seemed to work without errors. When I try the ePsiThermo model as follows it gives the error as displayed below:

ePsiThermo<pureMixture<constTransport<specieThermo <eConstThermo<perfectGas>>>>>;

Not implemented#0 Foam::error:rintStack(Foam::Ostream&) in "/home/bernarde/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/"
#1 Foam::error::abort() in "/home/bernarde/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/"
#2 Foam::basicThermo::h() in "/home/bernarde/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/"
#3 main in "/home/bernarde/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/rhoSimpleFoam"
#4 __libc_start_main in "/lib64/"
#5 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116

From function basicThermo::h()
in file basicThermo/basicThermo.C at line 260.

At this point I am really stuck and would appreciate any advice and help with this problem

OpenFoam 1.6
OpenSuse 11.1
WillemvdM is offline   Reply With Quote

Old   February 3, 2010, 07:44
New Member
Willem van der Meer
Join Date: Jan 2010
Posts: 3
Rep Power: 9
WillemvdM is on a distinguished road
It seems like changing my inlet conditions to 'n pressureInletVelocity helped with stability. It still have to nurture my relaxation factors.
WillemvdM is offline   Reply With Quote

Old   June 20, 2012, 03:54
Senior Member
Bernhard Linseisen
Join Date: May 2010
Location: Magdeburg/Geneva
Posts: 182
Blog Entries: 1
Rep Power: 9
Linse is on a distinguished road
I wonder what the results were in this case. Were they true to analytical solutions such as calculated from ISO 5167-2?

I am trying to do the same. While the simulations don't blow up, they definitely give results strongly differing from the ISO-solution.

For comparison reasons I use a case which at no point reaches Ma > 0.25, so incompressible treatment should be okay as well. Indeed simpleFoam does provide me with a result that does fit the ISO-solution quite well (12 % of difference for the pressure difference, qualitative very similar flow field).
Using the same case with rhoSimpleFoam nevertheless gives results which are quite off for the pressure difference (ca. 20 to 25%) and by no way does agree qualitatively with the flow field.

The mesh used was identical for both cases, simpleFoam and rhoSimpleFoam.

Attached you can see the velocity magnitude fields for both simulations.
Among other things the flow profile should become normal again after some way. It does that for simpleFoam, while rhoSimpleFoam keeps this strange kind of layering...

Any comments/suggestions are welcome!


Linse is offline   Reply With Quote


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
Flow past orifice aurel Main CFD Forum 0 January 27, 2009 07:46
fluid structure interaction taru agrawal FLUENT 4 September 10, 2007 03:12
How to apply negtive pressure to outlet bioman66 CFX 5 June 3, 2006 01:40
My Revised "Time Vs Energy" Article For Review Abhi Main CFD Forum 2 July 9, 2002 09:08
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 09:11

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