flow in the pipe
显示 更早的评论
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 个评论
Amit
2014-12-9
Which line you get error? Also can you use the '{} code' to post your complete code?
Image Analyst
2014-12-12
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?
采纳的回答
per isakson
2014-12-9
编辑：per isakson
2014-12-10
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 个评论
per isakson
2014-12-10
Do you assume an ideal model of fully developed laminar flow?
per isakson
2014-12-10
编辑：per isakson
2014-12-10
Your English is adequate. However, I didn't understand that you wanted to plot the displacement.
See my answer
a
2014-12-10
Arik Chattopadhyay
2021-3-31
can I get the entire code?
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Fluid Dynamics 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
