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

pressureTools postProcessing function objects in 2.2.0

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

Like Tree6Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   May 23, 2013, 17:16
Default pressureTools postProcessing function objects in 2.2.0
  #1
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18
immortality is on a distinguished road
here:
http://www.openfoam.org/version2.2.0...processing.php
its mentioned about pressureTools.how can use these function objects?
is there a tutorial or more guidance through?
new pressureTools - enables the calculation of pressure (from kinematic pressure), total pressure, pressure coefficient and total pressure coefficient.
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   May 23, 2013, 17:58
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,659
Blog Entries: 39
Rep Power: 99
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Hi Ehsan,

Quick answer:
  1. In the following page is a description of the base configurations for function objects: http://foam.sourceforge.net/docs/cpp/a00002.html
  2. In this page is the description on how to use the pressure tools: http://foam.sourceforge.net/docs/cpp...8.html#details
Best regards,
Bruno
Mojtaba.a, Julien and JR22 like this.
__________________
wyldckat is offline   Reply With Quote

Old   May 24, 2013, 06:27
Default
  #3
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18
immortality is on a distinguished road
thank you Bruno
can use a
Code:
operation       areaAverage;
with it?
totally can use other features of function objects like the one below with it?
and is there a totalTemperature too or isn't?
Code:
Average_left
    {
        type           faceSource;
        functionObjectLibs ("libfieldFunctionObjects.so");
        enabled on;
        log             yes;
        valueOutput     false;
        surfaceFormat   vtk;
        source          patch;
        sourceName      left;
        operation       areaAverage;
        fields
        (
            p
            //phi // surface fields not supported for sampledSurface
            U
            T
            rho
            h
            e
        );
        outputControl timeStep;//outputTime;//
        outputInterval 1; // 5000;
      }
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   May 25, 2013, 16:33
Default
  #4
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,659
Blog Entries: 39
Rep Power: 99
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Hi Ehsan,

A bit of a confusion of ideas going on here... OK, so to sort out by items:
  • Total temperature? I can find anything like that. And temperature in OpenFOAM is usually in Kelvin, so there isn't anything more "total temperature" than the unit that stands for the absolute temperature measurement!
  • The function object "pressureTools" must come first in the list of function objects, if you want to use the resulting fields elsewhere.
  • The function object "pressureTools" creates volume fields, depending on the calculation. There isn't a direct description, but from the source code, it can create the following field names, depending on the options:
    Code:
    total(p)
    static(p)
    total(p)_coeff
    static(p)_coeff
    These names can be used in the field list for the average function object.
Best regards,
Bruno
__________________
wyldckat is offline   Reply With Quote

Old   May 25, 2013, 16:43
Default
  #5
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18
immortality is on a distinguished road
you mean I can simply add total(p) to the function as below?without adding anything other?
Code:
Average_left
    {
        type           faceSource;
        functionObjectLibs ("libfieldFunctionObjects.so");
        enabled on;
        log             yes;
        valueOutput     false;
        surfaceFormat   vtk;
        source          patch;
        sourceName      left;
        operation       areaAverage;
        fields
        (
            p
            total(p)
            U
            T
            rho
            h
            e
        );
        outputControl timeStep;//outputTime;//
        outputInterval 1; // 5000;
      }
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   May 25, 2013, 16:51
Default
  #6
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,659
Blog Entries: 39
Rep Power: 99
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Where's the "pressureTools" entry?

I have not tested this, this is just copy-paste-adapt from the pages I posted links for:
Code:
    myFunctionObject    // user-defined name of function object entry
    {
        type            pressureTools;
        libs            ("libutilityFunctionObjects.so");
        //region          defaultRegion;
        enabled         yes;
        //timeStart       0;
        //timeEnd         10;
        outputControl timeStep;//outputTime;//
        outputInterval 1; // 5000;
        calcTotal   yes;
        calcCoeff   no;
    }

    Average_left
    {
        type           faceSource;
        functionObjectLibs ("libfieldFunctionObjects.so");
        enabled on;
        log             yes;
        valueOutput     false;
        surfaceFormat   vtk;
        source          patch;
        sourceName      left;
        operation       areaAverage;
        fields
        (
            p
            total(p)
            U
            T
            rho
            h
            e
        );
        outputControl timeStep;//outputTime;//
        outputInterval 1; // 5000;
      }
__________________

Last edited by wyldckat; May 25, 2013 at 16:52. Reason: changed "outputControl"
wyldckat is offline   Reply With Quote

Old   May 25, 2013, 17:43
Default
  #7
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18
immortality is on a distinguished road
Ok.thanks.is it true now?
Code:
totalP    // user-defined name of function object entry
    {
        type            pressureTools;
        libs            ("libutilityFunctionObjects.so");
        //region          defaultRegion;
        enabled         yes;
        //timeStart       0;
        //timeEnd         10;
        outputControl timeStep;//outputTime;//
        outputInterval 1; // 5000;
        calcTotal   yes;
        calcCoeff   no;
    }

//patch-left------------------------------------------------------------------
    flowInlet
    {
        type swakExpression;
        valueType patch;
        patchName left;
        expression "phi";
        accumulations (
            sum
        );
        verbose true; 
    }

   Average_left
    {
        type           faceSource;
        functionObjectLibs ("libfieldFunctionObjects.so");
        enabled on;
        log             yes;
        valueOutput     false;
        surfaceFormat   vtk;
        source          patch;
        sourceName      left;
        operation       areaAverage;
        fields
        (
            p
            //phi // surface fields not supported for sampledSurface
            U
            T
            rho
            h
            e
        );
        outputControl timeStep;//outputTime;//
        outputInterval 1; // 5000;
      }
      totalPressure_left
    {
        $Average_left;
        operation       weightedAverage;
        weighted rho;
        fields
        (
          total(p)
        );
        outputControl timeStep;//outputTime;//
        outputInterval 1; // 5000;
      }
what was the correct form of this part:
Code:
operation       weightedAverage;
        weighted rho;
because I have to use weightedAverage for totalPressure as the link you took me before in another thread.
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   May 25, 2013, 18:47
Default
  #8
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,659
Blog Entries: 39
Rep Power: 99
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
I think it's all OK. I can't see nothing wrong with it.
__________________
wyldckat is offline   Reply With Quote

Old   May 25, 2013, 19:04
Default
  #9
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18
immortality is on a distinguished road
hi Bruno
Really?i thought maybe those lines for weightedAverarage had another keywords as the case you sent to me before but i can't find now.please have a look into it again tomorrow.
Thanks.
Have a good night.
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   May 25, 2013, 19:12
Default
  #10
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,659
Blog Entries: 39
Rep Power: 99
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Hi Ehsan,

Sorry, I don't understand the question. The function objects code seems to be OK. I don't see any problem with it.
The reference page on how to use "faceSource" is this one: http://foam.sourceforge.net/docs/cpp...5.html#details

Wait... now I see the problem! It should be:
Code:
operation weightedAverage;
weightField rho;
Best regards,
Bruno
immortality likes this.
__________________
wyldckat is offline   Reply With Quote

Old   May 26, 2013, 04:30
Default
  #11
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18
immortality is on a distinguished road
Hi
if I want to use function object on a face,how should change these terms?
Code:
source          patch;
        sourceName      left;
found it:
Code:
source          sampledSurface;
        //sourceName      baffleFaces;
        sampledSurfaceDict
        {
          type cuttingPlane;
          planeType pointAndNormal;
          pointAndNormalDict
          {
             basePoint       (0.0345 0.002 0);
             normalVector    (1 0 0);
          }
          interpolate true;
        }
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   June 11, 2013, 17:40
Default
  #12
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18
immortality is on a distinguished road
in compressible case I should use pRef 0;
correct?
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   June 12, 2013, 10:29
Default
  #13
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18
immortality is on a distinguished road
it writes total(p) values in all time steps in the solution folder but not in the patch I have set:
Code:
# Source : patch left
# Faces  : 60
# Time	sum(magSf)	weightedAverage(total(p))
1.123595506e-08	1.2e-05
2.375195563e-08	1.2e-05
3.694449677e-08	1.2e-05
5.013703791e-08	1.2e-05
6.332957905e-08	1.2e-05
7.652212019e-08	1.2e-05
8.971466133e-08	1.2e-05
1.029072025e-07	1.2e-05
1.160997436e-07	1.2e-05
1.292922847e-07	1.2e-05
1.424848259e-07	1.2e-05
1.55677367e-07	1.2e-05
1.688699082e-07	1.2e-05
1.820624493e-07	1.2e-05
1.952549904e-07	1.2e-05
2.084475316e-07	1.2e-05
2.216400727e-07	1.2e-05
2.348326139e-07	1.2e-05
2.48025155e-07	1.2e-05
2.612176961e-07	1.2e-05
2.744102373e-07	1.2e-05
2.876027784e-07	1.2e-05
3.007953196e-07	1.2e-05
3.110777413e-07	1.2e-05
3.210621219e-07	1.2e-05
3.304918147e-07	1.2e-05
3.395392226e-07	1.2e-05
3.482294959e-07	1.2e-05
3.565855281e-07	1.2e-05
3.64730015e-07	1.2e-05
its the totalP function:
Code:
totalP    // user-defined name of function object entry
    {
        type            pressureTools;
        libs            ("libutilityFunctionObjects.so");
        //region          defaultRegion;
        enabled         yes;
        //timeStart       0;
        //timeEnd         10;
        outputControl timeStep;//outputTime;//
        outputInterval 1; // 5000;
        calcTotal   yes;
        calcCoeff   no;
        pRef 0;
    }
and I want on this patch:
Code:
totalPressure_left
      {
        $Average_left;
        operation       weightedAverage;
        weightField     rho;
        fields
        (
          total(p)
        );
      }
that its Average_left:
Code:
Average_left
    {
        type           faceSource;
        functionObjectLibs ("libfieldFunctionObjects.so");
        enabled on;
        log             yes;
        valueOutput     false;
        surfaceFormat   vtk;
        source          patch;
        sourceName      left;
        operation       areaAverage;
        fields
        (
            p
            //phi // surface fields not supported for sampledSurface
            T
            rho
            h
            U
        );
        outputControl timeStep;//outputTime;//
        outputInterval 1; // 5000;
      }
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   June 15, 2013, 15:54
Default
  #14
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,659
Blog Entries: 39
Rep Power: 99
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Hi Ehsan,

Took me a while to get to here, but here goes:
  1. You told me through an email that if you used swak4Foam, you only needed 1 function instead of 3.
  2. Guess what!? swak4Foam 0.2.4 is already available, at least through SVN and at my GitHub account https://github.com/wyldckat/swak4foam (look for the big ZIP button).
    (Side note: Bernhard was at the OFW8 these past days and I think he hasn't had the time to do the final announcements, but he has already provided these news over Twitter.
    ... edit: never mind, he's on top of things : swak4Foam - New release post #29)
  3. With this version of swak4Foam, which builds without problems with OpenFOAM 2.2.0, the following additional code for the function object should work as expected:
    Code:
    outputControlMode outputTime;
    outputInterval 1;


__________________


As for the original issue on the previous post, I'll quote what I've already sent you through emails - The following function objects worked well together:
Code:
functions
{
    writeMissingFields
    {
        type writeRegisteredObject;
        functionObjectLibs ( "libIOFunctionObjects.so" );
        objectNames ("phi");
        outputControl     outputTime;
    }
//................................................................................
    totalP    // user-defined name of function object entry
    {
        type            pressureTools;
        libs            ("libutilityFunctionObjects.so");
        //region          defaultRegion;
        enabled         yes;
        //timeStart       0;
        //timeEnd         10;
        outputControl timeStep;//outputTime;//
        outputInterval 1; // 5000;
        calcTotal   yes;
        calcCoeff   no;
        pRef 0;
    }


    Average_left
    {
        type           faceSource;
        functionObjectLibs ("libfieldFunctionObjects.so");
        enabled on;
        log             yes;
        valueOutput     false;
        surfaceFormat   vtk;
        source          patch;
        sourceName      left;
        operation       areaAverage;
        fields
        (
            p
            //phi // surface fields not supported for sampledSurface
            T
            rho
            h
            U
        );
        outputControl timeStep;//outputTime;//
        outputInterval 1; // 5000;
      }
      reloadTotalP
      {
          type        readFields;
          functionObjectLibs ("libfieldFunctionObjects.so");
          //region          defaultRegion;
          enabled         yes;
          //timeStart       0;
          //timeEnd         10;
          outputControl timeStep;//outputTime;//
          outputInterval 1; // 5000;
          fields
          (
              "total(p)"
          );
      }
      totalPressure_left
      {
        $Average_left;
        operation       weightedAverage;
        weightField     rho;
        
        fields
        (
          total(p)
        );
      }
}
Notes:
  • The function object "readFields" was needed for reloading back into memory the field "total(p)".
  • As for time folders that only have the field "total(p)": it's because of how the function object was created. It only calculates when it has to write the field to file and it is not registered into memory, which is why we have to use "readFields" for loading it back.
  • For this function object to always calculate the field and save it only occasionally, this would require to edit the original source code for the "pressureTools" function object.
Best regards,
Bruno
immortality likes this.
__________________

Last edited by wyldckat; June 15, 2013 at 16:36. Reason: see "edit:"
wyldckat is offline   Reply With Quote

Old   June 16, 2013, 10:36
Default
  #15
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18
immortality is on a distinguished road
dear Bruno
it wouldn't make any trouble to use it?because I'm in middle of the run and Bernhard has told maybe it will cause some inconsistency with before versions if I have understood correctly?
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   June 16, 2013, 11:09
Default
  #16
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,659
Blog Entries: 39
Rep Power: 99
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Quote:
Originally Posted by immortality View Post
Bernhard has told maybe it will cause some inconsistency with before versions if I have understood correctly?
Are you talking about the following?
Quote:
Originally Posted by gschaider View Post
- functionObjects now store their data in the sub-directory postProcessing according to the new convention OpenFOAM has (sorry: this may break existing scripts relying on the old location)
- slight change in the format written may break the restart of cases where swak-specific restart data was written with older swak-version (should not hit many people)
Mmm... the first detail isn't very relevant... it's just a matter of moving the folders from the previously generated swak-like folders into the "postProcessing" folder, if I'm not mistaken.
The second detail is likely related to groovyBC!?
__________________
wyldckat is offline   Reply With Quote

Old   June 16, 2013, 11:22
Default
  #17
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18
immortality is on a distinguished road
thanks.I'm compiling i now.it will replaced to before version automatically?
and I have to move old folders to postProcessing folder manually and no other issue will take place?
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   June 16, 2013, 11:39
Default
  #18
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,659
Blog Entries: 39
Rep Power: 99
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Quote:
Originally Posted by immortality View Post
thanks.I'm compiling i now.it will replaced to before version automatically?
Yes, I think so.

Quote:
Originally Posted by immortality View Post
and I have to move old folders to postProcessing folder manually and no other issue will take place?
Yes, move the ones in the case folder named "swak*". I haven't tested this yet myself, so I'm not 100% sure of this. But from what I briefly saw yesterday from the source code, I think it's all that is necessary.
immortality likes this.
__________________
wyldckat is offline   Reply With Quote

Old   July 31, 2014, 03:25
Default
  #19
New Member
 
Join Date: Mar 2014
Posts: 17
Rep Power: 4
BenJ is on a distinguished road
Hi Ehsan and Bruno,

I found your post very interesting, since I'm trying to use the pressureTools to compute the static pressure from an incompressible case that uses kinematic pressure.

I would have used the utility as follows:

{
type pressureTools;
libs ("libutilityFunctionObjects.so");
enabled yes;
outputControl timeStep;//outputTime;//
calcTotal no;
calcCoeff no;
}

It does yield a static(p) field, but the values are very low (1e-16).

Have I done something wrong ?

Best regards
BenJ
BenJ is offline   Reply With Quote

Old   August 15, 2014, 14:49
Default
  #20
Member
 
Brock Lee
Join Date: Sep 2012
Location: Midwest
Posts: 38
Rep Power: 6
GRAUPS is on a distinguished road
Quote:
It does yield a static(p) field, but the values are very low (1e-16).
Ben,

I just figured this out not too long ago. You have to specify a rho value for pressure tools to use to convert from kinematic to static pressure. Add these two lines to your pressure tools functionObject...

rhoName rhoInf;
rhoInf 1.225; // density value

______________________________________

Gentlemen,

I have an additional question related to this pressure tools function object. Is there a way to calculate static pressure for every time step and not have it write out the iteration folder every time step? I want to be able to monitor a mass-weighted static pressure at an outlet every time-step (for an in-compressible solver). I can't find a way to do this though without forcing pressure tools to write an iteration folder every time step. Here is my current pressure-tools setup...

Code:
    p_tools
    {
        type                pressureTools;
        functionObjectLibs  ("libutilityFunctionObjects.so");        
        enabled             yes;        
        outputControl       timeStep;
        outputInterval      1; 
        rhoName              rhoInf;
        rhoInf              831.1;
        calcTotal           no;
        calcCoeff           no;
    }
Thanks for any help you guys can provide!

Brock
GRAUPS is offline   Reply With Quote

Reply

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
compressible flow in turbocharger riesotto OpenFOAM 50 May 26, 2014 01:47
channelFoam for a 3D pipe AlmostSurelyRob OpenFOAM 3 June 24, 2011 13:06
ParaView for OF-1.6-ext Chrisi1984 OpenFOAM Installation 0 December 31, 2010 07:42
Droplet Evaporation Christian Main CFD Forum 2 February 27, 2007 07:27
REAL GAS UDF brian FLUENT 6 September 11, 2006 08:23


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