[sur4j]

Hello,

I have recently created a simulation which has a temperature dependant viscosity: (http://www.tfd.chalmers.se/~hani/kur...nFoam%20v2.pdf)

The solver and test case both worked fine and compiled okay, however when I tried to modify the simulation to apply this temperature dependant viscosity on the cavity tutorial geometry I had gotten an error:

--> FOAM FATAL IO ERROR: 
Unable to set reference cell for field p
    Please supply either pRefCell or pRefPoint


file: /home/sk/OpenFOAM/sk-2.2.1/run/tvCavity/system/fvSolution.SIMPLE from line 80 to line 87.

    From function void Foam::setRefCell
(
    const volScalarField&,
    const volScalarField&,
    const dictionary&,
    label& scalar&,
    bool
)
    in file cfdTools/general/findRefCell/findRefCell.C at line 125.

FOAM exiting

The code in fvSolution lines 80-87 is the following:

SIMPLE
{
    nNonOrthogonalCorrectors 0;

    residualControl
    {
	T 		1e-2;
        p               1e-2;
        U               1e-3;
        "(k|epsilon|omega)" 1e-3;
    }
}

Again, the tutorial case worked fine, this error only occured when I modified the blockMeshDict and boundary conditions to work with the cavity test case. This has really confused me and I would appreciate any help with this problem.

Thank you.

[Bernhard]

Add

    pRefCell        0;
    pRefValue       0;

To your system/fvSolution in the SIMPLE subdictionary.

You are changing more things, as you are trying to solve using a SIMPLE algorithm (apparently), which looks up different subdictionaries.