// Parameters
Re        20;   // Reynolds number
ncell_D    6;   // number of cells per cylinder diameter

Ubar      0.1;  // superficial velocity
rhoFluid  1000; // fluid density

D         0.2; // cylinder diameter
AR        2;   // cylinder aspect ratio
ncyl      16;  // number of cylinders

Cmin      0.2; // minimum domain coordinate 
Cmax      1.2; // maximum domain coordinate

// snappy settings
nCellsBetweenLevels  4; 
surfLevel            2; // surfaceRefinementLevel in snappyHexMesh; 
                        // increase to 4 for visually nicer mesh
gapLevel             3; 

// Calculations
r     #eval "$D/2"; // radius

sf    #calc "$ncyl*(M_PI*pow($r, 2)*$AR*$D)"; // solid fraction 
vf    #eval "1-$sf"; // void fraction  
Uint  #eval "$Ubar/$vf";  // interstitial velocity

nu    #eval "$Ubar*$D/$Re";  // fluid kinetic viscosity

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