https://linen.dev logo
Join Discord
Powered by
# function-objects
  • y
    and what do you get in your log file? You shoud see if the function is read and/or if you get an error/warning
  • b
    Yea, but... During the solution steps, nothing appears.
  • b
    I won't be able to see the beginning of the solution, because I didn't generate the .log. lol
  • y
    not a good habit 😂
  • y
    alright, for debugging try replacing ``writeControl outputTime;`` by ``writeControl timeStep;``, and show your controlDict file
  • b
    I'm lying, I managed to evaluate it, NOTHING appears. 😅
  • b
    Ok
  • b
    One minute
  • y
    then there is something wrong with the way you define / include your functions in your controlDict
  • b
    Copy code
    /--------------------------------- C++ -----------------------------------\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  9
     \\/     M anipulation  |
\---------------------------------------------------------------------------/
FoamFile
{
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

libs ();
     
application     multiphaseEulerFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         0.04; //~40000

deltaT          1e-6; 

writeControl    runTime;

writeInterval   4e-5; 

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable yes;

adjustTimeStep  no;

maxCo           0.5;

//maxAlphaCo      0.5;

maxDeltaT       0.01;

// ************************* //

cacheTemporaryObjects
();
    
functions
{
    //$surfaceFieldValueFunctions;   
    //#includeFunc Q(U.water)    
    
    #includeFunc faceZoneAverage(name=monitorZone, phi.water)
    #includeFunc faceZoneFlowRate(name=monitorZone)
    
    massFlowRate
    {
        type            surfaceFieldValue;
        functionObjectLibs ("libfieldFunctionObjects.so");
        writeControl    outputTime;
        enabled         true;
        log             true;
        writeFields     false;
        regionType      sampledSurface;
        name            monitorZone;
        operation       sum;
        fields          (phi.water phi.air);
        
        sampledSurfaceDict
        {
            type    faceZoneSet;
            setName monitorZone;
        }
    }
    
    fieldAverage
    {
        type            fieldAverage;
        functionObjectLibs ( "libfieldFunctionObjects.so" );
        writeControl    outputTime;        
        enabled         true;
        timeStart       0.01; //10000

        fields
        (
            U.air
            {
                 mean        on;
                 prime2Mean  on;
                 base        time;
            }

            U.water
            {
                 mean        on;
                 prime2Mean  on;
                 base        time;
            }

            alpha.air
            {
                 mean        on;
                 prime2Mean  off;
                 base        time;
            }

            alpha.water
            {
                 mean        on;
                 prime2Mean  off;
                 base        time;
            }
            
            nut.water
            {
                 mean        on;
                 prime2Mean  on;
                 base        time;
            }

            p
            {
                 mean        on;
                 prime2Mean  off;
                 base        time;
            }
            
            k.water
            {
                 mean        on;
                 prime2Mean  on;
                 base        time;
            }
        );
    }
    
}
    
surfaceFieldValueFunctions
{}

// ************************* //
  • b
  • b
    Copy code
    writeControl    outputTime;
    works...
  • y
    Aright, here is your problem:
    Copy code
    --> FOAM Warning : 
    From function bool Foam::functionObjectList::read()
    in file db/functionObjects/functionObjectList/functionObjectList.C at line 853
    Caught FatalError 
--> FOAM FATAL ERROR: 
Unknown sample type faceZoneSet

Valid sample types : 

8
(
cuttingPlane
distanceSurface
isoSurface
patch
patchInternalField
plane
thresholdCellFaces
triSurfaceMesh
)
  • y
    faceZoneAverage and faceZoneFlowRate seems to work just fine
    Copy code
    surfaceFieldValue faceZoneAverage(name=monitorZone,phi.water):
    total faces  = 2841
    total area   = 0.0366263


surfaceFieldValue faceZoneFlowRate(name=monitorZone):
    total faces  = 2841
    total area   = 0.0366263
  • b
    Yes, I saw this error, but it wasn't showing up before... the way I called the function was incorrect. That's why nothing was showing up. However, the options it offers me don't work for me. I created the region, but I don't know how to reference it.
  • b
    Yea
  • b
    In fact, one function would do the same as the other, it would be redundant.
  • y
    you can use ``foamInfo`` to get some information about functions, boundary conditions, etc. For instance: ``foamInfo surfaceFieldValue``
  • b
    Now, he's calculating the average... I don't know if he would have to calculate the sum.
  • y
    Since you have a faceZone, you should use:
    Copy code
    regionType      faceZone;
        name            monitorZone;
  • b
    I will try
  • y
    I'm surprised that faceZoneAverage on phi.water gives you the same result as faceZoneFlowRate
  • b
    It's because there was still air in there.
  • b
    I believe it will go wrong because
  • b
    Copy code
    Unknown sample type faceZoneSet

Valid sample types : 

8
(
cuttingPlane
distanceSurface
isoSurface
patch
patchInternalField
plane
thresholdCellFaces
triSurfaceMesh
)
  • y
    FYI, these functions (faceZoneAverage, faceZoneFlowRate, ...) are defined in ``$FOAM_ETC/caseDicts/postProcessing/``. It can be interesting to check how they are defined
  • y
    You get this error because you are using a ``sampledSurface`` as ``regionType`` rather than ``faceZone``
  • b
    OOOh... ok, thanks.
  • y
    ``regionType`` allows to define what kind of entity you want to use (faceZone, patch or sampledSurface). And depending on what you select there, you need to define either the name of the patch or faceZone, or a sampleSurfaceDict to use sampling functions
  • b
    Yes, I was not paying enough attention... 😅