CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Pre-Processing

How to handle it

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 20, 2008, 02:51
Default hi, i am doing a simulation a
  #1
Senior Member
 
weihong yao
Join Date: Mar 2009
Posts: 117
Rep Power: 8
ivanyao is on a distinguished road
hi,
i am doing a simulation about wind pressure on a 6m cube,i put the 6m*6m*6m cube into 192m*120m*60m box.when i extrudeMesh,it show follow message;
FOAM FATAL ERROR : Trying to specify a boundary face 4(0 3 7 4) on the face on cell 14 which is either an internal face or already belongs to some other patch. This is face 0 of patch 2 named front.#0 Foam::error::printStack(Foam:stream&)
#1 Foam::error::abort()
#2 Foam::polyMesh::polyMesh(Foam::IOobject const&, Foam::Field<foam::vector<double> > const&, Foam::List<foam::cellshape> const&, Foam::List<foam::list<foam::face> > const&, Foam::List<foam::word> const&, Foam::List<foam::word> const&, Foam::word const&, Foam::List<foam::word> const&)
#3 Foam::blockMesh::createTopology(Foam::IOdictionary &)
#4 Foam::blockMesh::blockMesh(Foam::IOdictionary&)
#5 main
#6 __libc_start_main
#7 __gxx_personality_v0 at /usr/src/packages/BUILD/glibc-2.3/csu/../sysdeps/i386/elf/start.S:122


From function polyMesh::polyMesh
(
const IOobject& io,
const pointField& points,
const cellShapeList& cellsAsShapes,
const faceListList& boundaryFaces,
const wordList& boundaryPatchTypes,
const wordList& boundaryPatchNames,
const word& defaultBoundaryPatchType
)
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 482.
FOAM aborting

my blockMesh file are:
convertToMeters 1;

vertices
(
(-3 -3 0)
(3 -3 0)
(3 3 0)
(-3 3 0)
(-3 -3 6)
(3 -3 6)
(3 3 6)
(-3 3 6)
(-39 -3 0)
(-39 3 0)
(-39 -3 6)
(-39 3 6)
(153 -3 0)
(153 3 0)
(153 -3 6)
(153 3 6)
(-3 -60 0)
(3 -60 0)
(-3 -60 6)
(3 -60 6)
(-3 60 0)
(3 60 0)
(-3 60 6)
(3 60 6)
(-39 -60 0)
(153 -60 0)
(153 60 0)
(-39 60 0)
(-39 -60 6)
(153 -60 6)
(153 60 6)
(-39 60 6)
(-39 -60 60)
(153 -60 60)
(153 60 60)
(-39 60 60)
(-39 -3 60)
(-39 3 60)
(-3 -60 60)
(3 -60 60)
(153 -3 60)
(153 3 60)
(-3 60 60)
(3 60 60)
(-3 -3 60)
(3 -3 60)
(3 3 60)
(-3 3 60)
);

blocks
(
hex (24 16 0 8 28 18 4 10) (36 57 6) simpleGrading (1 1 1)
hex (28 18 4 10 32 38 44 36) (36 57 54) simpleGrading (1 1 1)
hex (16 17 1 0 18 19 5 4) (6 57 6) simpleGrading (1 1 1)
hex (18 19 5 4 38 39 45 44) (6 57 54) simpleGrading (1 1 1)
hex (17 25 12 1 19 29 14 5) (150 57 6) simpleGrading (1 1 1)
hex (19 29 14 5 39 33 40 45) (150 57 54) simpleGrading (1 1 1)
hex (1 12 13 2 5 14 15 6) (150 6 6) simpleGrading (1 1 1)
hex (5 14 15 6 45 40 41 46) (150 6 54) simpleGrading (1 1 1)
hex (2 13 26 21 6 15 30 23) (150 57 6) simpleGrading (1 1 1)
hex (6 15 30 23 46 41 34 43) (150 57 54) simpleGrading (1 1 1)
hex (3 2 21 20 7 6 23 22) (6 57 6) simpleGrading (1 1 1)
hex (7 6 23 22 47 46 43 42) (6 57 54) simpleGrading (1 1 1)
hex (9 3 20 27 11 7 22 31) (36 57 6) simpleGrading (1 1 1)
hex (11 7 22 31 37 47 42 35) (36 57 54) simpleGrading (1 1 1)
hex (8 0 3 9 10 4 7 11) (36 6 6) simpleGrading (1 1 1)
hex (10 4 7 11 36 44 47 37) (36 6 54) simpleGrading (1 1 1)
hex (4 5 6 7 44 45 46 47) (6 6 54) simpleGrading (1 1 1)
hex (0 1 2 3 4 5 6 7) (6 6 6) edgeGrading (1 1 1 1 1 1 1 1 1 1 1 1)
);

edges
(
);

patches
(
patch inlet
(
(27 9 11 31)
(37 11 31 35)
(9 8 10 11)
(11 10 36 37)
(8 24 28 10)
(10 28 32 36)
)
patch outlet
(
(25 12 14 29)
(29 14 40 33)
(12 13 15 14)
(14 15 41 40)
(13 26 30 15)
(15 30 34 41)
)
wall front
(
(0 3 7 4)
)
wall <back>
(
(1 2 6 5)
)
wall <left>
(
(0 1 5 4)
)
wall <right>
(
(2 3 7 6)
)
wall <top>
(
(4 5 6 7)
)
empty <wall>
(
(24 16 18 28)
(28 18 38 32)
(16 17 19 18)
(18 19 39 38)
(17 25 29 19)
(19 29 33 39)
(26 21 23 30)
(30 23 43 34)
(21 20 22 23)
(23 22 42 43)
(20 27 31 22)
(22 31 35 42)
(32 38 44 36)
(44 36 37 47)
(47 37 35 42)
(38 39 45 44)
(44 45 46 47)
(47 46 43 42)
(39 33 40 45)
(40 45 46 41)
(41 46 43 34)
(16 24 8 0)
(0 8 9 3)
(3 9 27 20)
(17 16 0 1)
(2 3 20 21)
(25 17 1 12)
(12 1 2 13)
(13 2 21 26)
)
wall <bottom>
(
(0 1 2 3)
)
);

mergePatchPairs
(
);
would anyone help me?
ivanyao is offline   Reply With Quote

Old   August 20, 2008, 04:03
Default Run blockMesh with -blockTopol
  #2
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
Run blockMesh with -blockTopology option. This will write a .obj file which can be visualised with e.g. javaview or (with objToVTK) with paraview. Both have the option of displaying vertex numbers (with selecting and selection inspector in paraview). Very useful to debug blockMeshes.
mattijs is offline   Reply With Quote

Old   August 20, 2008, 05:39
Default I have written an application
  #3
Senior Member
 
Gavin Tabor
Join Date: Mar 2009
Posts: 181
Rep Power: 8
grtabor is on a distinguished road
I have written an application called simpleViewer based on blockMesh, which visualises the block structure of a blockMeshDict file for you. Uses FOX/OpenGL. It works for OF1.4.1; I'm guessing making it work for 1.5 would not be too hard. Email me at g.r.tabor@ex.ac.uk if you would like to try with this.

Gavin
grtabor is offline   Reply With Quote

Old   August 20, 2008, 22:36
Default hi Mattijs, it can't work.whe
  #4
Senior Member
 
weihong yao
Join Date: Mar 2009
Posts: 117
Rep Power: 8
ivanyao is on a distinguished road
hi Mattijs,
it can't work.when i type blockMesh <root> <case> -blockTopology.it show:
FOAM FATAL ERROR : Trying to specify a boundary face 4(0 3 7 4) on the face on cell 14 which is either an internal face or already belongs to some other patch. This is face 0 of patch 2 named front.#0 Foam::error::printStack(Foam:stream&)
#1 Foam::error::abort()
#2 Foam::polyMesh::polyMesh(Foam::IOobject const&, Foam::Field<foam::vector<double> > const&, Foam::List<foam::cellshape> const&, Foam::List<foam::list<foam::face> > const&, Foam::List<foam::word> const&, Foam::List<foam::word> const&, Foam::word const&, Foam::List<foam::word> const&)
#3 Foam::blockMesh::createTopology(Foam::IOdictionary &)
#4 Foam::blockMesh::blockMesh(Foam::IOdictionary&)
#5 main
#6 __libc_start_main
#7 __gxx_personality_v0 at /usr/src/packages/BUILD/glibc-2.3/csu/../sysdeps/i386/elf/start.S:122


From function polyMesh::polyMesh
(
const IOobject& io,
const pointField& points,
const cellShapeList& cellsAsShapes,
const faceListList& boundaryFaces,
const wordList& boundaryPatchTypes,
const wordList& boundaryPatchNames,
const word& defaultBoundaryPatchType
)
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 482.
FOAM aborting
so how do i fix it?
ivanyao is offline   Reply With Quote

Old   August 21, 2008, 16:21
Default I just went through your block
  #5
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
I just went through your blockMeshDict. Set all block sizes to (1 1 1) so there are no matching errors. Then commented out all patch specifications (front, back, left, right, top) that give errors. Now with forementioned -blockTopology start putting patches and block sizes back in.
mattijs is offline   Reply With Quote

Old   August 22, 2008, 21:35
Default hi, thank you!but now what ca
  #6
Senior Member
 
weihong yao
Join Date: Mar 2009
Posts: 117
Rep Power: 8
ivanyao is on a distinguished road
hi,
thank you!but now what can I do?
ivanyao is offline   Reply With Quote

Old   August 24, 2008, 21:27
Default hi Mattijs, i am poor in Engl
  #7
Senior Member
 
weihong yao
Join Date: Mar 2009
Posts: 117
Rep Power: 8
ivanyao is on a distinguished road
hi Mattijs,
i am poor in English.so i don't know what to do.
would you tell me the step in detail.please...
ivanyao 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
How to handle the gcc412 ivanyao OpenFOAM Bugs 1 March 28, 2015 23:35
Is UDF the only way to handle particle melting? marq FLUENT 3 February 23, 2008 02:48
Can CFX handle local hex refinement? JollyP CFX 2 October 31, 2006 10:00
Libstdcso6 cannot handle TLS data barryw OpenFOAM Installation 2 September 7, 2006 02:02
Can ALE handle air bubble breakups? Pei-Ying Hsieh Main CFD Forum 1 May 29, 2003 10:54


All times are GMT -4. The time now is 12:44.