CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Meshing & Mesh Conversion (https://www.cfd-online.com/Forums/openfoam-meshing/)
-   -   [mesh manipulation] Converting a 2Dmesh to axisymmetric (https://www.cfd-online.com/Forums/openfoam-meshing/61473-converting-2dmesh-axisymmetric.html)

gschaider April 20, 2005 15:57

Converting a 2Dmesh to axisymmetric
 
I just wrote a utility that takes a regular 2D-mesh (a mesh that is one cell thick) and converts it into a axi-symmetric mesh (as described in the UserGuide 6.2.2) by collapsing one (user-specified) boundary patch and adjusting the cordinates of the vertices.

If it is of general interest, I'd be happy to share it.

hsieh April 20, 2005 20:15

Hi, Bernhard, I am interest
 
Hi, Bernhard,

I am interested in this utility. Could you please email to me? Is the input using true 2D mesh (like shell element) or it takes a 3D mesh but, with only one cell thick in the Z direction? Thanks!

Pei
Peiying2003@yahoo.com

nils April 21, 2005 03:05

Hi, I am also interested! C
 
Hi,

I am also interested! Could you please also email it to me? Thanks a lot!


Nils :-)

gschaider April 21, 2005 03:51

Yes, the utilities uses a Open
 
Yes, the utilities uses a OpenFOAM-2D-mesh: a 3D-mesh that is one cell thick. Don't know if OF can handle true 2D-grids (I don't think so)

groens January 17, 2006 08:26

Hello Bernhard, I am starting
 
Hello Bernhard,
I am starting to do some work with FOAM.
I would be interested in your tool, too.

Thanks in Advance.
Thomas

gschaider January 17, 2006 10:01

Hello Thomas! The utility c
 
Hello Thomas!

The utility can be found at:

http://openfoamwiki.net/index.php/Contrib_MakeAxi alMesh

Let me know if it works for you (Judging from the date it was developed on 1.1, maybe something has to be changed in Make/options to get it to compile on 1.2)

engrstudentgmailcom February 1, 2008 11:08

Hey All, Im a mac user, for
 
Hey All,

Im a mac user, formerly a windows user, and I am not a linux programmer.

Im trying to work through a basic cfd tutorial for fluent, except in OpenFOAM.
http://courses.cit.cornell.edu/fluent/pipe1/index.htm

I made the mesh in gambit, and exported to fluent.
I found the utility to import the fluent mesh into OpenFOAM. Now I discover that one problem is that the new mesh isn't axisymmetric. I need to find a way to do that. (Im using CAELinux with its installation of OpenFOAM.)

I dont know how to compile with makefiles. A majority of my programming experience is MatLab.

How do i make the makeAxialMesh go from a .tgz to a working program in CAE Linux - Foam?

Thanks,
Mike Munroe

gschaider February 4, 2008 08:36

Hi Michael! Just untar the
 
Hi Michael!

Just untar the sources (the tgz is nothing but a zip-file for men). Propable you can do that from the GUI of your Linux-distro. Then go to the resulting directory and do as described in section 3.2 of the OF-Userguide (Compilation of applications). No programming should be required.

Bernhard

bonzodeb April 19, 2008 19:08

Hi all, I have a problem whe
 
Hi all,
I have a problem when I try to compile MakeAxialMesh .
I have:
1) In $HOME untar the file MakeAxialMesh1.4.1.tgz
2) cd ./makeAxialMesh
3)wmake

This is the output

SOURCE=makeAxialMesh.C ; g++ -m32 -Dlinux -DDP -Wall -Wno-strict-aliasing -Wextra -Wno-unused-parameter -Wold-style-cast -O3 -DNoRepository -ftemplate-depth-40 -I/home/openflu/OpenFOAM/OpenFOAM-1.4.1/src/dynamicMesh/lnInclude -I/home/openflu/OpenFOAM/OpenFOAM-1.4.1/src/cfdTools/lnInclude -I/home/openflu/OpenFOAM/OpenFOAM-1.4.1/src/meshTools/lnInclude -IlnInclude -I. -I/home/openflu/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/lnInclude -fPIC -pthread -c $SOURCE -o Make/linuxGccDPOpt/makeAxialMesh.o
In file included from /home/openflu/OpenFOAM/linux/gcc-4.2.1/bin/../lib/gcc/i686-pc-linux-gnu/4.2.1/.. /../../../include/c++/4.2.1/i686-pc-linux-gnu/bits/c++config.h:41,
from /home/openflu/OpenFOAM/linux/gcc-4.2.1/bin/../lib/gcc/i686-pc-linux-gnu/4.2.1/.. /../../../include/c++/4.2.1/string:45,
from /home/openflu/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/lnInclude/string.H:47,
from /home/openflu/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/lnInclude/stringList.H:36,
from /home/openflu/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/lnInclude/argList.H:36,
from makeAxialMesh.C:46:
/home/openflu/OpenFOAM/linux/gcc-4.2.1/bin/../lib/gcc/i686-pc-linux-gnu/4.2.1/.. /../../../include/c++/4.2.1/i686-pc-linux-gnu/bits/os_defines.h:44:22: error: features.h: No such file or directory
In file included from /home/openflu/OpenFOAM/linux/gcc-4.2.1/bin/../lib/gcc/i686-pc-linux-gnu/4.2.1/.. /../../../include/c++/4.2.1/bits/char_traits.h:45,
from /home/openflu/OpenFOAM/linux/gcc-4.2.1/bin/../lib/gcc/i686-pc-linux-gnu/4.2.1/.. /../../../include/c++/4.2.1/string:47,
from /home/openflu/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/lnInclude/string.H:47,
from /home/openflu/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/lnInclude/stringList.H:36,
from /home/openflu/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/lnInclude/argList.H:36,
from makeAxialMesh.C:46:
/home/openflu/OpenFOAM/linux/gcc-4.2.1/bin/../lib/gcc/i686-pc-linux-gnu/4.2.1/.. /../../../include/c++/4.2.1/cstring:52:20: error: string.h: No such file or directory
In file included from
...
(very long list of similar error messages)
...

/home/openflu/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/lnInclude/argList.H:36,
from makeAxialMesh.C:46:
/home/openflu/OpenFOAM/linux/gcc-4.2.1/bin/../lib/gcc/i686-pc-linux-gnu/4.2.1/.. /../../../include/c++/4.2.1/i686-pc-linux-gnu/bits/gthr.h:122: error: expected declaration before end of line
make: *** [Make/linuxGccDPOpt/makeAxialMesh.o] Error 1

I can't understand where I went wrong...
What's wrong??

Simone

bonzodeb April 21, 2008 06:27

I try also with the command s
 
I try also with the command
svn checkout https://openfoam-extend.svn.sourceforge.net/svnroot/openfoam-extend/trunk/Breede r/utilities/mesh/manipulation/MakeAxialMesh

but there is the same error.
It says "not such file or directory", but the file is there...


Solutions??

I try to write the file blockMeshDict by hand:
I want to simulate a pipe r=0.5 m
This is my blockMeshDict file:
/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.4 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

FoamFile
{
version 2.0;
format ascii;

root "";
case "";
instance "";
local "";

class dictionary;
object blockMeshDict;
}

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


convertToMeters 1;

vertices
(
( 0 0 0 ) // vertex 0
( 0 0.5 0 ) // vertex 1
( 0.04358 0.4981 0 ) // vertex 2
( 0 0 100 ) // vertex 3
( 0 0.5 100 ) // vertex 4
( 0.04358 0.4981 100 ) // vertex 5
);

blocks
(
hex (0 2 1 0 3 5 4 3) (1 10 1) simpleGrading (1 1 1)
);

edges
(
);

patches
(
patch inlet
(
(0 1 2 0 )
)

patch outlet
(
(3 5 4 3 )
)

patch wall
(
(1 4 5 2 )
)

wedge a
(
(4 1 0 3)
)

wedge b
(
( 0 2 5 3 )
)

empty collapsed
(
( 0 0 3 3 )
)
);

mergePatchPairs
(
);


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

This is the error
Reading block mesh description dictionary

Creating block mesh

Creating blockCorners

Creating curved edges

Creating blocks

Creating patches

Creating block mesh topology

Default patch type set to empty
#0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) in "/home/openflu/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/openflu/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0xb7f61420]
#3 Foam::wedgePolyPatch::initTransforms() in "/home/openflu/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 Foam::wedgePolyPatch::wedgePolyPatch(Foam::word const&, int, int, int, Foam::polyBoundaryMesh const&) in "/home/openflu/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#5 Foam::polyPatch::addwordConstructorToTable<foam::w edgepolypatch>::New(Foam::word const&, int, int, int, Foam::polyBoundaryMesh const&) in "/home/openflu/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#6 Foam::polyPatch::New(Foam::word const&, Foam::word const&, int, int, int, Foam::polyBoundaryMesh const&) in "/home/openflu/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#7 Foam::polyMesh::polyMesh(Foam::IOobject const&, Foam::Field<foam::vector<double> > const&, Foam::List<foam::cellshape> const&, Foam::List<foam::list<foam::face> > const&, Foam::List<foam::word> const&, Foam::List<foam::word> const&, Foam::word const&, Foam::List<foam::word> const&, bool) in "/home/openflu/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#8 Foam::blockMesh::createTopology(Foam::IOdictionary &) in "/home/openflu/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/blockMesh"
#9 Foam::blockMesh::blockMesh(Foam::IOdictionary&) in "/home/openflu/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/blockMesh"
#10 main in "/home/openflu/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/blockMesh"
#11 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#12 Foam::regIOobject::readIfModified() in "/home/openflu/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/blockMesh"
Floating point exception


Can somebody help me??

Simone

gschaider April 21, 2008 06:56

Hi Simone! Concerning your
 
Hi Simone!

Concerning your compilation problem from the 19th. Seems like some system headers are missing. I assume that
a) you're using the gcc that comes with OpenFOAM (and there is no other gcc installed)
b) you never compiled any applications on your machine before (otherwise the error would have happened there too)
The easiers way to get the necessary headers is by installing the vendor-gcc, which will pull in all the necessary dependencies (glibc-headers etc) (OF will keep on using the OF-gcc, but at least the headers will be there)

Bernhard

bonzodeb April 21, 2008 10:22

You are right! I have been so
 
You are right!
I have been solved:

$ aptitude install build-essential gcc

Thanks Bernhard.

Simone

mahendra January 13, 2009 05:03

hi everyone, I am the getting
 
hi everyone,
I am the getting the following error ,plz help me solve the problem.

Exec : makeAxialMesh Shaft_wall frontAndBackPlanes
Date : Jan 13 2009
Time : 15:43:00
Host : linux
PID : 10815
Case : /home/ms.wankhede/OpenFOAM/ms.wankhede-1.5/mrj/channel
nProcs : 1

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

Create mesh for time = 0

#0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) in "/home/ms.wankhede/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/lib OpenFOAM.so"
#1 Foam::sigSegv::sigSegvHandler(int) in "/home/ms.wankhede/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenF OAM.so"
#2 ?? in "/lib64/libc.so.6"
#3 std::basic_string<char,>, std::allocator<char> >::basic_string(std::string const&) in "/home/ms.wankhede/OpenFOAM/ThirdParty/gcc-4.3.1/platforms/linux64/lib64/libstdc ++.so.6"
#4 main in "/home/ms.wankhede/OpenFOAM/ms.wankhede-1.5/applications/bin/linux64GccDPOpt/mak eAxialMesh"
#5 __libc_start_main in "/lib64/libc.so.6"
#6 __gxx_personality_v0 in "/home/ms.wankhede/OpenFOAM/ms.wankhede-1.5/applications/bin/linux64GccDPOpt/mak eAx ialMesh"
Segmentation fault

gschaider January 13, 2009 06:57

Hi Mahend! Can't tell much
 
Hi Mahend!

Can't tell much from that error message. The strange thing is that it seems to happen during some string-handling operation. Have you got a Debug-enabled OF-version (see http://openfoamwiki.net/index.php/Main_FAQ#An_application_ends_with_a_segmentati on_fault._What_is_wrong.3F and the links going from it)? Could you provide us with a stack-trace from that one? Or the case?

Bernhard

PS: other question: do other applications that YOU compiled work? Have you recently updated your OF-installation?

mahendra January 13, 2009 07:57

Hello Bernhard... Yes I hav
 
Hello Bernhard...

Yes I have compiled the "ComputeTorque", "calcMassFlow" utilities and they work fine. Even the "makeAxialMesh" utility was compiled without any errors.

And I have not updated the OpenFOAM-1.5 to higher version, I am still using OpenFOAM-1.5.

Regards,
Mahendra

gschaider January 13, 2009 17:36

Hi Mahendra! makeAxialMesh
 
Hi Mahendra!

makeAxialMesh never touches the pressure-field (I hope). I'd need the mesh. At least the boundary file

On a second thought: the call you posted looks a bit strange to me (maybe the docu on the Wiki is a bit out of date): try
makeAxialMesh -axis Bottom_wall -wedge frontAndBackPlanes

Bernhard

mahendra January 13, 2009 23:51

Dear Bernhard hello, This i
 
Dear Bernhard hello,

This is what I am getting after using makeAxialMesh -axis Bottom_wall -wedge frontAndBackPlanes

lcfd@clust5:~/OpenFOAM/lcfd-1.5/Mahendra/Pipe> makeAxialMesh -axis Bottom_wall -wedge frontAndBackPlanes
/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.5 |
| \ / A nd | Web: http://www.OpenFOAM.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

Usage: makeAxialMesh <axis> <wedge> [-case dir] [-parallel] [-help] [-doc] [-srcDoc]


Invalid option: -axis
Invalid option: -wedge


FOAM exiting


I think the makeAxialMesh i have downloaded is for OpenFOAM-1.4.1 and I am using OpenFOAM-1.5. Is here the catch??

Regards,
Mahendra.

mahendra January 14, 2009 01:37

Dear Bernhard Hello ! I jus
 
Dear Bernhard Hello !

I just went through the User Guide of OpenFOAM-1.5 and found the thing regarding importing a 2D fluent Mesh to OpenFOAM-1.5 and converting it into axi-symmetric mesh. It says that "Creation of axi-symmetric mesh is currently not supported"

I just wanted to know is really the conversion not available in OpenFOAM-1.5?

Regards,
Mahendra

gschaider January 14, 2009 05:18

Hi Mahendra! @1.5: Definitl
 
Hi Mahendra!

@1.5: Definitly. Get the newest version (as described on the Wiki-page) before trying anything else. That one has the new Syntax.
@fluentMeshToFoam: That was never supported (meshes that were defined in Fluent as axisymetric). That is why makeAxialMesh got written in the first place

Bernhard

mahendra January 14, 2009 05:36

Dear Bernhard... I am unabl
 
Dear Bernhard...

I am unable to download the makeAxialMesh utility for OpenFOAM-1.5 from svn, Can you please help me how to download the utility?

Regards,
Mahendra

gschaider January 14, 2009 07:07

Go to http://openfoam-extend.
 
Go to
http://openfoam-extend.svn.sourceforge.net/viewvc/openfoam-extend/trunk/Breeder_ 1.5/utilities/mesh/manipulation/MakeAxialMesh/
and click on "Download GNU Traball"

Bernhard

mahendra January 15, 2009 06:14

Dear Bernhard, Thanks for t
 
Dear Bernhard,

Thanks for the makeAxialMesh Utility. I am now able to make a axial mesh in OpenFOAM-1.5.

mahendra January 17, 2009 08:54

Dear Bernhard ! One quick q
 
Dear Bernhard !

One quick question, how to specify axis to a grid?

mahendra January 19, 2009 00:50

Dear Bernhard, I am not abl
 
Dear Bernhard,

I am not able to use collapseEdges utility, It is giving me segmentation fault....

Can u shed some light on this?

Regards,
Mahendra

mahendra January 19, 2009 01:26

The error i am getting is like
 
The error i am getting is like this:

Cell:14998 uses faces:6(29649 30248 60346 60347 29349 29648) of which too many are marked for remov
al:
29649 30248 60346 60347 29349 29648
Cell:14999 uses faces:6(30249 30299 60348 60349 29350 29649) of which too many are marked for remov
al:
30249 30299 60348 60349 29350 29649
Morphing ...
--> FOAM Warning :
From function boundBox::boundBox(const pointField& points)
in file meshes/boundBox/boundBox.C at line 52
Cannot find bounding box for zero sized pointField, returning zero
--> FOAM Warning :
From function boundBox::boundBox(const pointField& points)
in file meshes/boundBox/boundBox.C at line 52
Cannot find bounding box for zero sized pointField, returning zero
Collapsing 0 small edges
Collapsing 0 in line edges
#0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) in "/home/lcfd/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOp t/libOpenFOAM.so"
#1 Foam::sigSegv::sigSegvHandler(int) in "/home/lcfd/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/lib OpenFOAM.so"
#2 ?? in "/lib64/libc.so.6"
#3 collapseHighAspectFaces(Foam::polyMesh const&, Foam::PackedList<1> const&, double, double, Foam ::edgeCollapser&) in "/home/lcfd/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/collapseEdg es"
#4 main in "/home/lcfd/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/collapseEdges "
#5 __libc_start_main in "/lib64/libc.so.6"
#6 __gxx_personality_v0 in "/home/lcfd/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/coll apseEdges"
Segmentation fault
lcfd@clust5:~/OpenFOAM/lcfd-1.5/Mahendra/Pipe_trial>

Regards,
Mahendra

gschaider January 19, 2009 11:31

Hi Mahendra! The axis is de
 
Hi Mahendra!

The axis is defined using the axis-option (it assumes that you have a planar patch of that name). With the offset-option you can move the axis away from that line. I'm not sure what happens if the patch is not planar

@the collapseEdge-crash: no idea. At first try to run checkMesh on your original (axial) mesh to see whether it is valid. It seems to me that the utilitiy is trying to remove too many edges. Decrease the tolerance parameter (what value are you using now?)

Bernhard

mahendra January 20, 2009 01:53

Dear Bernhard, I am using a
 
Dear Bernhard,

I am using a edge length as 3m and angle as 5 degree, what i found was when I used edge length as 0.001m, the collapseEdges utility worked fine saying that 0 edges and 0 faces modified.

I assumed that since my pipe length is 3m and wedge angle is 5 degrees I should use them while specifying the arguments to collapseEdges utilty and it resulted in segmentation fault.

Was I correct in using it or not???

Regards,
Mahendra

gschaider January 20, 2009 12:12

Hi Mahendra! The first para
 
Hi Mahendra!

The first parameter means "all edges shorter that this will be removed". Which would mean almost all edges. Vary the feature angle

Bernhard

sega January 20, 2009 15:27

Hello! I have just learned
 
Hello!

I have just learned that my meshes created with makeAxialMesh are not 1 cell but 5 cells thick (in the wedge direction).

Did I do something wrong?
Do you have any idea what may be the problem?

Greetings.
Sebastian.

gschaider January 21, 2009 15:33

Hi Sebastian! Have you chec
 
Hi Sebastian!

Have you checked the original (2D planar) mesh? This shouldn't happen unless the original is that thick, too

Bernhard

sega January 21, 2009 16:10

Yes I did. The mesh is just
 
Yes I did.

The mesh is just fine in both 2D and wedge.
Due to the big amount of cells and their small size this was just my mistake in seeing them correctly.

Sorry.

olwi January 21, 2009 16:34

Mahendra: Please have a look a
 
Mahendra: Please have a look at my post of September 9, 2008: http://www.cfd-online.com/OpenFOAM_D.../126/9114.html

I think your problem is in large differences in sizes of cells/faces in your mesh. Try the modified "collapseEdgesBetter" in my post. Use the option -areaFactor to specify an "areafactor" smaller than the default 1e-9. Maybe 1e-11?

/Ola

sega January 21, 2009 17:07

No, that is not the problem.
 
No, that is not the problem.
I was just too blind to see them correctly.
But thanks for your advice.!

mahendra January 21, 2009 23:48

Hi Bernhard, Still I am getti
 
Hi Bernhard,
Still I am getting the segmentation fault error with using collapseEdgesBetter.

Cell:49999 uses faces:6(100949 100999 201048 201049 97950 98949) of which too many are marked for removal:
100949 100999 201048 201049 97950 98949
Morphing ...
--> FOAM Warning :
From function boundBox::boundBox(const pointField& points)
in file meshes/boundBox/boundBox.C at line 52
Cannot find bounding box for zero sized pointField, returning zero
--> FOAM Warning :
From function boundBox::boundBox(const pointField& points)
in file meshes/boundBox/boundBox.C at line 52
Cannot find bounding box for zero sized pointField, returning zero
Collapsing 0 small edges
Collapsing 0 in line edges
#0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) in "/home/lcfd/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigSegv::sigSegvHandler(int) in "/home/lcfd/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib64/libc.so.6"
#3 collapseHighAspectFaces(Foam::polyMesh const&, Foam::PackedList<1> const&, double, double, Foam::edgeCollapser&) in "/home/lcfd/OpenFOAM/lcfd-1.5/applications/bin/linux64GccDPOpt/collapseEdgesBett er"
#4 main in "/home/lcfd/OpenFOAM/lcfd-1.5/applications/bin/linux64GccDPOpt/collapseEdgesBett er"
#5 __libc_start_main in "/lib64/libc.so.6"
#6 __gxx_personality_v0 in "/home/lcfd/OpenFOAM/lcfd-1.5/applications/bin/linux64GccDPOpt/collapseEdgesBett er"
Segmentation fault
lcfd@clust5:~/OpenFOAM/lcfd-1.5/Mahendra/Pipe_2D>

Regards,
Mahendra

gschaider January 22, 2009 10:23

Hi Mahendra! Try skipping t
 
Hi Mahendra!

Try skipping the collapseEdges-Stuff. Simply set all scalars to zeroGradient and the vectors to slip on the axis-patch (I think that are the settings that work) and try to use the mesh that way

Bernhard

squadron February 6, 2009 05:35

Hi, I have created a 3D mes
 
Hi,

I have created a 3D mesh with Gmsh, 1 element thick, and used gmshToFoam to convert it. checkMesh gave an "OK"

Now I want to use makeAxialMesh to convert my mesh into an axi-symmetrical one. This was the result:

Create time

Create mesh for time = 0

Plane of the grid: (0 0 1) (0.105253 0.0174027 0.0005)

The rotation-axis: ((-0.06 0 0.0005) (0.94 0 0.0005))

Creating wedge with an opening angel of 5 degrees

Radius to axis: min = 1.0842e-19 max 0.05
#0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) in "/home/age/OpenFoam/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigSegv::sigSegvHandler(int) in "/home/age/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0xb7ff3400]
#3 getDistance(Foam::polyPatch const&, Foam::line<foam::vector<double>, Foam::Vector<double> >&) in "/home/age/OpenFOAM/age-1.5/applications/bin/linuxGccDPOpt/makeAxialMesh"
#4 main in "/home/age/OpenFOAM/age-1.5/applications/bin/linuxGccDPOpt/makeAxialMesh"
#5 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#6 _start in "/home/age/OpenFOAM/age-1.5/applications/bin/linuxGccDPOpt/makeAxialMesh"
Segmentation fault

What did I do wrong?

FYI, I'm not only an OpenFoam newby, but also a Linux newby, so please explain the whole process step by step.

elorriaux February 6, 2009 08:30

Hello, i suppose you have f
 
Hello,

i suppose you have followed the wiki instructions and checked the different constraints : http://openfoamwiki.net/index.php/Contrib_MakeAxialMesh

As you have already your model in gmsh, you can directly revolve your mesh in gmsh :

- rotate -2.5° your plane from your main 2D plane
- extrude (revolve) 5°
- save mesh
- gmshToFoam

I've already done that with gmsh, it works fine.

You can also have another solution by revolving your mesh in OpenFoam with the extrudeMesh tool, but it's not the easier way to do, i think.

Good luck.

gschaider February 6, 2009 08:35

Hi Age! No idea (havn't see
 
Hi Age!

No idea (havn't seen that stacktrace). How big is the mesh? Could you pass the mesh to me and I'll have a look

Bernhard

squadron February 6, 2009 09:49

@ Etienne: So you don't use m
 
@ Etienne:
So you don't use makeAxialMesh at all? Just extrude using rotation in Gmsh?

@ Bernhard:
The mesh-file is very large. I will try some other things first, like simple geo's etc.

Thanks guys

elorriaux February 6, 2009 10:52

It depends on the case, if my
 
It depends on the case, if my geometry is done in gmsh, i revolve the mesh in gmsh. If i already have a mesh in gambit or a blockMeshDict, i use makeAxialMesh.

All those procedures should work, i've tried successfully many different ways to build axy cases in OpenFOAM.

If your msh file is too large, perhaps you can zip the .geo, it should be easier to transfer.


All times are GMT -4. The time now is 18:12.