[Paul_AERO]

Greetings,

I am a beginner of OF. I have trouble compiling user defined variable (pprime) in the system/controldict. (version: v2012)

functions
{
pprime
{

type coded;
libs (utilityFunctionObjects);
name pprime;

writeControl timeStep;
writeInterval 1;

codeWrite //(or codeExecute)
#{
const volScalarField& p = mesh().lookupObject<volScalarField>("p");
const volVectorField& U = mesh().lookupObject<volVectorField>("U");
volScalarField pprime
(
IOobject
(
"pprime",
mesh().time().timeName(),
p.mesh(),
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
p
);
pprime=p-101325;
#};

fields
(
pprime
{
base time;
}
);
}

}

I referred to Create registered object at runtime using a functionObject

But mine does not work, showing:
--> Could not load "...

Before running, is there any command to compile?
My current sequence to simulate is mpirun -np 16 rhoPimpleFoam -parallel

Thanks,

Paul

[clapointe]

Without seeing the full error it is hard to help. However, if we check out the codedFunctionObject.H file, we find the following template :

    {
        type coded;
        libs (utilityFunctionObjects);

        // Name of on-the-fly generated functionObject
        name writeMagU;
        codeWrite
        #{
            // Lookup U
            const volVectorField& U = mesh().lookupObject<volVectorField>("U");

            // Write
            mag(U)().write();
        #};
    }

So if we update this for your case, we get :

    pprime
    {
        type coded;
        libs (utilityFunctionObjects);

        // Name of on-the-fly generated functionObject
        name writePPrime;
        codeWrite
        #{
            // Lookup p, name pprime
            const volScalarField& pprime = mesh().lookupObject<volScalarField>("p");

            //Update pprime
            pprime -= dimensionedScalar("refP", pprime.dimensions(), 101325);

            // Write
            pprime.write();
        #};
    }

Note that I've not tested the above code, but if it doesn't work out of the box hopefully it'll point you in the right direction.

Caelan

[Paul_AERO]

Thanks Caelan!

I have resolved that issue

Paul