could anyone post a simple cylinder mesh
hi, how do I do a cylinder with blockmesh? Is it necessary do a quarter of a cylinder and stablish simmetry plane in the lateral faces of it or exists a mean of make a entire cylinder? My cylinder is totally wrong I think.
Could anyone posts a simple mesh to I see? |
Hi!
I don't have a blockMes
Hi!
I don't have a blockMesh Handy, but it's actualy quite easy once you have the general idea: you've got to compose the cylinder of _five_ blocks: one at the center and 4 blocks "bolted" to that (let's call them N,E,S and W according to the side of the center block they're glued to) For instance the N block shares his S-side with the centre-block, E-side with block E, W-side with block W and the N-side is part of the outer boundary of the cylinder |
Hi Guilherme,
I have a cylind
Hi Guilherme,
I have a cylinder-script file you could try. It does just what Bernhard suggests, and is easy to adapt to new number of cells/dimensions. To get a blockMeshDict file from it, do m4 cylinderMesh.m4 > blockMeshDict from a terminal window. If it doesn't work, you will need to install the m4-preprocessor. //Rasmus http://www.cfd-online.com/OpenFOAM_D...hment_icon.gif cylinderMesh.m4 |
Thanks Rasmus
The mesh work
Thanks Rasmus
The mesh works fine. Bye. |
Hi all,
I'm starting a new
Hi all,
I'm starting a new project which deals with a Large Eddy Simulation of an axisymmetric jet on a rotating and heated plate. The sketch of the computational domain is shown on the present figure. http://www.cfd-online.com/OpenFOAM_D...ges/1/4061.png The jet nozzle is circular and the jet is impinging on an heated rotor (in red on the attached picture). I have read some posts on cylinder mesh and was wondering about the best strategy for my configuration: - use the cylinderMesh.m4 script-file to decompose the cylinder with a square like on this picture: http://www.cfd-online.com/OpenFOAM_D...ges/1/4062.png - use a wedge type meshing strategy like on this picture: http://www.cfd-online.com/OpenFOAM_D...ges/1/4060.jpg For the later case, I was not able to find some useful info to make a 3D mesh with wedge type cells with blockmesh. Is it possible in OpenFOAM1.3 or this feature is only available for pseudo 3D domain (one cell in the third direction) ? Thanks for your help or suggestions. Francois |
Thanks Ville,
As you sugges
Thanks Ville,
As you suggested I will start with the with a 9 blocks mesh and play a bit with the cylinderMesh.m4 preprocessor script. But I don't realy understand why I'll loose in the resolution with wedge blocks. Anyway thanks a lot for your suggestions. I will post my blockMeshDict file when I'm done with it ... Have a nice day. Francois |
Hi, everyone
I am also inte
Hi, everyone
I am also interested in the wedge blocks,because when I use multiblock to simulate the spray process, it seems like the drops are allways spray from the corner of the "centre block". But, after I've read the userguider, I find that it seems like it's impossible to get the wedge blocks like the picture before using blockmesh? Could anyone can help me ? Thank you ~~~! Bobby |
Hello Cedric
Sorry to say,
Hello Cedric
Sorry to say, but that is not your only mistake... I have taken a look at your file and there is a lot wrong with it. Blocks are intersecting and not defined using a right-hand rule, you're arcs are also ill defined. Could you give an example (picture) of what you're geometry should look like? I've modified your file into the following, but I don't think this is the geometry you're looking for... /*-----------------------------------------------------*\ | ========= | | | \ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \ / O peration | Version: 1.3 | | \ / A nd | Web: http://www.openfoam.org | | \/ M anipulation | | \*-----------------------------------------------------*/ FoamFile { version 2.0; format ascii; root ""; case ""; instance ""; local ""; class dictionary; object blockMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * *// //meshGenApp blockMesh; convertToMeters 0.001; vertices ( ( 3.0 -3.0 0.0) ( 1.0 -1.0 0.0) ( 1.0 1.0 0.0) ( 3.0 3.0 0.0) ( -3.0 3.0 0.0) ( -1.0 1.0 0.0) ( -1.0 -1.0 0.0) ( -3.0 -3.0 0.0) ( 6.0 -6.0 0.0) ( 6.0 6.0 0.0) ( -6.0 6.0 0.0) //10 ( -6.0 -6.0 0.0) ( 6.0 -6.0 10.0) ( 3.0 -3.0 10.0) ( 3.0 3.0 10.0) ( 6.0 6.0 10.0) ( -6.0 6.0 10.0) ( -3.0 3.0 10.0) ( -3.0 -3.0 10.0) ( -6.0 -6.0 10.0) ); blocks ( hex (1 2 5 6 13 14 17 18) (6 6 10) simpleGrading (1 1 1) hex (2 3 4 5 14 15 16 17) (2 6 10) simpleGrading (1 1 1) hex (5 4 7 6 17 16 19 18) (2 6 10) simpleGrading (1 1 1) hex (0 1 6 7 12 13 18 19) (2 6 10) simpleGrading (1 1 1) hex (0 3 2 1 12 15 14 13) (6 2 10) simpleGrading (1 1 1) ); edges ( arc 5 2 (0.0 1.414214 0.0) arc 6 5 (-1.414214 0.0 0.0) arc 1 6 (0.0 -1.414214 0.0) arc 2 1 (1.414214 0.0 0.0) arc 17 14 (0.0 4.242641 10.0) arc 18 17 (-4.242641 0.0 10.0) arc 13 18 (0.0 -4.242641 10.0) arc 14 13 (4.242641 0.0 10.0) ); patches ( patch inlet ( (1 6 5 2) (3 2 5 4) (4 5 6 7) (7 6 1 0) (0 1 2 3) ) patch outlet ( (13 14 17 18) (14 15 16 17) (17 16 19 18) (12 13 18 19) (12 15 14 13) ) wall walls ( (3 4 16 15) (4 7 19 16) (7 0 12 19) (0 3 15 12) ) ); mergePatchPairs ( ); |
Hi Guido,
Thank you for your
Hi Guido,
Thank you for your reply. I hope it hadn't disturbed you too long ;o) this is my final geometry: http://www.cfd-online.com/OpenFOAM_D...ges/1/4361.jpg I made the left part (the diffuser) and it works correctly, but when I add the large dump downstream of the diffuser section, some mistake are coming. so what I gave you yesterday is only the dump with the outlet section of the diffuser. I also would like to use a buterfly mesh (see Francois Beaubert message upper) I give all also a piece of fortran code to do that easily where I can modifie all my parameters : http://www.cfd-online.com/OpenFOAM_D...s/mime_txt.gif CreateGridFoam_clausen.f.txt Cedric Thank you again for helping |
Hello Cedric,
I didn't take
Hello Cedric,
I didn't take me too long :-). Following your example, I have created a small diffuser with a dump, maybe this can point you in the right direction. I have made 7 blocks, 2 for the diffuser and 5 for the dump area. The dump area has been constructed as O-grid (You referred to it as a butterfly mesh, I think it is the same) The blocks should be created with the right-hand-rule. The arcs are created by the two corner points and a point on the line in between (in this case in the middle). The patches should be numbered in such order that the normal points outward. If you have anymore questions, let me know. :-) Guido /*-----------------------------------------------------*\ | ========= | | | \ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \ / O peration | Version: 1.3 | | \ / A nd | Web: http://www.openfoam.org | | \/ M anipulation | | \*-----------------------------------------------------*/ FoamFile { version 2.0; format ascii; root ""; case ""; instance ""; local ""; class dictionary; object blockMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * *// //meshGenApp blockMesh; convertToMeters 0.001; vertices ( ( 1.0 -1.0 -5.0) //0 ( 1.0 1.0 -5.0) ( -1.0 1.0 -5.0) ( -1.0 -1.0 -5.0) ( 1.0 -1.0 -2.0) //4 ( 1.0 1.0 -2.0) ( -1.0 1.0 -2.0) ( -1.0 -1.0 -2.0) ( 3.0 -3.0 0.0) //8 ( 3.0 3.0 0.0) ( -3.0 3.0 0.0) ( -3.0 -3.0 0.0) ( 6.0 -6.0 0.0) //12 ( 6.0 6.0 0.0) ( -6.0 6.0 0.0) ( -6.0 -6.0 0.0) ( 3.0 -3.0 10.0) //16 ( 3.0 3.0 10.0) ( -3.0 3.0 10.0) ( -3.0 -3.0 10.0) ( 6.0 -6.0 10.0) //20 ( 6.0 6.0 10.0) ( -6.0 6.0 10.0) ( -6.0 -6.0 10.0) ); blocks ( hex (0 1 2 3 4 5 6 7) (4 4 3) simpleGrading (1 1 1) hex (4 5 6 7 8 9 10 11) (4 4 2) simpleGrading (1 1 1) hex (8 9 10 11 16 17 18 19) (4 4 10) simpleGrading (1 1 1) hex (14 10 9 13 22 18 17 21) (6 4 10) simpleGrading (3 1 1) hex (15 11 10 14 23 19 18 22) (6 4 10) simpleGrading (3 1 1) hex (12 8 11 15 20 16 19 23) (6 4 10) simpleGrading (3 1 1) hex (13 9 8 12 21 17 16 20) (6 4 10) simpleGrading (3 1 1) ); edges ( arc 2 1 (0.0 1.414214 -5.0) arc 1 0 (1.414214 0.0 -5.0) arc 0 3 (0.0 -1.414214 -5.0) arc 3 2 (-1.414214 0.0 -5.0) arc 6 5 (0.0 1.414214 -2.0) arc 5 4 (1.414214 0.0 -2.0) arc 4 7 (0.0 -1.414214 -2.0) arc 7 6 (-1.414214 0.0 -2.0) arc 10 9 (0.0 4.242641 0.0) arc 9 8 (4.242641 0.0 0.0) arc 8 11 (0.0 -4.242641 0.0) arc 11 10 (-4.242641 0.0 0.0) arc 14 13 (0.0 8.485281 0.0) arc 13 12 (8.485281 0.0 0.0) arc 12 15 (0.0 -8.485281 0.0) arc 15 14 (-8.485281 0.0 0.0) arc 18 17 (0.0 4.242641 10.0) arc 17 16 (4.242641 0.0 10.0) arc 16 19 (0.0 -4.242641 10.0) arc 19 18 (-4.242641 0.0 10.0) arc 22 21 (0.0 8.485281 10.0) arc 21 20 (8.485281 0.0 10.0) arc 20 23 (0.0 -8.485281 10.0) arc 23 22 (-8.485281 0.0 10.0) ); patches ( patch inlet ( (0 3 2 1) ) patch outlet ( (16 17 18 19) (17 21 22 18) (18 22 23 19) (20 16 19 23) (20 21 17 16) ) ); mergePatchPairs ( ); |
Hello Guido,
Thank you very m
Hello Guido,
Thank you very much, it seems to work. I'm trying now to improve the mesh. If I have a problem, I won't hesitate to disturb you :-) Cedric |
Hello OpenFOAM users
I am n
Hello OpenFOAM users
I am new to OpenFOAM. With the help of forum knowledgebase I have been able to install and create meshes for some simple geometries. My interest is to simulate flow in a pipe with bend. Is it possible to generate geometry + mesh using the blockMesh tool. Looking forward to suggestions Kind Regards Jaswi |
You're misspelling fixedValue
You're misspelling fixedValue , capital V.
good luck |
Hello Rasmus or anyone else,
Hello Rasmus or anyone else,
I have tried to use the script provided by Rasmus. When I executed it, it worked fine. But to visualize it in paraFoam, when I tried, paraFoam window got open but on clicking "accept" in it, the window closed automatically showing the following FATAL ERROR message. [nikhil@localhost ~]$ paraFoam /home/nikhil/OpenFOAM/nikhil-1.4.1/run/tutorials/icoFoam cavityGrade --> FOAM FATAL IO ERROR : size 400 is not equal to the given value of 450 file: /home/nikhil/OpenFOAM/nikhil-1.4.1/run/tutorials/icoFoam/cavityGrade/0.8/p from line 25 to line 445. From function Field<type>::Field(const word& keyword, const dictionary& dict, const label s) in file lnInclude/Field.C at line 224. FOAM exiting I am pasting down the steps below. Somebody please tell me why I am not able to visualise in paraFoam and how I can make the necessary changes if required. Thank you, nikhil |
HI!!....perhaves you can visul
HI!!....perhaves you can visulize mesh by first clik on "constant" in paraFoam window and than "accept"..may be this error is because of differnt parameters.
regards Bhuvnesh Verma |
Hi Nikhil
The reason for yo
Hi Nikhil
The reason for your problems is seen in the following error: -> FOAM FATAL IO ERROR : size 400 is not equal to the given value of 450 It means that the number of elements in your data files, in this case in specific the pressure at t=0.8, does not correspond to the number of elements in your mesh. Best regards, Niels |
thanq bhuvanesh. ur suggestion
thanq bhuvanesh. ur suggestion really works.
thanq Niels. but i dint understand why itz not working when we click say "0.8" instead of "constant" in paraFoam. i mean, which parameters will change if we change the above options(constant and 0.8)? regards, nikhil |
The reason is, that in the 'co
The reason is, that in the 'constant' directory, the physical variables, i.e. pressure, velocity, etc is (probably) not loaded, thus a mismatch in the number of computational cells and the number of variables in the physical fields is not checked, thus you are allowed to see your mesh.
/ Niels |
HI!...reason is as Niels menti
HI!...reason is as Niels mentioned !..so for that you may try..'move to case->constant->..there you delet all neigh..,owne..ets previous stuff expect 'blockmeshdict'.than again form these stuff by blockmesh ./case.,than solve ur case. i think it must work now for anytime .
regards Bhuvnesh Verma |
thanq Niels n Bhuvanesh.
u
thanq Niels n Bhuvanesh.
u both were true. itz working when i deleted the old files, the problem was coz of new stuff not getting loaded in the presence of old files. regards, nikhil |
hi all,
I took your examples and try to adapt it to my case, where i only have to make a simple cylinder that i need to mesh inside as regularly as possible. the first solution proposed doesn't work on my version of openFOAM, but the others do but its a bit tricky, since i'm still learning. If someone has an example of what i need, please post it. Anyway thank you very much, i learned a lot from this discussion. :D |
it's ok, i did it ! :cool:
|
Quote:
I have the same problem! Did you earned the mesh?? I'm new with the OPenFOAM and for me it's difficult to understand how it works the mesh generation. Thanks! |
@libia87
OpenFOAM mesh how does it work? This is a broad question, and really depends on the problem you are solving. So, Plz discribe your problem. Note: OpenFOAM mesh making module has two major blocks. 1. blockMesh // this generate basic structured mesh for domain to be considered. 2. snappyHexMesh // This generates fine mesh for boundary ex a CAR or Cylinder. Also consider a cylinder placed in regular domain. If you want to mesh inside the cylnder domain Use following Patch in snappyHexMeshDict::castellatedMeshControls refinementSurfaces { myCylinder // user defined cylinder { level (2 2); faceZone cylinder; cellZone cylinder; zoneInside true; // inside or ZoneOutside for outer mesh } } . Also it is best you state your model. It will be much easier for us. regards CFDkid |
Hi,
blockmesh and snappyhexmesh are the two meshing programs that are included into OpemFoam. If you install Engrid you have another opensource mesher. With pyFoamDisplayBlockMesh.py an external program (you could visualize your blockmesh) => http://openfoamwiki.net/index.php/Co...m#Installation There are commercial mesher that can export directly as OF (ANSA, Hypermesh, Pointwise, CastNet) |
Hi,
has anyone tried to run script provided by Rasmus with the OF 2.0.1? I get an error when I try to run the blockMesh on blockMeshDict created after execution of the m4 script. blockMeshDict can be seen below. This is what I get after executing Rasmus's script: /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.0.1 | | \\ / A nd | Web: http://www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object blockMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // convertToMeters 0.001; //5 mm column diameter //10 cm length //Width of middle square section //how many cells in the square section //how many cells from square section to perimeter // how many cells from top to bottom vertices ( ( 1.25 0.0 1.25) // Vertex fiveoclocksqb = 0 (-1.25 0.0 1.25) // Vertex sevenoclocksqb = 1 (-1.25 0.0 -1.25) // Vertex elevenoclocksqb = 2 ( 1.25 0.0 -1.25) // Vertex oneoclocksqb = 3 ( 1.76776695455285 0.0 1.76776695137989) // Vertex fiveoclockcb = 4 (-1.76776695455285 0.0 1.76776695137989) // Vertex sevenoclockcb = 5 (-1.76776695455285 0.0 -1.76776695137989) // Vertex elevenoclockcb = 6 ( 1.76776695455285 0.0 -1.76776695137989) // Vertex oneoclockcb = 7 ( 1.25 100 1.25) // Vertex fiveoclocksqt = 8 (-1.25 100 1.25) // Vertex sevenoclocksqt = 9 (-1.25 100 -1.25) // Vertex elevenoclocksqt = 10 ( 1.25 100 -1.25) // Vertex oneoclocksqt = 11 ( 1.76776695455285 100 1.76776695137989) // Vertex fiveoclockct = 12 (-1.76776695455285 100 1.76776695137989) // Vertex sevenoclockct = 13 (-1.76776695455285 100 -1.76776695137989) // Vertex elevenoclockct = 14 ( 1.76776695455285 100 -1.76776695137989) // Vertex oneoclockct = 15 ); blocks ( //square block hex ( 1 0 3 2 9 8 11 10 ) (3 3 10) simpleGrading (1 1 1) //slice1 hex ( 5 4 0 1 13 12 8 9 ) (3 3 10) simpleGrading (1 1 1) //slice2 hex ( 1 2 6 5 9 10 14 13 ) (3 3 10) simpleGrading (1 1 1) //slice3 hex ( 2 3 7 6 10 11 15 14 ) (3 3 10) simpleGrading (1 1 1) //slice4 hex ( 3 0 4 7 11 8 12 15 ) (3 3 10) simpleGrading (1 1 1) ); //create the quarter circles edges ( arc 4 5 (0.0 0.0 2.5) arc 5 6 (-2.5 0.0 0.0) arc 6 7 (0.0 0.0 -2.5) arc 7 4 (2.5 0.0 0.0) arc 12 13 (0.0 100 2.5) arc 13 14 (-2.5 100 0.0) arc 14 15 (0.0 100 -2.5) arc 15 12 (2.5 100 0.0) ); patches ( patch outlet ( (2 4(0 3 2 1)) (2 4(0 4 7 3)) (2 4(4 0 1 5)) (2 4(1 2 6 5)) (2 4(3 7 6 2)) ) patch inlet ( (2 4(8 11 10 9)) (2 4(8 12 15 11)) (2 4(12 8 9 13)) (2 4(9 10 14 13)) (2 4(11 15 14 10)) ) wall walls ( (2 4(5 4 12 13)) (2 4(5 13 14 6)) (2 4(6 14 15 7)) (2 4(7 15 12 4)) ) ); mergePatchPairs ( ); |
1 Attachment(s)
PrzemekPL,
I believe it has to do with a change in the way that patches are now defined. Rename the attached version to .m4 and use the same method that hemph described: Quote:
-Smed |
Smed,
you were right. There was an issue with the patches definition, i.e. : (2 4(fiveoclocksqb oneoclocksqb elevenoclocksqb sevenoclocksqb)) Now, everything works, thanks! |
1 Attachment(s)
Hello
Here is my small contrubution to the forum.I modified the cylindermesher.m4 file to have the Z direction passing through the middle. I wanted to have arced inlet conditions for the turbinsiting tutorial so I modified the skript somebodyposted in this thread(Thanks!:)). So maybe somebody needs this version also.Just rename the file to cylindermesher.m4 and use it(You have to have the m4 preprocessor installed).I hope this helps. My Regards Burak |
Quote:
I'm trying to do basically the same mesh from this problem but using openfoam 2.1.1 and all the examples posted aren't working. Using the mesh generated from the file posted from Burak_1984 i get the blockMesh but when i try to see it on paraFoam the following error appears: --> FOAM FATAL IO ERROR: keyword outlet is undefined in dictionary "/home/lest/Downloads/C1/0/p::boundaryField" file: /home/lest/Downloads/C1/0/p::boundaryField from line 25 to line 35. From function dictionary::subDict(const word& keyword) const in file db/dictionary/dictionary.C at line 461. FOAM exiting Can anyone help me? |
1 Attachment(s)
Quote:
Quote:
|
create cylineder mesh by collapsing one edge
2 Attachment(s)
Dear all,
First make block. Collapse one side. Use arc for one side in block. Preapre mesh for half cylivder. see fig halfcyl. create blockmesh. use feature mirrorMesh by keeping and editing file mirrorMeshDict in system. See cylinder made in fig. fullcyl. Regards Dr Sachin Borse |
Quote:
Seems that some syntax has changed... Thanks |
Error in blockMeshDict
Hi Foamers,
I am very new to OpenFoam and I have modified the codes posted by the members in this thread as follows. Code:
/*-----------------------------------------------------*\ Quote:
|
Hi,
Have you visualized your mesh? There seems to be something wrong with your mesh if your maximum aspect ratio is 0. Check that the mesh is the type you want before proceeding further. Cheers, Antimony |
Hi Antimony,
Thanks for your reply, I got it. I have corrected the script as follows. Code:
/*-----------------------------------------------------*\ |
Update for OpenFOAM V5
He everyone there have been some updates to OpenFOAM since this was last posted I updated the .m4 file for the the version 5.
Paste in a text file and see under extension .m4: /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 5.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "constant/polyMesh"; object blockMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // changecom(//)changequote([,]) define(calc, [esyscmd(perl -e 'printf ($1)')]) define(VCOUNT, 0) define(vlabel, [[// ]Vertex $1 = VCOUNT define($1, VCOUNT)define([VCOUNT], incr(VCOUNT))]) meshGenApp blockMesh; convertToMeters 0.001; define(D, 10) //10 mm column diameter define(L, 500) //50 cm length define(PI, 3.14159265) define(R, calc(D/2)) define(CW, calc(D/4)) //Width of middle square section define(CX, calc(R*cos((PI/180)*45))) define(CZ, calc(R*sin((PI/180)*45))) define(NPS, 50) // how many cells in width of square section define(NPD, 50) // how many cells from square section edge to perimeter of cylinder define(NPY, 50) // how many cells from top to bottom vertices ( ( CW 0.0 CW) vlabel(fiveoclocksqb) (-CW 0.0 CW) vlabel(sevenoclocksqb) (-CW 0.0 -CW) vlabel(elevenoclocksqb) ( CW 0.0 -CW) vlabel(oneoclocksqb) ( CX 0.0 CZ) vlabel(fiveoclockcb) (-CX 0.0 CZ) vlabel(sevenoclockcb) (-CX 0.0 -CZ) vlabel(elevenoclockcb) ( CX 0.0 -CZ) vlabel(oneoclockcb) ( CW L CW) vlabel(fiveoclocksqt) (-CW L CW) vlabel(sevenoclocksqt) (-CW L -CW) vlabel(elevenoclocksqt) ( CW L -CW) vlabel(oneoclocksqt) ( CX L CZ) vlabel(fiveoclockct) (-CX L CZ) vlabel(sevenoclockct) (-CX L -CZ) vlabel(elevenoclockct) ( CX L -CZ) vlabel(oneoclockct) ); blocks ( //square block hex ( sevenoclocksqb fiveoclocksqb oneoclocksqb elevenoclocksqb sevenoclocksqt fiveoclocksqt oneoclocksqt elevenoclocksqt ) (NPS NPS NPY) simpleGrading (1 1 1) //slice1 hex ( sevenoclockcb fiveoclockcb fiveoclocksqb sevenoclocksqb sevenoclockct fiveoclockct fiveoclocksqt sevenoclocksqt ) (NPS NPD NPY) simpleGrading (1 1 1) //slice2 hex ( sevenoclocksqb elevenoclocksqb elevenoclockcb sevenoclockcb sevenoclocksqt elevenoclocksqt elevenoclockct sevenoclockct ) (NPS NPD NPY) simpleGrading (1 1 1) //slice3 hex ( elevenoclocksqb oneoclocksqb oneoclockcb elevenoclockcb elevenoclocksqt oneoclocksqt oneoclockct elevenoclockct ) (NPS NPD NPY) simpleGrading (1 1 1) //slice4 hex ( oneoclocksqb fiveoclocksqb fiveoclockcb oneoclockcb oneoclocksqt fiveoclocksqt fiveoclockct oneoclockct ) (NPS NPD NPY) simpleGrading (1 1 1) ); //create the quarter circles edges ( arc fiveoclockcb sevenoclockcb (0.0 0.0 R) arc sevenoclockcb elevenoclockcb (-R 0.0 0.0) arc elevenoclockcb oneoclockcb (0.0 0.0 -R) arc oneoclockcb fiveoclockcb (R 0.0 0.0) arc fiveoclockct sevenoclockct (0.0 L R) arc sevenoclockct elevenoclockct (-R L 0.0) arc elevenoclockct oneoclockct (0.0 L -R) arc oneoclockct fiveoclockct (R L 0.0) ); boundary ( outlet { type patch; faces ( (fiveoclocksqb oneoclocksqb elevenoclocksqb sevenoclocksqb) (fiveoclocksqb fiveoclockcb oneoclockcb oneoclocksqb) (fiveoclockcb fiveoclocksqb sevenoclocksqb sevenoclockcb) (sevenoclocksqb elevenoclocksqb elevenoclockcb sevenoclockcb) (oneoclocksqb oneoclockcb elevenoclockcb elevenoclocksqb) ); } inlet { type patch; faces ( (fiveoclocksqt oneoclocksqt elevenoclocksqt sevenoclocksqt) (fiveoclocksqt fiveoclockct oneoclockct oneoclocksqt) (fiveoclockct fiveoclocksqt sevenoclocksqt sevenoclockct) (sevenoclocksqt elevenoclocksqt elevenoclockct sevenoclockct) (oneoclocksqt oneoclockct elevenoclockct elevenoclocksqt) ); } walls { type wall; faces ( (sevenoclockcb fiveoclockcb fiveoclockct sevenoclockct) (sevenoclockcb sevenoclockct elevenoclockct elevenoclockcb) (elevenoclockcb elevenoclockct oneoclockct oneoclockcb) (oneoclockcb oneoclockct fiveoclockct fiveoclockcb) ); } ); // ************************************************** *********************** // |
Here is a pipe mesh done with OpenFOAM-dev and the #calc function. As well as the new methods to adress patches via their block and face number. Also naming the points etc is used. The values 0.278 and 0.35 are chosen in such a way, that the non orthogonality and skewness of the resulting mesh is optimal for the current mesh size. It also features an additional outer region for increased resolution of the boundary layer.
Code:
/*--------------------------------*- C++ -*----------------------------------*\ |
1/4 of pipe with blockmesh
hello,
i am fairly new with OpenFOAM and well, the preprocessing is a little bit complicated, OF it is super powerfull as software but a little bit complicated to get the hand on it.... i wanted to do a 1/4 of a cylinder with blockmesh, is it possible to do the blockmesh with polyhedral cells for the mesh? also how can i define a simple quarter of cylinder, i found some examples in the comments of a cylinder but it was difficult to follow them.... i wanted to finish with something where i can define the R and the lenght of the hight in Z with the 5 different patches defined... |
3 Attachment(s)
Hello everyone!
I am currently working with OpenFoam to create a cylinder. The cylinder I need should only be discretized in the h direction (Photo). Does anyone know how I can avoid the round basement being divided into so many small cells? (I have attached a photo). I´m quite new to OpenFoam and would be very happy if someone could help me! Regards! |
All times are GMT -4. The time now is 10:28. |