CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[blockMesh] merging two cylinders to create a t junction pipe

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 19, 2017, 10:31
Default merging two cylinders to create a t junction pipe
  #1
New Member
 
ozzy
Join Date: May 2017
Posts: 3
Rep Power: 8
decaf is on a distinguished road
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)
);


// ************************************************** *********************** //
decaf is offline   Reply With Quote

Old   July 20, 2017, 06:09
Default
  #2
New Member
 
ozzy
Join Date: May 2017
Posts: 3
Rep Power: 8
decaf is on a distinguished road
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.
decaf is offline   Reply With Quote

Old   July 27, 2017, 12:42
Default solution
  #3
New Member
 
ozzy
Join Date: May 2017
Posts: 3
Rep Power: 8
decaf is on a distinguished road
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
decaf is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


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


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