[maysmech]

Hi FOAMers;
i want solve lid driven cavity with LES. i used files in tutorial file pisoFoam/LES/.
After creating cavity mesh and adjusting boundary conditions in 0 and polymesh folders and typing pisoFoam leads to:

Create time

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi



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


file: /home/maysam/OpenFOAM/maysam-1.7.0/cavityRatio/10000LES/system/fvSolution::PISO from line 71 to line 72.

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

FOAM exiting

Please help me to find the problem solution.

Best,
Maysam

[rassilon]

Originally Posted by maysmech

Hi FOAMers;
i want solve lid driven cavity with LES. i used files in tutorial file pisoFoam/LES/.
After creating cavity mesh and adjusting boundary conditions in 0 and polymesh folders and typing pisoFoam leads to:

Create time

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi



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


file: /home/maysam/OpenFOAM/maysam-1.7.0/cavityRatio/10000LES/system/fvSolution::PISO from line 71 to line 72.

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

FOAM exiting

Please help me to find the problem solution.

Best,
Maysam

Did you work through this problem? I am encountering the same thing with interFoam now, and would be interested to know how to resolve it.


R

[maysmech]

Originally Posted by rassilon

Did you work through this problem? I am encountering the same thing with interFoam now, and would be interested to know how to resolve it.


R

hi,
i think its problem was for using system folder of another solver instead of pisoFoam and fvsolution and etc was not proper for pisoFoam.

[The King]

Regarding the error:

The solver needs to know what and where the reference pressure is. Add the following to your fvsolution file:

PISO
{
...
pRefPoint (-0.081 -0.0257 8.01);
pRefValue 1e5;
}

[anon_g]

hi everyone, this thread is probably "closed" because "solved",...

nevertheless i have a problem concerning just this pressure referencing, i hope somebody can help:

i'd like to sim a channel flow with pisoFoam (ras - if that maybe the cause of the problem, which i think is/should be unlikely), my checkMesh output gives

...Overall domain bounding box (0 0 -0.053) (0.99 0.031 0.053)...

for my referencing i have in fvSolution

PISO
{
...
pRefPoint (0.99 0.031 0);
pRefValue 8e5;
}


still, when i run the solver, i get the following error message:

--> FOAM FATAL IO ERROR:
Unable to set reference cell for field p
Reference point pRefPoint (0.99 0.031 0.053) found on 0 domains (should be one)

file: /home/users/<...>/run/multiphase/<...>/system/fvSolution::PISO from line 157 to line 167.
From function void Foam::setRefCell
(
const volScalarField&,
const volScalarField&,
const dictionary&,
label& scalar&,
bool
)
in file cfdTools/general/findRefCell/findRefCell.C at line 95.
FOAM exiting

***
any suggestions?
i should mention that i'm not quite an expert on code manipulation, if there might be a hint in the error output on what is actually going wrong , i might not have seen it. in that case, i would appreciate it if someone could shed some light...

[Bernhard]

It is a bit tricky, but the point that you use to define the reference cell, is on the boundary, and not inside the domain. Something like pRefPoint (0.899 0.0309 0); should select the desired cell.

[anon_g]

thx for the quick reply...

i thought about that too and tried 0.0309999 but without success,... but now i tried "every" parameter (x, y, z inside) and it worked...

anyone an idea why it doesnt work though? it works for interfoam?!

***EDIT***

is this (= not working if point is ON boundary) a bug or intended?

[Bernhard]

Probably z=0 is exactly on a cell face, so you should try a slightly bigger scalar there as well.