zhyang-dev/CFDEMcoupling-PFM
1
0
Fork 0
mirror of https://github.com/ParticulateFlow/CFDEMcoupling-PFM.git synced 2025-12-08 06:37:44 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

move blockMeshDict to system folder

Browse Source 
and remove m4 file
This commit is contained in:
danielque
2022-04-22 15:32:41 +02:00
parent c7532c919c
commit 5a1c957ba5
2 changed files with 50 additions and 232 deletions
Download Patch File Download Diff File Expand all files Collapse all files

153
tutorials/cfdemSolverRhoPimpleChem/PolydisperseFluidizedBed/R2_FB/CFD/constant/polyMesh/cylinderMesh.m4
View File

@ -1,153 +0,0 @@
// blockMesh : Block mesh description file
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.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 1.0;
define(D, 0.068) //68 mm column diameter
define(L, 0.150) //150 mm 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, 7) //how many cells in the square section
define(NPD, 3) //how many cells from square section to perimeter
define(NPY, 24) // 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)
);
patches
(
patch inlet
(
(fiveoclocksqb oneoclocksqb elevenoclocksqb sevenoclocksqb)
(fiveoclocksqb fiveoclockcb oneoclockcb oneoclocksqb)
(fiveoclockcb fiveoclocksqb sevenoclocksqb sevenoclockcb)
(sevenoclocksqb elevenoclocksqb elevenoclockcb sevenoclockcb)
(oneoclocksqb oneoclockcb elevenoclockcb elevenoclocksqb)
)
patch outlet
(
(fiveoclocksqt oneoclocksqt elevenoclocksqt sevenoclocksqt)
(fiveoclocksqt fiveoclockct oneoclockct oneoclocksqt)
(fiveoclockct fiveoclocksqt sevenoclocksqt sevenoclockct)
(sevenoclocksqt elevenoclocksqt elevenoclockct sevenoclockct)
(oneoclocksqt oneoclockct elevenoclockct elevenoclocksqt)
)
wall side-walls
(
(sevenoclockcb fiveoclockcb fiveoclockct sevenoclockct)
(sevenoclockcb sevenoclockct elevenoclockct elevenoclockcb)
(elevenoclockcb elevenoclockct oneoclockct oneoclockcb)
(oneoclockcb oneoclockct fiveoclockct fiveoclockcb)
)
);

63
tutorials/cfdemSolverRhoPimpleChem/PolydisperseFluidizedBed/R2_FB/CFD/constant/polyMesh/blockMeshDict → tutorials/cfdemSolverRhoPimpleChem/PolydisperseFluidizedBed/R2_FB/CFD/system/blockMeshDict
View File

@ -1,10 +1,9 @@
// blockMesh : Block mesh description file
/*--------------------------------*- C++ -*----------------------------------*\ /*--------------------------------*- C++ -*----------------------------------*\
| ========= | | ========= |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
| \\ / O peration | Version: 2.3.0 | \\ / O peration | Website: https://openfoam.org
| \\ / A nd | Web: www.OpenFOAM.org | \\ / A nd | Version: 6
| \\/ M anipulation | | \\/ M anipulation |
\*---------------------------------------------------------------------------*/ \*---------------------------------------------------------------------------*/
FoamFile FoamFile
{ {
@ -20,17 +19,16 @@ FoamFile
convertToMeters 1.0; convertToMeters 1.0;
// 68 mm column diameter // 68 mm column diameter
//150 mm length // 150 mm column length
// Width of middle square section 17 mm
//Width of middle square section // 7 cells in the square section
// 3 cells from square section to perimeter
//how many cells in the square section // 24 cells from top to bottom
//how many cells from square section to perimeter
// how many cells from top to bottom
vertices vertices
( (
// bottom
( 0.017 0.0 0.017) // Vertex fiveoclocksqb = 0 ( 0.017 0.0 0.017) // Vertex fiveoclocksqb = 0
(-0.017 0.0 0.017) // Vertex sevenoclocksqb = 1 (-0.017 0.0 0.017) // Vertex sevenoclocksqb = 1
(-0.017 0.0 -0.017) // Vertex elevenoclocksqb = 2 (-0.017 0.0 -0.017) // Vertex elevenoclocksqb = 2
@ -41,6 +39,7 @@ FoamFile
(-0.0240416305819187 0.0 -0.0240416305387665) // Vertex elevenoclockcb = 6 (-0.0240416305819187 0.0 -0.0240416305387665) // Vertex elevenoclockcb = 6
( 0.0240416305819187 0.0 -0.0240416305387665) // Vertex oneoclockcb = 7 ( 0.0240416305819187 0.0 -0.0240416305387665) // Vertex oneoclockcb = 7
// top
( 0.017 0.150 0.017) // Vertex fiveoclocksqt = 8 ( 0.017 0.150 0.017) // Vertex fiveoclocksqt = 8
(-0.017 0.150 0.017) // Vertex sevenoclocksqt = 9 (-0.017 0.150 0.017) // Vertex sevenoclocksqt = 9
(-0.017 0.150 -0.017) // Vertex elevenoclocksqt = 10 (-0.017 0.150 -0.017) // Vertex elevenoclocksqt = 10
@ -55,45 +54,19 @@ FoamFile
blocks blocks
( (
// square block // square block
hex ( hex ( 1 0 3 2 9 8 11 10) (7 7 24) simpleGrading (1 1 1)
1 0 3 2
9 8 11 10
)
(7 7 24)
simpleGrading (1 1 1)
// slice1 // slice1
hex ( hex ( 5 4 0 1 13 12 8 9) (7 3 24) simpleGrading (1 1 1)
5 4 0 1
13 12 8 9
)
(7 3 24)
simpleGrading (1 1 1)
//slice2 //slice2
hex ( hex ( 1 2 6 5 9 10 14 13) (7 3 24) simpleGrading (1 1 1)
1 2 6 5
9 10 14 13
)
(7 3 24)
simpleGrading (1 1 1)
// slice3 // slice3
hex ( hex ( 2 3 7 6 10 11 15 14) (7 3 24) simpleGrading (1 1 1)
2 3 7 6
10 11 15 14
)
(7 3 24)
simpleGrading (1 1 1)
//slice4 //slice4
hex ( hex ( 3 0 4 7 11 8 12 15) (7 3 24) simpleGrading (1 1 1)
3 0 4 7
11 8 12 15
)
(7 3 24)
simpleGrading (1 1 1)
); );
@ -109,7 +82,6 @@ FoamFile
arc 13 14 (-0.034 0.150 0.0) arc 13 14 (-0.034 0.150 0.0)
arc 14 15 (0.0 0.150 -0.034) arc 14 15 (0.0 0.150 -0.034)
arc 15 12 (0.034 0.150 0.0) arc 15 12 (0.034 0.150 0.0)
); );
patches patches
@ -139,6 +111,5 @@ FoamFile
(6 14 15 7) (6 14 15 7)
(7 15 12 4) (7 15 12 4)
) )
); );