CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Pressure stair-step behaviour using rhopimplefoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 11, 2014, 08:14
Default Pressure stair-step behaviour using rhopimplefoam
  #1
Senior Member
 
joegi
Join Date: Nov 2009
Location: genoa
Posts: 101
Rep Power: 16
joegi.geo is on a distinguished road
Hi,

I wonder if somebody else is experiencing this, I am having an strange stair-step behavior on the pressure when using rhopimplefoam and I don’t manage to pinpoint the reason of this.

The case consists of a simple straight pipe (Poiseuille case, so we have an analytical solution). I am computing the solution using icoFoam, pimpleFoam, and rhoPimpleFoam for the same mesh, physical properties, and same boundary conditions. The peak velocity is about 2.8 m/s, so I expect rhoPimpleFoam to yield to the same results as those of the incompressible solvers. While the velocity is pretty much the same, the pressure along the pipe axis of rhoPimpleFoam shows a stair-step (wiggles) behavior, whereas the pressure reported by the incompressible solvers has a smooth behavior (as expected and desired). At this point, I wonder if you have seen this behavior.

I have used many different numerical schemes and I always see the same behavior. Using finer meshes did not yield to better results.

By the way, I see the same behavior using the Heat transfer and buoyancy-driven flows solvers.


To reproduce the results:

Download the case file and go to rp5:

blockMesh -dict system/blockMeshDict

rhoPimpleFoam

sample -latestTime


Then go to pimple:

blockMesh

pimpleFoam

sample -latestTime


Then compare the solutions for both cases (pressure).


In figure 1 it is plot a comparison between icofoam and rhopimplefoam. In figure 4 it is plot a compariosn between fluent a openfoam. The numerical setup, mesh, bc and ic are identical in fluent.
Attached Images
File Type: jpg f1.jpg (19.8 KB, 14 views)
File Type: png f4.png (15.9 KB, 15 views)
Attached Files
File Type: gz pimple.tar.gz (6.4 KB, 7 views)
joegi.geo is offline   Reply With Quote

Old   December 12, 2014, 06:01
Default
  #2
Senior Member
 
Tom Fahner
Join Date: Mar 2009
Location: Breda, Netherlands
Posts: 634
Rep Power: 32
tomf will become famous soon enoughtomf will become famous soon enough
Send a message via MSN to tomf Send a message via Skype™ to tomf
Hi this is probably related to the setting for the writePrecision in your controlDict. If you change this to 12 instead of 6 you will loose the wiggles I assume.

Or you could use binary instead of ascii.

Regards,
Tom
tomf is offline   Reply With Quote

Old   December 12, 2014, 12:38
Default
  #3
Senior Member
 
joegi
Join Date: Nov 2009
Location: genoa
Posts: 101
Rep Power: 16
joegi.geo is on a distinguished road
Hi Tom,

Thank you for point that out.

This is the second time that this bites me with totally different solvers.

FYI, SnappyHexMesh is also sensitive to this issue.

Now I wonder why do they give the double precision version, and do not save the solution using all the significant digits (as default option).

Have a good one,

jg
joegi.geo is offline   Reply With Quote

Old   December 12, 2014, 13:10
Default
  #4
Senior Member
 
Tom Fahner
Join Date: Mar 2009
Location: Breda, Netherlands
Posts: 634
Rep Power: 32
tomf will become famous soon enoughtomf will become famous soon enough
Send a message via MSN to tomf Send a message via Skype™ to tomf
Hi

If you are referring to the tutorials as "default option" I would say that I never consider the tutorials as anything else than example cases showing the general working of a solver and the information that is needed. You should check any option they use and see if it fits your purpose and you would get a "best practice" for that solver, after all many tutorials use upwind for many terms in fvSchemes, which should not be used in production runs!

Kind regards,
Tom
tomf is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Setting up the pressure variation due to tornado in a duct(UDF)+animation guillaume1990 FLUENT 0 March 3, 2014 12:59
same geometry,structured and unstructured mesh,different behaviour. sharonyue OpenFOAM Running, Solving & CFD 13 January 2, 2013 23:40
Does star cd takes reference pressure? monica Siemens 1 April 19, 2007 12:26
what the result is negatif pressure at inlet chong chee nan FLUENT 0 December 29, 2001 06:13
Hydrostatic pressure in 2-phase flow modeling (long) DS & HB Main CFD Forum 0 January 8, 2000 16:00


All times are GMT -4. The time now is 01:23.