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

interFoam pressure miscalculation in 2.3 (wrt previous versions)

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 21, 2014, 08:37
Default interFoam pressure miscalculation in 2.3 (wrt previous versions)
  #1
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Singapore
Posts: 339
Rep Power: 8
Phicau is on a distinguished road
As I anticipated in Zagreb, I have been running a benchmark case to test each new version of OpenFOAM. It is just a simple dam break case in 3D for which pressure and free surface laboratory data are available:

http://dissertations.ub.rug.nl/facul...m.t.kleefsman/

Up to OpenFOAM version 2.2 the numerical results match the experimental ones with a high degree of accordance. However, in version 2.3 everything is off, and especially the pressure, having differences of more than 2000Pa (200% difference).

The bug has been reported: http://www.openfoam.org/mantisbt/view.php?id=1354

See the attached file. It includes a report and the cases to run in OpenFOAM 2.1.1 and 2.3.
Phicau is offline   Reply With Quote

Old   February 16, 2015, 10:52
Default any news on this bug?
  #2
New Member
 
Oystein Lande
Join Date: May 2014
Posts: 2
Rep Power: 0
oystelan is on a distinguished road
Hi Pablo,
Thank you for the release of IHFoam.
I read your bug report and the corresponding replies on the openfoam website.
The subject was marker as closed some time ago. Did you find any solution to the pressure problems in v.2.3 other than the explanation given in the bug report reply? (perhaps a way around the problem?)I assume that the same problem is present for v.2.3.1 (?)
oystelan is offline   Reply With Quote

Old   February 16, 2015, 11:34
Default
  #3
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Singapore
Posts: 339
Rep Power: 8
Phicau is on a distinguished road
Hi Řystein

welcome to IHFOAM!

Unfortunately the OpenFOAM response was childish and they did not investigate further because the techniques they rewrote for 2.3 work well for the single-phase solvers, or so they say. Therefore, the bug is still there, and I do not expect to see any further improvements soon.

It is a pity, but I suggest that you use version 2.2.2 or previous ones.

Best,

Pablo
Phicau is offline   Reply With Quote

Old   February 16, 2015, 11:54
Default
  #4
New Member
 
Oystein Lande
Join Date: May 2014
Posts: 2
Rep Power: 0
oystelan is on a distinguished road
Thanks, and thank you for the quick reply.
Unfortunately your answer is as i expected. I was hoping to get away with 2.3.0 (since this is was is currently installed on the hpc), but i will do as you suggest and install v 2.2.2.
Its a pitty....
oystelan is offline   Reply With Quote

Old   February 20, 2015, 08:52
Default
  #5
Member
 
Join Date: Dec 2009
Posts: 49
Rep Power: 8
katakgoreng is on a distinguished road
Hi Pablo,

Currently I'm using interFoam (OpenFoam 2.3) to calculate forces using funtionObject but I got a weird result. I do think this is due to pressure calculation error as you reported.
In the report, you said that you compiles bouyantPressure from OF 2.2.2 into OF 2.2.3. Would you mind elaborate a bit on the steps?

Kind regards,
katakgoreng
katakgoreng is offline   Reply With Quote

Old   February 20, 2015, 09:11
Default
  #6
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Singapore
Posts: 339
Rep Power: 8
Phicau is on a distinguished road
Hi katakgoreng,

yes, I copied buoyant pressure from version 2.2.2 and compiled it over in 2.3.0. I guess I needed to change something, but I don't remember it exactly. In any case It must have been minor stuff. I am not currently in my computer and will not be until mid next week, so I cannot provide the code now. If you are interested let me know, but note the following paragraph.

The important conclusion is that buoyantPressure yielded the same results as fixedFluxPressure, so the bug is definitely not in the boundary condition, but inside the solver. Probably in the new formulation of ddtPhiCorr and fvc::reconstruct.

Best,

Pablo
Phicau is offline   Reply With Quote

Old   November 25, 2015, 10:42
Default
  #7
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Singapore
Posts: 339
Rep Power: 8
Phicau is on a distinguished road
Dear all,

I have updates on this topic and the conclusion is that I am to blame for the mistake.

The setup of the cases I was comparing was different (initial value of turbulence). Once that was corrected I could check that both versions yield very similar results.

Sorry for any inconveniences caused.

Best regards,

Pablo
Phicau is offline   Reply With Quote

Reply

Tags
bugs, interfoam

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
how to set up Hydrostatic pressure distribution in interFoam? wes1204 OpenFOAM Pre-Processing 15 July 29, 2016 05:29
Pressure Outlet Guage pressure Mohsin FLUENT 36 April 29, 2016 17:16
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 06:27
InterPhaseChangeFoam ERROR shipman OpenFOAM Running, Solving & CFD 37 March 23, 2014 13:43
interFoam -- pressure loss FerdiFuchs OpenFOAM Running, Solving & CFD 1 February 12, 2014 10:16


All times are GMT -4. The time now is 14:35.