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

release on 2014-06-02_18-24-30

Browse Source
This commit is contained in:
goniva
2014-06-02 18:24:30 +02:00
parent 235449636c
commit e43cdce2ff
189 changed files with 1923 additions and 1199 deletions
Download Patch File Download Diff File Expand all files Collapse all files

18
tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/0/U
View File

@ -55,24 +55,6 @@ boundaryField
uniformValue table
(
(0.000 (0 0 0.002))
/*(0.010 (0 0 0.002))
(0.011 (0 0 0.004))
(0.020 (0 0 0.004))
(0.021 (0 0 0.006))
(0.030 (0 0 0.006))
(0.031 (0 0 0.008))
(0.040 (0 0 0.008))
(0.041 (0 0 0.010))
(0.050 (0 0 0.010))
(0.051 (0 0 0.012))
(0.060 (0 0 0.012))
(0.061 (0 0 0.014))
(0.070 (0 0 0.014))
(0.071 (0 0 0.016))
(0.080 (0 0 0.016))
(0.081 (0 0 0.018))
(0.090 (0 0 0.018))
(0.091 (0 0 0.020))*/
(0.100 (0 0 0.020))
);
}

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

@ -61,7 +61,8 @@ forceModels
//Archimedes
//volWeightedAverage
//totalMomentumExchange
particleCellVolume
//particleCellVolume
//fieldTimeAverage
);
momCoupleModels
@ -107,7 +108,7 @@ gradPForceProps
densityFieldName "rho";
voidfractionFieldName "voidfraction";
velocityFieldName "U";
//interpolation;
interpolation;
}
viscForceProps
@ -139,8 +140,12 @@ totalMomentumExchangeProps
}
GidaspowDragProps
{
verbose;
velFieldName "U";
densityFieldName "rho";
voidfractionFieldName "voidfraction";
interpolation;
phi 1;
}
DiFeliceDragProps
{
@ -155,6 +160,7 @@ KochHillDragProps
velFieldName "U";
densityFieldName "rho";
voidfractionFieldName "voidfraction";
interpolation;
}
BeetstraDragProps
@ -198,7 +204,7 @@ twoWayFilesProps
centreProps
{
alphaMin 0.10;
alphaMin 0.1;
}
engineProps
@ -222,7 +228,6 @@ dividedProps
twoWayMPIProps
{
//maxNumberOfParticles 10100;
liggghtsPath "../DEM/in.liggghts_resume";
}
twoWayM2MProps

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

@ -25,7 +25,7 @@ FoamFile
liggghtsCommandModels
(
runLiggghts
runLiggghts
);
// ************************************************************************* //

2
tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/octave/totalPressureDrop.m
View File

@ -30,7 +30,7 @@ phip = 1 % sphericity
epsilon = 0.451335 % void fraction
Ustart = 0.002
Uend = 0.02
timeStepSize = 0.001; % time interval of pressure data
timeStepSize = 0.0005; % time interval of pressure data
Tstart = 0;
Tend = t_sim(length(t_sim));
deltaU=(Uend-Ustart)/((Tend-Tstart)/timeStepSize);

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

@ -25,7 +25,7 @@ stopAt endTime;
endTime 0.1;
deltaT 0.001;
deltaT 0.0005;
writeControl adjustableRunTime;

0
tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/system/funkySetFieldsDict Executable file → Normal file
View File

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

@ -41,7 +41,7 @@ solvers
relTol 0;
}
"(voidfraction|Us|Ksl|dSmoothing)"
"(voidfraction|Us|Ksl|dSmoothing|UsNext|voidfractionNext)"
{
solver PCG;
preconditioner DIC;

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

@ -38,7 +38,7 @@ fix cylwalls all wall/gran model hertz tangential history primitive type 1 zcyli
#particle distributions and insertion
region bc cylinder z 0.0 0.0 0.012 0. 0.055 units box
fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 200 radius constant 0.0005
fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 2000 radius constant 0.0005
fix pdd1 all particledistribution/discrete 1. 1 pts1 1.0
fix ins all insert/pack seed 100001 distributiontemplate pdd1 vel constant 0. 0. -1. insert_every once overlapcheck yes all_in yes particles_in_region 10000 region bc

2
tutorials/cfdemSolverPiso/ErgunTestMPI/DEM/in.liggghts_resume
View File

@ -65,7 +65,7 @@ compute_modify thermo_temp dynamic yes
#insert the first particles so that dump is not empty
dump myDump all stl 1 post/dump_*.stl
#run 1
dump dmp all custom 5000 ../DEM/post/dump*.liggghts_restart id type type x y z vx vy vz fx fy fz f_dragforce[1] f_dragforce[2] f_dragforce[3] radius
dump dmp all custom 5000 ../DEM/post/dump*.liggghts_restart id type type x y z vx vy vz fx fy fz f_dragforce[1] f_dragforce[2] f_dragforce[3] radius
undump myDump
run 1

BIN
tutorials/cfdemSolverPiso/ErgunTestMPI/DEM/liggghts.restart
View File

Binary file not shown.

1
tutorials/cfdemSolverPiso/ErgunTestMPI/DEM/post/dummy
View File

@ -1 +0,0 @@
dummyfile

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

@ -21,7 +21,7 @@ logfileName="log_$headerText"
solverName="cfdemSolverPiso"
nrProcs="4"
machineFileName="none" # yourMachinefileName | none
debugMode="off" # on | off| strictdebugMode="off" # on | off| strict
debugMode="off" # on | off| strict # on | off| strict
testHarnessPath="$CFDEM_TEST_HARNESS_PATH"
runOctave="true"
postproc="false"
@ -82,22 +82,10 @@ fi
#- clean up case
echo "deleting data at: $casePath :\n"
rm -r $casePath/CFD/0.*
rm -r $casePath/CFD/callgrind.*
rm -r $casePath/CFD/*.out
rm -r $casePath/CFD/log.*
rm -r $casePath/CFD/octave/octave-core
rm -r $casePath/CFD/VTK
rm -r $casePath/CFD/processor*
rm -r $casePath/CFD/couplingFiles/*
rm -r $casePath/DEM/post/*
rm -r $casePath/DEM/log.*
rm -r $casePath/DEM/liggghts.restartCFDEM*
rm -r $casePath/CFD/probes
rm -r $casePath/CFD/postProcessing
rm -r $casePath/CFD/lagrangian
source $WM_PROJECT_DIR/bin/tools/CleanFunctions
cd $casePath/CFD
cleanCase
rm -r $casePath/CFD/clockData
rm -r $casePath/DEM/post/*
(cd $casePath/DEM/post && touch dummy)
echo "done"
#- preserve post directory
echo "dummyfile" >> $casePath/DEM/post/dummy

1
tutorials/cfdemSolverPiso/ErgunTestMPI_cgs/CFD/constant/couplingProperties
View File

@ -130,6 +130,7 @@ DiFeliceDragProps
velFieldName "U";
densityFieldName "rho";
voidfractionFieldName "voidfraction";
granVelFieldName "Us";
}
KochHillDragProps

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

@ -97,22 +97,12 @@ if [ $postproc == "true" ]
fi
#- clean up case
rm -rf $casePath/CFD/0.*
rm -r $casePath/CFD/callgrind.*
rm -r $casePath/CFD/*.out
echo "deleting data at: $casePath :\n"
source $WM_PROJECT_DIR/bin/tools/CleanFunctions
cd $casePath/CFD
cleanCase
rm -r $casePath/CFD/clockData
rm -rf $casePath/CFD/processor*
rm -r $casePath/CFD/VTK
rm -rf $casePath/CFD/patchAverage_pressureDrop
rm -rf $casePath/CFD/probes
rm -rf $casePath/CFD/postProcessing
rm -rf $casePath/CFD/particles
rm -r $casePath/CFD/log.*
rm -r $casePath/CFD/lagrangian
rm $casePath/log.liggghts
rm -r $casePath/DEM/post/*
rm $casePath/DEM/liggghts.restartCFDEM*
rm $casePath/DEM/post/dump.*
rm -r $casePath/DEM/log.*
#- preserve post directory
echo "dummyfile" >> $casePath/DEM/post/dummy
(cd $casePath/DEM/post && touch dummy)
echo "done"

6
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/CFD/0/U
View File

@ -55,7 +55,7 @@ boundaryField
uniformValue table
(
(0.000 (0 0 0.002))
(0.010 (0 0 0.002))
/*(0.010 (0 0 0.002))
(0.011 (0 0 0.004))
(0.020 (0 0 0.004))
(0.021 (0 0 0.006))
@ -72,9 +72,9 @@ boundaryField
(0.080 (0 0 0.016))
(0.081 (0 0 0.018))
(0.090 (0 0 0.018))
(0.091 (0 0 0.020))
(0.091 (0 0 0.020))*/
(0.100 (0 0 0.020))
);
);
}
outlet

3
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/CFD/0/p
View File

@ -35,8 +35,9 @@ boundaryField
outlet
{
//type zeroGradient;
type fixedValue;
value uniform 100000;
value uniform 1.0e5;
}
}

78
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/CFD/constant/couplingProperties_init
View File

@ -25,41 +25,43 @@ FoamFile
//===========================================================================//
// sub-models & settings
verbous;
modelType "A"; // A or B
couplingInterval 100;
voidFractionModel divided;
voidFractionModel divided;//centre;//
locateModel engine;
locateModel engine;//turboEngineM2M;//
meshMotionModel noMeshMotion;
regionModel allRegion;
IOModel "basicIO";
IOModel basicIO;
probeModel off;
dataExchangeModel twoWayMPI;//twoWayFiles;//oneWayVTK;//
dataExchangeModel twoWayMPI;//twoWayM2M;//twoWayFiles;//oneWayVTK;//
averagingModel dense;//dilute;//
clockModel off;
clockModel standardClock;//off;
smoothingModel off;
smoothingModel off;// localPSizeDiffSmoothing;// constDiffSmoothing; //
forceModels
(
//GidaspowDrag
DiFeliceDrag
//KochHillDrag
//BeetstraDrag
//DiFeliceDrag
KochHillDrag
gradPForce
viscForce
//Archimedes
//volWeightedAverage
//totalMomentumExchange
//particleCellVolume
);
momCoupleModels
@ -72,6 +74,21 @@ turbulenceModelType "RASProperties";//"LESProperties";//
//===========================================================================//
// sub-model properties
localPSizeDiffSmoothingProps
{
lowerLimit 0.1;
upperLimit 1e10;
dSmoothingLength 1.5e-3;
Csmoothing 1.0;
}
constDiffSmoothingProps
{
lowerLimit 0.1;
upperLimit 1e10;
smoothingLength 1.5e-3;
}
implicitCoupleProps
{
velFieldName "U";
@ -110,7 +127,7 @@ volWeightedAverageProps
);
upperThreshold 0.999;
lowerThreshold 0;
verbous;
verbose;
}
totalMomentumExchangeProps
{
@ -122,8 +139,11 @@ totalMomentumExchangeProps
}
GidaspowDragProps
{
verbose;
velFieldName "U";
densityFieldName "rho";
voidfractionFieldName "voidfraction";
phi 1;
}
DiFeliceDragProps
{
@ -134,8 +154,25 @@ DiFeliceDragProps
KochHillDragProps
{
//verbose;
velFieldName "U";
densityFieldName "rho";
voidfractionFieldName "voidfraction";
}
BeetstraDragProps
{
velFieldName "U";
densityFieldName "rho";
gravityFieldName "g";
rhoParticle 2000.;
voidfractionFieldName "voidfraction";
interpolation ;
useFilteredDragModel ;
useParcelSizeDependentFilteredDrag ;
k 0.05;
aLimit 0.0;
// verbose ;
}
virtualMassForceProps
@ -144,6 +181,13 @@ virtualMassForceProps
densityFieldName "rho";
}
particleCellVolumeProps
{
upperThreshold 0.999;
lowerThreshold 0.;
verbose;
}
oneWayVTKProps
{
couplingFilename "vtk_out%4.4d.vtk";
@ -165,13 +209,25 @@ engineProps
treeSearch true;
}
turboEngineM2MProps
{
turboEngineProps
{
treeSearch true;
}
}
dividedProps
{
alphaMin 0.1;
alphaMin 0.01;
scaleUpVol 1.0;
}
twoWayMPIProps
{
liggghtsPath "../DEM/in.liggghts_resume";
}
twoWayM2MProps
{
maxNumberOfParticles 10100;
liggghtsPath "../DEM/in.liggghts_resume";

74
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/CFD/constant/couplingProperties_restart
View File

@ -25,41 +25,43 @@ FoamFile
//===========================================================================//
// sub-models & settings
verbous;
modelType "A"; // A or B
couplingInterval 100;
voidFractionModel divided;
voidFractionModel divided;//centre;//
locateModel engine;
locateModel engine;//turboEngineM2M;//
meshMotionModel noMeshMotion;
regionModel allRegion;
IOModel "basicIO";
IOModel basicIO;
probeModel off;
dataExchangeModel twoWayMPI;//twoWayFiles;//oneWayVTK;//
dataExchangeModel twoWayMPI;//twoWayM2M;//twoWayFiles;//oneWayVTK;//
averagingModel dense;//dilute;//
clockModel off;
clockModel standardClock;//off;//
smoothingModel off;
smoothingModel off;// localPSizeDiffSmoothing;// constDiffSmoothing; //
forceModels
(
//GidaspowDrag
DiFeliceDrag
//KochHillDrag
//BeetstraDrag
//DiFeliceDrag
KochHillDrag
gradPForce
viscForce
//Archimedes
//volWeightedAverage
//totalMomentumExchange
//particleCellVolume
);
momCoupleModels
@ -72,6 +74,21 @@ turbulenceModelType "RASProperties";//"LESProperties";//
//===========================================================================//
// sub-model properties
localPSizeDiffSmoothingProps
{
lowerLimit 0.1;
upperLimit 1e10;
dSmoothingLength 1.5e-3;
Csmoothing 1.0;
}
constDiffSmoothingProps
{
lowerLimit 0.1;
upperLimit 1e10;
smoothingLength 1.5e-3;
}
implicitCoupleProps
{
velFieldName "U";
@ -110,7 +127,7 @@ volWeightedAverageProps
);
upperThreshold 0.999;
lowerThreshold 0;
verbous;
verbose;
}
totalMomentumExchangeProps
{
@ -122,8 +139,11 @@ totalMomentumExchangeProps
}
GidaspowDragProps
{
verbose;
velFieldName "U";
densityFieldName "rho";
voidfractionFieldName "voidfraction";
phi 1;
}
DiFeliceDragProps
{
@ -134,8 +154,25 @@ DiFeliceDragProps
KochHillDragProps
{
//verbose;
velFieldName "U";
densityFieldName "rho";
voidfractionFieldName "voidfraction";
}
BeetstraDragProps
{
velFieldName "U";
densityFieldName "rho";
gravityFieldName "g";
rhoParticle 2000.;
voidfractionFieldName "voidfraction";
interpolation ;
useFilteredDragModel ;
useParcelSizeDependentFilteredDrag ;
k 0.05;
aLimit 0.0;
// verbose ;
}
virtualMassForceProps
@ -144,6 +181,13 @@ virtualMassForceProps
densityFieldName "rho";
}
particleCellVolumeProps
{
upperThreshold 0.999;
lowerThreshold 0.;
verbose;
}
oneWayVTKProps
{
couplingFilename "vtk_out%4.4d.vtk";
@ -165,9 +209,17 @@ engineProps
treeSearch true;
}
turboEngineM2MProps
{
turboEngineProps
{
treeSearch true;
}
}
dividedProps
{
alphaMin 0.1;
alphaMin 0.01;
scaleUpVol 1.0;
}

4
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/CFD/constant/liggghtsCommands_init
View File

@ -29,6 +29,10 @@ liggghtsCommandModels
writeLiggghts
);
runLiggghtsProps
{
preNo false;
}
//- optional
writeLiggghtsProps

6
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/CFD/constant/liggghtsCommands_restart
View File

@ -27,5 +27,9 @@ liggghtsCommandModels
(
runLiggghts
);
// ************************************************************************* //
/*runLiggghtsProps
{
preNo false;
}*/

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

@ -23,13 +23,13 @@ startTime 0.05;
stopAt endTime;
endTime 0.1;//0.01;
endTime 0.1;
deltaT 0.001;
writeControl adjustableRunTime;
writeInterval 0.01;//0.01;
writeInterval 0.01;
purgeWrite 0;
@ -49,7 +49,7 @@ adjustTimeStep no;
maxCo 0.1;
//libs ( "libgroovyBC.so" );
//libs ( "libgroovyBC.so" "libfiniteVolumeCFDEM.so");
functions
(
@ -87,7 +87,7 @@ functions
);
// Fields to be probed
fields ( p U voidfraction volAverage_voidfraction);
fields ( p U voidfraction volAverage_voidfraction voidfractionNext voidfractionPrev);
// Write at same frequency as fields
outputControl timeStep;//outputTime;

8
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/CFD/system/controlDict_init
View File

@ -23,13 +23,13 @@ startTime 0;
stopAt endTime;
endTime 0.05;//0.01;
endTime 0.05;
deltaT 0.001;
writeControl adjustableRunTime;
writeInterval 0.01;//0.01;
writeInterval 0.01;
purgeWrite 0;
@ -49,7 +49,7 @@ adjustTimeStep no;
maxCo 0.1;
//libs ( "libgroovyBC.so" );
//libs ( "libgroovyBC.so" "libfiniteVolumeCFDEM.so");
functions
(
@ -87,7 +87,7 @@ functions
);
// Fields to be probed
fields ( p U voidfraction volAverage_voidfraction);
fields ( p U voidfraction volAverage_voidfraction voidfractionNext voidfractionPrev);
// Write at same frequency as fields
outputControl timeStep;//outputTime;

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

@ -23,13 +23,13 @@ startTime 0.05;
stopAt endTime;
endTime 0.1;//0.01;
endTime 0.1;
deltaT 0.001;
writeControl adjustableRunTime;
writeInterval 0.01;//0.01;
writeInterval 0.01;
purgeWrite 0;
@ -49,7 +49,7 @@ adjustTimeStep no;
maxCo 0.1;
//libs ( "libgroovyBC.so" );
//libs ( "libgroovyBC.so" "libfiniteVolumeCFDEM.so");
functions
(
@ -87,7 +87,7 @@ functions
);
// Fields to be probed
fields ( p U voidfraction volAverage_voidfraction);
fields ( p U voidfraction volAverage_voidfraction voidfractionNext voidfractionPrev);
// Write at same frequency as fields
outputControl timeStep;//outputTime;

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

@ -17,7 +17,7 @@ FoamFile
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
numberOfSubdomains 2;
numberOfSubdomains 4;
//- Keep owner and neighbour on same processor for faces in zones:
// preserveFaceZones (heater solid1 solid3);
@ -30,7 +30,7 @@ numberOfSubdomains 2;
simpleCoeffs
{
n (1 1 2);
n (2 2 1);
delta 0.001;
}

0
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/CFD/system/funkySetFieldsDict Executable file → Normal file
View File

34
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/CFD/system/fvSolution
View File

@ -17,7 +17,7 @@ FoamFile
solvers
{
p
"(p)"
{
solver PCG;
preconditioner DIC;
@ -33,7 +33,7 @@ solvers
relTol 0;
}
U
"(U|k|epsilon|R|nuTilda)"
{
solver PBiCG;
preconditioner DILU;
@ -41,34 +41,10 @@ solvers
relTol 0;
}
k
"(voidfraction|Us|Ksl|dSmoothing)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0;
}
epsilon
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0;
}
R
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0;
}
nuTilda
{
solver PBiCG;
preconditioner DILU;
solver PCG;
preconditioner DIC;
tolerance 1e-05;
relTol 0;
}

8
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/DEM/in.liggghts_init
View File

@ -1,5 +1,5 @@
# Pour granular particles into chute container, then induce flow
echo both
echo both
atom_style granular
atom_modify map array
communicate single vel yes
@ -36,12 +36,9 @@ fix zwalls1 all wall/gran model hertz tangential history primitive type 1 zplan
fix zwalls2 all wall/gran model hertz tangential history primitive type 1 zplane 0.0553
fix cylwalls all wall/gran model hertz tangential history primitive type 1 zcylinder 0.01385 0. 0.
#cfd coupling
#fix cfd all couple/cfd
#particle distributions and insertion
region bc cylinder z 0.0 0.0 0.012 0. 0.055 units box
fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 200 radius constant 0.0005
fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 2000 radius constant 0.0005
fix pdd1 all particledistribution/discrete 1. 1 pts1 1.0
fix ins all insert/pack seed 100001 distributiontemplate pdd1 vel constant 0. 0. -1. insert_every once overlapcheck yes all_in yes particles_in_region 10000 region bc
@ -49,7 +46,6 @@ fix ins all insert/pack seed 100001 distributiontemplate pdd1 vel constant 0. 0
#apply nve integration to all particles that are inserted as single particles
fix integr all nve/sphere
#screen output
compute 1 all erotate/sphere
thermo_style custom step atoms ke c_1 vol

13
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/DEM/in.liggghts_restart
View File

@ -1,5 +1,5 @@
# Pour granular particles into chute container, then induce flow
echo both
atom_style granular
atom_modify map array
communicate single vel yes
@ -8,10 +8,9 @@ boundary m m m
newton off
units si
processors 2 1 1
processors 2 2 1
#read the restart file
#read_restart ../DEM/liggghts.restartCFDEM
read_restart ../DEM/liggghts.restartCFDEM_0.050000
#do not do this here, the simulation box is in the restart file!
@ -36,8 +35,6 @@ pair_coeff * *
timestep 0.00001
fix gravi all gravity 9.81 vector 0.0 0.0 -1.0
#walls
#walls
fix zwalls1 all wall/gran model hertz tangential history primitive type 1 zplane 0.0
fix zwalls2 all wall/gran model hertz tangential history primitive type 1 zplane 0.0553
fix cylwalls all wall/gran model hertz tangential history primitive type 1 zcylinder 0.01385 0. 0.
@ -56,11 +53,11 @@ fix integr all nve/sphere
compute centerOfMass all com
#compute total dragforce
#compute dragtotal all reduce sum f_dragforce[1] f_dragforce[2] f_dragforce[3]
compute dragtotal all reduce sum f_dragforce[1] f_dragforce[2] f_dragforce[3]
#screen output
compute 1 all erotate/sphere
thermo_style custom step atoms ke c_1 vol c_centerOfMass[3] #c_dragtotal[1] c_dragtotal[2] c_dragtotal[3]
thermo_style custom step atoms ke c_1 vol c_centerOfMass[3] c_dragtotal[1] c_dragtotal[2] c_dragtotal[3]
thermo 10
thermo_modify lost ignore norm no
compute_modify thermo_temp dynamic yes
@ -68,7 +65,7 @@ compute_modify thermo_temp dynamic yes
#insert the first particles so that dump is not empty
dump myDump all stl 1 post/dump_*.stl
#run 1
dump dmp all custom 5000 ../DEM/post/dump.liggghts_restart id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
dump dmp all custom 5000 ../DEM/post/dump*.liggghts_restart id type type x y z vx vy vz fx fy fz f_dragforce[1] f_dragforce[2] f_dragforce[3] radius
undump myDump
run 1

9
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/DEM/in.liggghts_resume
View File

@ -8,7 +8,7 @@ boundary m m m
newton off
units si
processors 2 1 1
processors 2 2 1
#read the restart file
read_restart ../DEM/liggghts.restart
@ -35,7 +35,6 @@ pair_coeff * *
timestep 0.00001
fix gravi all gravity 9.81 vector 0.0 0.0 -1.0
#walls
fix zwalls1 all wall/gran model hertz tangential history primitive type 1 zplane 0.0
fix zwalls2 all wall/gran model hertz tangential history primitive type 1 zplane 0.0553
fix cylwalls all wall/gran model hertz tangential history primitive type 1 zcylinder 0.01385 0. 0.
@ -54,11 +53,11 @@ fix integr all nve/sphere
compute centerOfMass all com
#compute total dragforce
#compute dragtotal all reduce sum f_dragforce[1] f_dragforce[2] f_dragforce[3]
compute dragtotal all reduce sum f_dragforce[1] f_dragforce[2] f_dragforce[3]
#screen output
compute 1 all erotate/sphere
thermo_style custom step atoms ke c_1 vol c_centerOfMass[3] #c_dragtotal[1] c_dragtotal[2] c_dragtotal[3]
thermo_style custom step atoms ke c_1 vol c_centerOfMass[3] c_dragtotal[1] c_dragtotal[2] c_dragtotal[3]
thermo 10
thermo_modify lost ignore norm no
compute_modify thermo_temp dynamic yes
@ -66,7 +65,7 @@ compute_modify thermo_temp dynamic yes
#insert the first particles so that dump is not empty
dump myDump all stl 1 post/dump_*.stl
#run 1
dump dmp all custom 5000 ../DEM/post/dump.liggghts_restart id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
dump dmp all custom 5000 ../DEM/post/dump*.liggghts_restart id type type x y z vx vy vz fx fy fz f_dragforce[1] f_dragforce[2] f_dragforce[3] radius
undump myDump
run 1

BIN
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/DEM/liggghts.restart
View File

Binary file not shown.

60
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/DEM/post/dummy
View File

@ -1,60 +0,0 @@
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile
dummyfile

5
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/parCFDDEMrun.sh
View File

@ -19,12 +19,13 @@ logpath=$casePath
headerText="run_parallel_cfdemSolverPiso_ErgunTestMPI_CFDDEM"
logfileName="log_$headerText"
solverName="cfdemSolverPiso"
nrProcs="2"
nrProcs="4"
machineFileName="none" # yourMachinefileName | none
debugMode="off" # on | off| strict
reconstuctCase="true" # true | false
testHarnessPath="$CFDEM_TEST_HARNESS_PATH"
#--------------------------------------------------------------------------------#
#- call function to run a parallel CFD-DEM case
parCFDDEMrun $logpath $logfileName $casePath $headerText $solverName $nrProcs $machineFileName $debugMode
parCFDDEMrun $logpath $logfileName $casePath $headerText $solverName $nrProcs $machineFileName $debugMode $reconstuctCase

8
tutorials/cfdemSolverPiso/settlingTestMPI/CFD/0/Ksl
View File

@ -21,15 +21,11 @@ internalField uniform 0;
boundaryField
{
atmosphere
"(inlet|outlet|walls)"
{
type zeroGradient;
value uniform 0;
}
walls
{
type zeroGradient;
}
}

7
tutorials/cfdemSolverPiso/settlingTestMPI/CFD/0/U
View File

@ -21,7 +21,12 @@ internalField uniform (0 0 0);
boundaryField
{
atmosphere
inlet
{
type fixedValue;
value uniform (0 0 0);
}
outlet
{
type zeroGradient;
value uniform (0 0 0);

7
tutorials/cfdemSolverPiso/settlingTestMPI/CFD/0/Us
View File

@ -21,13 +21,10 @@ internalField uniform (0 0 0);
boundaryField
{
atmosphere
{
type zeroGradient;
}
walls
"(inlet|outlet|walls)"
{
type zeroGradient;
value uniform (0 0 0);
}
}

8
tutorials/cfdemSolverPiso/settlingTestMPI/CFD/0/k
View File

@ -26,12 +26,18 @@ boundaryField
value uniform 0;
}
atmosphere
inlet
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}
outlet
{
type outletInlet;
outletValue uniform 0;
value uniform 0;
}
}
// ************************************************************************* //

7
tutorials/cfdemSolverPiso/settlingTestMPI/CFD/0/nuSgs
View File

@ -20,12 +20,7 @@ internalField uniform 0;
boundaryField
{
walls
{
type zeroGradient;
}
atmosphere
".*"
{
type zeroGradient;
}

12
tutorials/cfdemSolverPiso/settlingTestMPI/CFD/0/p
View File

@ -25,18 +25,20 @@ boundaryField
type zeroGradient;
value uniform 0;
}
atmosphere
inlet
{
/* type totalPressure;
type zeroGradient;
value uniform 0;
}
outlet
{
type totalPressure;
p0 uniform 0;
U U;
phi phi;
rho rho;
psi none;
gamma 1;
value uniform 0;*/
type fixedValue;
value uniform 0;
}
}

10
tutorials/cfdemSolverPiso/settlingTestMPI/CFD/0/rho
View File

@ -20,16 +20,10 @@ internalField uniform 10;
boundaryField
{
walls
"(inlet|outlet|walls)"
{
type zeroGradient;
value uniform 1;
}
atmosphere
{
type zeroGradient;
value uniform 1;
value uniform 10;
}
}

13
tutorials/cfdemSolverPiso/settlingTestMPI/CFD/0/voidfraction
View File

@ -20,19 +20,10 @@ internalField uniform 1;
boundaryField
{
walls
{
type zeroGradient;
}
obstacle
{
type zeroGradient;
}
atmosphere
"(inlet|outlet|walls)"
{
type zeroGradient;
value uniform 1;
}
}

6
tutorials/cfdemSolverPiso/settlingTestMPI/CFD/constant/couplingProperties
View File

@ -46,9 +46,9 @@ dataExchangeModel twoWayMPI;//twoWayFiles;//oneWayVTK;//
averagingModel dense;//dilute;//
clockModel standardClock;//off; //
clockModel off; //standardClock;//
smoothingModel off; //constDiffSmoothing; //
smoothingModel constDiffSmoothing; //off;//
//useDDTvoidfraction;
@ -93,6 +93,8 @@ DiFeliceDragProps
velFieldName "U";
densityFieldName "rho";
voidfractionFieldName "voidfraction";
granVelFieldName "Us";
verbose;
}
SchillerNaumannDragProps

13
tutorials/cfdemSolverPiso/settlingTestMPI/CFD/constant/polyMesh/blockMeshDict
View File

@ -30,8 +30,10 @@ vertices
blocks
(
//hex (0 1 2 3 4 5 6 7) (5 25 5) simpleGrading (1 1 1) // regular grid
hex (0 1 2 3 4 5 6 7) (25 25 25) simpleGrading (1 1 1) // irrregular grid
hex (0 1 2 3 4 5 6 7) (5 25 5) simpleGrading (1 1 1) // regular grid
//hex (0 1 2 3 4 5 6 7) (25 125 25) simpleGrading (1 1 1) // irrregular fine grid
//hex (0 1 2 3 4 5 6 7) (25 25 25) simpleGrading (1 1 1) // irrregular coarse grid
//hex (0 1 2 3 4 5 6 7) (10 25 10) simpleGrading (1 1 1) // irrregular grid (ratio 2)
);
edges
@ -40,15 +42,18 @@ edges
patches
(
patch atmosphere
patch inlet
(
(3 7 6 2)
)
patch outlet
(
(1 5 4 0)
)
wall walls
(
(0 4 7 3)
(2 6 5 1)
(1 5 4 0)
(0 3 2 1)
(4 5 6 7)
)

41
tutorials/cfdemSolverPiso/settlingTestMPI/CFD/system/fvSolution
View File

@ -33,7 +33,7 @@ solvers
relTol 0;
}
U
"(U|k|epsilon|R|nuTilda)"
{
solver PBiCG;
preconditioner DILU;
@ -41,50 +41,13 @@ solvers
relTol 0;
}
k
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0;
}
epsilon
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0;
}
R
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0;
}
nuTilda
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0;
}
voidfraction
"(voidfraction|Ksl|UsNext|voidfractionNext)"
{
solver PCG;
preconditioner DIC;
tolerance 1e-09;
relTol 1e-06;
}
Ksl
{
$voidfraction
}
}
PISO

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

@ -29,7 +29,7 @@ postproc="false"
#--------------------------------------------------------------------------------#
#- call function to run a parallel CFD-DEM case
parCFDDEMrun $logpath $logfileName $casePath $headerText $solverName $nrProcs $machineFileName $debugMode
parCFDDEMrun $logpath $logfileName $casePath $headerText $solverName $nrProcs $machineFileName $debugMode "true"
if [ $runOctave == "true" ]
then