zhyang-dev/CFDEMcoupling-PFM
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

release on 2012-11-21_14-15-25

Browse Source
This commit is contained in:
cfdem
2012-11-21 14:15:25 +01:00
parent 5b13e00f72
commit 5d51595f96
24 changed files with 248 additions and 37 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@ -32,7 +32,7 @@ couplingInterval 100;
voidFractionModel divided;
locateModel engine;
locateModel engine;//
meshMotionModel noMeshMotion;
@ -49,8 +49,8 @@ clockModel standardClock;//off;
forceModels
(
//GidaspowDrag
//DiFeliceDrag
KochHillDrag
DiFeliceDrag
//KochHillDrag
gradPForce
viscForce
//Archimedes
@ -169,6 +169,11 @@ dividedProps
}
twoWayMPIProps
{
//maxNumberOfParticles 10100;
liggghtsPath "../DEM/in.liggghts_resume";
}
twoWayM2MProps
{
maxNumberOfParticles 10100;
liggghtsPath "../DEM/in.liggghts_resume";

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

@ -28,3 +28,8 @@ liggghtsCommandModels
runLiggghts
);
// ************************************************************************* //
runLiggghtsProps
{
preNo false;
}

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

@ -55,9 +55,20 @@ fprintf('so the result does not depend on density\n')
%==================================
rhoP = 2000 % particle density in kg/m3
g = 9.81 % gravity m/s2
Umf = dp^2*(rhoP-rhoG)*g/(150*muG)*(epsilon^3*phip^2)/(1-epsilon)
ReMF = Umf*dp*rhoG/muG % must be <20 !!!
%Umf = sqrt(phip*dp^2/1.75*(rhoP-rhoG)/rhoG*g*epsilon^3) % Re>1000
Umf = dp^2*(rhoP-rhoG)*g/(150*muG)*(epsilon^3*phip^2)/(1-epsilon);
ReMF = Umf*dp*rhoG/muG;
if(ReMF<20)
fprintf('applying eqn1 for Umf.\n')
elseif(ReMF>20 && ReMF<1000)
fprintf('applying eqn1 for Umf.\n')
elseif (ReMF>=1000)
fprintf('applying eqn2 for Umf.\n')
Umf = sqrt(dp*(rhoP-rhoG)*g/(1.75*rhoG)*epsilon^3*phip);
ReMF = Umf*dp*rhoG/muG;
end
Umf
ReMF
dpUmf= L * (
150*((1-epsilon)^2/epsilon^3)*((muG.*Umf)/(phip*dp)^2)

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

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

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

@ -8,7 +8,7 @@ boundary m m m
newton off
units si
processors 2 2 2
processors 2 2 1
#read the restart file
read_restart ../DEM/liggghts.restart

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

@ -19,7 +19,7 @@ logpath=$casePath
headerText="run_parallel_cfdemSolverPiso_ErgunTestMPI_CFDDEM"
logfileName="log_$headerText"
solverName="cfdemSolverPiso"
nrProcs="8"
nrProcs="4"
machineFileName="none" # yourMachinefileName | none
debugMode="off" # on | off | prof
testHarnessPath="$CFDEM_TEST_HARNESS_PATH"

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

@ -25,3 +25,4 @@ dummyfile
dummyfile
dummyfile
dummyfile
dummyfile

17
tutorials/cfdemSolverPisoMS/ErgunTestMPI/CFD/octave/totalPressureDrop.m
View File

@ -74,9 +74,20 @@ fprintf('so the result does not depend on pressure\n')
%==================================
%rhoP = 2000 % particle density in kg/m3
g = 9.81 % gravity m/s2
Umf = dp^2*(rhoP-rhoG)*g/(150*muG)*(epsilon^3*phip^2)/(1-epsilon)
ReMF = Umf*dp*rhoG/muG % must be <20 !!!
%Umf = sqrt(phip*dp^2/1.75*(rhoP-rhoG)/rhoG*g*epsilon^3) % Re>1000
Umf = dp^2*(rhoP-rhoG)*g/(150*muG)*(epsilon^3*phip^2)/(1-epsilon);
ReMF = Umf*dp*rhoG/muG;
if(ReMF<20)
fprintf('applying eqn1 for Umf.\n')
elseif(ReMF>20 && ReMF<1000)
fprintf('applying eqn1 for Umf.\n')
elseif (ReMF>=1000)
fprintf('applying eqn2 for Umf.\n')
Umf = sqrt(dp*(rhoP-rhoG)*g/(1.75*rhoG)*epsilon^3*phip);
ReMF = Umf*dp*rhoG/muG;
end
Umf
ReMF
dpUmf= L * (
150*((1-epsilon)^2/epsilon^3)*((muG.*Umf)/(phip*dp)^2)