CJESTMIR^.COM

#scilab -nwni -f filename.sci
-->exec 'code.sci'

function [eta0]=etafct(dp,dc,alpha,U,T,mu,rhop)
ap=dp/2;     //radius of the particle (m)
NR=dp/dc;	//Aspect Ratio
Dinf=(k*T)/(6*%pi*mu*ap);	//Diffusion coefficient, in an infinite medium
NPe=(U*dc)/(Dinf);	//Peclet NUmber, characterizing ratio of convective 
			//transport to diffusive transport
NvdW=A/(k*T);		//Van der Waals number, characterizing ratio of vdW 
			//interaction energy to the particle's thermal energy
NA=A/(12*%pi*mu*(ap^2)*U);	//Attraction number, represents combined 
	//influence of vdW attraction forces and fluid velocity on particle 
	//deposition rate due to interception
NG=(2/9)*(((ap^2)*(rhop-rhof)*g)/(mu*U));	//Gravity number, ratio of 
	//Stokes particle settling velocity to approach velocity of the fluid
Y=(1-f)^(1/3);	//Y-part of porosity-dependent parameter of Happel's model
AS=(2*(1-Y^5))/(2-3*Y+3*Y^5-2*Y^6);	//porosity-dependent parameter of 
	//Happel's model

eta0=(2.4*AS^(1/3)*NR^(-0.081)*NPe^(-0.715)*NvdW^(0.052))+(0.55*AS*NR^(1.675)
	*NA^(0.125))+(0.22*NR^(-0.24)*NG^(1.11)*NvdW^(0.053));	
	//long equation, glue the second line to the one above...
endfunction

eta0=etafct(dp,dc,alpha,U,T,mu,rhop);