Facing problem about plot the velocity profile

Debasis Roy (view profile)

on 11 Nov 2018
Latest activity Commented on by Souvick roy

Souvick roy (view profile)

on 4 Dec 2018
Respected sir i have attched one picture here , likethis picture i want to plot velocity(u) vs y plot
my formula for velocity profile is below whatever I have done in Matlab code I am writing below the formulae
U = 1./mu.*(t_o.*(y-delta)+4.*alpha.*sqrt(t_o.*t_e).*exp(- alpha./(2.*delta))+t_e.*alpha.*(1-exp(-y./alpha)))
% t_o=106;
t_e=156;
deltavals = 0.6 ;
alphavals = [-0.3,- 0.2 ,-0.1,0, 0.1,0.2,0.3 ];
yvals = 0.3:0.1:1;
[ alpha, y] = ndgrid( alphavals, yvals);
U = 1./mu.*(t_o.*(y-delta)+4.*alpha.*sqrt(t_o.*t_e).*exp(- alpha./(2.*delta))+t_e.*alpha.*(1-exp(-y./alpha)))
plot(y,U,'-r') end

Show 1 older comment
Debasis Roy

Debasis Roy (view profile)

on 11 Nov 2018
Respected sir value of variable t_o=100 t_e=150 mu=0.3 delta=0.6 alpha values are 0.1,0.2,0.3,0.3,0.4,-0.2,-0.5
Image Analyst

Image Analyst (view profile)

on 1 Dec 2018
OK, now, what is y?
t_o = 100;
t_e = 150
mu = 0.3
delta = 0.6
alpha = [0.1, 0.2, 0.3, 0.3, 0.4, -0.2, -0.5]
U = 1./mu.*(t_o.*(y-delta)+4.*alpha.*sqrt(t_o.*t_e).*exp(- alpha./(2.*delta))+t_e.*alpha.*(1-exp(-y./alpha)))
% t_o = 106;
t_e = 156;
deltavals = 0.6 ;
alphavals = [-0.3,- 0.2 ,-0.1,0, 0.1,0.2,0.3 ];
yvals = 0.3:0.1:1;
[ alpha, y] = ndgrid( alphavals, yvals);
U = 1./mu.*(t_o.*(y-delta)+4.*alpha.*sqrt(t_o.*t_e).*exp(- alpha./(2.*delta))+t_e.*alpha.*(1-exp(-y./alpha)))
plot(y,U,'-r')
If this is still a problem for you, give us everything we need to help you.
Souvick roy

Souvick roy (view profile)

on 4 Dec 2018
sir due to problem in value i have non dimensionalise my final equation. i have attached the picture. in wquation one is valid from y =0 to y=0.3 and equation two is valid for y=0.3 to 1
I want to divide y in 0.01 interval
balue of delta 0.3
value of k=100
value of beta 0.01