% Idelchik estimation

w = 0.1; % 2l/min
mu = 1.825E-5;%air 20degC
rho = 1.204;%air 20degC
L = 0.1; %[m]
d_gr = 1.5e-3;
phi = 0.76;
eps = 0.35;
B = 1.8; %1.8 bodies smooth surface, 4.0 bodies rough surface
g = 9.81; %[m/s^2]
T_in = 293;
T_ex = 293;


d_el = phi*d_gr;
Re = w*d_el*rho/mu;
lambda = 360*(1-eps)^2/(eps^3*Re)+B*(1-eps)/eps^3;
T_m = (T_in+T_ex)/2;
delta_psi_t = 2*(T_ex-T_in)/T_m;


kappa = inv(1/(w*mu)*(lambda*rho*w^2/(2*d_el)+delta_psi_t*rho*w^2/(2*L)));

disp('Re')
disp(Re)
disp('Permeability:')
disp(kappa)

Re
    7.5209

Permeability:
  6.9261e-010

Published with MATLAB® 7.8