[liguifan]

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

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

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

/*--------------------------------*- 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

/*--------------------------------*- 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

/*--------------------------------*- 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]

Originally Posted by liguifan

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]

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]

Originally Posted by liguifan

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"