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/)
-   -   [blockMesh] Meshing a wedge (https://www.cfd-online.com/Forums/openfoam-meshing/185336-meshing-wedge.html)

veronica2015 March 23, 2017 09:46

Meshing a wedge
 
I am creating a mesh of a simple wedge composed by two blocks. The case with another geometry runs well, but when I put my data I am getting the following error.

Create time

Creating block mesh from
"/home/veronica/Escritorio/new_simulations/KandM_H-D-10_2/constant/polyMesh/blockMeshDict"
Creating curved edges
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::arcEdge::calcAngle() in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libblockMesh.so"
#4 Foam::arcEdge::arcEdge(Foam::Field<Foam::Vector<do uble> > const&, Foam::Istream&) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libblockMesh.so"
#5 Foam::curvedEdge::addIstreamConstructorToTable<Foa m::arcEdge>::New(Foam::Field<Foam::Vector<double> > const&, Foam::Istream&) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libblockMesh.so"
#6 Foam::curvedEdge::New(Foam::Field<Foam::Vector<dou ble> > const&, Foam::Istream&) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libblockMesh.so"
#7 Foam::blockMesh::createTopology(Foam::IOdictionary const&, Foam::word const&) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libblockMesh.so"
#8 Foam::blockMesh::blockMesh(Foam::IOdictionary const&, Foam::word const&) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libblockMesh.so"
#9
in "/opt/openfoam230/platforms/linuxGccDPOpt/bin/blockMesh"
#10 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#11
in "/opt/openfoam230/platforms/linuxGccDPOpt/bin/blockMesh"
Excepción de coma flotante (`core' generado)

The blockMeshDict is:

convertToMeters 1;

vertices
(
(0 0 0) //0
(3.174516432e-3 0 5.541139044e-5) //1
(0.1904709859 0 3.324683426e-3) //2
(0 0.635 0) //3
(3.174516432e-3 0.635 5.541139044e-5) //4
(0.1904709859 0.635 3.324683426e-3) //5
(3.174516432e-3 0 -5.541139044e-5) //6
(0.1904709859 0 -3.324683426e-3) //7
(3.174516432e-3 0.635 -5.541139044e-5) //8
(0.1904709859 0.635 -3.324683426e-3) //9

);

blocks
(
hex (1 4 3 0 6 8 3 0) (34 112 1) simpleGrading (1 2.05671982525299e-4 1) //block 1
hex (2 5 4 1 7 9 8 6) (64 112 1) simpleGrading (150.0842276 2.05671982525299e-4 1) //block 2

);

edges
(
arc 2 7 (0.1905 0 0)
arc 5 12 (0.1905 0.635 0)
);

boundary
(
inlet
{
type patch;
faces
(
(3 4 8 3)

);
}
outlet
{
type patch;
faces
(
(2 7 9 5)

);
}
upperWall
{
type wall;
faces
(
(4 5 9 8)

);
}
lowerWall
{
type wall;
faces
(
(0 6 1 0)
(1 6 7 2)

);
}

wedge1
{
type wedge;
faces
(
(0 1 4 3)
(1 2 5 4)

);
}
wedge2
{
type wedge;
faces
(
(0 3 8 6)
(6 8 9 7)

);
}
);

mergePatchPairs
(
);

Any help please?

Thanks so much

piu58 March 23, 2017 12:11

One mistake is for sure, that you used nodes twice in hex element 1:
Code:

hex (1 4 3 0 6 8 3 0) (34 112 1) simpleGr

veronica2015 March 23, 2017 12:25

I have done that in other simulations and it has not been a problem. It corresponds to the first block that is a triangle.

piu58 March 23, 2017 13:06

Quote:

Originally Posted by veronica2015 (Post 642032)
I have done that in other simulations and it has not been a problem. It corresponds to the first block that is a triangle.

If you need a triangle take one: OF offers triangles.

If you mesh this way cell faces have zero area.

tareqkh March 23, 2017 17:04

What shape are you tying to apply wedge on? I could make that in gmsh quickly. Just provide your dimensions and it should be easy.

Cheers,


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