How to plot this velocity profile?

Respected sir,

sir I have derived to velocity profile that i have attached below .As I am new so please help me Sir please provide some hint to plot U2 velocity profile. here tau ,tua0, delta all are constant value except y. I am waiting for reply

[madhan ravi]

upload the code you have tried

Answer 1

[madhan ravi]

clear all.  %edited

close all

clc

mu=0.5;

t_o=100;

t_e=200;

n=15;

lambda=0.25;

delta = 0.25; 

dh=linspace(0.01,0.05,n);

y=5;

U = 1./(mu).*(t_o.*(y-delta)+4.*lambda.*sqrt(t_o.*t_e).*(exp(-y./(2.*lambda)))+t_e.*lambda.*(1-exp(-y./lambda)))

plot(dh,U,'-m')  % you have to mention which variables you want to plot

[madhan ravi]

Type U2 formula and paste it here I’ll help you further

[Debasis Roy]

Respected sir

As i have written math type that's why i am unable to paste here . For this reason i have attched picture file of this formula

[madhan ravi]

I suggest you to try typing something in matlab , i will correct it if it’s wrong

[Debasis Roy]

Respected sir

What ever i intially done or able do is

clear all
close all
clc
mu=0.5;
t_0=100;
t_e=200;
n=15;
lamda=0.25;
dh=linspace(0.01,0.05,n);
y=5;

Answer 2

[Debasis Roy]

Respected sir I am extremely sorry for above expression clearity as i am not specific about it *i want to See how velocity (U)is changing with the ratio of (lamda/delta)

following the value of constants

tau0= 100

tauE=250

lamda varies from case 1)0.1 to 0.5 case 2)lamda=0.3

delta varies from 0/025 to 0.87 what ever i am tried

 clear all
close all
clc
mu=0.5; t_0=100; t_e=200; n=15; lamda=0.25; dh=linspace(0.01,0.05,n); y=5;

i am unable to proceed please help

[Debasis Roy]

Respected Sir,

I want to see what is the change coming in velocity profile if I change the ratio (lamda /delta). if this ratio 0.11 . 0.01, 0.02 then how it vary with y

Answer 3

[Debasis Roy]

Sir if I write is it ok? but I am not getting the result. As I want different plot of u vs y corresponding to different values of lambda and delta

mu=0.6;
t_o=106;
t_e=156;
 deltavals = [0.6 1.1 1.6];
lambdavals = [0.5 0.6 0.7 0.8];
yvals = 0.3:0.1:1;
[delta, lambda, y] = ndgrid(deltavals, lambdavals, yvals);
U = 1./mu.*(t_o.*(y-delta)+4.*lambda.*sqrt(t_o.*t_e).*exp(-lambda./(2.*delta))+t_e.*lambda.*(1-exp(-y./lambda)))
plot(y,U,'-r')

Answer 4

[Aliya Altaf]

How I make this .

Answer 5

[Konde]

R=0.5;
L=10;
v=0.0028;
Pa=15;
Pb=30;
r=linspace(-R,R,20);
V=((Pb-Pa)*(R^2-r^.2)/(4vL));
figure;
plot(V,r,’g’);
Xlabel(‘velocity’);
Ylabel(‘radius’);
A=table(V’,r’,’variable names’,{‘velocity’,’radius’},824-velocity profile);
Disp(A);