Foam Error - Turbulence Model?
 

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

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?

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 :)

I don't see anything wrong, apart from this

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?

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

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...

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::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::DICPreconditioner::calcReciprocalD(Foam::Fie ld<double>&, Foam::lduMatrix const&) at ??:?
#4 Foam::DICPreconditioner::DICPreconditioner(Foam::l duMatrix::solver const&, Foam::dictionary const&) at ??:?
#5 Foam::lduMatrix::preconditioner::addsymMatrixConst ructorToTable<Foam::DICPreconditioner>::New(Foam:: lduMatrix::solver const&, Foam::dictionary const&) at ??:?
#6 Foam::lduMatrix::preconditioner::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

Error in OpenFOAM
 

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::printStack(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::primitiveEntry::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)

Hi all,
It would be good to share with us how you solved the problems you had...
Can you tell us how you solved them?

I am getting the following error:

--> FOAM FATAL ERROR:
Attempt to return primitive entry ITstream : /home/guillaume/OpenFOAM/Whirljet/WhirlJet7/0/epsilon.boundaryField.bc_outlet, line 37, IOstream: Version 2.0, format ASCII, line 0, OPENED, GOOD
primitiveEntry 'bc_outlet' comprises
on line 37 the word 'type'
on line 37 the word 'zeroGradient'
as a sub-dictionary

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

FOAM aborting

Thank you for your help,

Guillaume

Guillaume,
did you resolve your error ?
I've got the same in connection with snappyHexMesh
Best wishes
Peter