Re 60;
ncell_D    6;   // number of cells per equivalent sphere diameter
sf         0.2; // solid fraction 

nu        1;  // fluid kinematic viscosity
rhoFluid  1;  // fluid density

D         1;     // spherocylinder diameter
AR        4;     // spherocylinder aspect ratio
nscyl     200;   // number of spherocylinders

Cmin      0.; // minimum domain coordinate 

// global orientation vector
// horizontal: ex 1, ey 0, ez 0
// vertical:   ex 0, ey 0, ez 0 
ex       1;
ey       0; 
ez       0; 

// tangentiality components
Delta_mo 0.25; // Minimum boundary overlap
Delta_md 0.1;  // Minimum distance from the boundary

// snappy settings
nCellsBetweenLevels  4; 
surfLevel            3;
gapLevel             4; 

// time settings 
simulationEnd  10000;
simulationWrite 1000;

// searchableSurfaces parameters 
gapSize      0.3;  // minimum relative distance between particles

// Calculations
r       #eval "$D/2"; // radius
beta    #eval "0.5*(3*$AR-1)";
Vp      #calc "M_PI*$beta*pow($D,3)/6"; // spherocyl volume
Vsolid  #eval "$nscyl*$Vp";
Vdomain #eval "$Vsolid/$sf";
Ldomain #eval "cbrt($Vdomain)"; 
Cmax    #eval "$Cmin+$Ldomain"; 

// block mesh cell size
ncell #eval "ceil(($Cmax-$Cmin)/$D)*$ncell_D"; 

// volume-equivalent sphere
Deq  #eval "cbrt($beta)*$D";

// superficial velocity
Ubar   #eval "$Re*$nu/$Deq"; 

vf     #eval "1-$sf";    // void fraction  
Uint   #eval "$Ubar/$vf";  // interstitial velocity
