% Initialize some parameters
rhog  = 0.3;   % Density of gas
a     = 1e-6;  % Particle radius
mu    = 3e-5;  % Viscosity of gas
rhop  = 12700; % Density of particle
Omega = 5000;  % Flow rotation rate
Wg    = 1100;  % Axial flow velocity
Ro    = 0.15;  % Radius of chamber
Rs    = 0.14;  % Radius of skimmer
L     = 1.0;   % Length of chamber

mp   = rhop*4*pi/3*a^3; % mass particle

% Store the parameters in a vector for later use
p = [rhog, a, mu, Omega, Wg, mp];