[goguler]

Dear All,

I have been working on an interFoam based solver nowadays. I cannot understand why do we need pCorrTypes in interFoam. Because, we do not have pCorrTypes for example in pimpleDyMFoam.

Let me write the details.

In interFoam pcorr is defined in correctPhi.H as

volScalarField pcorr
    (
        IOobject
        (
            "pcorr",
            runTime.timeName(),
            mesh,
            IOobject::NO_READ,
            IOobject::NO_WRITE
        ),
        mesh,
        dimensionedScalar("pcorr", pd.dimensions(), 0),
        pcorrTypes
    );

Before that we have declaration of pCorrTypes as follows:

wordList pcorrTypes
    (
        pd.boundaryField().size(),
        zeroGradientFvPatchScalarField::typeName
    );

    forAll (pd.boundaryField(), i)
    {
        if (pd.boundaryField()[i].fixesValue())
        {
            pcorrTypes[i] = fixedValueFvPatchScalarField::typeName;
        }
    }

However, in pimpleDyMFoam pcorr is declared as follows:

volScalarField pcorr("pcorr", p);
pcorr *= 0;

I cannot get the idea why we cannot define pcorr in the same way with pimpleDyMFoam in interFoam.

It might be related to definition of pressure in interFoam vs pimpleDyMFoam. In the first one it is pd and in the second one it is p.

I would be happy to read your comments.