[Sylv]

Dear foamer,

I am using pimpleFoam to run some LES case. As you know, a key element of an LES is the filter width, which is known as delta in OpenFOAM (also in most of the litterature). This delta, is simply a volScalarField

For my simulation, I want to write this volScalarField delta.

In pimpleFoam, the turbulence model is loaded with

autoPtr<incompressible::turbulenceModel> turbulence
(
    incompressible::turbulenceModel::New(U, phi, laminarTransport)
);

The class turbulenceModel has no access to delta, but delta MUST exist somewhere. How can I access it??

Cheers,
Marcel

[dkxls]

Originally Posted by Sylv

The class turbulenceModel has no access to delta, but delta MUST exist somewhere. How can I access it??

The delta field is inside the LESModel class:
src/turbulenceModels/incompressible/LES/LESModel/LESModel.H

You cannot access the LESModel through the "turbulence pointer" (at least to my knowledge), but you could lookup the delta field (or the LESmodel for that matter) from the objectRegistry. Here is how you can lookup the delta field:

    volScalarField LESdelta
    (
        IOobject
        (
            "LESdelta",
            runTime.timeName(),
            mesh,
            IOobject::NO_READ,
            IOobject::NO_WRITE
        ),
        mesh,
        dimensionedScalar("LESdelta", dimLength, SMALL)
    );

    if (mesh.objectRegistry::foundObject<volScalarField>("delta"))
    {
        const volScalarField& delta =
            mesh.objectRegistry::lookupObject<volScalarField>("delta");

        LESdelta.internalField() = delta.internalField();
    }

Note that I haven't tested this code in particular, but I'm doing it in a similar way in some of my own classes.

[Sylv]

Thank for reply dkxls. your suggestion works perfectly.

Marcel

[pmarreiro]

Originally Posted by Sylv

Thank for reply dkxls. your suggestion works perfectly.

Marcel

Marcel, can you tell me where did you implemented the code?

Kind regards,
Pedro