CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Post-Processing

Reported phi values exactly twice as high as expected

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

LinkBack Thread Tools Search this Thread Display Modes
Old   November 28, 2016, 14:59
Question Reported phi values exactly twice as high as expected
New Member
Jay B.
Join Date: Nov 2016
Posts: 12
Rep Power: 6
New_Old is on a distinguished road
Hello everyone,

I am new to OpenFoam and today I was getting familiar with the run-time monitoring possibilities when I spotted some strange results.

I have spent quite some time trying to figure this out with no luck, but I believe I have now narrowed the problem down to the point where a more experienced user could easily tell me what's going on.

I attach here two small cases, ready to go on v1606+ (./Allrun).Those are somewhat modified versions of the multiRegionHeater tutorial case:
- case 1 is a transient run, very similar to the original tutorial provided with my version of OpenFOAM (v1606+ on Windows).
- case 2 is a steady-state analysis of the the same problem.

For both case I monitor what is going on at minX on the water side (inlet, fixed velocity), and obtain the following results:

Case 1:
- Total area:	4.00000000e-03
- average(U):	(1.00000000e-03 0.00000000e+00 0.00000000e+00)
- sum(phi):	-4.00000000e-03
Case 2:
- Total area:	4.0000000e-03
- average(U):	(1.0000000e-03 0.0000000e+00 0.0000000e+00)
- sum(phi):	-8.0000000e-06
The reported areas and velocities are correct. So is Phi for case one, assuming OpenFOAM reports a mass flow rate in this case.
However, I can't understand how Phi is computed in case two. A volumetric flow rate would explain the different order of magnitude but the reported flow rate would still be twice as high as I would expect.

I have tried to run case 2 on a single processor, which didn't change the results.

I would be extremely thankful if someone could have a look at this.

Attached Files
File Type: zip (150.5 KB, 1 views)
File Type: zip (150.9 KB, 2 views)
New_Old is offline   Reply With Quote

Old   December 2, 2016, 06:44
Lightbulb Update
New Member
Jay B.
Join Date: Nov 2016
Posts: 12
Rep Power: 6
New_Old is on a distinguished road
I believe I have identified what the problem was:
As a consequence of some shameless copy-pasting of the fvSolution files from a steady-state tutorial, and although I had defined the fluid density as constant (1000 kg/m3), those SIMPLE settings were limiting the density to 2 kg/m3.

Lesson learned...

    momentumPredictor on;
    nNonOrthogonalCorrectors 0;
    pRefCell        0;
    pRefValue       100000;
    rhoMin          0.2;
    rhoMax          2;
New_Old is offline   Reply With Quote


flow rate, flow rate weighting, phi

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
problem during mpi in server: expected Scalar, found on line 0 the word 'nan' muth OpenFOAM Running, Solving & CFD 3 August 27, 2018 04:18
The udf.h headers are unable to open- in VISUAL STUDIO 13 sanjeetlimbu Fluent UDF and Scheme Programming 4 May 2, 2016 05:38
How to install CGNS under windows xp? lzgwhy Main CFD Forum 1 January 11, 2011 18:44
RMS values too high! Usman CFX 1 January 22, 2008 16:53
Multicomponent fluid Andrea CFX 2 October 11, 2004 05:12

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