|
|
Re: analytical solution laminar flow in rectangular microchannel
Posted:
Dec 4, 2012 4:47 AM
|
|
> I expect the velocity field perpendicular to the
> channel entry (y positions) to the right and it still
> can be observed instantaneously along the channels(x
> position).
>
> It should be something like the velocity field in
> this figure link.
>
> http://ars.els-cdn.com/content/image/1-s2.0-S026322410
> 800078X-gr6.jpg
>
> Best Regards,
> Anne
x_array: Array of coordinates of x-position in the channel
y_array: Array of coordinates of height positions at the channel entry
velocity_array: Array of velocities at the height positions at the channel entry
(Arrays y_array and velocity_array have the same size).
Then you could proceed as follows:
n=length(x_array)
m=length(y_array)
x_quiver_array=zeros(m*n);
y_quiver_array=zeros(m*n);
x_velocity_quiver_array=zeros(m*n);
y_velocity_quiver_array=zeros(m*n);
For i=1:n
For j=1:m
x_quiver_array((i-1)*m+j)=x_array(i);
y_quiver_array((i-1)*m+j)=y_array(j);
x_velocity_quiver_array((i-1)*m+j)=velocity_array(j);
y_velocity_quiver_array((i-1)*m+j)=0.0;
end
end
quiver(x_quiver_array,y_quiver_array,x_velocity_quiver_array,y_velocity_quiver_array);
Best wishes
Torsten.
|
|
