diff --git a/tutorials/incompressible/MRFSimpleFoam/axialTurbine_ggi/Allrun b/tutorials/incompressible/MRFSimpleFoam/axialTurbine_ggi/Allrun index a420c6658..b3fb98fd0 100755 --- a/tutorials/incompressible/MRFSimpleFoam/axialTurbine_ggi/Allrun +++ b/tutorials/incompressible/MRFSimpleFoam/axialTurbine_ggi/Allrun @@ -3,16 +3,17 @@ # Source tutorial run functions . $WM_PROJECT_DIR/bin/tools/RunFunctions -application="MRFSimpleFoam" +application=`getApplication` #Create the mesh: m4 < constant/polyMesh/blockMeshDict.m4 > constant/polyMesh/blockMeshDict -blockMesh +runApplication blockMesh +transformPoints -scale "(1 20 1)" transformPoints -cylToCart "((0 0 0) (0 0 1) (1 0 0))" # Set 0-directory and create GGI set: cp -r 0_orig 0 -setSet -batch setBatchGgi -setsToZones -noFlipMap +runApplication setSet -batch setBatchGgi +runApplication setsToZones -noFlipMap runApplication $application diff --git a/tutorials/incompressible/MRFSimpleFoam/axialTurbine_ggi/constant/polyMesh/blockMeshDict.m4 b/tutorials/incompressible/MRFSimpleFoam/axialTurbine_ggi/constant/polyMesh/blockMeshDict.m4 index e6d01b412..af602bea2 100644 --- a/tutorials/incompressible/MRFSimpleFoam/axialTurbine_ggi/constant/polyMesh/blockMeshDict.m4 +++ b/tutorials/incompressible/MRFSimpleFoam/axialTurbine_ggi/constant/polyMesh/blockMeshDict.m4 @@ -20,7 +20,7 @@ define(calc, [esyscmd(perl -e 'printf ($1)')]) //define(calc, [esyscmd(echo $1 | bc | tr -d \\n)]) define(VCOUNT, 0) define(vlabel, [[// ]Vertex $1 = VCOUNT define($1, VCOUNT)define([VCOUNT], incr(VCOUNT))]) -define(pi, 3.14159265) +define(pi, calc(3.14159265/20)) define(hex2D, hex ($1b $2b $3b $4b $1t $2t $3t $4t)) define(quad2D, ($1b $2b $2t $1t)) @@ -200,7 +200,7 @@ blocks hex2D(GV1l, GV1r, GV2r, GV2l) (GVtc GVbac GVrc) - simpleGrading (1 0.2 1) + simpleGrading (1 1 1) hex2D(GV2l, GV2r, GV3r, GV3l) (GVtc GViac GVrc) @@ -215,7 +215,7 @@ blocks hex2D(RU1l, RU1r, RU2r, RU2l) rotor (RUtc RUbac RUrc) - simpleGrading (1 0.4 1) + simpleGrading (1 1 1) hex2D(RU2l, RU2r, RU3r, RU3l) rotor diff --git a/tutorials/incompressible/MRFSimpleFoam/axialTurbine_mixingPlane/Allrun b/tutorials/incompressible/MRFSimpleFoam/axialTurbine_mixingPlane/Allrun index a420c6658..b3fb98fd0 100755 --- a/tutorials/incompressible/MRFSimpleFoam/axialTurbine_mixingPlane/Allrun +++ b/tutorials/incompressible/MRFSimpleFoam/axialTurbine_mixingPlane/Allrun @@ -3,16 +3,17 @@ # Source tutorial run functions . $WM_PROJECT_DIR/bin/tools/RunFunctions -application="MRFSimpleFoam" +application=`getApplication` #Create the mesh: m4 < constant/polyMesh/blockMeshDict.m4 > constant/polyMesh/blockMeshDict -blockMesh +runApplication blockMesh +transformPoints -scale "(1 20 1)" transformPoints -cylToCart "((0 0 0) (0 0 1) (1 0 0))" # Set 0-directory and create GGI set: cp -r 0_orig 0 -setSet -batch setBatchGgi -setsToZones -noFlipMap +runApplication setSet -batch setBatchGgi +runApplication setsToZones -noFlipMap runApplication $application diff --git a/tutorials/incompressible/MRFSimpleFoam/axialTurbine_mixingPlane/constant/polyMesh/blockMeshDict.m4 b/tutorials/incompressible/MRFSimpleFoam/axialTurbine_mixingPlane/constant/polyMesh/blockMeshDict.m4 index e91a53476..b5d676b56 100644 --- a/tutorials/incompressible/MRFSimpleFoam/axialTurbine_mixingPlane/constant/polyMesh/blockMeshDict.m4 +++ b/tutorials/incompressible/MRFSimpleFoam/axialTurbine_mixingPlane/constant/polyMesh/blockMeshDict.m4 @@ -20,7 +20,7 @@ define(calc, [esyscmd(perl -e 'printf ($1)')]) //define(calc, [esyscmd(echo $1 | bc | tr -d \\n)]) define(VCOUNT, 0) define(vlabel, [[// ]Vertex $1 = VCOUNT define($1, VCOUNT)define([VCOUNT], incr(VCOUNT))]) -define(pi, 3.14159265) +define(pi, calc(3.14159265/20)) define(hex2D, hex ($1b $2b $3b $4b $1t $2t $3t $4t)) define(quad2D, ($1b $2b $2t $1t)) @@ -200,7 +200,7 @@ blocks hex2D(GV1l, GV1r, GV2r, GV2l) (GVtc GVbac GVrc) - simpleGrading (1 0.2 1) + simpleGrading (1 1 1) hex2D(GV2l, GV2r, GV3r, GV3l) (GVtc GViac GVrc) @@ -215,7 +215,7 @@ blocks hex2D(RU1l, RU1r, RU2r, RU2l) rotor (RUtc RUbac RUrc) - simpleGrading (1 0.4 1) + simpleGrading (1 1 1) hex2D(RU2l, RU2r, RU3r, RU3l) rotor diff --git a/tutorials/incompressible/pimpleDyMFoam/axialTurbine/0_orig/nut b/tutorials/incompressible/pimpleDyMFoam/axialTurbine/0_orig/nut new file mode 100644 index 000000000..68e3586ca --- /dev/null +++ b/tutorials/incompressible/pimpleDyMFoam/axialTurbine/0_orig/nut @@ -0,0 +1,151 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | foam-extend: Open Source CFD | +| \\ / O peration | Version: 3.1 | +| \\ / A nd | Web: http://www.extend-project.de | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object nut; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 2 -1 0 0 0 0]; + +internalField uniform 0; + +boundaryField +{ + GVINLET + { + type calculated; + value uniform 0; + } + GVOUTLET + { + type overlapGgi; + value uniform 0; + } + GVCYCLIC1 + { + type cyclicGgi; + value uniform 0; + } + GVCYCLIC2 + { + type cyclicGgi; + value uniform 0; + } + GVBLADE + { + type nutWallFunction; + Cmu 0.09; + kappa 0.41; + E 9.8; + value uniform 0; + } + GVHUB + { + type nutWallFunction; + Cmu 0.09; + kappa 0.41; + E 9.8; + value uniform 0; + } + GVSHROUD + { + type nutWallFunction; + Cmu 0.09; + kappa 0.41; + E 9.8; + value uniform 0; + } + RUINLET + { + type overlapGgi; + value uniform 0; + } + RUOUTLET + { + type overlapGgi; + value uniform 0; + } + RUCYCLIC1 + { + type cyclicGgi; + value uniform 0; + } + RUCYCLIC2 + { + type cyclicGgi; + value uniform 0; + } + RUBLADE + { + type nutWallFunction; + Cmu 0.09; + kappa 0.41; + E 9.8; + value uniform 0; + } + RUHUB + { + type nutWallFunction; + Cmu 0.09; + kappa 0.41; + E 9.8; + value uniform 0; + } + RUSHROUD + { + type nutWallFunction; + Cmu 0.09; + kappa 0.41; + E 9.8; + value uniform 0; + } + DTINLET + { + type overlapGgi; + value uniform 0; + } + DTOUTLET + { + type calculated; + value uniform 0; + } + DTCYCLIC1 + { + type cyclicGgi; + value uniform 0; + } + DTCYCLIC2 + { + type cyclicGgi; + value uniform 0; + } + DTHUB + { + type nutWallFunction; + Cmu 0.09; + kappa 0.41; + E 9.8; + value uniform 0; + } + DTSHROUD + { + type nutWallFunction; + Cmu 0.09; + kappa 0.41; + E 9.8; + value uniform 0; + } +} + + +// ************************************************************************* // diff --git a/tutorials/incompressible/pimpleDyMFoam/axialTurbine/Allrun b/tutorials/incompressible/pimpleDyMFoam/axialTurbine/Allrun index 348410f97..b3fb98fd0 100755 --- a/tutorials/incompressible/pimpleDyMFoam/axialTurbine/Allrun +++ b/tutorials/incompressible/pimpleDyMFoam/axialTurbine/Allrun @@ -3,16 +3,17 @@ # Source tutorial run functions . $WM_PROJECT_DIR/bin/tools/RunFunctions -application="pimpleDyMFoam" +application=`getApplication` #Create the mesh: m4 < constant/polyMesh/blockMeshDict.m4 > constant/polyMesh/blockMeshDict -blockMesh +runApplication blockMesh +transformPoints -scale "(1 20 1)" transformPoints -cylToCart "((0 0 0) (0 0 1) (1 0 0))" # Set 0-directory and create GGI set: cp -r 0_orig 0 -setSet -batch setBatchGgi -setsToZones -noFlipMap +runApplication setSet -batch setBatchGgi +runApplication setsToZones -noFlipMap runApplication $application diff --git a/tutorials/incompressible/pimpleDyMFoam/axialTurbine/constant/polyMesh/blockMeshDict.m4 b/tutorials/incompressible/pimpleDyMFoam/axialTurbine/constant/polyMesh/blockMeshDict.m4 index 517a737f9..9169949e7 100644 --- a/tutorials/incompressible/pimpleDyMFoam/axialTurbine/constant/polyMesh/blockMeshDict.m4 +++ b/tutorials/incompressible/pimpleDyMFoam/axialTurbine/constant/polyMesh/blockMeshDict.m4 @@ -20,7 +20,7 @@ define(calc, [esyscmd(perl -e 'printf ($1)')]) //define(calc, [esyscmd(echo $1 | bc | tr -d \\n)]) define(VCOUNT, 0) define(vlabel, [[// ]Vertex $1 = VCOUNT define($1, VCOUNT)define([VCOUNT], incr(VCOUNT))]) -define(pi, 3.14159265) +define(pi, calc(3.14159265/20)) define(hex2D, hex ($1b $2b $3b $4b $1t $2t $3t $4t)) define(quad2D, ($1b $2b $2t $1t)) @@ -200,7 +200,7 @@ blocks hex2D(GV1l, GV1r, GV2r, GV2l) (GVtc GVbac GVrc) - simpleGrading (1 0.2 1) + simpleGrading (1 1 1) hex2D(GV2l, GV2r, GV3r, GV3l) (GVtc GViac GVrc) @@ -215,7 +215,7 @@ blocks hex2D(RU1l, RU1r, RU2r, RU2l) rotor (RUtc RUbac RUrc) - simpleGrading (1 0.4 1) + simpleGrading (1 1 1) hex2D(RU2l, RU2r, RU3r, RU3l) rotor diff --git a/tutorials/incompressible/pimpleDyMFoam/axialTurbine/system/controlDict b/tutorials/incompressible/pimpleDyMFoam/axialTurbine/system/controlDict index f61c340dc..1ca2048f4 100644 --- a/tutorials/incompressible/pimpleDyMFoam/axialTurbine/system/controlDict +++ b/tutorials/incompressible/pimpleDyMFoam/axialTurbine/system/controlDict @@ -14,7 +14,7 @@ FoamFile } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -application simpleFoam; +application pimpleDyMFoam; startFrom startTime; diff --git a/tutorials/incompressible/simpleSRFFoam/axialTurbine/0_orig/Uabs b/tutorials/incompressible/simpleSRFFoam/axialTurbine/0_orig/Uabs new file mode 100644 index 000000000..65950e5dc --- /dev/null +++ b/tutorials/incompressible/simpleSRFFoam/axialTurbine/0_orig/Uabs @@ -0,0 +1,62 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | foam-extend: Open Source CFD | +| \\ / O peration | Version: 3.1 | +| \\ / A nd | Web: http://www.extend-project.de | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volVectorField; + location "0"; + object Uabs; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 0); + +boundaryField +{ + RUINLET + { + type calculated; + value uniform (0 0 0); + } + RUOUTLET + { + type calculated; + value uniform (0 0 0); + } + RUCYCLIC1 + { + type cyclicGgi; + value uniform (0 0 0); + } + RUCYCLIC2 + { + type cyclicGgi; + value uniform (0 0 0); + } + RUBLADE + { + type calculated; + value uniform (0 0 0); + } + RUHUB + { + type calculated; + value uniform (0 0 0); + } + RUSHROUD + { + type calculated; + value uniform (0 0 0); + } +} + + +// ************************************************************************* // diff --git a/tutorials/incompressible/simpleSRFFoam/axialTurbine/0_orig/Urel b/tutorials/incompressible/simpleSRFFoam/axialTurbine/0_orig/Urel new file mode 100644 index 000000000..0034d6bf4 --- /dev/null +++ b/tutorials/incompressible/simpleSRFFoam/axialTurbine/0_orig/Urel @@ -0,0 +1,61 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | foam-extend: Open Source CFD | +| \\ / O peration | Version: 3.1 | +| \\ / A nd | Web: http://www.extend-project.de | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volVectorField; + object Urel; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 -1); + +boundaryField +{ + RUINLET + { + type SRFVelocity; + inletValue uniform (0 0 -1); + relative no; + value uniform (0 0 0); + } + RUOUTLET + { + type zeroGradient; + } + RUCYCLIC1 + { + type cyclicGgi; + } + RUCYCLIC2 + { + type cyclicGgi; + } + RUBLADE + { + type fixedValue; + value uniform (0 0 0); + } + RUHUB + { + type fixedValue; + value uniform (0 0 0); + } + RUSHROUD + { + type SRFVelocity; + inletValue uniform (0 0 0); + relative yes; + value uniform (0 0 0); + } +} + +// ************************************************************************* // diff --git a/tutorials/incompressible/simpleSRFFoam/axialTurbine/0_orig/epsilon b/tutorials/incompressible/simpleSRFFoam/axialTurbine/0_orig/epsilon new file mode 100644 index 000000000..3352007f5 --- /dev/null +++ b/tutorials/incompressible/simpleSRFFoam/axialTurbine/0_orig/epsilon @@ -0,0 +1,75 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | foam-extend: Open Source CFD | +| \\ / O peration | Version: 3.1 | +| \\ / A nd | Web: http://www.extend-project.de | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object epsilon; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 2 -3 0 0 0 0]; + +internalField uniform 14.855; + +boundaryField +{ + RUINLET + { + type fixedValue; + value uniform 14.855; + } + RUOUTLET + { + type zeroGradient; + } + RUCYCLIC1 + { + type cyclicGgi; + value uniform 14.855; + } + RUCYCLIC2 + { + type cyclicGgi; + value uniform 14.855; + } + RUBLADE + { + type epsilonWallFunction; + U Urel; + refValue uniform 0; + value uniform 14.855; + Cmu 0.09; + kappa 0.41; + E 9.8; + } + RUHUB + { + type epsilonWallFunction; + U Urel; + refValue uniform 0; + value uniform 14.855; + Cmu 0.09; + kappa 0.41; + E 9.8; + } + RUSHROUD + { + type epsilonWallFunction; + U Urel; + refValue uniform 0; + value uniform 14.855; + Cmu 0.09; + kappa 0.41; + E 9.8; + } +} + +// ************************************************************************* // diff --git a/tutorials/incompressible/simpleSRFFoam/axialTurbine/0_orig/k b/tutorials/incompressible/simpleSRFFoam/axialTurbine/0_orig/k new file mode 100644 index 000000000..935906287 --- /dev/null +++ b/tutorials/incompressible/simpleSRFFoam/axialTurbine/0_orig/k @@ -0,0 +1,60 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | foam-extend: Open Source CFD | +| \\ / O peration | Version: 3.1 | +| \\ / A nd | Web: http://www.extend-project.de | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object k; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 2 -2 0 0 0 0]; + +internalField uniform 0.375; + +boundaryField +{ + RUINLET + { + type fixedValue; + value uniform 0.375; + } + RUOUTLET + { + type zeroGradient; + } + RUCYCLIC1 + { + type cyclicGgi; + value uniform 0.375; + } + RUCYCLIC2 + { + type cyclicGgi; + value uniform 0.375; + } + RUBLADE + { + type kqRWallFunction; + value uniform 0.375; + } + RUHUB + { + type kqRWallFunction; + value uniform 0.375; + } + RUSHROUD + { + type kqRWallFunction; + value uniform 0.375; + } +} + +// ************************************************************************* // diff --git a/tutorials/incompressible/simpleSRFFoam/axialTurbine/0_orig/nut b/tutorials/incompressible/simpleSRFFoam/axialTurbine/0_orig/nut new file mode 100644 index 000000000..a0ab9f0a7 --- /dev/null +++ b/tutorials/incompressible/simpleSRFFoam/axialTurbine/0_orig/nut @@ -0,0 +1,71 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | foam-extend: Open Source CFD | +| \\ / O peration | Version: 3.1 | +| \\ / A nd | Web: http://www.extend-project.de | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object nut; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 2 -1 0 0 0 0]; + +internalField uniform 0; + +boundaryField +{ + RUINLET + { + type calculated; + value uniform 0; + } + RUOUTLET + { + type calculated; + value uniform 0; + } + RUCYCLIC1 + { + type cyclicGgi; + value uniform 0; + } + RUCYCLIC2 + { + type cyclicGgi; + value uniform 0; + } + RUBLADE + { + type nutWallFunction; + Cmu 0.09; + kappa 0.41; + E 9.8; + value uniform 0; + } + RUHUB + { + type nutWallFunction; + Cmu 0.09; + kappa 0.41; + E 9.8; + value uniform 0; + } + RUSHROUD + { + type nutWallFunction; + Cmu 0.09; + kappa 0.41; + E 9.8; + value uniform 0; + } +} + + +// ************************************************************************* // diff --git a/tutorials/incompressible/simpleSRFFoam/axialTurbine/0_orig/p b/tutorials/incompressible/simpleSRFFoam/axialTurbine/0_orig/p new file mode 100644 index 000000000..150f9f242 --- /dev/null +++ b/tutorials/incompressible/simpleSRFFoam/axialTurbine/0_orig/p @@ -0,0 +1,54 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | foam-extend: Open Source CFD | +| \\ / O peration | Version: 3.1 | +| \\ / A nd | Web: http://www.extend-project.de | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + object p; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 2 -2 0 0 0 0]; + +internalField uniform 0; + +boundaryField +{ + RUINLET + { + type zeroGradient; + } + RUOUTLET + { + type fixedValue; + value uniform 0; + } + RUCYCLIC1 + { + type cyclicGgi; + } + RUCYCLIC2 + { + type cyclicGgi; + } + RUBLADE + { + type zeroGradient; + } + RUHUB + { + type zeroGradient; + } + RUSHROUD + { + type zeroGradient; + } +} + +// ************************************************************************* // diff --git a/tutorials/incompressible/simpleSRFFoam/axialTurbine/Allclean b/tutorials/incompressible/simpleSRFFoam/axialTurbine/Allclean new file mode 100755 index 000000000..c8c3c4e26 --- /dev/null +++ b/tutorials/incompressible/simpleSRFFoam/axialTurbine/Allclean @@ -0,0 +1,4 @@ +#!/bin/sh + +( cd constant/polyMesh && \rm -rf boundary* blockMeshDict *Zones* faces* neighbour* owner* points* sets ) +\rm -rf 0 [1-9]* VTK log* processor* diff --git a/tutorials/incompressible/simpleSRFFoam/axialTurbine/Allrun b/tutorials/incompressible/simpleSRFFoam/axialTurbine/Allrun new file mode 100755 index 000000000..b2d76583e --- /dev/null +++ b/tutorials/incompressible/simpleSRFFoam/axialTurbine/Allrun @@ -0,0 +1,21 @@ +#!/bin/bash +# Source tutorial run functions +. $WM_PROJECT_DIR/bin/tools/RunFunctions + +# Get application from system/controlDict +application=`getApplication` + +# make sure the application is compiled +#wmake ../$application + +m4 < constant/polyMesh/blockMeshDict.m4 > constant/polyMesh/blockMeshDict +runApplication blockMesh +transformPoints -scale "(1 20 1)" +transformPoints -cylToCart "((0 0 0) (0 0 1) (1 0 0))" + +# Set 0-directory and create GGI set: +cp -r 0_orig 0 +runApplication setSet -batch setBatchGgi +runApplication setsToZones -noFlipMap + +runApplication $application diff --git a/tutorials/incompressible/simpleSRFFoam/axialTurbine/allBlades.pvsm b/tutorials/incompressible/simpleSRFFoam/axialTurbine/allBlades.pvsm new file mode 100644 index 000000000..41675781b --- /dev/null +++ b/tutorials/incompressible/simpleSRFFoam/axialTurbine/allBlades.pvsm @@ -0,0 +1,8491 @@ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/tutorials/incompressible/simpleSRFFoam/axialTurbine/constant/RASProperties b/tutorials/incompressible/simpleSRFFoam/axialTurbine/constant/RASProperties new file mode 100644 index 000000000..b64f08334 --- /dev/null +++ b/tutorials/incompressible/simpleSRFFoam/axialTurbine/constant/RASProperties @@ -0,0 +1,191 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | foam-extend: Open Source CFD | +| \\ / O peration | Version: 3.1 | +| \\ / A nd | Web: http://www.extend-project.de | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object RASProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +RASModel RNGkEpsilon; + +turbulence on; + +printCoeffs on; + +laminarCoeffs +{ +} + +kEpsilonCoeffs +{ + Cmu 0.09; + C1 1.44; + C2 1.92; + alphaEps 0.76923; +} + +RNGkEpsilonCoeffs +{ + Cmu 0.0845; + C1 1.42; + C2 1.68; + alphak 1.39; + alphaEps 1.39; + eta0 4.38; + beta 0.012; +} + +kOmegaSSTCoeffs +{ + alphaK1 0.85034; + alphaK2 1.0; + alphaOmega1 0.5; + alphaOmega2 0.85616; + gamma1 0.5532; + gamma2 0.4403; + beta1 0.0750; + beta2 0.0828; + betaStar 0.09; + a1 0.31; + c1 10; + + Cmu 0.09; +} + +NonlinearKEShihCoeffs +{ + Cmu 0.09; + C1 1.44; + C2 1.92; + alphak 1; + alphaEps 0.76932; + A1 1.25; + A2 1000; + Ctau1 -4; + Ctau2 13; + Ctau3 -2; + alphaKsi 0.9; +} + +LienCubicKECoeffs +{ + C1 1.44; + C2 1.92; + alphak 1; + alphaEps 0.76923; + A1 1.25; + A2 1000; + Ctau1 -4; + Ctau2 13; + Ctau3 -2; + alphaKsi 0.9; +} + +QZetaCoeffs +{ + Cmu 0.09; + C1 1.44; + C2 1.92; + alphaZeta 0.76923; + anisotropic no; +} + +LaunderSharmaKECoeffs +{ + Cmu 0.09; + C1 1.44; + C2 1.92; + alphaEps 0.76923; +} + +LamBremhorstKECoeffs +{ + Cmu 0.09; + C1 1.44; + C2 1.92; + alphaEps 0.76923; +} + +LienCubicKELowReCoeffs +{ + Cmu 0.09; + C1 1.44; + C2 1.92; + alphak 1; + alphaEps 0.76923; + A1 1.25; + A2 1000; + Ctau1 -4; + Ctau2 13; + Ctau3 -2; + alphaKsi 0.9; + Am 0.016; + Aepsilon 0.263; + Amu 0.00222; +} + +LienLeschzinerLowReCoeffs +{ + Cmu 0.09; + C1 1.44; + C2 1.92; + alphak 1; + alphaEps 0.76923; + Am 0.016; + Aepsilon 0.263; + Amu 0.00222; +} + +LRRCoeffs +{ + Cmu 0.09; + Clrr1 1.8; + Clrr2 0.6; + C1 1.44; + C2 1.92; + Cs 0.25; + Ceps 0.15; + alphaEps 0.76923; +} + +LaunderGibsonRSTMCoeffs +{ + Cmu 0.09; + Clg1 1.8; + Clg2 0.6; + C1 1.44; + C2 1.92; + C1Ref 0.5; + C2Ref 0.3; + Cs 0.25; + Ceps 0.15; + alphaEps 0.76923; + alphaR 1.22; +} + +SpalartAllmarasCoeffs +{ + alphaNut 1.5; + Cb1 0.1355; + Cb2 0.622; + Cw2 0.3; + Cw3 2; + Cv1 7.1; + Cv2 5.0; +} + +wallFunctionCoeffs +{ + kappa 0.4187; + E 9; +} + +// ************************************************************************* // diff --git a/tutorials/incompressible/simpleSRFFoam/axialTurbine/constant/SRFProperties b/tutorials/incompressible/simpleSRFFoam/axialTurbine/constant/SRFProperties new file mode 100644 index 000000000..f4fa26190 --- /dev/null +++ b/tutorials/incompressible/simpleSRFFoam/axialTurbine/constant/SRFProperties @@ -0,0 +1,26 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | foam-extend: Open Source CFD | +| \\ / O peration | Version: 3.1 | +| \\ / A nd | Web: http://www.extend-project.de | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object SRFProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +SRFModel rpm; + +axis (0 0 1); + +rpmCoeffs +{ + rpm -95.49; //-10 rad/s +} + +// ************************************************************************* // diff --git a/tutorials/incompressible/simpleSRFFoam/axialTurbine/constant/polyMesh/blockMeshDict.m4 b/tutorials/incompressible/simpleSRFFoam/axialTurbine/constant/polyMesh/blockMeshDict.m4 new file mode 100644 index 000000000..fa8078273 --- /dev/null +++ b/tutorials/incompressible/simpleSRFFoam/axialTurbine/constant/polyMesh/blockMeshDict.m4 @@ -0,0 +1,244 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | foam-extend: Open Source CFD | +| \\ / O peration | Version: 3.1 | +| \\ / A nd | Web: http://www.extend-project.de | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object blockMeshDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // +// General macros to create 2D/extruded-2D meshes + +changecom(//)changequote([,]) +define(calc, [esyscmd(perl -e 'printf ($1)')]) +//define(calc, [esyscmd(echo $1 | bc | tr -d \\n)]) +define(VCOUNT, 0) +define(vlabel, [[// ]Vertex $1 = VCOUNT define($1, VCOUNT)define([VCOUNT], incr(VCOUNT))]) +define(pi, calc(3.14159265/20)) + +define(hex2D, hex ($1b $2b $3b $4b $1t $2t $3t $4t)) +define(quad2D, ($1b $2b $2t $1t)) +define(frontQuad, ($1t $2t $3t $4t)) +define(backQuad, ($1b $4b $3b $2b)) + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +convertToMeters 1; + +// HUB AND SHROUD RADIUS +// Hub radius (m) +define(hr, 0.05) +// Shroud radius (m) +define(sr, 0.1) + +// RUNNER REGION GEOMETRY AND MESH PROPERTIES +// Runner inlet axial length (m) +define(RUial, 0.02) +// Runner axial length (m) +define(RUal, 0.1) +// Runner outlet axial length (m) +define(RUoal, 0.02) +// Number of runner blades per 360 degrees (integer!) +define(RUnb, 5) +// Number of cells in radial direction in runner +define(RUrc, 10) +// Number of cells in tangential direction between runner blades +define(RUtc, 10) +// Number of cells in axial direction at runner inlet +define(RUiac, 2) +// Number of cells in axial direction between runner blades +define(RUbac, 10) +// Number of cells in axial direction at runner outlet +define(RUoac, 2) + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +// TANGENTIAL PITCHES (RADIANS) +// Runner region +define(RUp, calc(2*pi/RUnb)) + +// TANGENTIAL SHIFTS BETWEEN AXIAL LEVELS (BOTTOM-UP) +// Runner region +// Tangential shift from level RU0 to RU1 +define(RUts01, calc(-1/10*RUp)) +// Tangential shift from level RU1 to RU2 +define(RUts12, calc(-4/5*RUp)) +// Tangential shift from level RU2 to RU3 +define(RUts23, calc(-1/10*RUp)) + +// AXIAL/TANGENTIAL BASE POINTS FOR EACH LEVEL (BOTTOM-UP): +// (CENTER OF RUNNER SET TO THETA=0, Z=0) +// Runner: +define(RUa0, calc(-RUoal-0.5*RUal)) //Center runner +define(RUt0, calc(-0.5*RUp-(0.5*RUts12))) //Center runner +define(RUt1, calc(RUt0+RUts01)) +define(RUt2, calc(RUt1+RUts12)) +define(RUt3, calc(RUt2+RUts23)) + +vertices //(radial [m], tangential [radians], axial [m]) +( +//Runner hub: + (hr RUt0 RUa0) vlabel(RU0lb) + (hr calc(RUt0+RUp) RUa0) vlabel(RU0rb) + (hr RUt1 calc(RUa0+RUoal)) vlabel(RU1lb) + (hr calc(RUt1+RUp) calc(RUa0+RUoal)) vlabel(RU1rb) + (hr RUt2 calc(RUa0+RUoal+RUal)) vlabel(RU2lb) + (hr calc(RUt2+RUp) calc(RUa0+RUoal+RUal)) vlabel(RU2rb) + (hr RUt3 calc(RUa0+RUoal+RUal+RUial)) vlabel(RU3lb) + (hr calc(RUt3+RUp) calc(RUa0+RUoal+RUal+RUial)) vlabel(RU3rb) + +//Runner shroud: + (sr RUt0 RUa0) vlabel(RU0lt) + (sr calc(RUt0+RUp) RUa0) vlabel(RU0rt) + (sr RUt1 calc(RUa0+RUoal)) vlabel(RU1lt) + (sr calc(RUt1+RUp) calc(RUa0+RUoal)) vlabel(RU1rt) + (sr RUt2 calc(RUa0+RUoal+RUal)) vlabel(RU2lt) + (sr calc(RUt2+RUp) calc(RUa0+RUoal+RUal)) vlabel(RU2rt) + (sr RUt3 calc(RUa0+RUoal+RUal+RUial)) vlabel(RU3lt) + (sr calc(RUt3+RUp) calc(RUa0+RUoal+RUal+RUial)) vlabel(RU3rt) +); + +blocks +( +//Runner: + hex2D(RU0l, RU0r, RU1r, RU1l) + rotor + (RUtc RUoac RUrc) + simpleGrading (1 1 1) + + hex2D(RU1l, RU1r, RU2r, RU2l) + rotor + (RUtc RUbac RUrc) + simpleGrading (1 1 1) + + hex2D(RU2l, RU2r, RU3r, RU3l) + rotor + (RUtc RUiac RUrc) + simpleGrading (1 1 1) +); + +edges // Inappropriate with arc due to coordinate conversion +( +//Runner + spline RU1lt RU2lt + ( + (sr calc(RUt1+0.65*(RUt2-(RUt1))) calc(RUa0+RUoal+0.5*RUal)) + ) + spline RU1lb RU2lb + ( + (hr calc(RUt1+0.65*(RUt2-(RUt1))) calc(RUa0+RUoal+0.5*RUal)) + ) + spline RU1rt RU2rt + ( + (sr calc(RUt1+RUp+0.75*(RUt2-(RUt1))) calc(RUa0+RUoal+0.5*RUal)) + ) + spline RU1rb RU2rb + ( + (hr calc(RUt1+RUp+0.75*(RUt2-(RUt1))) calc(RUa0+RUoal+0.5*RUal)) + ) +); + +boundary +( + RUINLET + { + type patch; + faces + ( + quad2D(RU3r, RU3l) + ); + } + + RUOUTLET + { + type patch; + faces + ( + quad2D(RU0l, RU0r) + ); + } + + RUCYCLIC1 + { + type cyclicGgi; + shadowPatch RUCYCLIC2; + zone RUCYCLIC1Zone; + bridgeOverlap false; + rotationAxis (0 0 1); + rotationAngle 72; + separationOffset (0 0 0); + faces + ( + quad2D(RU1l, RU0l) + quad2D(RU3l, RU2l) + ); + } + + RUCYCLIC2 + { + type cyclicGgi; + shadowPatch RUCYCLIC1; + zone RUCYCLIC2Zone; + bridgeOverlap false; + rotationAxis (0 0 1); + rotationAngle -72; + separationOffset (0 0 0); + faces + ( + quad2D(RU0r, RU1r) + quad2D(RU2r, RU3r) + ); + } + + //RUCYCLIC + //{ + // type cyclic; + // faces + // ( + // quad2D(RU1l, RU0l) + // quad2D(RU3l, RU2l) + // quad2D(RU0r, RU1r) + // quad2D(RU2r, RU3r) + // ); + //} + + RUBLADE + { + type wall; + faces + ( + quad2D(RU2l, RU1l) + quad2D(RU1r, RU2r) + ); + } + + RUHUB + { + type wall; + faces + ( + backQuad(RU0l, RU0r, RU1r, RU1l) + backQuad(RU1l, RU1r, RU2r, RU2l) + backQuad(RU2l, RU2r, RU3r, RU3l) + ); + } + + RUSHROUD + { + type wall; + faces + ( + frontQuad(RU0l, RU0r, RU1r, RU1l) + frontQuad(RU1l, RU1r, RU2r, RU2l) + frontQuad(RU2l, RU2r, RU3r, RU3l) + ); + } +); + +// ************************************************************************* // diff --git a/tutorials/incompressible/simpleSRFFoam/axialTurbine/constant/transportProperties b/tutorials/incompressible/simpleSRFFoam/axialTurbine/constant/transportProperties new file mode 100644 index 000000000..74e582a5c --- /dev/null +++ b/tutorials/incompressible/simpleSRFFoam/axialTurbine/constant/transportProperties @@ -0,0 +1,37 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | foam-extend: Open Source CFD | +| \\ / O peration | Version: 3.1 | +| \\ / A nd | Web: http://www.extend-project.de | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object transportProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +transportModel Newtonian; + +nu nu [0 2 -1 0 0 0 0] 1.5e-05; + +CrossPowerLawCoeffs +{ + nu0 nu0 [0 2 -1 0 0 0 0] 1e-06; + nuInf nuInf [0 2 -1 0 0 0 0] 1e-06; + m m [0 0 1 0 0 0 0] 1; + n n [0 0 0 0 0 0 0] 1; +} + +BirdCarreauCoeffs +{ + nu0 nu0 [0 2 -1 0 0 0 0] 1e-06; + nuInf nuInf [0 2 -1 0 0 0 0] 1e-06; + k k [0 0 1 0 0 0 0] 0; + n n [0 0 0 0 0 0 0] 1; +} + +// ************************************************************************* // diff --git a/tutorials/incompressible/simpleSRFFoam/axialTurbine/setBatchGgi b/tutorials/incompressible/simpleSRFFoam/axialTurbine/setBatchGgi new file mode 100644 index 000000000..a6f500cba --- /dev/null +++ b/tutorials/incompressible/simpleSRFFoam/axialTurbine/setBatchGgi @@ -0,0 +1,3 @@ +faceSet RUCYCLIC1Zone new patchToFace RUCYCLIC1 +faceSet RUCYCLIC2Zone new patchToFace RUCYCLIC2 +quit diff --git a/tutorials/incompressible/simpleSRFFoam/axialTurbine/system/controlDict b/tutorials/incompressible/simpleSRFFoam/axialTurbine/system/controlDict new file mode 100644 index 000000000..02638f43c --- /dev/null +++ b/tutorials/incompressible/simpleSRFFoam/axialTurbine/system/controlDict @@ -0,0 +1,62 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | foam-extend: Open Source CFD | +| \\ / O peration | Version: 3.1 | +| \\ / A nd | Web: http://www.extend-project.de | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object controlDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +application simpleSRFFoam; + +startFrom startTime; + +startTime 0; + +stopAt endTime; + +endTime 100; + +deltaT 1; + +writeControl timeStep; + +writeInterval 20; + +purgeWrite 0; + +writeFormat ascii; + +writePrecision 6; + +writeCompression compressed; + +timeFormat general; + +timePrecision 6; + +runTimeModifiable yes; + +// Compute the flux value on each side of a GGI interface +functions +( + ggiCheck + { + // Type of functionObject + type ggiCheck; + + phi phi; + + // Where to load it from (if not already in solver) + functionObjectLibs ("libcheckFunctionObjects.so"); + } +); + +// ************************************************************************* // diff --git a/tutorials/incompressible/simpleSRFFoam/axialTurbine/system/decomposeParDict b/tutorials/incompressible/simpleSRFFoam/axialTurbine/system/decomposeParDict new file mode 100644 index 000000000..97f4c4ff7 --- /dev/null +++ b/tutorials/incompressible/simpleSRFFoam/axialTurbine/system/decomposeParDict @@ -0,0 +1,67 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | foam-extend: Open Source CFD | +| \\ / O peration | Version: 3.1 | +| \\ / A nd | Web: http://www.extend-project.de | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object decomposeParDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +numberOfSubdomains 8; + +method metis; + +globalFaceZones +( + RUCYCLIC1Zone + RUCYCLIC2Zone +); + +simpleCoeffs +{ + n (2 2 1); + delta 0.001; +} + +hierarchicalCoeffs +{ + n (1 1 1); + delta 0.001; + order xyz; +} + +metisCoeffs +{ + processorWeights + ( + 1 + 1 + 1 + 1 + 1 + 1 + 1 + 1 + ); +} + +manualCoeffs +{ + dataFile ""; +} + +distributed no; + +roots +( +); + +// ************************************************************************* // + diff --git a/tutorials/incompressible/simpleSRFFoam/axialTurbine/system/fvSchemes b/tutorials/incompressible/simpleSRFFoam/axialTurbine/system/fvSchemes new file mode 100644 index 000000000..e7128dae6 --- /dev/null +++ b/tutorials/incompressible/simpleSRFFoam/axialTurbine/system/fvSchemes @@ -0,0 +1,71 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | foam-extend: Open Source CFD | +| \\ / O peration | Version: 3.1 | +| \\ / A nd | Web: http://www.extend-project.de | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object fvSchemes; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +ddtSchemes +{ + default steadyState; +} + +gradSchemes +{ + default Gauss linear; + grad(p) Gauss linear; + grad(Urel) Gauss linear; +} + +divSchemes +{ + default none; + div(phi,Urel) Gauss linearUpwind Gauss linear; + div(phi,k) Gauss upwind; + div(phi,epsilon) Gauss upwind; + div(phi,omega) Gauss upwind; + div(phi,R) Gauss upwind; + div(R) Gauss linear; + div(phi,nuTilda) Gauss upwind; + div((nuEff*dev(grad(Urel).T()))) Gauss linear; +} + +laplacianSchemes +{ + default none; + laplacian(nuEff,Urel) Gauss linear corrected; + laplacian((1|A(Urel)),p) Gauss linear corrected; + laplacian(DkEff,k) Gauss linear corrected; + laplacian(DepsilonEff,epsilon) Gauss linear corrected; + laplacian(DomegaEff,omega) Gauss linear corrected; + laplacian(DREff,R) Gauss linear corrected; + laplacian(DnuTildaEff,nuTilda) Gauss linear corrected; +} + +interpolationSchemes +{ + default linear; + interpolate(Urel) linear; +} + +snGradSchemes +{ + default corrected; +} + +fluxRequired +{ + default no; + p; +} + +// ************************************************************************* // diff --git a/tutorials/incompressible/simpleSRFFoam/axialTurbine/system/fvSolution b/tutorials/incompressible/simpleSRFFoam/axialTurbine/system/fvSolution new file mode 100644 index 000000000..b1cb8bf06 --- /dev/null +++ b/tutorials/incompressible/simpleSRFFoam/axialTurbine/system/fvSolution @@ -0,0 +1,86 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | foam-extend: Open Source CFD | +| \\ / O peration | Version: 3.1 | +| \\ / A nd | Web: http://www.extend-project.de | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object fvSolution; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solvers +{ + p + { + solver PCG; + preconditioner DIC; + tolerance 1e-06; + relTol 0.01; + }; + Urel + { + solver PBiCG; + preconditioner DILU; + tolerance 1e-05; + relTol 0.1; + }; + k + { + solver PBiCG; + preconditioner DILU; + tolerance 1e-05; + relTol 0.1; + }; + epsilon + { + solver PBiCG; + preconditioner DILU; + tolerance 1e-05; + relTol 0.1; + }; + omega + { + solver PBiCG; + preconditioner DILU; + tolerance 1e-05; + relTol 0.1; + }; + R + { + solver PBiCG; + preconditioner DILU; + tolerance 1e-05; + relTol 0.1; + }; + nuTilda + { + solver PBiCG; + preconditioner DILU; + tolerance 1e-05; + relTol 0.1; + }; +} + +SIMPLE +{ + nNonOrthogonalCorrectors 0; +} + +relaxationFactors +{ + p 0.3; + Urel 0.7; + k 0.7; + epsilon 0.7; + omega 0.7; + R 0.7; + nuTilda 0.7; +} + +// ************************************************************************* //