exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Modifications to run with either ggi or mixingPlane

Browse source
This commit is contained in:
Håkan Nilsson 2014-05-25 20:00:13 +02:00
parent 4fd91ff21d
commit 62bc3564d1
9 changed files with 366 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

1
tutorials/incompressible/MRFSimpleFoam/axialTurbine/Allclean
View file

@ -1,4 +1,5 @@
#!/bin/sh
( cd constant/polyMesh && rm -r boundary* blockMeshDict *Zones* faces* neighbour* owner* points* sets )
( cd system && rm changeDictionaryDict fvSolution )
rm -r 0 [1-9]* VTK log* processor*

10
tutorials/incompressible/MRFSimpleFoam/axialTurbine/Allrun → tutorials/incompressible/MRFSimpleFoam/axialTurbine/Allrun_ggi
View file

@ -9,12 +9,16 @@ application="MRFSimpleFoam"
#Create the mesh:
./makeMesh
#Fix constant/boundary:
changeDictionary
# Set 0-directory and create GGI set:
cp -r 0_orig 0
setSet -batch setBatchGgi
setsToZones -noFlipMap
#Fix constant/boundary:
cp system/changeDictionaryDict_ggi system/changeDictionaryDict
changeDictionary
#Pick settings:
cp system/fvSolution_ggi system/fvSolution
runApplication $application

24
tutorials/incompressible/MRFSimpleFoam/axialTurbine/Allrun_mixingPlane Executable file
View file

@ -0,0 +1,24 @@
#!/bin/sh
# Source tutorial run functions
. $WM_PROJECT_DIR/bin/tools/RunFunctions
application="MRFSimpleFoam"
#application="simpleFoam"
#Create the mesh:
./makeMesh
# Set 0-directory and create GGI set:
cp -r 0_orig 0
setSet -batch setBatchGgi
setsToZones -noFlipMap
#Fix constant/boundary:
cp system/changeDictionaryDict_mixingPlane system/changeDictionaryDict
changeDictionary
#Pick settings:
cp system/fvSolution_mixingPlane system/fvSolution
runApplication $application

0
tutorials/incompressible/MRFSimpleFoam/axialTurbine/system/changeDictionaryDict → tutorials/incompressible/MRFSimpleFoam/axialTurbine/system/changeDictionaryDict_ggi
View file

260
tutorials/incompressible/MRFSimpleFoam/axialTurbine/system/changeDictionaryDict_mixingPlane Normal file
View file

@ -0,0 +1,260 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 3.0 |
| \\ / A nd | Web: http://www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object changeDictionaryDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dictionaryReplacement
{
epsilon
{
boundaryField
{
GVOUTLET
{
type mixingPlane;
}
RUINLET
{
type mixingPlane;
}
RUOUTLET
{
type mixingPlane;
}
DTINLET
{
type mixingPlane;
}
}
}
k
{
boundaryField
{
GVOUTLET
{
type mixingPlane;
}
RUINLET
{
type mixingPlane;
}
RUOUTLET
{
type mixingPlane;
}
DTINLET
{
type mixingPlane;
}
}
}
U
{
boundaryField
{
GVOUTLET
{
type mixingPlane;
}
RUINLET
{
type mixingPlane;
}
RUOUTLET
{
type mixingPlane;
}
DTINLET
{
type mixingPlane;
}
}
}
p
{
boundaryField
{
GVOUTLET
{
type mixingPlane;
}
RUINLET
{
type mixingPlane;
}
RUOUTLET
{
type mixingPlane;
}
DTINLET
{
type mixingPlane;
}
}
}
boundary
{
GVINLET
{
type patch;
}
GVCYCLIC1
{
type cyclicGgi;
shadowPatch GVCYCLIC2;
zone GVCYCLIC1Zone;
bridgeOverlap false;
rotationAxis (0 0 1);
rotationAngle 72;
separationOffset (0 0 0);
}
GVCYCLIC2
{
type cyclicGgi;
shadowPatch GVCYCLIC1;
zone GVCYCLIC2Zone;
bridgeOverlap false;
rotationAxis (0 0 1);
rotationAngle -72;
separationOffset (0 0 0);
}
RUCYCLIC1
{
type cyclicGgi;
shadowPatch RUCYCLIC2;
zone RUCYCLIC1Zone;
bridgeOverlap false;
rotationAxis (0 0 1);
rotationAngle 72;
separationOffset (0 0 0);
}
RUCYCLIC2
{
type cyclicGgi;
shadowPatch RUCYCLIC1;
zone RUCYCLIC2Zone;
bridgeOverlap false;
rotationAxis (0 0 1);
rotationAngle -72;
separationOffset (0 0 0);
}
DTCYCLIC1
{
type cyclicGgi;
shadowPatch DTCYCLIC2;
zone DTCYCLIC1Zone;
bridgeOverlap false;
rotationAxis (0 0 1);
rotationAngle 72;
separationOffset (0 0 0);
}
DTCYCLIC2
{
type cyclicGgi;
shadowPatch DTCYCLIC1;
zone DTCYCLIC2Zone;
bridgeOverlap false;
rotationAxis (0 0 1);
rotationAngle -72;
separationOffset (0 0 0);
}
GVOUTLET
{
type mixingPlane;
shadowPatch RUINLET;
zone GVOUTLETZone;
coordinateSystem
{
type cylindrical;
//name mixingCS;
origin (0 0 0);
e1 (1 0 0);
e3 (0 0 1);
}
ribbonPatch
{
sweepAxis Theta;
stackAxis R;
discretisation bothPatches;
}
}
RUINLET
{
type mixingPlane;
shadowPatch GVOUTLET;
zone RUINLETZone;
coordinateSystem
{
type cylindrical;
//name mixingCS;
origin (0 0 0);
e1 (1 0 0);
e3 (0 0 1);
}
ribbonPatch
{
sweepAxis Theta;
stackAxis R;
discretisation bothPatches;
}
}
RUOUTLET
{
type mixingPlane;
shadowPatch DTINLET;
zone RUOUTLETZone;
coordinateSystem
{
type cylindrical;
//name mixingCS;
origin (0 0 0);
e1 (1 0 0);
e3 (0 0 1);
}
ribbonPatch
{
sweepAxis Theta;
stackAxis R;
discretisation bothPatches;
}
}
DTINLET
{
type mixingPlane;
shadowPatch RUOUTLET;
zone DTINLETZone;
coordinateSystem
{
type cylindrical;
//name mixingCS;
origin (0 0 0);
e1 (1 0 0);
e3 (0 0 1);
}
ribbonPatch
{
sweepAxis Theta;
stackAxis R;
discretisation bothPatches;
}
}
DTOUTLET
{
type patch;
}
}
}
// ************************************************************************* //

2
tutorials/incompressible/MRFSimpleFoam/axialTurbine/system/controlDict
View file

@ -14,7 +14,7 @@ FoamFile
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application simpleFoam;
application MRFSimpleFoam;
startFrom startTime;

7
tutorials/incompressible/MRFSimpleFoam/axialTurbine/system/fvSchemes
View file

@ -61,4 +61,11 @@ fluxRequired
p;
}
mixingPlane
{
default areaAveraging;
//U fluxAveragingAdjustMassFlow;
//p zeroGradientAreaAveragingMix;
}
// ************************************************************************* //

0
tutorials/incompressible/MRFSimpleFoam/axialTurbine/system/fvSolution → tutorials/incompressible/MRFSimpleFoam/axialTurbine/system/fvSolution_ggi
View file

66
tutorials/incompressible/MRFSimpleFoam/axialTurbine/system/fvSolution_mixingPlane Normal file
View file

@ -0,0 +1,66 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 3.0 |
| \\ / A nd | Web: http://www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
// For meshes with
// disparities in angle span across the mixingPlane interface, a
// solver for asymmetrical matrices is required, like BiCGStab. When
// in doubt, just select a solver for asymmetrical matrices for p.
p
{
solver BiCGStab;
preconditioner DILU;
tolerance 1e-08;
relTol 0.05;
};
U
{
solver BiCGStab;
preconditioner DILU;
tolerance 1e-07;
relTol 0.1;
};
k
{
solver BiCGStab;
preconditioner DILU;
tolerance 1e-07;
relTol 0.1;
};
epsilon
{
solver BiCGStab;
preconditioner DILU;
tolerance 1e-07;
relTol 0.1;
};
}
SIMPLE
{
nNonOrthogonalCorrectors 0;
}
relaxationFactors
{
p 0.3;
U 0.7;
k 0.7;
epsilon 0.7;
}
// ************************************************************************* //