ICODE-ADC/CFDEMcoupling-PUBLIC
4
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

release on 2017-11-30_15-39-15

Browse Source
This commit is contained in:
Christoph Goniva
2017-11-30 15:39:15 +01:00
parent d03549ba88
commit 5cb9190126
266 changed files with 2325 additions and 6827 deletions
Download Patch File Download Diff File Expand all files Collapse all files

0
tutorials/ParScale/cfdemSolverPisoSTM/packedBedTempParScale/CFD/constant/polyMesh/blockMeshDict → tutorials/ParScale/cfdemSolverPisoSTM/packedBedTempParScale/CFD/system/blockMeshDict
View File

168
tutorials/ParScale/cfdemSolverPisoSTM/packedBedTempParScale/CFD/system/controlDict
View File

@ -52,168 +52,6 @@ maxCo 0.1;
maxDeltaT 1;
//libs ( "libgroovyBC.so" );
functions
(
probes
{
type probes;
// Where to load it from
functionObjectLibs ( "libsampling.so" );
// Name of the directory for probe data
name probes;
probeLocations
(
(0.05 0.05 -0.499)
(0.05 0.05 1.499)
);
// Fields to be probed
fields ( p T );
// Write at same frequency as fields
outputControl timeStep;//outputTime;
outputInterval 1;
}
volFlow_inlet
{
type faceSource;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;//timeStep;
log true;
writeArea true; //write for OF 3.0 to be consistent
valueOutput true;
source patch;
sourceName inlet;
operation sum;
fields
(
phi
);
surfaceFormat null;
}
volFlow_outlet
{
type faceSource;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;//timeStep;
log true;
writeArea true; //write for OF 3.0 to be consistent
valueOutput true;
source patch;
sourceName outlet;
operation sum;
fields
(
phi
);
surfaceFormat null;
}
volFlow_wall
{
type faceSource;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;//timeStep;
log true;
writeArea true; //write for OF 3.0 to be consistent
valueOutput true;
source patch;
sourceName wall;
operation sum;
fields
(
phi
);
surfaceFormat null;
}
/*
// simpleFunctionObjects
heatFlux
{
type patchHeatFlux;
functionObjectLibs
(
"libsimpleFunctionObjects.so"
);
verbose true;
patches
(
inlet
outlet
);
fields // name of temp field
(
T
);
cp 1007; // cp in [J/(kg*K)]
factor 1.188; // density for incomp solvers!
}
massFlux
{
type patchMassFlow;
functionObjectLibs
(
"libsimpleFunctionObjects.so"
);
verbose true;
patches
(
inlet
outlet
);
factor 1.188; // density for incomp solvers!
}
pressureDrop
{
type patchAverage;
functionObjectLibs
(
"libsimpleFunctionObjects.so"
);
verbose true;
patches
(
inlet
outlet
);
fields
(
p
);
factor 1;
}
T
{
type patchAverage;
functionObjectLibs
(
"libsimpleFunctionObjects.so"
);
verbose true;
patches
(
inlet
outlet
);
fields // name of temp field
(
T
);
factor 1;
}*/
);
// ************************************************************************* //
//#include "libsAndFunctionsOF3x"
//#include "libsAndFunctionsOF4x"
#include "libsAndFunctionsOF5x"

164
tutorials/ParScale/cfdemSolverPisoSTM/packedBedTempParScale/CFD/system/libsAndFunctionsOF3x Normal file
View File

@ -0,0 +1,164 @@
functions
(
probes
{
type probes;
// Where to load it from
functionObjectLibs ( "libsampling.so" );
// Name of the directory for probe data
name probes;
probeLocations
(
(0.05 0.05 -0.499)
(0.05 0.05 1.499)
);
// Fields to be probed
fields ( p T );
// Write at same frequency as fields
outputControl timeStep;//outputTime;
outputInterval 1;
}
volFlow_inlet
{
type faceSource;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;//timeStep;
log true;
writeArea true; //write for OF 3.0 to be consistent
valueOutput true;
source patch;
sourceName inlet;
operation sum;
fields
(
phi
);
surfaceFormat null;
}
volFlow_outlet
{
type faceSource;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;//timeStep;
log true;
writeArea true; //write for OF 3.0 to be consistent
valueOutput true;
source patch;
sourceName outlet;
operation sum;
fields
(
phi
);
surfaceFormat null;
}
volFlow_wall
{
type faceSource;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;//timeStep;
log true;
writeArea true; //write for OF 3.0 to be consistent
valueOutput true;
source patch;
sourceName wall;
operation sum;
fields
(
phi
);
surfaceFormat null;
}
/*
// simpleFunctionObjects
heatFlux
{
type patchHeatFlux;
functionObjectLibs
(
"libsimpleFunctionObjects.so"
);
verbose true;
patches
(
inlet
outlet
);
fields // name of temp field
(
T
);
cp 1007; // cp in [J/(kg*K)]
factor 1.188; // density for incomp solvers!
}
massFlux
{
type patchMassFlow;
functionObjectLibs
(
"libsimpleFunctionObjects.so"
);
verbose true;
patches
(
inlet
outlet
);
factor 1.188; // density for incomp solvers!
}
pressureDrop
{
type patchAverage;
functionObjectLibs
(
"libsimpleFunctionObjects.so"
);
verbose true;
patches
(
inlet
outlet
);
fields
(
p
);
factor 1;
}
T
{
type patchAverage;
functionObjectLibs
(
"libsimpleFunctionObjects.so"
);
verbose true;
patches
(
inlet
outlet
);
fields // name of temp field
(
T
);
factor 1;
}*/
);
// ************************************************************************* //

170
tutorials/ParScale/cfdemSolverPisoSTM/packedBedTempParScale/CFD/system/libsAndFunctionsOF4x Normal file
View File

@ -0,0 +1,170 @@
functions
{
probes
{
libs ( "libsampling.so" );
type probes;
name probes;
writeControl timeStep; //outputTime;
writeInterval 1;
fields ( p T );
probeLocations
(
(0.05 0.05 -0.499)
(0.05 0.05 1.499)
);
}
volFlow_inlet
{
type surfaceRegion;
functionObjectLibs ("libfieldFunctionObjects.so");
writeControl writeTime;
log yes;
writeFields true; //write for OF 3.0 to be consistent
regionType patch;
name inlet;
operation sum;
valueOutput true;
fields
(
phi
);
surfaceFormat null;
}
volFlow_outlet
{
type surfaceRegion;
functionObjectLibs ("libfieldFunctionObjects.so");
writeControl writeTime;
log yes;
writeFields true; //write for OF 3.0 to be consistent
regionType patch;
name outlet;
operation sum;
valueOutput true;
fields
(
phi
);
surfaceFormat null;
}
volFlow_wall
{
type surfaceRegion;
functionObjectLibs ("libfieldFunctionObjects.so");
writeControl writeTime;
log yes;
writeFields true; //write for OF 3.0 to be consistent
regionType patch;
name wall;
operation sum;
valueOutput true;
fields
(
phi
);
surfaceFormat null;
}
// simpleFunctionObjects
/*heatFlux
{
type patchHeatFlux;
functionObjectLibs
(
"libsimpleFunctionObjects.so"
);
verbose true;
outputControlMode timeStep;
outputInterval 1;
writeStartTime no;
patches
(
inlet
outlet
);
fields // name of temp field
(
T
);
cp 1007; // cp in [J/(kg*K)]
factor 1.188; // density for incomp solvers!
}*/
massFlux
{
type patchMassFlow;
functionObjectLibs
(
"libsimpleFunctionObjects.so"
);
verbose true;
outputControlMode timeStep;
outputInterval 1;
writeStartTime no;
patches
(
inlet
outlet
);
factor 1.188; // density for incomp solvers!
}
pressureDrop
{
type patchAverage;
functionObjectLibs
(
"libsimpleFunctionObjects.so"
);
verbose true;
outputControlMode timeStep;
outputInterval 1;
writeStartTime no;
patches
(
inlet
outlet
);
fields
(
p
);
factor 1;
}
T
{
type patchAverage;
functionObjectLibs
(
"libsimpleFunctionObjects.so"
);
verbose true;
outputControlMode timeStep;
outputInterval 1;
writeStartTime no;
patches
(
inlet
outlet
);
fields // name of temp field
(
T
);
factor 1;
}
};
// ************************************************************************* //

170
tutorials/ParScale/cfdemSolverPisoSTM/packedBedTempParScale/CFD/system/libsAndFunctionsOF5x Normal file
View File

@ -0,0 +1,170 @@
functions
{
probes
{
libs ( "libsampling.so" );
type probes;
name probes;
writeControl timeStep; //outputTime;
writeInterval 1;
fields ( p T );
probeLocations
(
(0.05 0.05 -0.499)
(0.05 0.05 1.499)
);
}
volFlow_inlet
{
type surfaceFieldValue;
functionObjectLibs ("libfieldFunctionObjects.so");
writeControl writeTime;
log yes;
writeFields true; //write for OF 3.0 to be consistent
regionType patch;
name inlet;
operation sum;
valueOutput true;
fields
(
phi
);
surfaceFormat null;
}
volFlow_outlet
{
type surfaceFieldValue;
functionObjectLibs ("libfieldFunctionObjects.so");
writeControl writeTime;
log yes;
writeFields true; //write for OF 3.0 to be consistent
regionType patch;
name outlet;
operation sum;
valueOutput true;
fields
(
phi
);
surfaceFormat null;
}
volFlow_wall
{
type surfaceFieldValue;
functionObjectLibs ("libfieldFunctionObjects.so");
writeControl writeTime;
log yes;
writeFields true; //write for OF 3.0 to be consistent
regionType patch;
name wall;
operation sum;
valueOutput true;
fields
(
phi
);
surfaceFormat null;
}
// simpleFunctionObjects
/*heatFlux
{
type patchHeatFlux;
functionObjectLibs
(
"libsimpleFunctionObjects.so"
);
verbose true;
outputControlMode timeStep;
outputInterval 1;
writeStartTime no;
patches
(
inlet
outlet
);
fields // name of temp field
(
T
);
cp 1007; // cp in [J/(kg*K)]
factor 1.188; // density for incomp solvers!
}*/
massFlux
{
type patchMassFlow;
functionObjectLibs
(
"libsimpleFunctionObjects.so"
);
verbose true;
outputControlMode timeStep;
outputInterval 1;
writeStartTime no;
patches
(
inlet
outlet
);
factor 1.188; // density for incomp solvers!
}
pressureDrop
{
type patchAverage;
functionObjectLibs
(
"libsimpleFunctionObjects.so"
);
verbose true;
outputControlMode timeStep;
outputInterval 1;
writeStartTime no;
patches
(
inlet
outlet
);
fields
(
p
);
factor 1;
}
T
{
type patchAverage;
functionObjectLibs
(
"libsimpleFunctionObjects.so"
);
verbose true;
outputControlMode timeStep;
outputInterval 1;
writeStartTime no;
patches
(
inlet
outlet
);
fields // name of temp field
(
T
);
factor 1;
}
};
// ************************************************************************* //

2
tutorials/ParScale/cfdemSolverPisoSTM/packedBedTempParScale/parCFDDEMrun.sh
View File

@ -42,7 +42,7 @@ if [ $runOctave == "true" ]
rm *.png
#- run octave
octave totalPressureDropAndNusselt.m
octave --no-gui totalPressureDropAndNusselt.m
#- show plots
eog cfdemSolverPisoSTM_Nusselt.png &

0
tutorials/cfdemPostproc/fillCylinder/CFD/constant/polyMesh/blockMeshDict → tutorials/cfdemPostproc/fillCylinder/CFD/system/blockMeshDict
View File

12
tutorials/cfdemPostproc/fillCylinder/CFD/system/controlDict
View File

@ -52,7 +52,7 @@ maxCo 0.1;
// ************************************************************************* //
functions
(
{
probes
{
@ -63,10 +63,10 @@ functions
name probes;
probeLocations
(
(0.00003 0 0.0001)
(0.00003 0 0.01)
(0.00003 0 0.02)
(0.00003 0 0.0529)
(0.00003 1e-4 0.0001)
(0.00003 1e-4 0.01)
(0.00003 1e-4 0.02)
(0.00003 1e-4 0.0529)
);
// Fields to be probed
@ -76,4 +76,4 @@ functions
outputControl timeStep;//outputTime;
outputInterval 1;
}
);
}

0
tutorials/cfdemSolverIB/twoSpheresGlowinskiMPI/CFD/constant/polyMesh/blockMeshDict → tutorials/cfdemSolverIB/twoSpheresGlowinskiMPI/CFD/system/blockMeshDict
View File

6
tutorials/cfdemSolverIB/twoSpheresGlowinskiMPI/CFD/system/controlDict
View File

@ -56,14 +56,14 @@ maxDeltaT 1;
functions
(
{
probes
{
type probes;
// Where to load it from
functionObjectLibs ( "libsampling.so" );
// Name of the directory for probe data
name drag;
name probes;
probeLocations
(
(0.5001 0.5001 0.0001)
@ -79,6 +79,6 @@ functions
outputControl timeStep;//outputTime;
outputInterval 1;
}
);
}
// ************************************************************************* //

2
tutorials/cfdemSolverIB/twoSpheresGlowinskiMPI/parCFDDEMrun.sh
View File

@ -34,7 +34,7 @@ if [ $runOctave == "true" ]
then
cd $casePath/CFD/octave
octave postproc.m
octave --no-gui postproc.m
evince pos_y_two_part_rec_glow.eps
evince vel_y_two_part_rec_glow.eps
#display pos_y_two_part_rec_glow.png &

15
tutorials/cfdemSolverPiso/ErgunTestMPI/Allrun.sh
View File

@ -27,6 +27,7 @@ fi
# keep old couplingProperties
cp $casePath/CFD/constant/couplingProperties $casePath/CFD/constant/couplingProperties_backup
cp $casePath/CFD/0/p $casePath/CFD/0/p_backup
# change to M2M + type B
#changeDictionary -constant -dict changeDicts/changeDictionaryDict_1 -case $casePath/CFD
@ -43,9 +44,21 @@ cp $casePath/CFD/constant/couplingProperties $casePath/CFD/constant/couplingProp
# change to MPI + type Bfull
#changeDictionary -constant -dict changeDicts/changeDictionaryDict_5 -case $casePath/CFD
# change to cfdemSolverPimpleImExFace + MPI + type Bfull + explicit gradP
#changeDictionary -constant -dict changeDicts/changeDictionaryDict_6 -case $casePath/CFD
# change to M2M + type B + fixedFluxPressure
#changeDictionary -dict changeDicts/changeDictionaryDict_0 -case $casePath/CFD
#changeDictionary -constant -dict changeDicts/changeDictionaryDict_1 -case $casePath/CFD
# change to MPI + type A + fixedFluxPressure
#changeDictionary -dict changeDicts/changeDictionaryDict_0 -case $casePath/CFD
#changeDictionary -constant -dict changeDicts/changeDictionaryDict_2 -case $casePath/CFD
#- run parallel CFD-DEM in new terminal
#gnome-terminal --title='cfdemSolverPiso ErgunTestMPI CFD' -e "bash $casePath/parCFDDEMrun.sh"
. $casePath/parCFDDEMrun.sh
# restore old couplingProperties
mv $casePath/CFD/constant/couplingProperties_backup $casePath/CFD/constant/couplingProperties
mv $casePath/CFD/constant/couplingProperties_backup $casePath/CFD/constant/couplingProperties
mv $casePath/CFD/0/p_backup $casePath/CFD/0/p

16
tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/0/k
View File

@ -17,27 +17,27 @@ FoamFile
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0.001;
internalField uniform 1e-6;
boundaryField
{
wall
{
type kqRWallFunction;
value uniform 0.001;
value $internalField;
}
outlet
{
type inletOutlet;
inletValue uniform 0.001;
value uniform 0.001;
inletValue $internalField;
value $internalField;
}
inlet
{
/*type turbulentIntensityKineticEnergyInlet;
intensity 0.05;
value uniform 0.001;*/
type zeroGradient;
type turbulentIntensityKineticEnergyInlet;
intensity 0.01;
value $internalField;
//type zeroGradient;
}
}

11
tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/0/nut
View File

@ -23,21 +23,18 @@ boundaryField
{
wall
{
type nutWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0;
type nutkWallFunction;
value $internalField;
}
outlet
{
type calculated;
value uniform 0;
value $internalField;
}
inlet
{
type calculated;
value uniform 0;
value $internalField;
}
}

2
tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/constant/couplingProperties
View File

@ -257,7 +257,7 @@ centreProps
alphaMin 0.1;
}
"(engineProps|turboEngineM2MProps)"
"(engineProps|turboEngineProps|turboEngineM2MProps)"
{
treeSearch true;
}

11
tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/constant/turbulenceProperties
View File

@ -16,7 +16,16 @@ FoamFile
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
//simulationType RASModel; //OFversion24x
simulationType laminar; //OFversion30x
simulationType laminar;//RAS; //OFversion30x
RAS
{
RASModel kEpsilon;//kOmegaSST;
turbulence on;
printCoeffs on;
}
// ************************************************************************* //

0
tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/constant/polyMesh/blockMeshDict → tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/system/blockMeshDict
View File

75
tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/system/controlDict
View File

@ -57,8 +57,7 @@ maxCo 0.1;
);*/
functions
(
{
probes
{
type probes;
@ -73,91 +72,65 @@ functions
);
// Fields to be probed
fields ( p U voidfraction volAverage_voidfraction voidfractionNext voidfractionPrev);
fields ( p U voidfraction volAverage_voidfraction);
// Write at same frequency as fields
outputControl timeStep;//outputTime;
outputInterval 1;
}
/*pressureDrop
{
type patchAverage;
functionObjectLibs
(
"libsimpleFunctionObjects.so"
);
verbose true;
patches
(
inlet
outlet
);
fields
(
p
);
factor 1;
}*/
volFlow_inlet
{
type faceSource;
type surfaceFieldValue;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;//timeStep;
outputControl timeStep;
outputInterval 1;
log true;
writeArea true; //write for OF 3.0 to be consistent
valueOutput true;
source patch;
sourceName inlet;
writeFields false;
valueOutput false;
regionType patch;
name inlet;
operation sum;
fields
(
phi
);
surfaceFormat null;
}
volFlow_outlet
{
type faceSource;
type surfaceFieldValue;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;//timeStep;
outputControl timeStep;
outputInterval 1;
log true;
writeArea true; //write for OF 3.0 to be consistent
valueOutput true;
source patch;
sourceName outlet;
writeFields false;
valueOutput false;
regionType patch;
name outlet;
operation sum;
fields
(
phi
);
surfaceFormat null;
}
volFlow_wall
{
type faceSource;
type surfaceFieldValue;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;//timeStep;
outputControl timeStep;
outputInterval 1;
log true;
writeArea true; //write for OF 3.0 to be consistent
valueOutput true;
source patch;
sourceName wall;
writeFields false;
valueOutput false;
regionType patch;
name wall;
operation sum;
fields
(
phi
);
surfaceFormat null;
}
);
}
// ************************************************************************* //

6
tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/system/fvSchemes
View File

@ -52,5 +52,11 @@ fluxRequired
p ;
}
// kOmegaSST
wallDist
{
method meshWave;
}
// ************************************************************************* //

4
tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/system/fvSolution
View File

@ -31,7 +31,7 @@ solvers
relTol 0;
}
"(U|k|epsilon|R|nuTilda)"
"(U|k|epsilon|omega|R|nuTilda)"
{
solver PBiCG;
preconditioner DILU;
@ -39,7 +39,7 @@ solvers
relTol 0.1;
}
"(U|k|epsilon|R|nuTilda)Final"
"(U|k|epsilon|omega|R|nuTilda)Final"
{
$U;
relTol 0;

2
tutorials/cfdemSolverPiso/ErgunTestMPI/parCFDDEMrun.sh
View File

@ -43,7 +43,7 @@ if [ $runOctave == "true" ]
rm cfdemSolverPiso_ErgunTestMPI.png
#- run octave
octave totalPressureDrop.m
octave --no-gui totalPressureDrop.m
#- show plot
eog cfdemSolverPiso_ErgunTestMPI.png

0
tutorials/cfdemSolverPiso/ErgunTestMPI_cgs/CFD/constant/polyMesh/blockMeshDict → tutorials/cfdemSolverPiso/ErgunTestMPI_cgs/CFD/system/blockMeshDict
View File

4
tutorials/cfdemSolverPiso/ErgunTestMPI_cgs/CFD/system/controlDict
View File

@ -149,7 +149,7 @@ DimensionedConstants
}
functions
(
{
probes
{
@ -210,5 +210,5 @@ functions
);
factor 1;
}*/
);
}
// ************************************************************************* //

2
tutorials/cfdemSolverPiso/ErgunTestMPI_cgs/DEM/in.liggghts_init
View File

@ -17,7 +17,7 @@ neigh_modify delay 0
# Material properties required for granular pair styles
soft_particles yes
fix m1 all property/global youngsModulus peratomtype 5.e6
fix m2 all property/global poissonsRatio peratomtype 0.45
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.3

2
tutorials/cfdemSolverPiso/ErgunTestMPI_cgs/DEM/in.liggghts_run
View File

@ -19,7 +19,7 @@ neighbor 0.1 bin
neigh_modify delay 0
# Material properties required for granular pair styles
soft_particles yes
fix m1 all property/global youngsModulus peratomtype 5.e6
fix m2 all property/global poissonsRatio peratomtype 0.45
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.3

2
tutorials/cfdemSolverPiso/ErgunTestMPI_cgs/parCFDDEMrun.sh
View File

@ -43,7 +43,7 @@ if [ $runOctave == "true" ]
rm cfdemSolverPiso_ErgunTestMPI.eps
#- run octave
octave totalPressureDrop.m
octave --no-gui totalPressureDrop.m
#- show plot
evince cfdemSolverPiso_ErgunTestMPI.eps

2
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/Allrun.sh
View File

@ -69,7 +69,7 @@ if [ $runOctave == "true" ]
rm cfdemSolverPiso_ErgunTestMPI.eps
#- run octave
octave totalPressureDrop.m
octave --no-gui totalPressureDrop.m
#- show plot
evince cfdemSolverPiso_ErgunTestMPI.eps

0
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/CFD/constant/polyMesh/blockMeshDict → tutorials/cfdemSolverPiso/ErgunTestMPI_restart/CFD/system/blockMeshDict
View File

4
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/CFD/system/controlDict
View File

@ -52,7 +52,7 @@ maxCo 0.1;
//libs ( "libgroovyBC.so" "libfiniteVolumeCFDEM.so");
functions
(
{
probes
{
@ -113,5 +113,5 @@ functions
);
factor 1;
}*/
);
}
// ************************************************************************* //

4
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/CFD/system/controlDict_restart
View File

@ -52,7 +52,7 @@ maxCo 0.1;
//libs ( "libgroovyBC.so" "libfiniteVolumeCFDEM.so");
functions
(
{
probes
{
@ -113,5 +113,5 @@ functions
);
factor 1;
}*/
);
}
// ************************************************************************* //

4
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/CFD/system/controlDict_run
View File

@ -52,7 +52,7 @@ maxCo 0.1;
//libs ( "libgroovyBC.so" "libfiniteVolumeCFDEM.so");
functions
(
{
probes
{
@ -113,5 +113,5 @@ functions
);
factor 1;
}*/
);
}
// ************************************************************************* //

14
tutorials/cfdemSolverPiso/periodicChannel/CFD/octave/checkVolFlow.m
View File

@ -13,7 +13,7 @@ Ainlet=1.000000e-02; % cross sectional area in m2
Utarget=10; % target superficial velocity in m/s
Vtarget=Utarget*Ainlet; % target flowrate in m3/s
filename = 'faceSource.dat';
filename = 'surfaceFieldValue.dat';
timedirectory = '0';
%Daten einlesen
@ -21,21 +21,21 @@ timedirectory = '0';
data = transpose(load(strcat(path,'/',timedirectory,'/',filename))); %data
[x,y]=size(data);
time=data(1,1+skip:y)
volFlow_inlet=data(3,1+skip:y)
volFlow_inlet=data(2,1+skip:y)
%Daten einlesen
path = '../../CFD/postProcessing/volFlow_outlet';
data = transpose(load(strcat(path,'/',timedirectory,'/',filename))); %data
[x,y]=size(data)
time=data(1,1+skip:y);
volFlow_outlet=data(3,1+skip:y);
volFlow_outlet=data(2,1+skip:y);
%Daten einlesen
path = '../../CFD/postProcessing/volFlow_wall';
data = transpose(load(strcat(path,'/',timedirectory,'/',filename))); %data
[x,y]=size(data)
time=data(1,1+skip:y);
volFlow_wall=data(3,1+skip:y);
volFlow_wall=data(2,1+skip:y);
%- time integrate data to vol entered per TS
deltatT=time(2)-time(1)
@ -65,7 +65,7 @@ figure1 = figure('PaperPositionMode','manual','PaperUnits','centimeters',...
% Create axes
axes1 = axes('Parent',figure1,'YGrid','on','XGrid','on','LineWidth',1,...
'FontWeight','light','FontSize',11,'FontName','Helvetica-Narrow');
'FontWeight','normal','FontSize',11,'FontName','Helvetica-Narrow');
box(axes1,'on');
hold(axes1,'all');
@ -97,7 +97,7 @@ figure2 = figure('PaperPositionMode','manual','PaperUnits','centimeters',...
% Create axes
axes2 = axes('Parent',figure2,'YGrid','on','XGrid','on','LineWidth',1,...
'FontWeight','light','FontSize',11,'FontName','Helvetica-Narrow');
'FontWeight','normal','FontSize',11,'FontName','Helvetica-Narrow');
box(axes2,'on');
hold(axes2,'all');
@ -129,7 +129,7 @@ figure3 = figure('PaperPositionMode','manual','PaperUnits','centimeters',...
% Create axes
axes3 = axes('Parent',figure3,'YGrid','on','XGrid','on','LineWidth',1,...
'FontWeight','light','FontSize',11,'FontName','Helvetica-Narrow');
'FontWeight','normal','FontSize',11,'FontName','Helvetica-Narrow');
box(axes3,'on');
hold(axes3,'all');

0
tutorials/cfdemSolverPiso/periodicChannel/CFD/constant/polyMesh/blockMeshDict → tutorials/cfdemSolverPiso/periodicChannel/CFD/system/blockMeshDict
View File

63
tutorials/cfdemSolverPiso/periodicChannel/CFD/system/controlDict
View File

@ -52,7 +52,7 @@ maxCo 0.1;
libs ( "libfvOptionsCFDEM.so" );
functions
(
{
probes
{
@ -63,8 +63,8 @@ functions
name probes;
probeLocations
(
(0.00003 0.05 0.05)
(0.39999 0.05 0.05)
(0.00003 0.050001 0.050001)
(0.39999 0.050001 0.050001)
);
// Fields to be probed
@ -75,65 +75,58 @@ functions
outputInterval 1;
}
volFlow_inlet
{
type faceSource;
type surfaceFieldValue;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;//timeStep;
log true;
writeArea true; //write for OF 3.0 to be consistent
valueOutput true;
source patch;
sourceName inlet;
outputControl timeStep;
outputInterval 1;
log false;
writeFields false;
valueOutput false;
regionType patch;
name inlet;
operation sum;
fields
(
phi
);
surfaceFormat null;
}
volFlow_outlet
{
type faceSource;
type surfaceFieldValue;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;//timeStep;
log true;
writeArea true; //write for OF 3.0 to be consistent
valueOutput true;
source patch;
sourceName outlet;
outputControl timeStep;
outputInterval 1;
log false;
writeFields false;
valueOutput false;
regionType patch;
name outlet;
operation sum;
fields
(
phi
);
surfaceFormat null;
}
volFlow_wall
{
type faceSource;
type surfaceFieldValue;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;//timeStep;
log true;
writeArea true; //write for OF 3.0 to be consistent
valueOutput true;
source patch;
sourceName wall;
outputControl timeStep;
outputInterval 1;
log false;
writeFields false;
valueOutput false;
regionType patch;
name wall;
operation sum;
fields
(
phi
);
surfaceFormat null;
}
);
}
// ************************************************************************* //

4
tutorials/cfdemSolverPiso/periodicChannel/CFD/system/fvOptions
View File

@ -24,9 +24,11 @@ momentumSource
{
selectionMode all; // apply force to all cells
fieldNames (U); // name of velocity field
fields (U); // name of velocity field
Ubar (10 0 0);// desired mean velocity
relaxation 0.9; // optional relaxation factor
twoPhase false;
coupled true;
voidfractionField "voidfraction";
}
}

2
tutorials/cfdemSolverPiso/periodicChannel/parCFDDEMrun.sh
View File

@ -44,7 +44,7 @@ if [ $runOctave == "true" ]
rm *.png
#- run octave
octave checkVolFlow.m
octave --no-gui checkVolFlow.m
#- show plot
eog volflow.png

30
tutorials/cfdemSolverPiso/settlingTestMPI/CFD/0/f
View File

@ -1,30 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object f;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -2 -2 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
"(inlet|outlet|walls)"
{
type zeroGradient;
value uniform (0 0 0);
}
}
// ************************************************************************* //

0
tutorials/cfdemSolverPiso/settlingTestMPI/CFD/constant/polyMesh/blockMeshDict → tutorials/cfdemSolverPiso/settlingTestMPI/CFD/system/blockMeshDict
View File

2
tutorials/cfdemSolverPiso/settlingTestMPI/parCFDDEMrun.sh
View File

@ -40,7 +40,7 @@ if [ $runOctave == "true" ]
rm cfdemSolverPiso_settlingTestMPI.png
#- run octave
octave settlingVelocity.m
octave --no-gui settlingVelocity.m
#- show plot
eog cfdemSolverPiso_settlingTestMPI.png

0
tutorials/cfdemSolverPiso/voidfractionTest/CFD/constant/polyMesh/blockMeshDict → tutorials/cfdemSolverPiso/voidfractionTest/CFD/system/blockMeshDict
View File

4
tutorials/cfdemSolverPiso/voidfractionTest/CFD/system/controlDict
View File

@ -52,7 +52,7 @@ maxCo 1;
maxDeltaT 1;
functions
(
{
probes
{
type probes;
@ -72,6 +72,6 @@ functions
outputControl timeStep;//outputTime;
outputInterval 1;
}
);
}
// ************************************************************************* //

2
tutorials/cfdemSolverPiso/voidfractionTest/parCFDDEMrun.sh
View File

@ -41,7 +41,7 @@ if [ $runOctave == "true" ]
rm averageDDTvoidfraction.png
#- run octave
octave particleCellVolume.m
octave --no-gui particleCellVolume.m
#- show plot
eog averageDDTvoidfraction.png

0
tutorials/cfdemSolverPisoSTM/packedBedTemp/CFD/constant/polyMesh/blockMeshDict → tutorials/cfdemSolverPisoSTM/packedBedTemp/CFD/system/blockMeshDict
View File

58
tutorials/cfdemSolverPisoSTM/packedBedTemp/CFD/system/controlDict
View File

@ -54,7 +54,7 @@ maxDeltaT 1;
//libs ( "libgroovyBC.so" );
functions
(
{
probes
{
@ -79,62 +79,56 @@ functions
volFlow_inlet
{
type faceSource;
type surfaceFieldValue;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;//timeStep;
log true;
writeArea true; //write for OF 3.0 to be consistent
valueOutput true;
source patch;
sourceName inlet;
outputControl timeStep;
outputInterval 1;
log false;
writeFields false;
valueOutput false;
regionType patch;
name inlet;
operation sum;
fields
(
phi
);
surfaceFormat null;
}
volFlow_outlet
{
type faceSource;
type surfaceFieldValue;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;//timeStep;
log true;
writeArea true; //write for OF 3.0 to be consistent
valueOutput true;
source patch;
sourceName outlet;
outputControl timeStep;
outputInterval 1;
log false;
writeFields false;
valueOutput false;
regionType patch;
name outlet;
operation sum;
fields
(
phi
);
surfaceFormat null;
}
volFlow_wall
{
type faceSource;
type surfaceFieldValue;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;//timeStep;
log true;
writeArea true; //write for OF 3.0 to be consistent
valueOutput true;
source patch;
sourceName wall;
outputControl timeStep;
outputInterval 1;
log false;
writeFields false;
valueOutput false;
regionType patch;
name wall;
operation sum;
fields
(
phi
);
surfaceFormat null;
}
/*
@ -215,5 +209,5 @@ functions
);
factor 1;
}*/
);
}
// ************************************************************************* //

2
tutorials/cfdemSolverPisoSTM/packedBedTemp/parCFDDEMrun.sh
View File

@ -42,7 +42,7 @@ if [ $runOctave == "true" ]
rm *.png
#- run octave
octave totalPressureDropAndNusselt.m
octave --no-gui totalPressureDropAndNusselt.m
#- show plots
eog cfdemSolverPisoSTM_Nusselt.png &

0
tutorials/cfdemSolverPisoScalar/packedBedTemp/CFD/constant/polyMesh/blockMeshDict → tutorials/cfdemSolverPisoScalar/packedBedTemp/CFD/system/blockMeshDict
View File

58
tutorials/cfdemSolverPisoScalar/packedBedTemp/CFD/system/controlDict
View File

@ -54,7 +54,7 @@ maxDeltaT 1;
//libs ( "libgroovyBC.so" );
functions
(
{
probes
{
@ -79,62 +79,56 @@ functions
volFlow_inlet
{
type faceSource;
type surfaceFieldValue;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;//timeStep;
log true;
writeArea true; //write for OF 3.0 to be consistent
valueOutput true;
source patch;
sourceName inlet;
outputControl timeStep;
outputInterval 1;
log false;
writeFields false;
valueOutput false;
regionType patch;
name inlet;
operation sum;
fields
(
phi
);
surfaceFormat null;
}
volFlow_outlet
{
type faceSource;
type surfaceFieldValue;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;//timeStep;
log true;
writeArea true; //write for OF 3.0 to be consistent
valueOutput true;
source patch;
sourceName outlet;
outputControl timeStep;
outputInterval 1;
log false;
writeFields false;
valueOutput false;
regionType patch;
name outlet;
operation sum;
fields
(
phi
);
surfaceFormat null;
}
volFlow_wall
{
type faceSource;
type surfaceFieldValue;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;//timeStep;
log true;
writeArea true; //write for OF 3.0 to be consistent
valueOutput true;
source patch;
sourceName wall;
outputControl timeStep;
outputInterval 1;
log false;
writeFields false;
valueOutput false;
regionType patch;
name wall;
operation sum;
fields
(
phi
);
surfaceFormat null;
}
/*
@ -215,5 +209,5 @@ functions
);
factor 1;
}*/
);
}
// ************************************************************************* //

2
tutorials/cfdemSolverPisoScalar/packedBedTemp/parCFDDEMrun.sh
View File

@ -42,7 +42,7 @@ if [ $runOctave == "true" ]
rm *.png
#- run octave
octave totalPressureDropAndNusselt.m
octave --no-gui totalPressureDropAndNusselt.m
#- show plots
eog cfdemSolverPisoScalar_Nusselt.png &