CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Solving problem in "cavity with pisoFoam LES"

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   August 10, 2010, 15:40
Question Solving problem in "cavity with pisoFoam LES"
  #1
Senior Member
 
maysmech's Avatar
 
Join Date: Jan 2010
Posts: 347
Blog Entries: 2
Rep Power: 8
maysmech is on a distinguished road
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
maysmech is offline   Reply With Quote

Old   September 30, 2010, 00:37
Default
  #2
Member
 
Join Date: Mar 2009
Location: Sydney, New South Wales, Australia
Posts: 42
Rep Power: 8
rassilon is on a distinguished road
Quote:
Originally Posted by maysmech View Post
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
rassilon is offline   Reply With Quote

Old   October 1, 2010, 17:15
Default
  #3
Senior Member
 
maysmech's Avatar
 
Join Date: Jan 2010
Posts: 347
Blog Entries: 2
Rep Power: 8
maysmech is on a distinguished road
Quote:
Originally Posted by rassilon View Post
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.
maysmech is offline   Reply With Quote

Old   May 2, 2011, 17:11
Default
  #4
New Member
 
Arnout
Join Date: Nov 2010
Posts: 23
Rep Power: 6
The King is on a distinguished road
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;
}
The King is offline   Reply With Quote

Old   July 24, 2012, 05:12
Default
  #5
New Member
 
Join Date: Sep 2011
Posts: 13
Rep Power: 5
anon_g is on a distinguished road
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...
anon_g is offline   Reply With Quote

Old   July 24, 2012, 05:24
Default
  #6
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 12
Bernhard is on a distinguished road
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.
Bernhard is offline   Reply With Quote

Old   July 24, 2012, 06:13
Default
  #7
New Member
 
Join Date: Sep 2011
Posts: 13
Rep Power: 5
anon_g is on a distinguished road
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?

Last edited by anon_g; July 24, 2012 at 06:45.
anon_g is offline   Reply With Quote

Old   July 24, 2012, 06:38
Default
  #8
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 12
Bernhard is on a distinguished road
Probably z=0 is exactly on a cell face, so you should try a slightly bigger scalar there as well.
Bernhard is offline   Reply With Quote

Reply

Tags
cavity, les, pisofoam

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33
lift and drag on ship superstructures vaina74 OpenFOAM Running, Solving & CFD 3 June 8, 2010 12:30
MRFSimpleFOAM goes divergenced! renyun0511 OpenFOAM Running, Solving & CFD 0 November 19, 2009 03:11
Negative value of k causing simulation to stop velan OpenFOAM Running, Solving & CFD 1 October 17, 2008 05:36
Modeling in micron scale using icoFoam m9819348 OpenFOAM Running, Solving & CFD 7 October 27, 2007 00:36


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