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
|