[polo]

Dear all,

OpenFoam Version: 2.2.1 (newest Release) and 2.2.x
swak4FOAM Version 0.2.4 Install from this post #8 http://www.cfd-online.com/Forums/openfoam-installation/118194-swak4foam-0-2-3-of2-2-x-installation-error.html#post434217

the zip Archive ...

swak is working the way it should I think, I can calculate all fields I want, but the following problem I have:

I calculate div(U*U) via the new codedFunctionObject in OpenFOAM and the same I do with swak, but the Results are different!? I don't know why??

I split the Problem and what I can figure out is, that the div operation seems to be make the difference...

so calculate U*U via codedFunctionObject gives the same values as swak
also mag(U*U), grad(U*U), curl(U*U) gives same values for both calculation methods.

BUT if i do div(U*U) and what I want finally to do div(div(U*U)) the values are different between codedFunctionObject and swak?

Have anybody an explanation? What could be the Problem? I've read the swak Reference guide and also the README, but I can't figure out what I do wrong?

One other question, is swak searching for div Scheme in fvScheme?? I'm not absolutly sure but I recognize that if I calculate div(U*U) via codedFunctionObject and have in my fvScheme div(U*U) Gauss linear; it changed if I place there div(U*U) Gauss cubic;

Finally question: which values should I trust? values (field) I get from codedFunctionObject or from swak?

this is my swak expression

DivdivUU
           {
                type expressionField;
                outputControl   timeStep; // timeStep or outputTime
                outputInterval 1;
                fieldName DivdivUU_swak;

                expression "div(div(U*U))";
                autowrite  true;

           }

Thanks for your assistance!

Best Regrads!

Ps: I prefere to use swak, because its easier for me in use!