interFoam pressure miscalculation in 2.3 (wrt previous versions)

Phicau · July 21, 2014, 08:37

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

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

July 21, 2014, 08:37	interFoam pressure miscalculation in 2.3 (wrt previous versions)	#1
Phicau Senior Member Pablo Higuera Join Date: Jan 2011 Location: Auckland Posts: 627 Rep Power: 19	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.

February 16, 2015, 09:52	any news on this bug?	#2
oystelan New Member Oystein Lande Join Date: May 2014 Posts: 2 Rep Power: 0	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 (?)

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
how to set up Hydrostatic pressure distribution in interFoam?	wes1204	OpenFOAM Pre-Processing	24	January 19, 2024 10:07
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 12:43
interFoam -- pressure loss	FerdiFuchs	OpenFOAM Running, Solving & CFD	1	February 12, 2014 09:16

February 16, 2015, 10:34		#3
Phicau Senior Member Pablo Higuera Join Date: Jan 2011 Location: Auckland Posts: 627 Rep Power: 19	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

February 16, 2015, 10:54		#4
oystelan New Member Oystein Lande Join Date: May 2014 Posts: 2 Rep Power: 0	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....

February 20, 2015, 07:52		#5
katakgoreng Member Join Date: Dec 2009 Posts: 49 Rep Power: 16	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

February 20, 2015, 08:11		#6
Phicau Senior Member Pablo Higuera Join Date: Jan 2011 Location: Auckland Posts: 627 Rep Power: 19	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

November 25, 2015, 09:42		#7
Phicau Senior Member Pablo Higuera Join Date: Jan 2011 Location: Auckland Posts: 627 Rep Power: 19	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