release on 2013-03-07_14-30-36

tutorials/cfdemPostproc/fillCylinder/CFD/constant/couplingProperties

@@ -38,7 +38,7 @@ meshMotionModel noMeshMotion;
 
 regionModel allRegion;
 
-IOModel off;
+IOModel basicIO;
 
 dataExchangeModel oneWayVTK;
 

tutorials/cfdemPostproc/fillCylinder/DEM/post/dummy

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

tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/octave/totalPressureDrop.m

@@ -5,6 +5,7 @@ clc;
 %====================================%
 % simulation data 1
 %====================================%
+rhoG = 10			% density in kg/m3
 %path = '../probes/0/p';
 path = '../probes/0/p';
 columns=22;
@@ -12,7 +13,7 @@ headerlines=4;
 data = loaddata(path,columns,headerlines);
 data=transpose(data);
 [x,y]=size(data)
-dp_sim = (data(:,2)-data(:,y))/10000;
+dp_sim = (data(:,2)-data(:,y))*rhoG; %conversion to Pa!
 t_sim = data(:,1);
 %fprintf('final pressureDrop of sim = %f Pa\n',dp_sim(length(dp_sim)) )
 
@@ -36,14 +37,13 @@ deltaU=(Uend-Ustart)/((Tend-Tstart)/timeStepSize);
 U = Ustart+deltaU:deltaU:Uend;  % velocity over time
 Ua = U / epsilon;		% physical velocity
 L = 0.0156			% length of bed
-rhoG = 10			% density in kg/m3
 nuG = 1.5*10^-4			% kinemat Visk in m2/s
 muG = nuG*rhoG			% dynam visc in Pa s
 
 dpErgun= L * (
                 150*((1-epsilon)^2/epsilon^3)*((muG.*U)/(phip*dp)^2) 
               +1.75*((1-epsilon)/epsilon^3)*((rhoG.*U.^2)/(phip*dp))
-        )/10000/rhoG;
+        );
 
 fprintf('NOTE: this pressure is divided by density (according to CFD solver)\n')
 fprintf('so the result does not depend on density\n')
@@ -71,20 +71,18 @@ end
 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
+        );
+%dpUmf2=(L*(1-epsilon)*(rhoP-rhoG)*g+pHydr)
 %====================================%
 % plot data
 %====================================%
-length(U)
-length(dp_sim)
 figure(2)
 plot(U,dp_sim)
 title("Ergun pressure drop vs. simulation")
 a=strcat("analytical (Ergun), Umf=",num2str(Umf),", dpUmf=",num2str(dpUmf));
 legend(a,"simulation")
 xlabel("velocity in [m/s]")
-ylabel("pressure drop [bar]")
+ylabel("pressure drop [Pa]")
 axis([0,Uend,0,dpErgun(length(dpErgun))])
 
 figure(1)
@@ -93,7 +91,7 @@ title("Ergun pressure drop vs. simulation")
 a=strcat("analytical (Ergun), Umf=",num2str(Umf),", dpUmf=",num2str(dpUmf));
 legend(a,"simulation","analyt. deltaP at Umf")
 xlabel("velocity in [m/s]")
-ylabel("pressure drop [bar]")
+ylabel("pressure drop [Pa]")
 axis([0,Uend,0,dpErgun(length(dpErgun))])
 
 %print('cfdemSolverPiso_settlingTest.eps','-deps2')

tutorials/cfdemSolverPiso/ErgunTestMPI_restart/CFD/octave/totalPressureDrop.m

@@ -5,6 +5,7 @@ clc;
 %====================================%
 % simulation data 1
 %====================================%
+rhoG = 10			% density in kg/m3
 %path = '../probes/0/p';
 path = '../probes/0/p';
 columns=22;
@@ -12,7 +13,7 @@ headerlines=4;
 data = loaddata(path,columns,headerlines);
 data=transpose(data);
 [x,y]=size(data)
-dp_sim = (data(:,2)-data(:,y))/10000;
+dp_sim = (data(:,2)-data(:,y))*rhoG; % conversion to Pa
 t_sim = data(:,1);
 %fprintf('final pressureDrop of sim = %f Pa\n',dp_sim(length(dp_sim)) )
 
@@ -26,7 +27,7 @@ headerlines=4;
 data = loaddata(path,columns,headerlines);
 data=transpose(data);
 [x,y]=size(data)
-dp_sim_2 = (data(:,2)-data(:,y))/10000;
+dp_sim_2 = (data(:,2)-data(:,y))*rhoG; % conversion to Pa
 t_sim_2 = data(:,1);
 %fprintf('final pressureDrop of sim = %f Pa\n',dp_sim(length(dp_sim)) )
 
@@ -52,14 +53,13 @@ deltaU=(Uend-Ustart)/((Tend-Tstart)/timeStepSize);
 U = Ustart+deltaU:deltaU:Uend;  % velocity over time
 Ua = U / epsilon;		% physical velocity
 L = 0.0156			% length of bed
-rhoG = 10			% density in kg/m3
 nuG = 1.5*10^-4			% kinemat Visk in m2/s
 muG = nuG*rhoG			% dynam visc in Pa s
 
 dpErgun= L * (
                 150*((1-epsilon)^2/epsilon^3)*((muG.*U)/(phip*dp)^2) 
               +1.75*((1-epsilon)/epsilon^3)*((rhoG.*U.^2)/(phip*dp))
-        )/10000/rhoG;
+        );
 
 fprintf('NOTE: this pressure is divided by density (according to CFD solver)\n')
 fprintf('so the result does not depend on density\n')
@@ -86,7 +86,7 @@ end
 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;
+        );
 
 %====================================%
 % plot data
@@ -98,7 +98,7 @@ dpUmf= L * (
 %a=strcat("analytical (Ergun), Umf=",num2str(Umf),", dpUmf=",num2str(dpUmf));
 %legend(a,"simulation")
 %xlabel("velocity in [m/s]")
-%ylabel("pressure drop [bar]")
+%ylabel("pressure drop [Pa]")
 %axis([0,Uend,0,dpErgun(length(dpErgun))])
 
 figure(1)
@@ -107,7 +107,7 @@ title("Ergun pressure drop vs. simulation")
 a=strcat("analytical (Ergun), Umf=",num2str(Umf),", dpUmf=",num2str(dpUmf));
 legend(a,"simulation")
 xlabel("velocity in [m/s]")
-ylabel("pressure drop [bar]")
+ylabel("pressure drop [Pa]")
 axis([0,Uend,0,dpErgun(length(dpErgun))])
 
 %print('cfdemSolverPiso_settlingTest.eps','-deps2')

tutorials/cfdemSolverPiso/ErgunTestMPI_restart/CFD/system/controlDict

@@ -19,11 +19,11 @@ application     pisoFoam;
 
 startFrom       startTime;
 
-startTime       0;
+startTime       0.05;
 
 stopAt          endTime;
 
-endTime         0.05;//0.01;
+endTime         0.1;//0.01;
 
 deltaT          0.001;
 

tutorials/cfdemSolverPiso/ErgunTestMPI_restart/DEM/post/dummy

@@ -35,3 +35,5 @@ dummyfile
 dummyfile
 dummyfile
 dummyfile
+dummyfile
+dummyfile