[samiam1000]

Dear Foamers,

I am writing a function object and I do have some issues to be solved. I hope someone can help:

I am attaching a very simple case. Could you have a look?

I case: Enter the case directory and launch: "scalarTransportFoam"
Everything is ok. In the 10 dir I have the T*2 field.

II case: Enter the case directory and launch: "scalarTransportFoam -noFunctionObjects" and as a post processing run "execFlowFunctionObjects -noZero -latestTime" I get this error:

samuele@pc-zampini:~/Documenti/PhD/02tesiDomenico/workingDirLocal/dummyCase/passiveScalarBaseGeometry$ execFlowFunctionObjects -noZero -latestTime
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.3.0-f5222ca19ce6
Exec   : execFlowFunctionObjects -noZero -latestTime
Date   : Jul 31 2014
Time   : 10:53:20
Host   : "pc-zampini"
PID    : 21188
Case   : /home/samuele/Documenti/PhD/02tesiDomenico/workingDirLocal/dummyCase/passiveScalarBaseGeometry
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 10

Using dynamicCode for functionObject difference at line 57 in "/home/samuele/Documenti/PhD/02tesiDomenico/workingDirLocal/dummyCase/passiveScalarBaseGeometry/system/controlDict.functions.difference"
Time = 10
    Reading phi
    Reading U
    Reading p
--> FOAM Warning : 
--> FOAM FATAL IO ERROR: 
cannot find file

file: /home/samuele/Documenti/PhD/02tesiDomenico/workingDirLocal/dummyCase/passiveScalarBaseGeometry/10/p at line 0.

    From function regIOobject::readStream()
    in file db/regIOobject/regIOobjectRead.C at line 73.

samuele@pc-zampini:~/Documenti/PhD/02tesiDomenico/workingDirLocal/dummyCase/passiveScalarBaseGeometry$

Hence, I copy a "p" file in the 10 directory. Nevertheless, I get this:

samuele@pc-zampini:~/Documenti/PhD/02tesiDomenico/workingDirLocal/dummyCase/passiveScalarBaseGeometry$ execFlowFunctionObjects -noZero -latestTime
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.3.0-f5222ca19ce6
Exec   : execFlowFunctionObjects -noZero -latestTime
Date   : Jul 31 2014
Time   : 10:52:07
Host   : "pc-zampini"
PID    : 21126
Case   : /home/samuele/Documenti/PhD/02tesiDomenico/workingDirLocal/dummyCase/passiveScalarBaseGeometry
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 10

Using dynamicCode for functionObject difference at line 57 in "/home/samuele/Documenti/PhD/02tesiDomenico/workingDirLocal/dummyCase/passiveScalarBaseGeometry/system/controlDict.functions.difference"
Time = 10
    Reading phi
    Reading U
    Reading p
Looking up field T



--> FOAM FATAL ERROR: 

    request for volScalarField T from objectRegistry region0 failed
    available objects of type volScalarField are
1(p)

    From function objectRegistry::lookupObject<Type>(const word&) const
    in file /home/samuele/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 198.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) in "/home/samuele/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::error::abort() in "/home/samuele/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const& Foam::objectRegistry::lookupObject<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> >(Foam::word const&) const in "/home/samuele/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#3  Foam::mixingEvaluationFunctionObject::write() at /usr/include/c++/4.7/bits/basic_string.h:536
#4  Foam::OutputFilterFunctionObject<Foam::mixingEvaluationFunctionObject>::execute(bool) at ~/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/OutputFilterFunctionObject.C:173
#5  Foam::functionObjectList::execute(bool) in "/home/samuele/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6  
 in "/home/samuele/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/execFlowFunctionObjects"
#7  
 in "/home/samuele/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/execFlowFunctionObjects"
#8  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9  
 in "/home/samuele/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/execFlowFunctionObjects"
Annullato
samuele@pc-zampini:~/Documenti/PhD/02tesiDomenico/workingDirLocal/dummyCase/passiveScalarBaseGeometry$

Could anyone help?

Also, one more question: how abaout if I want to get an average on each section?

Thanks a lot,
Samuele

[Bernhard]

You also need to -noFlow flag for execFlowFunctionObjects. All variables will then be loaded.

[samiam1000]

Hi Bernard and thanks for helping.

Also, about my secondo issue, I was thinking about exploiting the sample function and then to average it.

The point is that I would insert a loop in the sampleDict. Is that possible?

I wrote something like:

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

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

interpolationScheme cellPoint; 
   
  surfaceFormat       foamFile; 
points = 100;
L = 3000;

dx = $points/$L;

for i = 1:$points
{   
  surfaces 
  ( 
      plane%d, i
      { 
		  type	plane;
          basePoint (0.0001 -0.0016+($dx*%i)  0.0001); 
          normalVector (0 1 0); 
      }
 
  ); 
}
end
   
fields          ( U );

Any idea?

Thanks again.

[samiam1000]

I tryed the following solution:

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

application     scalarTransportFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         10; 		// T = 7.5e-3

deltaT          1;

writeControl    adjustableRunTime;

writeInterval   1;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

adjustTimeStep	yes;

maxCo			0.9;


functions
{
new_planes
{
    type            faceSource;
    functionObjectLibs ("libfieldFunctionObjects.so");
    enabled         true;
    outputControl   timeStep;
    log             true;
    valueOutput     true;
    source          sampledSurface;
    surfaceFormat   vtk;

	sampledSurfaceDict
	{
		type            plane;
		basePoint       (0.0001 -0.00115 0.0001);
		normalVector    (0 1 0);
	}

	sampledSurfaceDict
	{
		type            plane;
		basePoint       (0.0001 -0.00065 0.0001);
		normalVector    (0 1 0);
	}

    operation       areaAverage;
	fields
	(
		T
	);
}

The point is that it does not allow me to write more than one value of T. How can I do this?