CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Native Meshers: blockMesh

Problems meshing wedge type blocks

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 27, 2009, 20:05
Default Problems meshing wedge type blocks
  #1
New Member
 
Alan Harrland
Join Date: Mar 2009
Posts: 21
Rep Power: 8
Alan is on a distinguished road
Hi all,

I am trying to mesh a wedge shape in blockMesh, and I am having trouble getting it to work properly. I am getting the warning:

Quote:
Creating block mesh topology
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 134
zero or negative pyramid volume: -0.0075 for face 1
so I assume I have the ordering of some of points wrong, but I can't for the life of me figure out which ones are wrong. My blockMesh file is:

Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.5                                   |
|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 1;

vertices        
(
   (0 0 -0.05)
   (1 0 -0.05)
   (1.5 0 -0.05)
   (3 0 -0.05)
   (1 0.85 -0.05)
   (1.5 0.85 -0.05)
   (1 1 -0.05)
   (1 1.15 -0.05)
   (1.5 1.15 -0.05)
   (0 2 -0.05)
   (1 2 -0.05)
   (1.5 2 -0.05)
   (3 2 -0.05)
   (3 0.85 -0.05)
   (3 1.15 -0.05)
   (0 0 0.05)
   (1 0 0.05)
   (1.5 0 0.05)
   (3 0 0.05)
   (1 0.85 0.05)
   (1.5 0.85 0.05)
   (1 1 0.05)
   (1 1.15 0.05)
   (1.5 1.15 0.05)
   (0 2 0.05)
   (1 2 0.05)
   (1.5 2 0.05)
   (3 2 0.05)
   (3 0.85 0.05)
   (3 1.15 0.05)
);

blocks          
(
    hex (0 1 10 9 15 16 25 24)   (100 200 1) simpleGrading (1 1 1)
    hex (1 2 5 4 16 17 20 19)    (50 85 1) simpleGrading (1 1 1)
    hex (4 5 6 6 19 20 21 21)    (50 15 1) simpleGrading (1 1 1)
    hex (8 7 6 6 23 22 21 21)    (50 15 1) simpleGrading (1 1 1)
    hex (7 8 11 10 22 23 26 25)  (50 85 1) simpleGrading (1 1 1)
    hex (2 3 14 5 17 18 29 20)   (150 85 1) simpleGrading (1 1 1)
    hex (5 14 13 8 20 29 28 23)  (150 30 1) simpleGrading (1 1 1)
    hex (8 13 12 11 23 28 27 26) (150 85 1) simpleGrading (1 1 1)
);

edges           
(
);

patches         
(
    patch inlet 
    (
        (0 15 24 9)
    )
    patch outlet 
    (
        (3 18 29 14)
        (14 29 28 13)
        (13 28 27 12)
    )
    symmetryPlane bottom 
    (
        (0 1 16 15)
        (1 2 17 16)
        (2 3 18 17)
    )
    symmetryPlane top 
    (
        (9 10 25 24)
        (10 11 26 25)
    (11 12 27 26)
    )
    patch obstacle
    (
        (6 8 23 21)
        (6 21 20 5)
    (5 20 23 8)
    )
);

mergePatchPairs
(
);

// ************************************************************************* //
I assume the problem occurs around blocks 2 and 3 as they are the wedge shapes. If I play around with the ordering it defines the wedges as inside out (getting negative volumes) so I can't see why I have done it wrong.

Can anyone point me in the right direction here? Thanks very much in advance for your help.

Alan.
Alan 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
Compile problem ivanyao OpenFOAM Running, Solving & CFD 1 October 12, 2012 09:31
reconstructParMesh not working with an axisymetric case francesco OpenFOAM Bugs 4 May 8, 2009 05:49
RasInterFoam STRANGE RESULTS AT BOUNDARY kumar2 OpenFOAM Running, Solving & CFD 8 March 24, 2008 19:38
meshing problems with CFX 10 Cedric PERE CFX 1 March 15, 2007 19:09
ICEM CFD Meshing Problems: Hexa Dominant Evan CFX 6 July 3, 2006 12:37


All times are GMT -4. The time now is 07:51.