flow in the pipe
69 views (last 30 days)
Show older comments
Hi i want simulation of flow in the pipe.
but my script shows something wrong. T_T
I want to show second images like this Image.
how can i do more?
dbstop if error
Ps=10;
Pf=50;
mu=8.94e-4;
R=0.10;
l=1.00;
if Ps==Pf
disp(' Same pressure dont show Velocity in pipe.');
end
N=100;
r=linspace(-R,R,N);
S=pi*R^2;
Vmax=(abs(Ps-Pf))*(R^2)/(4*mu*l);
Vmean=Vmax/2;
V=Vmax*(1-(r.^2)/R^2);
Q=Vmean*S;
subplot(2, 1, 1)
plot(V,r,'b.');
grid on, xlabel(' Velocity (m/s)'), ylabel(' Radius (m)');
T=100*ones(100,R*2)
for iter = 1:100
for t=1:100
va=V*t
end
subplot(2, 1, 2)
imagesc(va')
drawnow
end
2 Comments
Amit
on 9 Dec 2014
Which line you get error? Also can you use the '{} code' to post your complete code?
Image Analyst
on 12 Dec 2014
a's "Flag" of this comment moved here, and the flag removed:
well... my english is poor. I want to show flow's displacement using V in second plot. like animation second plot means pipe. starting point velocity in pipe is 0. but t=0,1...,100 , flow is moving. velocity is side velocity slower than the center. velocity slows the farther away from the center. Do you understand?
Accepted Answer
per isakson
on 9 Dec 2014
Edited: per isakson
on 10 Dec 2014
The inner for-loop overwrites va one hundred times. The outer for-loop draws the "same" image one hundred times in the same place.
In plain English, what do you want to show in the second plot?
 
 
Replace
T=100*ones(100,R*2);
...
% end of script
by
% T=100*ones(100,R*2);
subplot(2, 1, 2)
plot( 0, r, 'r.' );
axh = gca;
set( axh, 'XLim', [0,1e3] )
hold all
for t = 1:6
displacement = V*t;
plot( axh, displacement, r, 'r.' );
pause( 1 )
end
and run the script. That will take you a bit closer to the goal.
5 Comments
More Answers (0)
See Also
Categories
Find more on Fluid Dynamics in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!