flow rate calculation in the microfluidic channels using data obtained from PIV using the Purday approximation.

 

%Program uses the maximum velocity obtained from PIV to calculate the flow
%rate
%Define time step (us)
time step=90;

%Channel width (um)
w=1500;
%Channel height (um)
h=212.5;
%Read in the image vectors
v=loadvec(’*.vc7’);
%For each time step/file:
%Access the x-component of velocity data, vx and find the maximum velocity
%Average the velocity data along each column to represent a global average
%for the divided pixel down the image window. Store in an array
maximum velocity=zeros(1,numel(v));
for i = 1 : numel(v)
%Average each column and store in a row each
average velocity along y=mean(v(i).vx,1);
%account for negatives (oscillatory flow)
avg=mean(average velocity along y);
if avg¿0
maximum velocity(i)=max(average velocity along y);
else
maximum velocity(i)=min(average velocity along y);
end
end
%Comes out currently in units of pixels
%Multiply by pixel to micron scaling factor for 10x magnification lens.
%Number of microns/pixels
pixel to micron=550/839.75;
%So now in units of microns
maximum velocity=maximum velocity*pixel to micron;
%Convert to m
maximum velocity=maximum velocity*10ˆ-6;

%Divide by time step to get velocity in m/s
maximum velocity=maximum velocity/(time step*10ˆ-6);
%Calculate channel area (in mˆ2)
A=w*h*(10ˆ-6)ˆ2;
%Calculate flow rate through Purday Approximation - Calculate n and m as
%well
alpha=h/w;
%Calculate n
if alpha ¿= (1/3);
n=2+0.3*(alpha-(1/3));
else
n=2;
end
%Calculate m
m=1.7+0.5*alphaˆ-1.4;
%Flow rate in mˆ3/s
Q=(m/(m+1))*(n/(n+1))*A*maximum velocity;
%Convert to ml/min
Q=Q*60*10ˆ3*10ˆ3;