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 |
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 |
I have done that in other simulations and it has not been a problem. It corresponds to the first block that is a triangle.
|
Quote:
If you mesh this way cell faces have zero area. |
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. |