CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Solving problem in "cavity with pisoFoam LES" (http://www.cfd-online.com/Forums/openfoam-solving/79046-solving-problem-cavity-pisofoam-les.html)

maysmech August 10, 2010 15:40

Solving problem in "cavity with pisoFoam LES"
 
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:


HTML Code:


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 September 30, 2010 00:37

Quote:

Originally Posted by maysmech (Post 271030)
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:


HTML Code:


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 October 1, 2010 17:15

Quote:

Originally Posted by rassilon (Post 277204)
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 May 2, 2011 17:11

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;
}

rabbit July 24, 2012 05:12

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 July 24, 2012 05:24

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.

rabbit July 24, 2012 06:13

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 July 24, 2012 06:38

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


All times are GMT -4. The time now is 20:05.