Last active
August 23, 2018 10:30
-
-
Save Foadsf/22edbd7a7bdefefc24005753c8f89d40 to your computer and use it in GitHub Desktop.
axisymmetric step OpenFOAM
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
//following this: https://openfoamwiki.net/index.php/Main_ContribExamples/AxiSymmetric | |
// /cygdrive/c/Users/sojoodif/surfdrive/sync/MyPhD/simulation/OpenFOAM/20180131_2 | |
/*--------------------------------*- C++ -*----------------------------------*\ | |
| ========= | | | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | |
| \\ / O peration | Version: 5 | | |
| \\ / A nd | Web: www.OpenFOAM.org | | |
| \\/ M anipulation | | | |
\*---------------------------------------------------------------------------*/ | |
FoamFile | |
{ | |
version 2.0; | |
format ascii; | |
root ""; | |
case ""; | |
instance ""; | |
local ""; | |
class dictionary; | |
object blockMeshDict; | |
} | |
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // | |
convertToMeters 0.001; | |
wa 5.0*constant::mathematical::pi/180; // wedge angle must be smaller than 5 degrees | |
cr 10.0; // cylinder radius | |
pr 8.0; //piston radius | |
cl 50.0; // length of cylinder | |
px 20.0; // position of piston | |
pl 10.0; // length of piston | |
px2 #calc "$px+$pl"; // end of the piston | |
ms 50; //meshsize | |
vertices | |
( | |
(0 0 0) //00 | |
(#calc "$pr*cos($wa/2)" #calc "$pr*sin($wa/2)" 0) //01 | |
(#calc "$pr*cos($wa/2)" #calc "-$pr*sin($wa/2)" 0) //02 | |
(#calc "$pr*cos($wa/2)" #calc "-$pr*sin($wa/2)" $px) //03 | |
(#calc "$pr*cos($wa/2)" #calc "$pr*sin($wa/2)" $px) //04 | |
(0 0 $px) //05 | |
(#calc "$cr*cos($wa/2)" #calc "$cr*sin($wa/2)" $px) //06 | |
(#calc "$cr*cos($wa/2)" #calc "$cr*sin($wa/2)" 0) //07 | |
(#calc "$cr*cos($wa/2)" #calc "-$cr*sin($wa/2)" 0) //08 | |
(#calc "$cr*cos($wa/2)" #calc "-$cr*sin($wa/2)" $px) //09 | |
(#calc "$cr*cos($wa/2)" #calc "-$cr*sin($wa/2)" $px2) //10 | |
(#calc "$cr*cos($wa/2)" #calc "$cr*sin($wa/2)" $px2) //11 | |
(#calc "$pr*cos($wa/2)" #calc "$pr*sin($wa/2)" $px2) //12 | |
(#calc "$pr*cos($wa/2)" #calc "-$pr*sin($wa/2)" $px2) //13 | |
(#calc "$pr*cos($wa/2)" #calc "-$pr*sin($wa/2)" $cl) //14 | |
(#calc "$pr*cos($wa/2)" #calc "$pr*sin($wa/2)" $cl) //15 | |
(#calc "$cr*cos($wa/2)" #calc "$cr*sin($wa/2)" $cl) //16 | |
(#calc "$cr*cos($wa/2)" #calc "-$cr*sin($wa/2)" $cl) //17 | |
(0 0 $cl) //18 | |
(0 0 $px2) //19 | |
); | |
blocks | |
( | |
hex (0 2 1 0 5 3 4 5) ($ms 1 $ms) simpleGrading (1 1 1) | |
hex (1 2 8 7 4 3 9 6) (1 $ms $ms) simpleGrading (1 1 1) | |
hex (4 3 9 6 12 13 10 11) (1 $ms $ms) simpleGrading (1 1 1) | |
hex (12 13 10 11 15 14 17 16) (1 $ms $ms) simpleGrading (1 1 1) | |
hex (19 13 12 19 18 14 15 18) ($ms 1 $ms) simpleGrading (1 1 1) | |
); | |
edges | |
( | |
); | |
boundary | |
( | |
inlet | |
{ | |
type patch; | |
faces | |
( | |
(0 2 1 0) | |
(1 2 8 7) | |
); | |
} | |
outlet | |
{ | |
type patch; | |
faces | |
( | |
(18 15 14 18) | |
(15 16 17 14) | |
); | |
} | |
walls | |
{ | |
type wall; | |
faces | |
( | |
(7 8 9 6) | |
(6 9 10 11) | |
(11 10 17 16) | |
(5 4 3 5) | |
(12 13 3 4) | |
(19 13 12 19) | |
); | |
} | |
wedgeFront | |
{ | |
type wedge; | |
faces | |
( | |
(0 1 4 5) | |
(1 7 6 4) | |
(4 6 11 12) | |
(12 11 16 15) | |
(19 12 15 18) | |
); | |
} | |
wedgeBack | |
{ | |
type wedge; | |
faces | |
( | |
(0 2 3 5) | |
(2 8 9 3) | |
(3 9 10 13) | |
(13 10 17 14) | |
(19 13 14 18) | |
); | |
} | |
axis | |
{ | |
type empty; | |
faces | |
( | |
(0 5 5 0) | |
(19 18 18 19) | |
); | |
} | |
); | |
mergePatchPairs | |
( | |
); | |
// ************************************************************************* // |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment