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

Error of using MakeAxialMesh ultility

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 14, 2013, 21:23
Default Error of using MakeAxialMesh ultility
  #1
Member
 
Guifan Li
Join Date: Apr 2011
Location: New York City, U.S.
Posts: 96
Rep Power: 6
liguifan is on a distinguished road
Hello everyone,

I am trying to generate a very simple 2D wedge which is part of a cylinder water tank. I started from the begging by follwoing the example here
http://openfoamwiki.net/index.php/Contrib_MakeAxialMesh
However, I met the following errors
Quote:
liguifan@liguifan-laptop:~/Desktop/damBreak2$ makeAxialMesh -axis axi -wedge frontAndBackPlanes
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.1-221db2718bbb
Exec : makeAxialMesh -axis axi -wedge frontAndBackPlanes
Date : Apr 14 2013
Time : 21:14:05
Host : "liguifan-laptop"
PID : 2821
Case : /home/liguifan/Desktop/damBreak2
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 mesh for time = 0

Using old mode. Dictionary not used
Plane of the grid: (0 0 1) (0.02 0.04 0.005)

The rotation-axis: ((0 0 0.005) (0.04 0 0.005))

Creating wedge with an opening angle of 5 degrees

Projecting nodes
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Uninterpreted:
#3
in "/home/liguifan/OpenFOAM/liguifan-2.1.1/platforms/linuxGccDPOpt/bin/makeAxialMesh"
#4
in "/home/liguifan/OpenFOAM/liguifan-2.1.1/platforms/linuxGccDPOpt/bin/makeAxialMesh"
#5 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#6
in "/home/liguifan/OpenFOAM/liguifan-2.1.1/platforms/linuxGccDPOpt/bin/makeAxialMesh"
Floating point exception
Can anyone help me with that?

Here are my files under folder 0 and my Gmsh file
Gmsh
Quote:
Point(1) = {0, 0, 0};
Point(2) = {0.04, 0, 0};
Line(1) = {1, 2};
Extrude {0, 0.08, 0} {
Line{1};Layers{10};Recombine;
}

Extrude {0, 0, 0.01} {
Surface{5};Layers{1};Recombine;
}
Transfinite Surface "*";
Recombine Surface "*";
Transfinite Volume "*";

Physical Surface("frontAndBackPlanes") = {27,5};
Physical Surface("outerWall") = {22};
Physical Surface("axi") = {14};
Physical Surface("atmosphere") = {26};
Physical Surface("lowerWall") = {18};
Physical Volume("internal") = {1};

alpha1.org
Quote:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object alpha;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 0 0 0 0];

internalField uniform 0;

boundaryField
{

outerWall
{
type zeroGradient;
}

lowerWall
{
type zeroGradient;
}
frontAndBackPlanes
{
type wedge;
}

atmosphere
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}

axi
{
type wedge;
}
}

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

dimensions [1 -1 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
outerWall
{
type buoyantPressure;
value uniform 0;
}
frontAndBackPlanes
{
type wedge;
//value uniform 0;
}

lowerWall
{
type buoyantPressure;
value uniform 0;
}


atmosphere
{
type totalPressure;
p0 uniform 0;
U U;
phi phi;
rho rho;
psi none;
gamma 1;
value uniform 0;
}

axi
{
type patch;
}
}

// ************************************************** *********************** //

constant/Mesh/boundary
Quote:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/polyMesh";
object boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

5
(
frontAndBackPlanes
{
type wedge;
nFaces 1280;
startFace 6816;
}
axi
{
type wedge;
nFaces 64;
startFace 8096;
}
lowerWall
{
type patch;
nFaces 160;
startFace 8160;
}
outerWall
{
type patch;
nFaces 64;
startFace 8320;
}
atmosphere
{
type patch;
nFaces 160;
startFace 8384;
}
)

// ************************************************** *********************** //
Thanks in advance
liguifan is offline   Reply With Quote

Old   April 15, 2013, 17:14
Default
  #2
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,912
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by liguifan View Post
Hello everyone,

I am trying to generate a very simple 2D wedge which is part of a cylinder water tank. I started from the begging by follwoing the example here
http://openfoamwiki.net/index.php/Contrib_MakeAxialMesh
However, I met the following errors


Can anyone help me with that?

Here are my files under folder 0 and my Gmsh file
Gmsh


alpha1.org


p_rgh


constant/Mesh/boundary


Thanks in advance
Not quite sure. I'm not too familiar with gmsh, but it might be that the problem is that in order to work correctly the axis should be in the z=0-plane (and the mesh centered around z=0 ... z should go from -0.005 to 0.005). But I'm not completely sure. Haven't touched the source for MAM for quite some time
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   April 15, 2013, 17:23
Default
  #3
Member
 
Guifan Li
Join Date: Apr 2011
Location: New York City, U.S.
Posts: 96
Rep Power: 6
liguifan is on a distinguished road
Thanks for the reply.

There is one thing I am confused.
In tutorial of NozzelFlow2D, its BC for axis is empty. However in my case after I execuated makeAxialMesh command. The BC for my axis became symmetryPlane type. Do I have to change this into empty in the constant/polyMesh/boundary file and all the other BC's file in 0/ folder ?

Kind regards,
Guifan
liguifan is offline   Reply With Quote

Old   April 16, 2013, 05:39
Default
  #4
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,912
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by liguifan View Post
Thanks for the reply.

There is one thing I am confused.
In tutorial of NozzelFlow2D, its BC for axis is empty. However in my case after I execuated makeAxialMesh command. The BC for my axis became symmetryPlane type. Do I have to change this into empty in the constant/polyMesh/boundary file and all the other BC's file in 0/ folder ?

Kind regards,
Guifan
See what works best for you. If I remember correctly the reason why MAM sets it to symmetry is that this works better if the area of the axis-patch is not 0 and it is therefor not "collapsed away"
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider 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 make a pipe using makeAxialMesh tomloh OpenFOAM Mesh Utilities 15 September 3, 2013 05:52
makeAxialMesh and axisymmetric flow yassepid OpenFOAM Mesh Utilities 4 April 24, 2013 20:01
makeAxialMesh - installation issues simulator4711 OpenFOAM Mesh Utilities 3 March 21, 2012 08:59
collapseEdges (after makeAxialMesh) not worked mechy OpenFOAM 2 November 17, 2011 16:09
makeAxialMesh > collapseEdges!! billynoe OpenFOAM 17 February 21, 2011 05:29


All times are GMT -4. The time now is 20:28.