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

Problem with Mesh Volume

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 29, 2012, 11:01
Default Problem with Mesh Volume
  #1
Member
 
Join Date: Aug 2011
Posts: 33
Rep Power: 6
seboxx is on a distinguished road
Dear Foamers,

I am struggling with a Problem for several weeks already and starting to get really desperate. I have posted it in the main forum already, but I figured posting a new thread here as I have narrowed the problem down to blockMesh, I guess so at least.

The Problem is I want to simulate a thin layer with a bump, I have modelled the bump with a blockMeshDict.
What is wrong here is that the volume is not recognized correctly. You can see in checkMesh the total volume of 1256.62. But this cannot be correct, since a flat disc with height one and radius 20 would have exactly that volume (Pi*r^2)

But the graphical output in paraview is correct. Whats wrong here?

Now when I do checkMesh, this is the output
Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.0.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.0.0-a317a4e7cd55
Exec   : checkMesh
Date   : Jan 29 2012
Time   : 15:48:06
Host   : ernie
PID    : 7269
Case   : /home/boie/OpenFOAM/OpenFOAM-reloaded/boie-2.0.0/simulations/tests/justbump
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
    points:           80802
    internal points:  0
    faces:            160400
    internal faces:   79600
    cells:            40000
    boundary patches: 3
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     40000
    prisms:        0
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    0
    polyhedra:     0

Checking topology...
    Boundary definition OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces ...
    Patch               Faces    Points   Surface topology                  
    frontWall           40000    40401    ok (non-closed singly connected)  
    backWall            40000    40401    ok (non-closed singly connected)  
    circular            800      1600     ok (non-closed singly connected)  

Checking geometry...
    Overall domain bounding box (-20 -20 -0.0312941) (20 20 21)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (8.25305e-21 7.15265e-20 -3.39808e-19) OK.
    Max cell openness = 2.24782e-16 OK.
    Max aspect ratio = 321.104 OK.
    Minumum face area = 0.000702619. Maximum face area = 0.283752.  Face area magnitudes OK.
    Min volume = 0.000692074. Max volume = 0.0803204.  Total volume = 1256.62.  Cell volumes OK.
    Mesh non-orthogonality Max: 87.5911 average: 38.3172
   *Number of severely non-orthogonal faces: 624.
    Non-orthogonality check OK.
  <<Writing 624 non-orthogonal faces to set nonOrthoFaces
    Face pyramids OK.
    Max skewness = 2.81207 OK.

Mesh OK.

End
Attached Images
File Type: jpg bump.jpg (11.4 KB, 14 views)
Attached Files
File Type: txt blockMeshDict.txt (4.0 KB, 5 views)
seboxx is offline   Reply With Quote

Old   January 30, 2012, 04:17
Default
  #2
Senior Member
 
BastiL
Join Date: Mar 2009
Posts: 471
Rep Power: 11
bastil is on a distinguished road
I don't see you problem. Your formula is right and I get:

3.14*20^2=3.14*400=1256.

So for me everything seems to be fine.
bastil is offline   Reply With Quote

Old   January 30, 2012, 04:21
Default
  #3
Member
 
Join Date: Aug 2011
Posts: 33
Rep Power: 6
seboxx is on a distinguished road
Hi bastil,

well the problem is the volume must be larger than this.
See the attached picture, its a bump and not a disc anymore. The Problem is similiar to taking a cube and a cube with a large dent in it but with the same length of the edges. The volume should be different
seboxx is offline   Reply With Quote

Old   January 30, 2012, 09:19
Default
  #4
Member
 
Join Date: Aug 2011
Posts: 33
Rep Power: 6
seboxx is on a distinguished road
To illustrate my problem I have made 2 more screenshots. In the first picture you look from on top on the geometry.
In the second picture is the same geometry with a different angle of observation. checkMesh reports no errors and a total volume of:
Quote:
Total volume = 400.
But that is the volume like if i project the geometry into the xy-plane.
On the second picture you see the real geometry and clearly the volume must be larger...

I really have no idea what is going wrong here
Attached Images
File Type: jpg projection.jpg (11.0 KB, 12 views)
File Type: jpg slanted.jpg (12.5 KB, 15 views)
seboxx 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
Anyone using GID mesher? problem with .igs file and volume mesh afo3 Main CFD Forum 1 November 24, 2010 14:31
Ansys 12 Workbench problem when generating volume mesh Sveink ANSYS 5 January 6, 2010 04:17
problem on volume mesh wan CFX 0 January 29, 2008 03:41
Volume mesh problem in STARDESIGN Jaedal CD-adapco 2 October 17, 2007 09:23
Help! i need help Mesh volume problem Sergi CFX 1 February 1, 2006 17:29


All times are GMT -4. The time now is 03:30.