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

Foam Error - Turbulence Model?

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

Like Tree2Likes
  • 1 Post By student666
  • 1 Post By student666

Reply
 
LinkBack Thread Tools Display Modes
Old   July 29, 2014, 10:05
Question Foam Error - Turbulence Model?
  #1
Member
 
Gowain
Join Date: May 2014
Location: Ireland
Posts: 86
Rep Power: 5
eaglemckenna is on a distinguished road
Hi
I am trying to do a thermal evaluation on the human body (see attached image for idea of geometry) - I am using buoyantBoussinesqueSimpleFoam solver and receive the following error - I wish to use kepsilon model for the human body - I am relatively new to openFoam and would appreciate some guidance and insight into the following error message:

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Creating turbulence model

Selecting RAS turbulence model kEpsilon


--> FOAM FATAL IO ERROR:
Cannot find patchField entry for body

file: /home/gowain/OpenFOAM/OpenFOAM-2.3.0/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/humanmodel/0/k.boundaryField from line 29 to line 30.

From function GeometricField<Type, PatchField, GeoMesh>::GeometricBoundaryField::readField(const DimensionedField<Type, GeoMesh>&, const dictionary&)
in file /home/gowain/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/GeometricBoundaryField.C at line 209.

FOAM exiting

many thanks
g
Attached Images
File Type: jpg sim1prisvoldivide.jpg (74.0 KB, 21 views)
eaglemckenna is offline   Reply With Quote

Old   July 29, 2014, 10:16
Default
  #2
Senior Member
 
M. C.
Join Date: May 2013
Location: Italy
Posts: 260
Rep Power: 10
student666 is on a distinguished road
I think you defined a "body" patch in constant/polymesh/boundary, but you're missing the same on 0 folder in k file...read carfeully OF always says where you're wrong...at least for the 99% of the cases .

did you obtain the mesh by using Engrid?
eaglemckenna likes this.
student666 is offline   Reply With Quote

Old   July 29, 2014, 11:34
Default
  #3
Member
 
Gowain
Join Date: May 2014
Location: Ireland
Posts: 86
Rep Power: 5
eaglemckenna is on a distinguished road
Hi - thanks for the response
I created the model on blender and meshed the domain and body with engrid (see attached)
I thought that was the problem also but it appears not, unless my patch is wrongly defined, here are my k and boundary files:

FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/polyMesh";
object boundary;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
9
(
back
{
type empty;
nFaces 1725;
startFace 3014960;
}
body
{
type wall;
nFaces 84480;
startFace 3016685;
}
bottom
{
type empty;
nFaces 1719;
startFace 3101165;
}
front
{
type empty;
nFaces 1771;
startFace 3102884;
}
inlet
{
type patch;
nFaces 96;
startFace 3104655;
}
left
{
type empty;
nFaces 1595;
startFace 3104751;
}
outlet
{
type patch;
nFaces 100;
startFace 3106346;
}
right
{
type empty;
nFaces 1831;
startFace 3106446;
}
top
{
type empty;
nFaces 1659;
startFace 3108277;
}
)

and

FoamFile
{
version 2.0;
format binary;
class volScalarField;
location "0";
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0.1;

boundaryField
{
boundaryField
{

body
{
type kqRwallFunction;
value uniform 0.1;
}


}

many thanks
Attached Images
File Type: jpg StandingBodySIM1.jpg (38.8 KB, 10 views)
eaglemckenna is offline   Reply With Quote

Old   July 29, 2014, 12:24
Default
  #4
Senior Member
 
M. C.
Join Date: May 2013
Location: Italy
Posts: 260
Rep Power: 10
student666 is on a distinguished road
I don't see anything wrong, apart from this

Quote:

boundaryField
{
boundaryField
{
2 times. but probably you copy & paste, so this is correct into your file.

The only suggestion I have is to copy one working file from another case.

Furthermore, are you sure you have to define into polyMesh/boundary: top, bottom, left, right, front, back as empty?

If I rembember well, empty is used to reduce the dimension of the problem (cavity case), why don't you use periodic, cyclic or other boundary type?
eaglemckenna likes this.
student666 is offline   Reply With Quote

Old   July 30, 2014, 08:30
Default
  #5
Member
 
Gowain
Join Date: May 2014
Location: Ireland
Posts: 86
Rep Power: 5
eaglemckenna is on a distinguished road
Thanks for the response - I managed to get the error sorted but alas I have encountered another error on the nut file that has me somewhat perplexed - I would be grateful if you could shed some light on it if possible:
--> FOAM FATAL ERROR:
Attempt to return primitive entry ITstream : /home/gowain/OpenFOAM/OpenFOAM-2.3.0/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/humanmodel/0/nut.boundaryField.outlet, line 60, IOstream: Version 2.0, format ASCII, line 0, OPENED, GOOD
primitiveEntry 'outlet' comprises
on line 60 the word 'type'
on line 60 the word 'calculated'
as a sub-dictionary

From function const dictionary& primitiveEntry::dict() const
in file db/dictionary/primitiveEntry/primitiveEntry.C at line 191.

FOAM aborting

thanks so much
apologies again
g
eaglemckenna is offline   Reply With Quote

Old   July 30, 2014, 12:34
Default
  #6
Senior Member
 
M. C.
Join Date: May 2013
Location: Italy
Posts: 260
Rep Power: 10
student666 is on a distinguished road
how did you managed the problem before?

About new one, I think there's a mismatching between calculated and boundary type.

I face problem like this one when I made following error:
on boundary I defined the type patch eg: inlet = patch, while in the 0 folder I set inlet = empty;

basic type are in contrast.

I suggest to check this way...let me know...
student666 is offline   Reply With Quote

Old   July 31, 2014, 16:47
Default
  #7
Member
 
Gowain
Join Date: May 2014
Location: Ireland
Posts: 86
Rep Power: 5
eaglemckenna is on a distinguished road
HI -
I managed to fix the above two errors - it was simple mistake; however I now get a running error when I run the solver buoyantBoussinesqSimpleFoam - I am still relatively new to OF and I am sure the error lies in case files - maybe. Or it is due to my limited computer resources - I intend to run the simulation from starttime 0 to endtime 100 whereby delta t is 0.25 at a write interval of 4 - I am sure this corresponds to a courant number <1
I receive the following error at 6 seconds into the sim - always at 6 seconds:

Time = 6

DILUPBiCG: Solving for T, Initial residual = 0.900926, Final residual = 0.011827, No Iterations 1
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam:ICPreconditioner::calcReciprocalD(Foam::Fie ld<double>&, Foam::lduMatrix const&) at ??:?
#4 Foam:ICPreconditioner:ICPreconditioner(Foam::l duMatrix::solver const&, Foam::dictionary const&) at ??:?
#5 Foam::lduMatrix:reconditioner::addsymMatrixConst ructorToTable<Foam:ICPreconditioner>::New(Foam:: lduMatrix::solver const&, Foam::dictionary const&) at ??:?
#6 Foam::lduMatrix:reconditioner::New(Foam::lduMatr ix::solver const&, Foam::dictionary const&) at ??:?
#7 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#8 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#9 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#10
at ??:?
#11
at ??:?
#12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13
at ??:?
Floating point exception (core dumped)

- Any pointers into what is wrong?
kind regards
g
eaglemckenna is offline   Reply With Quote

Old   December 14, 2016, 13:52
Default Error in OpenFOAM
  #8
New Member
 
Pedram
Join Date: Feb 2016
Location: United States
Posts: 11
Rep Power: 3
pedramtx is on a distinguished road
Dear Foamers,
I have been trying to do a simulation using pisoFoam solver. However, I face the following error while doing the simulation.
Any help is greatly appreciated.

Pedram


Create time

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model kOmega


--> FOAM FATAL ERROR:
Attempt to return primitive entry ITstream : /scratch/user//constant/RASProperties.kOmegaCoeffs, line 24, IOstream: Version 2.0, format ASCII, line 0, OPENED, GOOD
primitiveEntry 'kOmegaCoeffs' comprises
on line 24 the word 'printCoeffs'
on line 24 the word 'on'
as a sub-dictionary

From function const dictionary& primitiveEntry::dict() const
in file db/dictionary/primitiveEntry/primitiveEntry.C at line 191.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) in "/software/easybuild/software/OpenFOAM/2.3.1-intel-2015B/OpenFOAM-2.3.1/platforms/linux64IccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/software/easybuild/software/OpenFOAM/2.3.1-intel-2015B/OpenFOAM-2.3.1/platforms/linux64IccDPOpt/lib/libOpenFOAM.so"
#2 Foam:rimitiveEntry::dict() const in "/software/easybuild/software/OpenFOAM/2.3.1-intel-2015B/OpenFOAM-2.3.1/platforms/linux64IccDPOpt/lib/libOpenFOAM.so"
#3 Foam::dictionary::subOrEmptyDict(Foam::word const&, bool) const in "/software/easybuild/software/OpenFOAM/2.3.1-intel-2015B/OpenFOAM-2.3.1/platforms/linux64IccDPOpt/lib/libOpenFOAM.so"
#4 Foam::incompressible::RASModel::RASModel(Foam::wor d const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/software/easybuild/software/OpenFOAM/2.3.1-intel-2015B/OpenFOAM-2.3.1/platforms/linux64IccDPOpt/lib/libincompressibleRASModels.so"
#5 Foam::incompressible::RASModels::kOmega::kOmega(Fo am::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&, Foam::word const&) in "/software/easybuild/software/OpenFOAM/2.3.1-intel-2015B/OpenFOAM-2.3.1/platforms/linux64IccDPOpt/lib/libincompressibleRASModels.so"
#6 Foam::incompressible::RASModel::adddictionaryConst ructorToTable<Foam::incompressible::RASModels::kOm ega>::New(Foam::GeometricField<Foam::Vector<double >, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/software/easybuild/software/OpenFOAM/2.3.1-intel-2015B/OpenFOAM-2.3.1/platforms/linux64IccDPOpt/lib/libincompressibleRASModels.so"
#7 Foam::incompressible::RASModel::New(Foam::Geometri cField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/software/easybuild/software/OpenFOAM/2.3.1-intel-2015B/OpenFOAM-2.3.1/platforms/linux64IccDPOpt/lib/libincompressibleRASModels.so"
#8 Foam::incompressible::turbulenceModel::addturbulen ceModelConstructorToTable<Foam::incompressible::RA SModel>::NewturbulenceModel(Foam::GeometricField<F oam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/software/easybuild/software/OpenFOAM/2.3.1-intel-2015B/OpenFOAM-2.3.1/platforms/linux64IccDPOpt/lib/libincompressibleRASModels.so"
#9 Foam::incompressible::turbulenceModel::New(Foam::G eometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/software/easybuild/software/OpenFOAM/2.3.1-intel-2015B/OpenFOAM-2.3.1/platforms/linux64IccDPOpt/lib/libincompressibleTurbulenceModel.so"
#10
in "/software/easybuild/software/OpenFOAM/2.3.1-intel-2015B/OpenFOAM-2.3.1/platforms/linux64IccDPOpt/bin/pisoFoam"
#11 __libc_start_main in "/lib64/libc.so.6"
#12
in "/software/easybuild/software/OpenFOAM/2.3.1-intel-2015B/OpenFOAM-2.3.1/platforms/linux64IccDPOpt/bin/pisoFoam"
Aborted (core dumped)
pedramtx is offline   Reply With Quote

Reply

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
Take derivative of mean velocity in paraFoam hiuluom OpenFOAM Paraview & paraFoam 13 April 26, 2016 06:44
How to define to right point for locationInMesh Mirage12 OpenFOAM Native Meshers: snappyHexMesh and Others 7 March 13, 2016 15:07
simpleFoam parallel AndrewMortimer OpenFOAM Running, Solving & CFD 12 August 7, 2015 18:45
Using starToFoam clo OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 33 September 26, 2012 04:04
StarToFoam error Kart OpenFOAM Meshing & Mesh Conversion 1 February 4, 2010 05:38


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