zhyang-dev/CFDEMcoupling-PFM
1
0
Fork 0
mirror of https://github.com/ParticulateFlow/CFDEMcoupling-PFM.git synced 2025-12-08 06:37:44 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

release on 2012-11-23_14-05-54

Browse Source
This commit is contained in:
cfdem
2012-11-23 14:05:54 +01:00
parent 552fcb8700
commit 80403ff260
53 changed files with 3579 additions and 94 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
tutorials/cfdemSolverIB/twoSpheresGlowinskiMPI/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/cfdemSolverIB/twoSpheresGlowinskiMPI/DEM/in.liggghts_init
View File

@ -13,7 +13,7 @@ boundary f f f
newton off
units si
processors 2 2 2
processors 2 2 1
region reg block 0 1 0 1 0 4 units box
create_box 1 reg

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

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

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

@ -67,14 +67,12 @@ elseif (ReMF>=1000)
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)
+1.75*((1-epsilon)/epsilon^3)*((rhoG.*Umf.^2)/(phip*dp))
)/10000/rhoG;
%dpUmf2=(L*(1-epsilon)*(rhoP-rhoG)*g+pHydr)/10000
%====================================%
% plot data
%====================================%

14
tutorials/cfdemSolverPiso/ErgunTestMPI_restart/CFD/octave/totalPressureDrop.m
View File

@ -71,9 +71,17 @@ 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
dpUmf= L * (
150*((1-epsilon)^2/epsilon^3)*((muG.*Umf)/(phip*dp)^2)