CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Mesh Utilities (http://www.cfd-online.com/Forums/openfoam-meshing-utilities/)
-   -   Error of using MakeAxialMesh ultility (http://www.cfd-online.com/Forums/openfoam-meshing-utilities/116179-error-using-makeaxialmesh-ultility.html)

liguifan April 14, 2013 21:23

Error of using MakeAxialMesh ultility
 
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::printStack(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

gschaider April 15, 2013 17:14

Quote:

Originally Posted by liguifan (Post 420541)
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

liguifan April 15, 2013 17:23

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

gschaider April 16, 2013 05:39

Quote:

Originally Posted by liguifan (Post 420776)
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"


All times are GMT -4. The time now is 00:57.