[alf12]

Hi,

I would like to use the Pe tool, I use pisoFoam solver with kOmegaSST turbulence model.

When I run the Pe tool from my case directory, I get the following error message :

    Reading phi
    Calculating Pe
Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kOmegaSST


--> FOAM FATAL ERROR: 

    lookup of RASProperties from objectRegistry region0 successful
    but it is not a RASModel, it is a dictionary
...

What is going wrong ?

Thanks

[alf12]

I still haven't found what is wrong. When I am trying to use 'Pe' tool with cavity tutorial (incompressible/pisoFoam/ras/cavity), I get the same error. Should that work ?

[chegdan]

What version of openfoam are you using?

[alf12]

I was using version 1.6 and now I have switched to 1.7.1 but the issue is still present.

I have skirted the issue by commenting line 103 in the Pe.C utility (actually I have done another tool based on the Pe tool to deal only with RAS). I hope I have not messed up calculations.

IOdictionary RASProperties(RASPropertiesHeader);

[chegdan]

Originally Posted by alf12

I was using version 1.6 and now I have switched to 1.7.1 but the issue is still present.

I have skirted the issue by commenting line 103 in the Pe.C utility (actually I have done another tool based on the Pe tool to deal only with RAS). I hope I have not messed up calculations.

IOdictionary RASProperties(RASPropertiesHeader);

It seems that there has been a change in teh structure of the RASProperties file since the Pe tool has been made. The RASProperties file looks different from 1.5-dev to 1.6.x and 1.7.x and without the turbulence model coefficients being defined in the RASProperties file in later versions. It seems to be that your change doesn't affect the calculation since the conculation actually uses a laminartransport model to determine the nueff in the Pe calculation...rather than one calculated from a turbulence model.

this is seen in the lines

                IOdictionary RASProperties(RASPropertiesHeader);

                singlePhaseTransportModel laminarTransport(U, phi);

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

                PePtr.set
                (
                    new surfaceScalarField
                    (
                        IOobject
                        (
                            "Pe",
                            runTime.timeName(),
                            mesh,
                            IOobject::NO_READ
                        ),
                        mag(phi)
                       /(
                            mesh.magSf()
                          * mesh.surfaceInterpolation::deltaCoeffs()
                          * fvc::interpolate(RASModel->nuEff())
                         )
                    )
                );

where the RASModel uses the laminarTransport transport property. This seems to be ok as long as there is a nu in the transport properties file to read into the Pe tool. hope that helps.

Dan

[maysmech]

Hi,
Is it any solution found using Pe utility?

[alf12]

Is it any solution found using Pe utility?

I'm still using version 1.7.1, so nothing has changed for me. Solution is to comment line 103 and to recompile Pe tool (likely as a separate tool)