// Parameters

// solid fraction (below 0.1)
sf       0.02;    

// interstitial velocity
Uint        1;

// domain size 
Ldomain     2; 
// number of cells per sphere diameter (blockmesh)
ncell_D     6;

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

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

// residual control
re_control      #calc "1e-6";

// Calculations
Vdomain  #eval "$Ldomain*$Ldomain*$Ldomain"; 
Vsphere  #eval "$sf*$Vdomain"; 
r3       #calc "0.75*$Vsphere/M_PI"; 

// sphere radius
r        #eval "cbrt($r3)"; 
D        #eval "2*$r"; 

Cmin    #eval "-0.5*$Ldomain";
Cmax    #eval "0.5*$Ldomain"; 

ncell   #eval "ceil($Ldomain/$D)*$ncell_D";  

// reference point
xref #eval "$Cmax"; 
yref #eval "$Cmax";
zref #eval "$Cmax";
