Mixing plane tutorial update

This commit is contained in:
Hrvoje Jasak 2013-10-09 11:00:52 +01:00
parent 74630e76e2
commit d571209d75
58 changed files with 6882 additions and 717 deletions

File diff suppressed because it is too large Load diff

View file

@ -0,0 +1,74 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: 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 0.1;
boundaryField
{
outflow
{
type zeroGradient;
}
inflow
{
type fixedValue;
value uniform 0.1;
}
upstreamPerio1
{
type cyclicGgi;
value uniform 0.1;
}
upstreamPerio2
{
type cyclicGgi;
value uniform 0.1;
}
upstreamMixingPlanePatch
{
type mixingPlane;
value uniform 0.1;
}
downstreamMixingPlanePatch
{
type mixingPlane;
value uniform 0.1;
}
downstreamWall
{
type symmetryPlane;
}
upstreamWall
{
type symmetryPlane;
}
downstreamPerio1
{
type cyclicGgi;
value uniform 0.1;
}
downstreamPerio2
{
type cyclicGgi;
value uniform 0.1;
}
}
// ************************************************************************* //

View file

@ -0,0 +1,69 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: 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 0.1;
boundaryField
{
downstreamMixingPlanePatch
{
type mixingPlane;
}
upstreamMixingPlanePatch
{
type mixingPlane;
}
downstreamWall
{
type symmetryPlane;
}
upstreamWall
{
type symmetryPlane;
}
downstreamPerio1
{
type cyclicGgi;
}
downstreamPerio2
{
type cyclicGgi;
}
upstreamPerio1
{
type cyclicGgi;
}
upstreamPerio2
{
type cyclicGgi;
}
inflow
{
type fixedValue;
value uniform 0.1;
}
outflow
{
type zeroGradient;
}
}
// ************************************************************************* //

View file

@ -0,0 +1,74 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: 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.1;
boundaryField
{
outflow
{
type zeroGradient;
}
inflow
{
type fixedValue;
value uniform 0.01;
}
upstreamPerio1
{
type cyclicGgi;
value uniform 0.1;
}
upstreamPerio2
{
type cyclicGgi;
value uniform 0.1;
}
upstreamMixingPlanePatch
{
type mixingPlane;
value uniform 0.1;
}
downstreamMixingPlanePatch
{
type mixingPlane;
value uniform 0.1;
}
downstreamWall
{
type symmetryPlane;
}
upstreamWall
{
type symmetryPlane;
}
downstreamPerio1
{
type cyclicGgi;
value uniform 0.1;
}
downstreamPerio2
{
type cyclicGgi;
value uniform 0.1;
}
}
// ************************************************************************* //

View file

@ -0,0 +1,70 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: 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.1;
boundaryField
{
downstreamPerio1
{
type cyclicGgi;
}
downstreamPerio2
{
type cyclicGgi;
}
upstreamPerio1
{
type cyclicGgi;
}
upstreamPerio2
{
type cyclicGgi;
}
downstreamWall
{
type symmetryPlane;
}
upstreamWall
{
type symmetryPlane;
}
downstreamMixingPlanePatch
{
type mixingPlane;
}
upstreamMixingPlanePatch
{
type mixingPlane;
}
inflow
{
type fixedValue;
value uniform 0.01;
}
outflow
{
type zeroGradient;
}
}
// ************************************************************************* //

View file

@ -0,0 +1,75 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: 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
{
outflow
{
type calculated;
value uniform 0;
}
inflow
{
type calculated;
value uniform 0;
}
upstreamPerio1
{
type cyclicGgi;
value uniform 0;
}
upstreamPerio2
{
type cyclicGgi;
value uniform 0;
}
upstreamMixingPlanePatch
{
type mixingPlane;
value uniform 0;
}
downstreamMixingPlanePatch
{
type mixingPlane;
value uniform 0;
}
downstreamWall
{
type symmetryPlane;
}
upstreamWall
{
type symmetryPlane;
}
downstreamPerio1
{
type cyclicGgi;
value uniform 0;
}
downstreamPerio2
{
type cyclicGgi;
value uniform 0;
}
}
// ************************************************************************* //

View file

@ -0,0 +1,69 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5-dev |
| \\ / A nd | Revision: 1326 |
| \\/ M anipulation | Web: http://www.OpenFOAM.org |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [ 0 2 -2 0 0 0 0 ];
internalField uniform 0;
boundaryField
{
downstreamMixingPlanePatch
{
type mixingPlane;
}
upstreamMixingPlanePatch
{
type mixingPlane;
}
downstreamWall
{
type symmetryPlane;
}
upstreamWall
{
type symmetryPlane;
}
upstreamPerio1
{
type cyclicGgi;
}
upstreamPerio2
{
type cyclicGgi;
}
downstreamPerio1
{
type cyclicGgi;
}
downstreamPerio2
{
type cyclicGgi;
}
inflow
{
type zeroGradient;
}
outflow
{
type fixedValue;
value uniform 0;
}
}
// ************************************************************************* //

File diff suppressed because it is too large Load diff

View file

@ -0,0 +1,72 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
upstreamMixingPlanePatch
{
type mixingPlane;
}
downstreamMixingPlanePatch
{
type mixingPlane;
}
upstreamPerio1
{
type cyclicGgi;
}
upstreamPerio2
{
type cyclicGgi;
}
downstreamPerio1
{
type cyclicGgi;
}
downstreamPerio2
{
type cyclicGgi;
}
downstreamWall
{
type symmetryPlane;
}
upstreamWall
{
type symmetryPlane;
}
inflow
{
type surfaceNormalFixedValue;
refValue uniform -10;
value uniform (0 0 0);
}
outflow
{
type inletOutlet;
inletValue uniform (0 0 0);
value uniform (0 0 0);
}
}
// ************************************************************************* //

View file

@ -0,0 +1,69 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: 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 0.1;
boundaryField
{
downstreamMixingPlanePatch
{
type mixingPlane;
}
upstreamMixingPlanePatch
{
type mixingPlane;
}
downstreamWall
{
type symmetryPlane;
}
upstreamWall
{
type symmetryPlane;
}
downstreamPerio1
{
type cyclicGgi;
}
downstreamPerio2
{
type cyclicGgi;
}
upstreamPerio1
{
type cyclicGgi;
}
upstreamPerio2
{
type cyclicGgi;
}
inflow
{
type fixedValue;
value uniform 0.1;
}
outflow
{
type zeroGradient;
}
}
// ************************************************************************* //

View file

@ -0,0 +1,70 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: 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.01;
boundaryField
{
downstreamPerio1
{
type cyclicGgi;
}
downstreamPerio2
{
type cyclicGgi;
}
upstreamPerio1
{
type cyclicGgi;
}
upstreamPerio2
{
type cyclicGgi;
}
downstreamWall
{
type symmetryPlane;
}
upstreamWall
{
type symmetryPlane;
}
downstreamMixingPlanePatch
{
type mixingPlane;
}
upstreamMixingPlanePatch
{
type mixingPlane;
}
inflow
{
type fixedValue;
value uniform 0.01;
}
outflow
{
type zeroGradient;
}
}
// ************************************************************************* //

View file

@ -0,0 +1,69 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5-dev |
| \\ / A nd | Revision: 1326 |
| \\/ M anipulation | Web: http://www.OpenFOAM.org |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [ 0 2 -2 0 0 0 0 ];
internalField uniform 0;
boundaryField
{
downstreamMixingPlanePatch
{
type mixingPlane;
}
upstreamMixingPlanePatch
{
type mixingPlane;
}
downstreamWall
{
type symmetryPlane;
}
upstreamWall
{
type symmetryPlane;
}
upstreamPerio1
{
type cyclicGgi;
}
upstreamPerio2
{
type cyclicGgi;
}
downstreamPerio1
{
type cyclicGgi;
}
downstreamPerio2
{
type cyclicGgi;
}
inflow
{
type zeroGradient;
}
outflow
{
type fixedValue;
value uniform 0;
}
}
// ************************************************************************* //

View file

@ -1,74 +0,0 @@
// -*- C++ -*-
// File generated by PyFoam - sorry for the ugliness
FoamFile
{
version 2.0;
format ascii;
root "";
case "dom1_dom2_GGI";
instance "0";
local "";
class volVectorField;
object U;
}
dimensions [ 0 1 -1 0 0 0 0 ];
internalField uniform (0 0 0);
boundaryField
{
upstreamMixingPlanePatch
{
type mixingPlane;
}
downstreamMixingPlanePatch
{
type mixingPlane;
}
upstreamPerio1
{
type cyclicGgi;
}
upstreamPerio2
{
type cyclicGgi;
}
downstreamPerio1
{
type cyclicGgi;
}
downstreamPerio2
{
type cyclicGgi;
}
downstreamWall
{
type symmetryPlane;
}
upstreamWall
{
type symmetryPlane;
}
inflow
{
type profile1DfixedValue;
fileName "absolu_cr.csv";
fileFormat "turboCSV";
interpolateCoord "Z";
fieldName "Velocity";
}
outflow
{
type inletOutlet;
inletValue uniform (0 0 0);
value uniform (0 0 0);
}
} // ************************************************************************* //

View file

@ -1,68 +0,0 @@
// -*- C++ -*-
// File generated by PyFoam - sorry for the ugliness
FoamFile
{
version 2.0;
format ascii;
root "/home/page/OpenFOAM/bk3421-1.3/run/calculs/simpleFoam/MFR";
case "secteur360";
instance "" 0 "" ;
local "";
class volScalarField;
object alpha;
}
dimensions [ 0 0 0 0 0 0 0 ];
internalField uniform 1;
boundaryField
{
downstreamMixingPlanePatch
{
type mixingPlane;
}
upstreamMixingPlanePatch
{
type mixingPlane;
}
upstreamPerio1
{
type cyclicGgi;
}
upstreamPerio2
{
type cyclicGgi;
}
downstreamPerio1
{
type cyclicGgi;
}
downstreamPerio2
{
type cyclicGgi;
}
downstreamWall
{
type symmetryPlane;
}
upstreamWall
{
type symmetryPlane;
}
inflow
{
type zeroGradient;
}
outflow
{
type zeroGradient;
}
}

View file

@ -1,71 +0,0 @@
// -*- C++ -*-
// File generated by PyFoam - sorry for the ugliness
FoamFile
{
version 2.0;
format ascii;
root ".";
case "dom3_sFoam";
instance "0";
local "";
class volScalarField;
object epsilon;
}
dimensions [ 0 2 -3 0 0 0 0 ];
internalField uniform 0.1;
boundaryField
{
downstreamMixingPlanePatch
{
type mixingPlane;
}
upstreamMixingPlanePatch
{
type mixingPlane;
}
downstreamWall
{
type symmetryPlane;
}
upstreamWall
{
type symmetryPlane;
}
downstreamPerio1
{
type cyclicGgi;
}
downstreamPerio2
{
type cyclicGgi;
}
upstreamPerio1
{
type cyclicGgi;
}
upstreamPerio2
{
type cyclicGgi;
}
inflow
{
type profile1DfixedValue;
fileName "absolu_cr.csv";
fileFormat "turboCSV";
interpolateCoord "Z";
fieldName "Epsilon";
}
outflow
{
type zeroGradient;
}
} // ************************************************************************* //

View file

@ -1,71 +0,0 @@
// -*- C++ -*-
// File generated by PyFoam - sorry for the ugliness
FoamFile
{
version 2.0;
format ascii;
root ".";
case "dom3_sFoam";
instance "0";
local "";
class volScalarField;
object k;
}
dimensions [ 0 2 -2 0 0 0 0 ];
internalField uniform 0.1;
boundaryField
{
downstreamPerio1
{
type cyclicGgi;
}
downstreamPerio2
{
type cyclicGgi;
}
upstreamPerio1
{
type cyclicGgi;
}
upstreamPerio2
{
type cyclicGgi;
}
downstreamWall
{
type symmetryPlane;
}
upstreamWall
{
type symmetryPlane;
}
downstreamMixingPlanePatch
{
type mixingPlane;
}
upstreamMixingPlanePatch
{
type mixingPlane;
}
inflow
{
type profile1DfixedValue;
fileName "absolu_cr.csv";
fileFormat "turboCSV";
interpolateCoord "Z";
fieldName "K";
}
outflow
{
type zeroGradient;
}
} // ************************************************************************* //

View file

@ -1,71 +0,0 @@
// -*- C++ -*-
// File generated by PyFoam - sorry for the ugliness
FoamFile
{
version 2.0;
format ascii;
root "/home/math/OpenFOAM/test_MATH/test_GGI/simpleFoam_MFR_stitch/MFR/360deg_1on1";
case "dom1_dom2_GGI";
instance "0";
local "";
class volScalarField;
object p;
}
dimensions [ 0 2 -2 0 0 0 0 ];
internalField uniform 0;
boundaryField
{
downstreamMixingPlanePatch
{
type mixingPlane;
}
upstreamMixingPlanePatch
{
type mixingPlane;
}
downstreamWall
{
type symmetryPlane;
}
upstreamWall
{
type symmetryPlane;
}
upstreamPerio1
{
type cyclicGgi;
}
upstreamPerio2
{
type cyclicGgi;
}
downstreamPerio1
{
type cyclicGgi;
}
downstreamPerio2
{
type cyclicGgi;
}
inflow
{
type zeroGradient;
}
outflow
{
type fixedValue;
value uniform 0;
}
} // ************************************************************************* //

View file

@ -1,69 +0,0 @@
// -*- C++ -*-
// File generated by PyFoam - sorry for the ugliness
FoamFile
{
version 2.0;
format ascii;
root "/home/math/OpenFOAM/test_MATH/test_GGI/simpleFoam_MFR_stitch/MFR/360deg_1on1";
case "dom1_dom2_GGI";
instance "0";
local "";
class surfaceScalarField;
object phi;
}
dimensions [ 0 3 -1 0 0 0 0 ];
internalField uniform 0;
boundaryField
{
downstreamMixingPlanePatch
{
type mixingPlane;
}
upstreamMixingPlanePatch
{
type mixingPlane;
}
downstreamWall
{
type symmetryPlane;
}
upstreamWall
{
type symmetryPlane;
}
upstreamPerio1
{
type cyclicGgi;
}
upstreamPerio2
{
type cyclicGgi;
}
downstreamPerio1
{
type cyclicGgi;
}
downstreamPerio2
{
type cyclicGgi;
}
inflow
{
type fixedValue;
value uniform 0;
}
outflow
{
type fixedValue;
value uniform 0;
}
} // ************************************************************************* //

View file

@ -1,50 +1,7 @@
#!/bin/bash
#
# Source clean functions
#!/bin/sh
# Source tutorial clean functions
. $WM_PROJECT_DIR/bin/tools/CleanFunctions
#set -x
cleanSampling ()
{
rm -rf sets > /dev/null 2>&1
}
cleanPostProcessing ()
{
rm -rf postProcessing/*.png > /dev/null 2>&1
rm -rf postProcessing/log.sample > /dev/null 2>&1
rm -rf postProcessing/log.foamLog > /dev/null 2>&1
rm -rf postProcessing/log.gnuplotPlotResiduals > /dev/null 2>&1
rm -rf postProcessing/log.gnuplotCompareAll > /dev/null 2>&1
rm -rf postProcessing/logs > /dev/null 2>&1
}
cleanFilesFromM4 ()
{
rm {constant/polyMesh/boundary,constant/polyMesh/blockMeshDict} > /dev/null 2>&1
}
cleanTimeZero ()
{
rm -rf 0 > /dev/null 2>&1
}
cleanPyFoam ()
{
rm PlyParser_FoamFileParser_parsetab.* > /dev/null 2>&1
}
cleanVTK ()
{
rm -rf VTK > /dev/null 2>&1
}
# Cleaning up the case
cleanSampling
cleanPostProcessing
cleanCase
cleanFilesFromM4
cleanTimeZero
cleanPyFoam
cleanVTK
\rm -rf *.py

View file

@ -1,22 +1,11 @@
#!/bin/bash
#
#set -x
# Source tutorial run functions
. $WM_PROJECT_DIR/bin/tools/RunFunctions
# Load additional RunFunctions
. ./RunFunctionsSupplemental
export caseName=`basename $PWD`
#We stop this script when we encounter a problem
#trap "exit -1" ERR
echo "$caseName: Creating the file blockMeshDict using m4"
(cd constant/polyMesh; m4 -P blockMeshDict.m4 > blockMeshDict)
echo "$caseName: Running blockMesh"
runApplicationAndReportOnError blockMesh
runApplication blockMesh
echo "$caseName: Updating the boundary file"
# Let's see if pyFoamChangeMixingPlaneBoundary.py is available
@ -78,22 +67,15 @@ fi
echo "$caseName: Creating the starting time directory"
cp -r 0_orig 0
cp -r save 0
echo "$caseName: Creating zones for the mixingPlane"
runApplicationAndReportOnError setSet -batch setBatchMixingPlane
runApplicationAndReportOnError setsToZones -noFlipMap
mv log.setSet log.setSet_MixingPlane
mv log.setsToZones log.setsToZones_MixingPlane
runApplicationAndReportOnError setSet -batch setBatchGGIZone
runApplicationAndReportOnError setsToZones -noFlipMap
mv log.setSet log.setSet_GGIZone
mv log.setsToZones log.setsToZones_GGIZone
runApplication setSet -batch setBatch.batch
runApplication setsToZones -noFlipMap
echo "$caseName: Running potentialFoam"
runApplicationAndReportOnError potentialFoam
runApplication potentialFoam
echo "$caseName: Running simpleFoam"
runApplicationAndReportOnError simpleFoam
runApplication simpleFoam

View file

@ -1,66 +0,0 @@
#---------------------------------*- sh -*-------------------------------------
# ========= |
# \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
# \\ / O peration |
# \\ / A nd | Copyright held by original author
# \\/ M anipulation |
#------------------------------------------------------------------------------
# License
# This file is part of OpenFOAM.
#
# OpenFOAM is free software; you can redistribute it and/or modify it
# under the terms of the GNU General Public License as published by the
# Free Software Foundation; either version 2 of the License, or (at your
# option) any later version.
#
# OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
# FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
# for more details.
#
# You should have received a copy of the GNU General Public License
# along with OpenFOAM; if not, write to the Free Software Foundation,
# Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
#
# Script
# RunFunctionsSupplemental
#
# Description
# Additionnal run functions, based on the original RunFunctions
#
# Author:
# Martin Beaudoin, Hydro-Quebec. All rights reserved
#
#------------------------------------------------------------------------------
# Exception handler: dump the tail of the log file on error
verbose_report_on_runApplication_error()
{
if [ "$reportOnErrorOnlyOnce" ] ; then
echo "Error running $APP_RUN..."
if [ -f log.$APP_RUN ] ; then
tail log.$APP_RUN
fi
fi
# Do not report again when exiting
unset reportOnErrorOnlyOnce
}
# Variation of runApplication with exception handling
runApplicationAndReportOnError ()
{
trap 'verbose_report_on_runApplication_error' ERR
APP_RUN=$1; shift
reportOnErrorOnlyOnce=1
if [ -f log.$APP_RUN ] ; then
echo "$APP_RUN already run on $PWD: remove log file to run"
else
echo "Running $APP_RUN on $PWD"
$APP_RUN $* > log.$APP_RUN 2>&1
fi
}

View file

@ -0,0 +1,161 @@
//process this file using: m4 -P blockMeshDict.m4 > blockMeshDict
//m4 definitions -----------------------------
//Geometry -----------------------------------
// 2 planes levels
// Angle span for inner block
// Angle span for outer block
// Radial dimensions
// Mesh parameters
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
convertToMeters 1;
vertices
(
//Plane A:
//Bottom of curved block
(0.499999999999998 0.86602540378444 0.0) // Vertex A0 = 0
(0.999999999999996 1.73205080756888 0.0) // Vertex A1 = 1
(-0.866025403784436 0.500000000000004 0.0) // Vertex A2 = 2
(-1.73205080756887 1.00000000000001 0.0) // Vertex A3 = 3
//Plane B:
//Top of curved block
(0.499999999999998 0.86602540378444 0.1) // Vertex B0 = 4
(0.999999999999996 1.73205080756888 0.1) // Vertex B1 = 5
(-0.866025403784436 0.500000000000004 0.1) // Vertex B2 = 6
(-1.73205080756887 1.00000000000001 0.1) // Vertex B3 = 7
//Plane A: Bottom of straight block
(1.99238939618349 0.174311485495316 0.0) // Vertex B4 = 8
(2.98858409427524 0.261467228242975 0.0) // Vertex B5 = 9
(1.50941916044554 1.31211805798101 0.0) // Vertex B6 = 10
(2.26412874066832 1.96817708697152 0.0) // Vertex B7 = 11
//Plane B: Top of straight block
(1.99238939618349 0.174311485495316 0.1) // Vertex C0 = 12
(2.98858409427524 0.261467228242975 0.1) // Vertex C1 = 13
(1.50941916044554 1.31211805798101 0.1) // Vertex C2 = 14
(2.26412874066832 1.96817708697152 0.1) // Vertex C3 = 15
);
blocks
(
hex ( 0 1 3 2 4 5 7 6 ) (25 17 1) simpleGrading (1 1 1.0)
hex ( 8 9 11 10 12 13 15 14 ) (25 27 1) simpleGrading (1 1 1.0)
);
edges
(
// --- PLANE A: Bottom of curved block
arc 0 2 (-0.25881904510252 0.965925826289069 0.0)
arc 1 3 (-0.517638090205039 1.93185165257814 0.0)
// --- PLANE B: Top of curved block
arc 4 6 (-0.25881904510252 0.965925826289069 0.1)
arc 5 7 (-0.517638090205039 1.93185165257814 0.1)
// --- PLANE A: Bottom of straight block
arc 8 10 (1.84100970690488 0.781462256978547 0.0)
arc 9 11 (2.76151456035732 1.17219338546782 0.0)
// --- PLANE B: Top of straight block
arc 12 14 (1.84100970690488 0.781462256978547 0.1)
arc 13 15 (2.76151456035732 1.17219338546782 0.1)
);
patches
(
patch outflow
(
( 9 11 15 13 )
)
patch inflow
(
( 0 4 6 2 )
)
cyclicGgi upstreamPerio1
(
(8 9 13 12)
)
cyclicGgi upstreamPerio2
(
(10 14 15 11)
)
mixingPlane upstreamMixingPlanePatch
(
( 8 12 14 10 )
)
mixingPlane downstreamMixingPlanePatch
(
( 1 3 7 5 )
)
symmetryPlane downstreamWall
(
( 0 2 3 1 )
( 4 5 7 6 )
)
symmetryPlane upstreamWall
(
( 13 15 14 12)
(8 10 11 9)
)
cyclicGgi downstreamPerio1
(
(0 1 5 4)
)
cyclicGgi downstreamPerio2
(
(2 6 7 3)
)
);
mergePatchPairs
(
);

View file

@ -34,13 +34,8 @@ m4_define(grading, 1.0)
FoamFile
{
version 2.0; format ascii;
root "";
case "";
instance "";
local "";
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
@ -108,7 +103,7 @@ patches
(
patch outflow
(
( B5 B7 C3 C1 )
( B5 B7 C3 C1 )
)
patch inflow
(
@ -128,7 +123,7 @@ patches
mixingPlane upstreamMixingPlanePatch
(
( B4 C0 C2 B6 )
( B4 C0 C2 B6 )
)
mixingPlane downstreamMixingPlanePatch
@ -139,14 +134,14 @@ patches
symmetryPlane downstreamWall
(
( A0 A2 A3 A1 )
( A0 A2 A3 A1 )
( B0 B1 B3 B2 )
)
symmetryPlane upstreamWall
(
( C1 C3 C2 C0)
(B4 B6 B7 B5)
( C1 C3 C2 C0)
(B4 B6 B7 B5)
)
cyclicGgi downstreamPerio1
@ -156,7 +151,7 @@ patches
cyclicGgi downstreamPerio2
(
(A2 B2 B3 A3)
(A2 B2 B3 A3)
)
);

View file

@ -0,0 +1,153 @@
// -*- C++ -*-
// File generated by PyFoam - sorry for the ugliness
FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/polyMesh";
object boundary;
}
(
outflow
{
type patch;
nFaces 27;
startFace 2106;
}
inflow
{
type patch;
nFaces 17;
startFace 2133;
}
upstreamPerio1
{
type cyclicGgi;
nFaces 25;
startFace 2150;
shadowPatch upstreamPerio2;
zone upstreamPerio1Zone;
bridgeOverlap false;
rotationAxis (0 0 1);
rotationAngle 36;
separationOffset (0 0 0);
}
upstreamPerio2
{
type cyclicGgi;
nFaces 25;
startFace 2175;
shadowPatch upstreamPerio1;
zone upstreamPerio2Zone;
bridgeOverlap false;
rotationAxis (0 0 1);
rotationAngle -36;
separationOffset (0 0 0);
}
upstreamMixingPlanePatch
{
type mixingPlane;
nFaces 27;
startFace 2200;
shadowPatch downstreamMixingPlanePatch;
zone upstreamMixingPlaneZone;
coordinateSystem
{
name mixingCS;
type cylindrical;
origin (0 0 0);
e1 (1 0 0);
e3 (0 0 1);
}
ribbonPatch
{
sweepAxis Theta;
stackAxis Z;
discretisation bothPatches;
}
}
downstreamMixingPlanePatch
{
type mixingPlane;
nFaces 17;
startFace 2227;
shadowPatch upstreamMixingPlanePatch; /* empty */
zone downstreamMixingPlaneZone;
coordinateSystem
{
name mixingCS;
type cylindrical;
origin (0 0 0);
e1 (1 0 0);
e3 (0 0 1);
}
ribbonPatch
{
sweepAxis Theta;
stackAxis Z;
discretisation bothPatches;
}
}
downstreamWall
{
type symmetryPlane;
nFaces 850;
startFace 2244;
}
upstreamWall
{
type symmetryPlane;
nFaces 1350;
startFace 3094;
}
downstreamPerio1
{
type cyclicGgi;
nFaces 25;
startFace 4444;
shadowPatch downstreamPerio2;
zone downstreamPerio1Zone;
bridgeOverlap false;
rotationAxis (0 0 1);
rotationAngle 90;
separationOffset (0 0 0);
}
downstreamPerio2
{
type cyclicGgi;
nFaces 25;
startFace 4469;
shadowPatch downstreamPerio1;
zone downstreamPerio2Zone;
bridgeOverlap false;
rotationAxis (0 0 1);
rotationAngle -90;
separationOffset (0 0 0);
}
)
// ************************************************************************* //

View file

@ -1,136 +1,134 @@
// -*- C++ -*-
// File generated by PyFoam - sorry for the ugliness
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/polyMesh";
object boundary;
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/polyMesh";
object boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
10
(
outflow
{
type patch;
nFaces 27;
startFace 2106;
}
inflow
{
type patch;
nFaces 17;
startFace 2133;
}
upstreamPerio1
{
type cyclicGgi;
nFaces 25;
startFace 2150;
shadowPatch upstreamPerio2;
zone upstreamPerio1Zone;
bridgeOverlap false;
rotationAxis (0 0 1);
rotationAngle 36;
separationOffset (0 0 0);
}
upstreamPerio2
{
type cyclicGgi;
nFaces 25;
startFace 2175;
shadowPatch upstreamPerio1;
zone upstreamPerio2Zone;
bridgeOverlap false;
rotationAxis (0 0 1);
rotationAngle -36;
separationOffset (0 0 0);
}
upstreamMixingPlanePatch
{
type mixingPlane;
nFaces 27;
startFace 2200;
shadowPatch downstreamMixingPlanePatch;
coordinateSystem
outflow
{
type cylindrical;
origin (0 0 0);
axis (0 0 1);
direction (1 1 0);
type patch;
nFaces 27;
startFace 2106;
}
assembly both;
orientation dirZ_spanTheta;
}
downstreamMixingPlanePatch
{
type mixingPlane;
nFaces 17;
startFace 2227;
shadowPatch upstreamMixingPlanePatch;
coordinateSystem
inflow
{
type cylindrical;
origin (0 0 0);
axis (0 0 1);
direction (1 1 0);
type patch;
nFaces 17;
startFace 2133;
}
assembly both;
orientation dirZ_spanTheta;
}
downstreamWall
upstreamPerio1
{
type cyclicGgi;
nFaces 25;
startFace 2150;
shadowPatch upstreamPerio2;
zone upstreamPerio1Zone;
bridgeOverlap false;
rotationAxis (0 0 1);
rotationAngle 36;
separationOffset (0 0 0);
}
upstreamPerio2
{
type cyclicGgi;
nFaces 25;
startFace 2175;
shadowPatch upstreamPerio1;
zone upstreamPerio2Zone;
bridgeOverlap false;
rotationAxis (0 0 1);
rotationAngle -36;
separationOffset (0 0 0);
}
upstreamMixingPlanePatch
{
type mixingPlane;
nFaces 27;
startFace 2200;
shadowPatch downstreamMixingPlanePatch;
zone upstreamMixingPlaneZone;
ribbonPatch
{
discretisation bothPatches;
sweepAxis Theta;
stackAxis Z;
}
{
type symmetryPlane;
nFaces 850;
startFace 2244;
}
upstreamWall
{
type symmetryPlane;
nFaces 1350;
startFace 3094;
}
downstreamPerio1
{
type cyclicGgi;
nFaces 25;
startFace 4444;
shadowPatch downstreamPerio2;
zone downstreamPerio1Zone;
bridgeOverlap false;
rotationAxis (0 0 1);
rotationAngle 90;
separationOffset (0 0 0);
}
downstreamPerio2
{
type cyclicGgi;
nFaces 25;
startFace 4469;
shadowPatch downstreamPerio1;
zone downstreamPerio2Zone;
bridgeOverlap false;
rotationAxis (0 0 1);
rotationAngle -90;
separationOffset (0 0 0);
}
coordinateSystem
{
type cylindrical;
origin (0 0 0);
axis (0 0 1);
direction (1 1 0);
}
}
downstreamMixingPlanePatch
{
type mixingPlane;
nFaces 17;
startFace 2227;
shadowPatch upstreamMixingPlanePatch;
zone downstreamMixingPlaneZone;
coordinateSystem
{
type cylindrical;
origin (0 0 0);
axis (0 0 1);
direction (1 1 0);
}
assembly both;
orientation dirZ_spanTheta;
}
downstreamWall
{
type symmetryPlane;
nFaces 850;
startFace 2244;
}
upstreamWall
{
type symmetryPlane;
nFaces 1350;
startFace 3094;
}
downstreamPerio1
{
type cyclicGgi;
nFaces 25;
startFace 4444;
shadowPatch downstreamPerio2;
zone downstreamPerio1Zone;
bridgeOverlap false;
rotationAxis (0 0 1);
rotationAngle 90;
separationOffset (0 0 0);
}
downstreamPerio2
{
type cyclicGgi;
nFaces 25;
startFace 4469;
shadowPatch downstreamPerio1;
zone downstreamPerio2Zone;
bridgeOverlap false;
rotationAxis (0 0 1);
rotationAngle -90;
separationOffset (0 0 0);
}
)
// ************************************************************************* //

View file

@ -0,0 +1,72 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
upstreamMixingPlanePatch
{
type mixingPlane;
}
downstreamMixingPlanePatch
{
type mixingPlane;
}
upstreamPerio1
{
type cyclicGgi;
}
upstreamPerio2
{
type cyclicGgi;
}
downstreamPerio1
{
type cyclicGgi;
}
downstreamPerio2
{
type cyclicGgi;
}
downstreamWall
{
type symmetryPlane;
}
upstreamWall
{
type symmetryPlane;
}
inflow
{
type surfaceNormalFixedValue;
refValue uniform -10;
value uniform (0 0 0);
}
outflow
{
type inletOutlet;
inletValue uniform (0 0 0);
value uniform (0 0 0);
}
}
// ************************************************************************* //

View file

@ -0,0 +1,69 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: 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 0.1;
boundaryField
{
downstreamMixingPlanePatch
{
type mixingPlane;
}
upstreamMixingPlanePatch
{
type mixingPlane;
}
downstreamWall
{
type symmetryPlane;
}
upstreamWall
{
type symmetryPlane;
}
downstreamPerio1
{
type cyclicGgi;
}
downstreamPerio2
{
type cyclicGgi;
}
upstreamPerio1
{
type cyclicGgi;
}
upstreamPerio2
{
type cyclicGgi;
}
inflow
{
type fixedValue;
value uniform 0.1;
}
outflow
{
type zeroGradient;
}
}
// ************************************************************************* //

View file

@ -0,0 +1,70 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: 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.01;
boundaryField
{
downstreamPerio1
{
type cyclicGgi;
}
downstreamPerio2
{
type cyclicGgi;
}
upstreamPerio1
{
type cyclicGgi;
}
upstreamPerio2
{
type cyclicGgi;
}
downstreamWall
{
type symmetryPlane;
}
upstreamWall
{
type symmetryPlane;
}
downstreamMixingPlanePatch
{
type mixingPlane;
}
upstreamMixingPlanePatch
{
type mixingPlane;
}
inflow
{
type fixedValue;
value uniform 0.01;
}
outflow
{
type zeroGradient;
}
}
// ************************************************************************* //

View file

@ -0,0 +1,69 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5-dev |
| \\ / A nd | Revision: 1326 |
| \\/ M anipulation | Web: http://www.OpenFOAM.org |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [ 0 2 -2 0 0 0 0 ];
internalField uniform 0;
boundaryField
{
downstreamMixingPlanePatch
{
type mixingPlane;
}
upstreamMixingPlanePatch
{
type mixingPlane;
}
downstreamWall
{
type symmetryPlane;
}
upstreamWall
{
type symmetryPlane;
}
upstreamPerio1
{
type cyclicGgi;
}
upstreamPerio2
{
type cyclicGgi;
}
downstreamPerio1
{
type cyclicGgi;
}
downstreamPerio2
{
type cyclicGgi;
}
inflow
{
type zeroGradient;
}
outflow
{
type fixedValue;
value uniform 0;
}
}
// ************************************************************************* //

View file

@ -2,4 +2,6 @@ faceSet upstreamPerio1Zone new patchToFace upstreamPerio1
faceSet upstreamPerio2Zone new patchToFace upstreamPerio2
faceSet downstreamPerio1Zone new patchToFace downstreamPerio1
faceSet downstreamPerio2Zone new patchToFace downstreamPerio2
faceSet upstreamMixingPlaneZone new patchToFace upstreamMixingPlanePatch
faceSet downstreamMixingPlaneZone new patchToFace downstreamMixingPlanePatch
quit

View file

@ -1,3 +0,0 @@
faceSet upstreamMixingPlaneZone new patchToFace upstreamMixingPlanePatch
faceSet downstreamMixingPlaneZone new patchToFace downstreamMixingPlanePatch
quit

View file

@ -1,4 +1,4 @@
/*---------------------------------------------------------------------------*\
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.3 |
@ -12,19 +12,13 @@ FoamFile
{
version 2.0;
format ascii;
root "/home/page/OpenFOAM/bk3421-1.3/run/calculs/simpleFoam/MFR/sector_1on1";
case "dom1Dom2";
instance "system";
local "";
class dictionary;
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application potentialFoam;
application simpleFoam;
startFrom startTime;
@ -38,7 +32,7 @@ deltaT 1;
writeControl timeStep;
writeInterval 1;
writeInterval 10;
purgeWrite 0;
@ -56,9 +50,6 @@ graphFormat raw;
runTimeModifiable yes;
libs ( "libOpenFoamTurbo.so" );
// ************************************************************************* //
// Compute the flux value on each side of a GGI interface
functions
(
@ -82,9 +73,11 @@ functions
masterPatchName upstreamMixingPlanePatch;
masterPatchScaleFactor 10.0;
shadowPatchScaleFactor 4.0;
phi phi;
phi phi;
// Where to load it from (if not already in solver)
functionObjectLibs ("libcheckFunctionObjects.so");
}
);
// ************************************************************************* //

View file

@ -71,7 +71,7 @@ mixingPlane
{
default areaAveraging;
p areaAveraging;
U fluxAveraging;
U areaAveraging;
k fluxAveraging;
epsilon fluxAveraging;

View file

@ -0,0 +1,51 @@
# git-ls-files --others --exclude-from=.git/info/exclude
# Lines that start with '#' are comments.
# editor and misc backup files - anywhere
*~
.*~
*.bak
*.bak[0-9][0-9]
*.orig
*.orig[0-9][0-9]
\#*\#
# file-browser settings - anywhere
.directory
# CVS recovered versions - anywhere
.#*
# SVN directories - anywhere
.svn/
# OpenFOAM results
[0-9]*/
!/0/
processor*
!*/polyMesh/blockMeshDict
cellToRegion*
log*
# packages - anywhere
*.tar.bz2
*.tar.gz
*.tar
*.tgz
*.gtgz
# Pictures and movies
*.png
*.jpg
*.jpeg
*.bmp
*.png
*.avi
*.mp4
*.mpg
#end-of-file

View file

@ -0,0 +1,194 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "2000";
object UFluxAvg;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
bladeDomB
{
type fixedValue;
value uniform (0 0 0);
}
bladeDomA
{
type fixedValue;
value uniform (0 0 0);
}
outletDomB
{
type inletOutlet;
inletValue uniform (0 0 0);
value uniform (0 0 0);
}
inletDomA
{
type fixedValue;
value nonuniform List<vector>
100
(
(2.27476 -1.38401 0)
(2.27329 -1.38642 0)
(2.27182 -1.38882 0)
(2.27035 -1.39123 0)
(2.26887 -1.39363 0)
(2.2674 -1.39604 0)
(2.26592 -1.39844 0)
(2.26443 -1.40083 0)
(2.26295 -1.40323 0)
(2.26146 -1.40563 0)
(2.25997 -1.40802 0)
(2.25848 -1.41041 0)
(2.25699 -1.4128 0)
(2.25549 -1.41519 0)
(2.25399 -1.41758 0)
(2.25249 -1.41997 0)
(2.25098 -1.42235 0)
(2.24947 -1.42474 0)
(2.24796 -1.42712 0)
(2.24645 -1.4295 0)
(2.24494 -1.43187 0)
(2.24342 -1.43425 0)
(2.2419 -1.43663 0)
(2.24037 -1.439 0)
(2.23885 -1.44137 0)
(2.23732 -1.44374 0)
(2.23579 -1.44611 0)
(2.23426 -1.44848 0)
(2.23272 -1.45084 0)
(2.23119 -1.45321 0)
(2.22965 -1.45557 0)
(2.2281 -1.45793 0)
(2.22656 -1.46029 0)
(2.22501 -1.46264 0)
(2.22346 -1.465 0)
(2.22191 -1.46735 0)
(2.22035 -1.46971 0)
(2.2188 -1.47206 0)
(2.21723 -1.4744 0)
(2.21567 -1.47675 0)
(2.21415 -1.47903 0)
(2.21267 -1.48125 0)
(2.21119 -1.48346 0)
(2.2097 -1.48567 0)
(2.20822 -1.48788 0)
(2.20673 -1.49008 0)
(2.20524 -1.49229 0)
(2.20374 -1.4945 0)
(2.20225 -1.4967 0)
(2.20075 -1.4989 0)
(2.19925 -1.5011 0)
(2.19775 -1.5033 0)
(2.19624 -1.5055 0)
(2.19474 -1.50769 0)
(2.19323 -1.50988 0)
(2.19172 -1.51208 0)
(2.1902 -1.51427 0)
(2.18869 -1.51646 0)
(2.18717 -1.51865 0)
(2.18565 -1.52083 0)
(2.18408 -1.52308 0)
(2.18247 -1.52539 0)
(2.18085 -1.5277 0)
(2.17923 -1.53001 0)
(2.17761 -1.53232 0)
(2.17599 -1.53462 0)
(2.17436 -1.53693 0)
(2.17273 -1.53923 0)
(2.1711 -1.54153 0)
(2.16947 -1.54383 0)
(2.16783 -1.54613 0)
(2.16619 -1.54842 0)
(2.16455 -1.55071 0)
(2.16291 -1.553 0)
(2.16126 -1.55529 0)
(2.15962 -1.55758 0)
(2.15796 -1.55987 0)
(2.15631 -1.56215 0)
(2.15466 -1.56444 0)
(2.153 -1.56672 0)
(2.15134 -1.569 0)
(2.14967 -1.57127 0)
(2.14801 -1.57355 0)
(2.14634 -1.57582 0)
(2.14467 -1.57809 0)
(2.143 -1.58036 0)
(2.14133 -1.58263 0)
(2.13965 -1.5849 0)
(2.13797 -1.58716 0)
(2.13629 -1.58943 0)
(2.1346 -1.59169 0)
(2.13292 -1.59395 0)
(2.13123 -1.59621 0)
(2.12953 -1.59846 0)
(2.12784 -1.60072 0)
(2.12614 -1.60297 0)
(2.12445 -1.60522 0)
(2.12274 -1.60747 0)
(2.12104 -1.60972 0)
(2.11933 -1.61196 0)
)
;
}
outletDomA
{
type mixingPlane;
fluxAveraging true;
surfaceAveraging false;
value uniform (0 0 0);
}
inletDomB
{
type mixingPlane;
fluxAveraging false;
surfaceAveraging false;
value uniform (0 0 0);
}
ggi1DomB
{
type ggi;
}
ggi1DomA
{
type ggi;
}
ggi2DomB
{
type ggi;
}
ggi2DomA
{
type ggi;
}
topAndBottom
{
type empty;
}
cyclicDomA
{
type cyclic;
}
cyclicDomB
{
type cyclic;
}
}
// ************************************************************************* //

View file

@ -0,0 +1,86 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "2000";
object pFluxAvg;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
bladeDomB
{
type zeroGradient;
}
bladeDomA
{
type zeroGradient;
}
outletDomB
{
type fixedValue;
value uniform 0;
}
inletDomA
{
type zeroGradient;
}
outletDomA
{
type mixingPlane;
fluxAveraging false;
surfaceAveraging false;
value uniform 0;
}
inletDomB
{
type mixingPlane;
fluxAveraging false;
surfaceAveraging false;
value uniform 0;
}
ggi1DomB
{
type ggi;
}
ggi1DomA
{
type ggi;
}
ggi2DomB
{
type ggi;
}
ggi2DomA
{
type ggi;
}
topAndBottom
{
type empty;
}
cyclicDomA
{
type cyclic;
}
cyclicDomB
{
type cyclic;
}
}
// ************************************************************************* //

View file

@ -0,0 +1,10 @@
#!/bin/sh
# Source tutorial clean functions
. $WM_PROJECT_DIR/bin/tools/CleanFunctions
cleanTimeDirectories
\rm -rf 0
\cp -r save 0
\rm -rf constant/polyMesh/sets constant/polyMesh/pointZones* constant/polyMesh/faceZones* constant/polyMesh/cellZones*
\rm -rf VTK

View file

@ -0,0 +1,13 @@
#!/bin/bash
# Source tutorial run functions
. $WM_PROJECT_DIR/bin/tools/RunFunctions
\cp -f boundary constant/polyMesh/boundary
\rm -rf 0
\cp -r save 0
runApplication setSet -batch setBatch.batch
runApplication setsToZones -noFlipMap
runApplication potentialFoam
runApplication simpleFoam

View file

@ -0,0 +1,141 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext_mp |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "2/polyMesh";
object boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
13
(
bladeDomB
{
type wall;
nFaces 176;
startFace 49680;
}
bladeDomA
{
type wall;
nFaces 176;
startFace 49856;
}
outletDomB
{
type patch;
nFaces 100;
startFace 50032;
}
inletDomA
{
type patch;
nFaces 100;
startFace 50132;
}
outletDomA
{
type mixingPlane;
nFaces 100;
startFace 50232;
shadowPatch inletDomB;
zone outletDomAZone;
coordinateSystem
{
type cylindrical;
name mixingCS;
origin (0 0 0);
e1 (1 0 0);
e3 (0 0 1);
}
ribbonPatch
{
sweepAxis Theta;
stackAxis Z;
discretisation bothPatches;
}
}
inletDomB
{
type mixingPlane;
nFaces 100;
startFace 50332;
shadowPatch outletDomA;
zone inletDomBZone;
}
ggi1DomB
{
type ggi;
nFaces 128;
startFace 50432;
shadowPatch ggi2DomB;
zone ggi1DomBZone;
bridgeOverlap false;
}
ggi1DomA
{
type ggi;
nFaces 128;
startFace 50560;
shadowPatch ggi2DomA;
zone ggi1DomAZone;
bridgeOverlap false;
}
ggi2DomB
{
type ggi;
nFaces 176;
startFace 50688;
shadowPatch ggi1DomB;
zone ggi2DomBZone;
bridgeOverlap false;
}
ggi2DomA
{
type ggi;
nFaces 176;
startFace 50864;
shadowPatch ggi1DomA;
zone ggi2DomAZone;
bridgeOverlap false;
}
topAndBottom
{
type empty;
nFaces 50560;
startFace 51040;
}
cyclicDomA
{
type cyclic;
nFaces 212;
startFace 101600;
featureCos 0.9;
transform rotational;
rotationAxis (0 0 1);
rotationCentre (0 0 0);
rotationAngle 6;
}
cyclicDomB
{
type cyclic;
nFaces 188;
startFace 101812;
featureCos 0.9;
transform rotational;
rotationAxis (0 0 1);
rotationCentre (0 0 0);
rotationAngle 6;
}
)
// ************************************************************************* //

View file

@ -0,0 +1,200 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5.x |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object RASProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
RASModel laminar; //kEpsilon;
turbulence off; //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;
}
realizableKECoeffs
{
Cmu 0.09;
A0 4.0;
C2 1.9;
alphak 1;
alphaEps 0.833333;
}
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;
}
// ************************************************************************* //

View file

@ -0,0 +1,141 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext_mp |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "2/polyMesh";
object boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
13
(
bladeDomB
{
type wall;
nFaces 176;
startFace 49680;
}
bladeDomA
{
type wall;
nFaces 176;
startFace 49856;
}
outletDomB
{
type patch;
nFaces 100;
startFace 50032;
}
inletDomA
{
type patch;
nFaces 100;
startFace 50132;
}
outletDomA
{
type mixingPlane;
nFaces 100;
startFace 50232;
shadowPatch inletDomB;
zone outletDomAZone;
coordinateSystem
{
type cylindrical;
name mixingCS;
origin (0 0 0);
e1 (1 0 0);
e3 (0 0 1);
}
ribbonPatch
{
sweepAxis Theta;
stackAxis Z;
discretisation bothPatches;
}
}
inletDomB
{
type mixingPlane;
nFaces 100;
startFace 50332;
shadowPatch outletDomA;
zone inletDomBZone;
}
ggi1DomB
{
type ggi;
nFaces 128;
startFace 50432;
shadowPatch ggi2DomB;
zone ggi1DomBZone;
bridgeOverlap false;
}
ggi1DomA
{
type ggi;
nFaces 128;
startFace 50560;
shadowPatch ggi2DomA;
zone ggi1DomAZone;
bridgeOverlap false;
}
ggi2DomB
{
type ggi;
nFaces 176;
startFace 50688;
shadowPatch ggi1DomB;
zone ggi2DomBZone;
bridgeOverlap false;
}
ggi2DomA
{
type ggi;
nFaces 176;
startFace 50864;
shadowPatch ggi1DomA;
zone ggi2DomAZone;
bridgeOverlap false;
}
topAndBottom
{
type empty;
nFaces 50560;
startFace 51040;
}
cyclicDomA
{
type cyclic;
nFaces 212;
startFace 101600;
featureCos 0.9;
transform rotational;
rotationAxis (0 0 1);
rotationCentre (0 0 0);
rotationAngle 6;
}
cyclicDomB
{
type cyclic;
nFaces 188;
startFace 101812;
featureCos 0.9;
transform rotational;
rotationAxis (0 0 1);
rotationCentre (0 0 0);
rotationAngle 6;
}
)
// ************************************************************************* //

View file

@ -0,0 +1,12 @@
#
[Name]
Profil azimutal RZ
[Spatial Fields]
Z
[Data]
Z [ m ], Velocity Axial [ m s^-1 ], Velocity Radial [ m s^-1 ], Velocity Circumferential [ m s^-1 ], Pressure [ Pa ], Turbulence Kinetic Energy [ m^2 s^-2 ], Turbulence Eddy Dissipation [ m^2 s^-3 ]
-1.0 , 0.0000, -1.5 , -2.2 , 144.57, 0.000265875 , 0.0000712773
1.0 , 0.0000, -1.5 , -2.2 , 144.57, 0.000265875 , 0.0000712773
1 #
2
3 [Name]
4 Profil azimutal RZ
5 [Spatial Fields]
6 Z
7 [Data]
8 Z [ m ], Velocity Axial [ m s^-1 ], Velocity Radial [ m s^-1 ], Velocity Circumferential [ m s^-1 ], Pressure [ Pa ], Turbulence Kinetic Energy [ m^2 s^-2 ], Turbulence Eddy Dissipation [ m^2 s^-3 ]
9 -1.0 , 0.0000, -1.5 , -2.2 , 144.57, 0.000265875 , 0.0000712773
10 1.0 , 0.0000, -1.5 , -2.2 , 144.57, 0.000265875 , 0.0000712773

View file

@ -0,0 +1,37 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
transportModel Newtonian;
nu nu [0 2 -1 0 0 0 0] 0.01;
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;
}
// ************************************************************************* //

View file

@ -0,0 +1,194 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "2000";
object UFluxAvg;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
bladeDomB
{
type fixedValue;
value uniform (0 0 0);
}
bladeDomA
{
type fixedValue;
value uniform (0 0 0);
}
outletDomB
{
type inletOutlet;
inletValue uniform (0 0 0);
value uniform (0 0 0);
}
inletDomA
{
type fixedValue;
value nonuniform List<vector>
100
(
(2.27476 -1.38401 0)
(2.27329 -1.38642 0)
(2.27182 -1.38882 0)
(2.27035 -1.39123 0)
(2.26887 -1.39363 0)
(2.2674 -1.39604 0)
(2.26592 -1.39844 0)
(2.26443 -1.40083 0)
(2.26295 -1.40323 0)
(2.26146 -1.40563 0)
(2.25997 -1.40802 0)
(2.25848 -1.41041 0)
(2.25699 -1.4128 0)
(2.25549 -1.41519 0)
(2.25399 -1.41758 0)
(2.25249 -1.41997 0)
(2.25098 -1.42235 0)
(2.24947 -1.42474 0)
(2.24796 -1.42712 0)
(2.24645 -1.4295 0)
(2.24494 -1.43187 0)
(2.24342 -1.43425 0)
(2.2419 -1.43663 0)
(2.24037 -1.439 0)
(2.23885 -1.44137 0)
(2.23732 -1.44374 0)
(2.23579 -1.44611 0)
(2.23426 -1.44848 0)
(2.23272 -1.45084 0)
(2.23119 -1.45321 0)
(2.22965 -1.45557 0)
(2.2281 -1.45793 0)
(2.22656 -1.46029 0)
(2.22501 -1.46264 0)
(2.22346 -1.465 0)
(2.22191 -1.46735 0)
(2.22035 -1.46971 0)
(2.2188 -1.47206 0)
(2.21723 -1.4744 0)
(2.21567 -1.47675 0)
(2.21415 -1.47903 0)
(2.21267 -1.48125 0)
(2.21119 -1.48346 0)
(2.2097 -1.48567 0)
(2.20822 -1.48788 0)
(2.20673 -1.49008 0)
(2.20524 -1.49229 0)
(2.20374 -1.4945 0)
(2.20225 -1.4967 0)
(2.20075 -1.4989 0)
(2.19925 -1.5011 0)
(2.19775 -1.5033 0)
(2.19624 -1.5055 0)
(2.19474 -1.50769 0)
(2.19323 -1.50988 0)
(2.19172 -1.51208 0)
(2.1902 -1.51427 0)
(2.18869 -1.51646 0)
(2.18717 -1.51865 0)
(2.18565 -1.52083 0)
(2.18408 -1.52308 0)
(2.18247 -1.52539 0)
(2.18085 -1.5277 0)
(2.17923 -1.53001 0)
(2.17761 -1.53232 0)
(2.17599 -1.53462 0)
(2.17436 -1.53693 0)
(2.17273 -1.53923 0)
(2.1711 -1.54153 0)
(2.16947 -1.54383 0)
(2.16783 -1.54613 0)
(2.16619 -1.54842 0)
(2.16455 -1.55071 0)
(2.16291 -1.553 0)
(2.16126 -1.55529 0)
(2.15962 -1.55758 0)
(2.15796 -1.55987 0)
(2.15631 -1.56215 0)
(2.15466 -1.56444 0)
(2.153 -1.56672 0)
(2.15134 -1.569 0)
(2.14967 -1.57127 0)
(2.14801 -1.57355 0)
(2.14634 -1.57582 0)
(2.14467 -1.57809 0)
(2.143 -1.58036 0)
(2.14133 -1.58263 0)
(2.13965 -1.5849 0)
(2.13797 -1.58716 0)
(2.13629 -1.58943 0)
(2.1346 -1.59169 0)
(2.13292 -1.59395 0)
(2.13123 -1.59621 0)
(2.12953 -1.59846 0)
(2.12784 -1.60072 0)
(2.12614 -1.60297 0)
(2.12445 -1.60522 0)
(2.12274 -1.60747 0)
(2.12104 -1.60972 0)
(2.11933 -1.61196 0)
)
;
}
outletDomA
{
type mixingPlane;
fluxAveraging true;
surfaceAveraging false;
value uniform (0 0 0);
}
inletDomB
{
type mixingPlane;
fluxAveraging false;
surfaceAveraging false;
value uniform (0 0 0);
}
ggi1DomB
{
type ggi;
}
ggi1DomA
{
type ggi;
}
ggi2DomB
{
type ggi;
}
ggi2DomA
{
type ggi;
}
topAndBottom
{
type empty;
}
cyclicDomA
{
type cyclic;
}
cyclicDomB
{
type cyclic;
}
}
// ************************************************************************* //

View file

@ -0,0 +1,86 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "2000";
object pFluxAvg;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
bladeDomB
{
type zeroGradient;
}
bladeDomA
{
type zeroGradient;
}
outletDomB
{
type fixedValue;
value uniform 0;
}
inletDomA
{
type zeroGradient;
}
outletDomA
{
type mixingPlane;
fluxAveraging false;
surfaceAveraging false;
value uniform 0;
}
inletDomB
{
type mixingPlane;
fluxAveraging false;
surfaceAveraging false;
value uniform 0;
}
ggi1DomB
{
type ggi;
}
ggi1DomA
{
type ggi;
}
ggi2DomB
{
type ggi;
}
ggi2DomA
{
type ggi;
}
topAndBottom
{
type empty;
}
cyclicDomA
{
type cyclic;
}
cyclicDomB
{
type cyclic;
}
}
// ************************************************************************* //

View file

@ -0,0 +1,7 @@
faceSet ggi1DomBZone new patchToFace ggi1DomB
faceSet ggi1DomAZone new patchToFace ggi1DomA
faceSet ggi2DomBZone new patchToFace ggi2DomB
faceSet ggi2DomAZone new patchToFace ggi2DomA
faceSet outletDomAZone new patchToFace outletDomA
faceSet inletDomBZone new patchToFace inletDomB
quit

View file

@ -0,0 +1,65 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / 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;
class dictionary;
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application simpleFoam;
startFrom startTime;
startTime 0;
stopAt endTime;
endTime 2000;
deltaT 1;
writeControl timeStep;
writeInterval 100;
cycleWrite 0;
writeFormat ascii;
writePrecision 6;
writeCompression compressed;
timeFormat general;
timePrecision 6;
runTimeModifiable yes;
functions
(
mixingPlaneCheck
{
// Type of functionObject
type mixingPlaneCheck;
phi phi;
// Where to load it from (if not already in solver)
functionObjectLibs ("libcheckFunctionObjects.so");
}
);
// ************************************************************************* //

View file

@ -0,0 +1,111 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object createPatchDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// This application/dictionary controls:
// - optional: create new patches from boundary faces (either given as
// a set of patches or as a faceSet)
// - always: order faces on coupled patches such that they are opposite. This
// is done for all coupled faces, not just for any patches created.
// - optional: synchronise points on coupled patches.
// 1. Create cyclic:
// - specify where the faces should come from
// - specify the type of cyclic. If a rotational specify the rotationAxis
// and centre to make matching easier
// - pointSync true to guarantee points to line up.
// 2. Correct incorrect cyclic:
// This will usually fail upon loading:
// "face 0 area does not match neighbour 2 by 0.0100005%"
// " -- possible face ordering problem."
// - change patch type from 'cyclic' to 'patch' in the polyMesh/boundary file.
// - loosen match tolerance to get case to load
// - regenerate cyclic as above
// Tolerance used in matching faces. Absolute tolerance is span of
// face times this factor. To load incorrectly matches meshes set this
// to a higher value.
matchTolerance 1E-3;
// Do a synchronisation of coupled points after creation of any patches.
// Note: this does not work with points that are on multiple coupled patches
// with transformations.
pointSync false;
// Patches to create.
patchInfo
(
{
// Name of new patch
name cyclicDomA;
// Type of new patch
dictionary
{
type cyclic;
// Optional: explicitly set transformation tensor.
// Used when matching and synchronising points.
//transform translational;
//separationVector (-2289 0 0);
transform rotational;
rotationAngle 6.0;
rotationAxis (0 0 1);
rotationCentre (0 0 0);
}
// How to construct: either from 'patches' or 'set'
constructFrom patches;
// If constructFrom = patches : names of patches. Wildcards allowed.
patches (per1DomA per2DomA);
// If constructFrom = set : name of faceSet
set f0;
}
{
// Name of new patch
name cyclicDomB;
// Type of new patch
dictionary
{
type cyclic;
// Optional: explicitly set transformation tensor.
// Used when matching and synchronising points.
//transform translational;
//separationVector (-2289 0 0);
transform rotational;
rotationAngle 6.0;
rotationAxis (0 0 1);
rotationCentre (0 0 0);
}
// How to construct: either from 'patches' or 'set'
constructFrom patches;
// If constructFrom = patches : names of patches. Wildcards allowed.
patches (per1DomB per2DomB);
// If constructFrom = set : name of faceSet
set f0;
}
);
// ************************************************************************* //

View file

@ -0,0 +1,39 @@
// -*- C++ -*-
// File generated by PyFoam - sorry for the ugliness
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object decomposeParDict;
}
numberOfSubdomains 4;
method manual; //metis; //simple;
globalFaceZones (
ggi1DomBZone
ggi1DomAZone
ggi2DomBZone
ggi2DomAZone
);
manualCoeffs
{
dataFile "manualDecomp.txt";
}
simpleCoeffs
{
n (1 4 1);
delta 0.001;
}
distributed no;
roots
(
);

View file

@ -0,0 +1,37 @@
// -*- C++ -*-
// File generated by PyFoam - sorry for the ugliness
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object decomposeParDict;
}
numberOfSubdomains 4;
method manual; //metis; //simple;
globalFaceZones (
ggi1DomAZone
ggi2DomAZone
);
manualCoeffs
{
dataFile "manualDecomp.txt";
}
simpleCoeffs
{
n (1 4 1);
delta 0.001;
}
distributed no;
roots
(
);

View file

@ -0,0 +1,35 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object faceSetDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// Name of set to operate on
name outletDomAZone;
// One of clear/new/invert/add/delete|subset/list
action new;
// Actions to apply to pointSet. These are all the topoSetSource's ending
// in ..ToFace (see the meshTools library).
topoSetSources
(
// All faces of patch
patchToFace
{
name outletDomA; // Name of patch, regular expressions allowed
}
);
// ************************************************************************* //

View file

@ -0,0 +1,79 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / 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;
class dictionary;
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default steadyState; //Euler;
}
gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
}
divSchemes
{
default none;
div(phi,U) Gauss linearUpwind Gauss; //upwind;
div((nuEff*dev(grad(U).T()))) Gauss linear;
div(phi,epsilon) Gauss upwind;
div(phi,k) Gauss upwind;
}
laplacianSchemes
{
default none;
laplacian(nu,U) Gauss linear corrected;
laplacian(rAU,pcorr) Gauss linear corrected;
laplacian(rAU,p) Gauss linear corrected;
laplacian(nuEff,U) Gauss linear corrected;
laplacian((1|A(U)),p) Gauss linear corrected;
laplacian(DepsilonEff,epsilon) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(1,p) Gauss linear corrected;
}
interpolationSchemes
{
default linear;
interpolate(HbyA) linear;
interpolate(1|A) linear;
}
snGradSchemes
{
default corrected;
}
fluxRequired
{
default no;
pcorr;
p;
}
mixingPlane
{
default areaAveraging;
U areaAveraging;
p areaAveraging;
}
// ************************************************************************* //

View file

@ -0,0 +1,82 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / 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;
class dictionary;
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
p
{
solver CG;
preconditioner Cholesky;
minIter 0;
maxIter 1000;
tolerance 1e-07;
relTol 0.01;
}
U
{
solver BiCGStab;
preconditioner DILU;
minIter 0;
maxIter 1000;
tolerance 1e-07;
relTol 0;
}
epsilon
{
solver BiCGStab;
preconditioner DILU;
minIter 0;
maxIter 1000;
tolerance 1e-07;
relTol 0;
}
k
{
solver BiCGStab;
preconditioner DILU;
minIter 0;
maxIter 1000;
tolerance 1e-07;
relTol 0;
}
}
SIMPLE
{
nNonOrthogonalCorrectors 1;
pRefCell 0;
pRefValue 0;
convergence 1e-5;
}
relaxationFactors
{
p 0.3;
U 0.7;
k 0.7;
epsilon 0.7;
}
// ************************************************************************* //