MATLAB: Flow in the pipe

flowpipeviscous fluid

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

Best Answer

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?
&nbsp
&nbsp
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.
Related Question