|
|
Re: analytical solution laminar flow in rectangular microchannel
Posted:
Nov 16, 2012 7:12 AM
|
|
> Hello!
> I got another code for the previous question. It
> works like this.
> v=0.0003;
> w=100e-6;
> h=100e-6;
> y=0;
> z=0;
> i=0;
> while(i<100)
> n=1;
> temp=0;
> ans=0;
> while (n<10000)
> a=sin(n*pi*z/h);
> b=cosh(n*pi*y/h);
> c=cosh(n*pi*w/(2*h));
> d=(1/(n^3));
> e=tanh(n*pi*w/(2*h));
> f=(192*h)/((n^5)*(pi^5)*w);
> g=(48*v)/(pi^3);
> value=(g*d*(1-(b/c))*a)/(1-(f*e));
> ans=ans+value;
> n=n+2;
> end
> result((i+1),1)=z;
> result((i+1),2)=ans;
> i=i+1;
> z=h/999*i;
> end
> xlswrite('Result',result);
>
> I want to try to use the previous code to generate
> same result as above. But I can't get similar result.
Can you give me the link from where you got
the formula for the axial velocity as you use
it in your code above ?
Best wishes
Torsten.
|
|
