|
|
Re: analytical solution laminar flow in rectangular microchannel
Posted:
Oct 29, 2012 3:58 AM
|
|
> I'm new to matlab. I need to assign and analysis
> velocity profile for the Poiseuille flow in a
> rectangular microchannel. I need help on the
> equations below.
>
> a= sigma_(n,odd to inf)(1/n^3) [1-(cosh(n*pi*
> y/h))/(cosh(n*pi* w/h))] sin(n*pi* z/h)
> b= 1-sigma_(n,odd to inf)(192h/(n^5* pi^5* w)
> *tanh(n*pi*w/h))
>
> u_x (y,z)=48*Q/(pi^3*h*w)*a/b
>
>
> Message was edited by: Anne
Something like
n=1;
term1=1.0;
term2=1.0;
eps=1e-8;
n=1;
a=0;
b=0;
While (term1 > eps) | (term2 > eps)
term1 = 1/n^3 [1-(cosh(n*pi* y/h))/(cosh(n*pi* w/h))]* sin(n*pi* z/h);
term2 = 192*h/(n^5* pi^5* w)*tanh(n*pi*w/h)
a=a+term1;
b=b+term2;
n=n+2;
end
b=1-b;
ux=48*Q/(pi^3*h*w)*a/b;
Best wishes
Torsten.
|
|
