CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Bugs (https://www.cfd-online.com/Forums/openfoam-bugs/)
-   -   interFoam pressure miscalculation in 2.3 (wrt previous versions) (https://www.cfd-online.com/Forums/openfoam-bugs/139286-interfoam-pressure-miscalculation-2-3-wrt-previous-versions.html)

Phicau July 21, 2014 08:37

interFoam pressure miscalculation in 2.3 (wrt previous versions)
 
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.

oystelan February 16, 2015 09:52

any news on this bug?
 
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 (?)

Phicau February 16, 2015 10:34

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

oystelan February 16, 2015 10:54

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....

katakgoreng February 20, 2015 07:52

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

Phicau February 20, 2015 08:11

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 November 25, 2015 09:42

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


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