CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Post-Processing (http://www.cfd-online.com/Forums/openfoam-post-processing/)
-   -   Missing output of sampledSet (http://www.cfd-online.com/Forums/openfoam-post-processing/125418-missing-output-sampledset.html)

letzel October 25, 2013 03:31

Missing output of sampledSet
 
With the controlDict of tutorials/incompressible/icoFoam/cavity modified as follows,
Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.2.2                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application    icoFoam;

startFrom      startTime;

startTime      0;

stopAt          endTime;

endTime        0.5;

deltaT          0.005;

writeControl    timeStep;

writeInterval  20;

purgeWrite      0;

writeFormat    ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision  6;

runTimeModifiable true;


functions
(
  profile01
  {
    type sets;
    enabled true;
    outputMode timeStep;
    outputInterval 1;
    verbose false;
    set (
      type uniform;
      axis y;
      start ( 0.5 0 0.05 );
      end  ( 0.5 1 0.05 );
      nPoints ( 100 );
        );
    setFormat xmgr;
    fields ( U );
  }

);

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

the expected run-time output of the sampledSet is not written to disc, but the screen log shows no error:
Code:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading transportProperties

Reading field p

Reading field U

Reading/calculating face flux field phi


Starting time loop

Time = 0.005

Courant Number mean: 0 max: 0
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 2.96338e-06, No Iterations 8
DILUPBiCG:  Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG:  Solving for p, Initial residual = 1, Final residual = 7.55402e-07, No Iterations 35
time step continuity errors : sum local = 5.03809e-09, global = 4.40722e-19, cumulative = 4.40722e-19
DICPCG:  Solving for p, Initial residual = 0.523591, Final residual = 9.72352e-07, No Iterations 34
time step continuity errors : sum local = 1.07766e-08, global = 2.93897e-19, cumulative = 7.34619e-19
ExecutionTime = 0.01 s  ClockTime = 0 s

Time = 0.01

Courant Number mean: 0.0976805 max: 0.585722
DILUPBiCG:  Solving for Ux, Initial residual = 0.148584, Final residual = 7.15711e-06, No Iterations 6
DILUPBiCG:  Solving for Uy, Initial residual = 0.256618, Final residual = 8.94127e-06, No Iterations 6
DICPCG:  Solving for p, Initial residual = 0.379232, Final residual = 3.38648e-07, No Iterations 34
time step continuity errors : sum local = 3.15698e-09, global = 1.94222e-19, cumulative = 9.28841e-19
DICPCG:  Solving for p, Initial residual = 0.286937, Final residual = 5.99637e-07, No Iterations 33
time step continuity errors : sum local = 6.08774e-09, global = 5.80474e-19, cumulative = 1.50932e-18
ExecutionTime = 0.01 s  ClockTime = 0 s

Time = 0.015

How can I achieve the expected output?

letzel October 25, 2013 10:49

While I could not solve the problem initially reported for the test case above, I did solve it for a production case. Formatting was wrong.

The following is a working example just for reference. Note that not all libs may be required for this excerpt.
Code:

libs (
        "libOpenFOAM.so"
        "libsampling.so"
        "libsimpleSwakFunctionObjects.so"
        "libswakFunctionObjects.so"
    );
functions
(
  profiles
  {
    type sets;
    enabled true;
    verbose true;
    interpolationScheme cellPoint;
    outputControlMode outputTime;
    setFormat xmgr;
    fields ( wuav wsusav wvav wsvsav Iav );
    sets
    (
      profile01
      {
        type uniform;
        axis y;
        start ( 20  0 56 );
        end  ( 20 50 56 );
        nPoints 100;
      }
      profile02
      {
        type uniform;
        axis y;
        start ( 40  0 56 );
        end  ( 40 50 56 );
        nPoints 100;
      }
    );
  }
);



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