|
[Sponsors] |
[blockMesh] merging two cylinders to create a t junction pipe |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
July 19, 2017, 10:31 |
merging two cylinders to create a t junction pipe
|
#1 |
New Member
ozzy
Join Date: May 2017
Posts: 3
Rep Power: 8 |
Hello
I am trying to create in block mesh a t juntion using two cylinders. I managed to run blockmesh without a problem but when i run the the code it complains that FOAM FATAL ERROR: Number of poly-patches = 4 in blockMeshDict, are not equal to the number of patch models = 3, de$ if you look at my blockmesh and my drawing(http://imgur.com/a/83nhB), i already know i made a patch called mergewall and a masterpatch and these are not defined in the boundarydict because there is nothing to define. a part of the wall merge wall should merge with the masterpatch and it should create an internal boundary that i dont have to define. I don't know how to do this. I have seen some related posts but so far couldn't do it. thanks for any help and i ll post a solution if i figure it out. /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: dev | | \\ / A nd | Web: www.OpenFOAM.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object blockMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // convertToMeters 1e-2; vertices ( ( 0.100000 -0.100000 0.000000) ( -0.100000 -0.100000 0.000000) ( -0.100000 0.100000 0.000000) ( 0.100000 0.100000 0.000000) ( 0.141421 -0.141421 0.000000) ( -0.141421 -0.141421 0.000000) ( -0.141421 0.141421 0.000000) ( 0.141421 0.141421 0.000000) ( 0.100000 -0.100000 1.000000) ( -0.100000 -0.100000 1.000000) ( -0.100000 0.100000 1.000000) ( 0.100000 0.100000 1.000000) ( 0.141421 -0.141421 1.000000) ( -0.141421 -0.141421 1.000000) ( -0.141421 0.141421 1.000000) ( 0.141421 0.141421 1.000000) ( 0.050000 -0.183333 0.450000) ( -0.050000 -0.183333 0.450000) ( -0.050000 -0.183333 0.550000) ( 0.050000 -0.183333 0.550000) ( 0.070711 -0.183333 0.429289) ( -0.070711 -0.183333 0.429289) ( -0.070711 -0.183333 0.570711) ( 0.070711 -0.183333 0.570711) ( 0.050000 -1.183333 0.450000) ( -0.050000 -1.183333 0.450000) ( -0.050000 -1.183333 0.550000) ( 0.050000 -1.183333 0.550000) ( 0.070711 -1.183333 0.429289) ( -0.070711 -1.183333 0.429289) ( -0.070711 -1.183333 0.570711) ( 0.070711 -1.183333 0.570711) ); blocks ( //block0 hex (1 0 3 2 9 8 11 10) square (5 5 50) simpleGrading (1 1 1) //block1 hex (0 4 7 3 8 12 15 11) innerCircle (5 5 50) simpleGrading (1 1 1) //block2 hex (3 7 6 2 11 15 14 10) innerCircle (5 5 50) simpleGrading (1 1 1) //block3 hex (2 6 5 1 10 14 13 9) innerCircle (5 5 50) simpleGrading (1 1 1) //block4 hex (1 5 4 0 9 13 12 8) innerCircle (5 5 50) simpleGrading (1 1 1) //block5 hex ( 17 16 19 18 25 24 27 26 ) square (5 5 10) simpleGrading (1 1 1) //block5 hex ( 16 20 23 19 24 28 31 27 ) innerCircle (5 5 10) simpleGrading (1 1 1) //block7 hex ( 19 23 22 18 27 31 30 26 ) innerCircle (5 5 10) simpleGrading (1 1 1) //block8 hex ( 18 22 21 17 26 30 29 25 ) innerCircle (5 5 10) simpleGrading (1 1 1) //block9 hex ( 17 21 20 16 25 29 28 24 ) innerCircle (5 5 10) simpleGrading (1 1 1) ); edges ( arc 7 4 ( 0.200000 0.000000 0.000000) arc 4 5 ( 0.000000 -0.200000 0.000000) arc 5 6 ( -0.200000 0.000000 0.000000) arc 6 7 ( 0.000000 0.200000 0.000000) arc 15 12 ( 0.200000 0.000000 1.000000) arc 12 13 ( 0.000000 -0.200000 1.000000) arc 13 14 ( -0.200000 0.000000 1.000000) arc 14 15 ( 0.000000 0.200000 1.000000) arc 3 0 ( 0.110000 0.000000 0.000000) arc 0 1 ( 0.000000 -0.110000 0.000000) arc 1 2 ( -0.110000 0.000000 0.000000) arc 2 3 ( 0.000000 0.110000 0.000000) arc 11 8 ( 0.110000 0.000000 1.000000) arc 8 9 ( 0.000000 -0.110000 1.000000) arc 9 10 ( -0.110000 0.000000 1.000000) arc 10 11 ( 0.000000 0.110000 1.000000) arc 23 20 ( 0.100000 -0.183333 0.500000) arc 20 21 ( 0.000000 -0.183333 0.400000) arc 21 22 ( -0.100000 -0.183333 0.500000) arc 22 23 ( 0.000000 -0.183333 0.600000) arc 31 28 ( 0.100000 -1.183333 0.500000) arc 28 29 ( 0.000000 -1.183333 0.400000) arc 29 30 ( -0.100000 -1.183333 0.500000) arc 30 31 ( 0.000000 -1.183333 0.600000) arc 19 16 ( 0.060000 -0.183333 0.500000) arc 16 17 ( 0.000000 -0.183333 0.440000) arc 17 18 ( -0.060000 -0.183333 0.500000) arc 18 19 ( 0.000000 -0.183333 0.560000) arc 27 24 ( 0.060000 -1.183333 0.500000) arc 24 25 ( 0.000000 -1.183333 0.440000) arc 25 26 ( -0.060000 -1.183333 0.500000) arc 26 27 ( 0.000000 -1.183333 0.560000) ); patches ( wall walls ( (3 0 1 2) (3 7 4 0) (2 6 7 3) (1 5 6 2) (0 4 5 1) (4 7 15 12) (6 5 13 14) (7 6 14 15) ( 20 23 31 28 ) ( 21 20 28 29 ) ( 22 21 29 30 ) ( 23 22 30 31 ) ) patch inlet ( ( 27 26 25 24 ) ( 27 24 28 31 ) ( 26 27 31 30 ) ( 25 26 30 29 ) ( 24 25 29 28 ) ) patch outlet ( (11 10 9 8) (11 8 12 15) (10 11 15 14) (9 10 14 13) (8 9 13 12) ) patch masterPatch ( ( 19 16 17 18 ) ( 19 23 20 16 ) ( 18 22 23 19 ) ( 17 21 22 18 ) ( 16 20 21 17 ) ) wall mergewall ( (5 4 12 13) ) ); mergePatchPairs ( (masterPatch mergewall) ); // ************************************************** *********************** // |
|
July 20, 2017, 06:09 |
|
#2 |
New Member
ozzy
Join Date: May 2017
Posts: 3
Rep Power: 8 |
I was able to run the code however it quits quickly after. Here is what i did to be able run it :
(maybe the solution i am looking for is somewhere here Remove internal patch(es)) i run the command blockMesh then i run the command stitchMesh masterPatch slavePatch this creates a 2e-06 folder with polymesh folder inside same structure as the one in constant folder but there is no blockMeshDict so i replace the files in constant with these cp ../../2e-06/polyMesh/* . and i do the following changes in the boundary file i correct the folder location to constant/polyMesh instead of 2e-06/polyMesh then i correct the number of patches from 5 to 4 then i remove the masterZone which has 0 faces then i go back and run the usual ./Allrun command with the blockMesh command removed so the code starts to run but there is this error in the log file --> FOAM FATAL ERROR: No base point for face 23517, 4(7652 7653 7658 7657), produces a decomposition that has a minimum volume greater than tolerance. Last edited by decaf; July 20, 2017 at 09:18. |
|
July 27, 2017, 12:42 |
solution
|
#3 |
New Member
ozzy
Join Date: May 2017
Posts: 3
Rep Power: 8 |
The solution is forget about blockmesh.
I am using os X so what i ended up doing is : i installed autocad which has some free acadamic version to make geometries and i created easily using only line, polyline, rectangle circle and sweep, merge commands complicated geometries i needed. then i used inventor fusion, again free, and installable on os X to convert a file to step file. then i installed virtualbox to install ubuntu to install salome there is this really nice video that describes hot to mesh starting from a step file https://www.youtube.com/watch?v=1zQbU-E4k1U salome runs absolutely fine on my laptop and i easily created the mesh one little thing i had to was to scale my geometry because for some reason i started in mm ended up in cm but scaling is trivial in salome. then as described in the video convert the mesh to files for openfoam using ideasUnvToFoam yourfile.unv i also did checkMesh which came out all ok but checkmesh -allgeometry gave 2 fails maybe its just a warning. the code (dsmcFoamPlus) runs fine, now i need to see if it makes any sense physicswise. my gemorty is not a t juntion anymore since i can do any geometry i want : ) http://imgur.com/a/BINKZ |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[ANSYS Meshing] t junction pipe. | aqil1 | ANSYS Meshing & Geometry | 0 | November 28, 2018 04:14 |
need help about double pipe heat exchanger with chtMultiRegionSimpleFoam | wuyangzhen | OpenFOAM | 10 | December 12, 2017 00:19 |
How do I create vector plots in the cross section of a pipe? | Raff94 | FLUENT | 1 | February 16, 2016 00:55 |
create plane for 3D pipe | Tolgahan | FLUENT | 2 | May 7, 2007 14:14 |
Terrible Mistake In Fluid Dynamics History | Abhi | Main CFD Forum | 12 | July 8, 2002 09:11 |