CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   Problem in running ICEM grid in Openfoam (http://www.cfd-online.com/Forums/openfoam/80939-problem-running-icem-grid-openfoam.html)

Tarak October 11, 2010 17:08

Problem in running ICEM grid in Openfoam
 
Hii,

I am having problem in running OF case for a simple rectangular channel, the grid being generated in ICEM and converted to OF format by fluentMeshToFoam utility (the mesh being first converted to .msh from ICEM).

The error I'm getting is:

Create time

Create mesh for time = 0


Reading transportProperties

Reading field p



--> FOAM FATAL ERROR:
Attempt to cast type patch to type lduInterface

From function refCast<To>(From&)
in file /home/44792/OpenFOAM/OpenFOAM-1.7.1/src/OpenFOAM/lnInclude/typeInfo.H at line 114.

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&) in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 Foam::Ostream& Foam::operator<< <Foam::error>(Foam::Ostream&, Foam::errorManip<Foam::error>) in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/channelFoam"
#3 Foam::lduInterface const& Foam::refCast<Foam::lduInterface const, Foam::fvPatch const>(Foam::fvPatch const&) in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libfiniteVolume.so"
#4 Foam::cyclicFvPatchField<double>::cyclicFvPatchFie ld(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libfiniteVolume.so"
#5 Foam::fvPatchField<double>::adddictionaryConstruct orToTable<Foam::cyclicFvPatchField<double> >::New(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libfiniteVolume.so"
#6 Foam::fvPatchField<double>::New(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/channelFoam"
#7 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::GeometricB oundaryField(Foam::fvBoundaryMesh const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/channelFoam"
#8 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readField(Foam::dictionary const&) in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/channelFoam"
#9 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readField(Foam::Istream&) in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/channelFoam"
#10 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&) in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/channelFoam"
#11 main in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/channelFoam"
#12 __libc_start_main in "/lib64/libc.so.6"
#13 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/channelFoam"
Aborted


Can someone please offer me any help regarding this?

My aim is to run ICEM generated grids with OF. Can someone doing the same help me with the steps?


Thanking You,
Tarak

kalle October 12, 2010 02:39

Hi Tarak,

please check that the settings in your fields (U, p, ...) in the start directory (0) corresponds to the polyMesh/boundary file. If a boundary patch is given one type in the bouindary file, and a non-corresponding boundary condition is applied for the fields you get this type of errors...

Possibly during mesh conversion, the boundary file is written in a way which does not suit your problem...

Good luck!

K

c_blake October 12, 2010 16:53

What format did you output from ICEM to then convert to OpenFOAM format?
I output an ICEM mesh as Fluent v6 format then used fluentMeshToFoam utility to convert the mesh and it worked fine. Though I had to change some of the patch types in the boundary file within the polymesh directory as kalle also stated.

Tarak October 12, 2010 22:11

@Karl and Chris
 
@karl,

I did change the BC names is the files you mentioned, and after doing that only this error comes.

@Chris

I set the output solver as Fluent V6 in ICEM and converted it using fluentMeshToFoam utility. Can you please tell me what modifications did you made in the files other than changing the BC names. Also please let me know which BC"s to be kept as patch.

Thanking You,
Tarak

c_blake October 13, 2010 15:43

It depends on what your boundary conditions are in the U/p 0 directory are. They just have to be the same. For example you couuldn't have a boundary as a patch in the boundary file (polymesh directory) and then set up boundary conditions for a wall in the 0 directory. have a look in the user guide, which sets out which patch types corrospond with which boundary conditions. (at least I'm fairly sure it does, haven't looked at it in a while)

My suggestion would be to start off with a very simple mesh with simple boundary conditions, for example just modelling flow through a rectangular domain with free slip walls around the sides. This is very simple to do, then once you get this working, you can slowly modify it to what you want to model. This is a good way to find out where the problems occur, rather than just jumping right into the final problem.

Tarak November 5, 2010 16:01

Hii,

After checking consistencies for the BC's and mesh, I am still getting this error message. can you please give me a hint about any way to solve this?



Create time

Create mesh for time = 0


Reading transportProperties

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting LES turbulence model mdesDynOneEqEddy
#0 Foam::error::printStack(Foam::Ostream&) in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 __restore_rt at sigaction.c:0
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libincompressibleLESModels.so"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libincompressibleLESModels.so"
#6 Foam::incompressible::LESModels::mdesDynOneEqEddy: :F1(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/tnandi/OpenFOAM/tnandi-1.7.1/lib/linux64GccDPOpt/libmymdeslesIncompressibleLESModels.so"
#7 Foam::incompressible::LESModels::mdesDynOneEqEddy: :mdesDynOneEqEddy(Foam::GeometricField<Foam::Vecto r<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/home/tnandi/OpenFOAM/tnandi-1.7.1/lib/linux64GccDPOpt/libmymdeslesIncompressibleLESModels.so"
#8 Foam::incompressible::LESModel::adddictionaryConst ructorToTable<Foam::incompressible::LESModels::mde sDynOneEqEddy>::New(Foam::GeometricField<Foam::Vec tor<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/home/tnandi/OpenFOAM/tnandi-1.7.1/lib/linux64GccDPOpt/libmymdeslesIncompressibleLESModels.so"
#9 Foam::incompressible::LESModel::New(Foam::Geometri cField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libincompressibleLESModels.so"
#10 main in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/channelFoam"
#11 __libc_start_main in "/lib64/libc.so.6"
#12 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/storage/software/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/channelFoam"
Floating point exception

Thanks,
Tarak

desmoge September 9, 2011 17:51

Hi there I'm new here and I have almost the same problem Tarak has got. I've been trying to export in OF a mesh which has been edited by me just for the node distribution using ICEM. The previous setting works just fine but once I increase the number of nodes and the their distribution (spline), I get this error by using 'fluentMeshtoFoam':


Dimension of grid: 2
Number of points: 1052150
Reading points
Number of cells: 1048856
Other readCellGroupData: 20 1 c13f8 1 3
Reading uniform cells
Other readCellGroupData: 21 c13f9 100118 1 3
Reading uniform cells
number of faces: 2101006
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Read zone1:32 name:FLUID patchTypeID:fluid
Reading zone data
Read zone1:33 name:DRZONE patchTypeID:fluid
Reading zone data
Read zone1:34 name:int_FLUID patchTypeID:interior
Reading zone data
Read zone1:35 name:int_DRZONE patchTypeID:interior
Reading zone data
Read zone1:36 name:CRVS_AF patchTypeID:wall
Reading zone data
Read zone1:37 name:P_OUT patchTypeID:pressure-outlet
Reading zone data
Read zone1:38 name:V_IN patchTypeID:velocity-inlet
Reading zone data
Read zone1:39 name:HUMP01 patchTypeID:wall
Reading zone data
Read zone1:40 name:HUMP02 patchTypeID:wall
Reading zone data
Read zone1:41 name:HUMP03 patchTypeID:wall
Reading zone data
Read zone1:42 name:HUMP04 patchTypeID:wall
Reading zone data
Read zone1:43 name:HUMP05 patchTypeID:wall
Reading zone data
Read zone1:44 name:HUMP06 patchTypeID:wall
Reading zone data
Read zone1:45 name:HUMP07 patchTypeID:wall
Reading zone data
Read zone1:46 name:HUMP08 patchTypeID:wall
Reading zone data
Read zone1:47 name:HUMP09 patchTypeID:wall
Reading zone data
Read zone1:48 name:HUMP10 patchTypeID:wall
Reading zone data
Read zone1:49 name:HUMP11 patchTypeID:wall
Reading zone data
Read zone1:50 name:HUMP12 patchTypeID:wall
Reading zone data
Read zone1:51 name:HUMP13 patchTypeID:wall
Reading zone data
Read zone1:52 name:HUMP14 patchTypeID:wall
Reading zone data
Read zone1:53 name:HUMP15 patchTypeID:wall
Reading zone data
Read zone1:54 name:HUMP17 patchTypeID:wall
Reading zone data
Read zone1:55 name:DRZONE patchTypeID:fluid
Reading zone data
Read zone1:56 name:HUMP16 patchTypeID:wall
Reading zone data


FINISHED LEXING


dimension of grid: 2
Grid is 2-D. Extruding in z-direction by: 0.512242
Creating shapes for 2-D cells
Building patch-less mesh...--> FOAM Warning :
From function polyMesh::polyMesh(... construct from shapes...)
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 614
Found 2104300 undefined faces in mesh; adding to default patch.
done.

Building boundary and internal patches.
Creating patch 0 for zone: 34 start: 1 end: 1580198 type: interior name: int_FLUID
Creating patch 1 for zone: 35 start: 1580199 end: 2093730 type: interior name: int_DRZONE
Creating patch 2 for zone: 36 start: 2093731 end: 2095833 type: wall name: CRVS_AF
Creating patch 3 for zone: 37 start: 2095834 end: 2096521 type: pressure-outlet name: P_OUT
Creating patch 4 for zone: 38 start: 2096522 end: 2099570 type: velocity-inlet name: V_IN
Creating patch 5 for zone: 39 start: 2099571 end: 2099614 type: wall name: HUMP01
Creating patch 6 for zone: 40 start: 2099615 end: 2099658 type: wall name: HUMP02
Creating patch 7 for zone: 41 start: 2099659 end: 2099702 type: wall name: HUMP03
Creating patch 8 for zone: 42 start: 2099703 end: 2099746 type: wall name: HUMP04
Creating patch 9 for zone: 43 start: 2099747 end: 2099790 type: wall name: HUMP05
Creating patch 10 for zone: 44 start: 2099791 end: 2099834 type: wall name: HUMP06
Creating patch 11 for zone: 45 start: 2099835 end: 2099878 type: wall name: HUMP07
Creating patch 12 for zone: 46 start: 2099879 end: 2099922 type: wall name: HUMP08
Creating patch 13 for zone: 47 start: 2099923 end: 2099966 type: wall name: HUMP09
Creating patch 14 for zone: 48 start: 2099967 end: 2100010 type: wall name: HUMP10
Creating patch 15 for zone: 49 start: 2100011 end: 2100054 type: wall name: HUMP11
Creating patch 16 for zone: 50 start: 2100055 end: 2100098 type: wall name: HUMP12
Creating patch 17 for zone: 51 start: 2100099 end: 2100142 type: wall name: HUMP13
Creating patch 18 for zone: 52 start: 2100143 end: 2100186 type: wall name: HUMP14
Creating patch 19 for zone: 53 start: 2100187 end: 2100230 type: wall name: HUMP15
Creating patch 20 for zone: 54 start: 2100231 end: 2100274 type: wall name: HUMP17
Creating patch 21 for zone: 55 start: 2100275 end: 2100962 type: fluid name: DRZONE


--> FOAM FATAL ERROR:
fluent patch type fluid not recognised.

From function fluentToFoam::main(int argc, char *argv[])
in file fluentMeshToFoam.L at line 1344.

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&) in "/home/gennaro/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/gennaro/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2
in "/home/gennaro/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/bin/fluentMeshToFoam"
#3 __libc_start_main in "/lib/libc.so.6"
#4
in "/home/gennaro/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/bin/fluentMeshToFoam"
Aborted


Can somebody help me with this? I seriously don't get the problem. Apparently it does not recognise patches but how can this happen just changing nodes number and distribution? (Output solver used in ICEM is always Fluent v6). Thanks guys


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