[matthias]

Hi folks,

I'm using the coded function object to create and calculate the field epsilon

   epsilonAverage
    {
        functionObjectLibs ("libutilityFunctionObjects.so");
        type coded;
        redirectType average1;
        outputControl timeStep;
        code
        #{
            volScalarField epsilon
            (
                IOobject
                (
                    "epsilon",
                    mesh().time().timeName(),
                    mesh(),
                    IOobject::NO_READ,
                    IOobject::AUTO_WRITE
                ),
                2*mesh().lookupObject<volScalarField>("nu")*magSqr(symm(fvc::grad(mesh().lookupObject<volVectorField>("U"))))
            );

            //epsilon.write();
        #};

    }

When I call the fieldAverage function object and apply the averaging to epsilon, it isn't known since the field epsilon is not registered globally but only within the code directive.

So is it possible to use the fieldAverage function object for this purpose?

(Of course I can implement the code into the solver directly but using the coded function object I must not change the solver)


Best regards

Matthias

[wyldckat]

Greetings Matthias,

I'm not certain I understand the problem you're describing, but my understanding is that you're using something like the utility execFlowFunctionObjects and using that function object to calculate the average of epsilon field!? The part I'm not certain is why you're using a coded function object...

OK, there are two possible answers, from what I can figure out from your question:

1. Use the option "-noFlow" when calling execFlowFunctionObjects: http://openfoamwiki.net/index.php/Ex...unctionObjects
2. Or use another function object to load up the fields you need (before the one you're using):
   

       readFields
       {
           functionObjectLibs ( "libfieldFunctionObjects.so" );
           type readFields;
           fields (p U epsilon);
       }

Best regards,
Bruno

[matthias]

Hi Bruno,

actually I don't use the execFlowFunctionObjects. My intention was to calculate the dissipation field (epsilon) on-the-fly without modifying the solver.

I can calculate the epsilon field in system/controlDict using the coded function object but how does the averaging work? After the coded function object has been finished the epsilon field can be written to disk but it is not registered! If you call the averaging function object after coded function object it won't find a field named epsilon.

I will have a look at execFlowFunctionObjects if it will do the job.

Best regards

Matthias

[fabian_roesler]

Hi Matthias

Have a look into this thread:

http://www.cfd-online.com/Forums/ope...ionobject.html

In post #12 (http://www.cfd-online.com/Forums/ope...tml#post387459) Eelco explains the solution to his problem using functionobjects.

Cheers

Fabian