halfWidth 0.02[m] "Channel half-width"
height 10[m] "Channel height"
nelemWidth 20 "Number of mesh elements along half-width"
nelemHeight 100 "Number of mesh elements along height"
va 1.65[m/s] "Average fluid velocity"
rhof 1.3[kg/m^3] "Air density"
nuf 15.7e-6[m^2/s] "Air kinematic viscosity"
muf nuf*rhof "Air dynamic viscosity"
ustar_exp 0.11775[m/s] "Expected wall velocity"
rhop 769*rhof "Particle density"
St 125 "Particle Stokes number"
taup St*nuf/ustar_exp^2 "Particle relaxation time"
t1 675*nuf/ustar_exp^2 "Time for t_plus = 675"
t2 1125*nuf/ustar_exp^2 "Time for t_plus = 1125"
dp sqrt(18*taup*nuf*rhof/rhop) "Particle diameter"
Np 10000 "Number of particles"