[Selvinus]

Hi there!

I'm currently developing my own LES Dynamic Smagorinsky solver applying Finite Elements. For verifying my own implementation I apply OpenFOAM and compare different values obtained from different fields.

When applying this dynamic smagorinsky method the scalar fields LijMij and MijMij are computed and used in the lagrangian averaging PDEs together with a time scale variable invT. I would like to write these fields together with the standard fields U, p, nuSgs and so on. Any help on how to do this?

I have tried adding

    MM
    (
        IOobject
        (
            "MM",
            runTime_.timeName(),
            mesh_,
            IOobject::MUST_READ,
            IOobject::AUTO_WRITE
        ),
        mesh_
    ),
    invT
    (
        IOobject
        (
            "invT",
            runTime_.timeName(),
            mesh_,
            IOobject::MUST_READ,
            IOobject::AUTO_WRITE
        ),
        mesh_
    ),

Into my dynLagrangian.C together with

    volScalarField MM;
    volScalarField invT;

into dynLagrangian.H.

When compiling these files and running pisoFoam the fields are written, however there are no values! They are zero all over the place.. Any help?

Thanks!

Selvinus