# How to plot axial velocity profiles

30 views (last 30 days)
Bernard Hippolyte on 31 Jul 2019
Commented: infinity on 2 Aug 2019
Hello,
I am studying flow conditions over a flat plate. In my analysis I'd like to make a graph representing the velocity profile at different locations along the x-axis.
I have 3 velocity vectors for 3 different locations along the x-axis, whose values change according to the height with respect to the flow surface. I have created a first plot with that code :
z = [0,9.99999953400000e-05,0.000219999991200000,0.000363999986200000,0.000536799977600000,0.000744159964100000,0.000992991944000000,0.00129159031500000,0.00164990835500000,0.00207988999700000,0.00259586796200000,0.00321504151600000,0.00395804977800000,0.00484965968200000,0.00591959154100000,0.00720350974200000,0.00874421157500000,0.0105930537600000,0.0120000000000000];
v1 = [0,0.00325397021600000,0.00701125750100000,0.0113216226500000,0.0162260455500000,0.0217406444100000,0.0278199545200000,0.0343086455900000,0.0409495839500000,0.0474558756800000,0.0535523554300000,0.0589996739700000,0.0636245292600000,0.0673341282000000,0.0701542241100000,0.0722371526700000,0.0737607573100000,0.0748035608600000,0.0750795450500000];
v2 = [0,0.00366377801700000,0.00791285352500000,0.0128013963500000,0.0183614737800000,0.0245786185200000,0.0313532627000000,0.0384678052200000,0.0456212766800000,0.0525096049400000,0.0588643367300000,0.0644851997700000,0.0692671336100000,0.0731802814600000,0.0762524020600000,0.0785669430900000,0.0802213325000000,0.0812515896100000,0.0814175550000000] ;
v3 =[0,0.00520340445400000,0.0111829468300000,0.0179754812400000,0.0255750928700000,0.0338966479500000,0.0427261762400000,0.0517216408400000,0.0604887040700000,0.0686496402800000,0.0758974471700000,0.0820465381000000,0.0870401606400000,0.0909297100700000,0.0938583875400000,0.0960021018900000,0.0974763949700000,0.0983292911300000,0.0984241062400000] ;
plot(v1,z,v2,z,v3,z)
legend('x = 0 mm','x = 50 mm','x = 100mm')
xlabel('velocity [m/s]')
ylabel('height [m]')
But contrary to what my matlab code does, I would like to create a new axis (x-axis) representing the location of each velocity profile as shown in the following figure :
Despite my research on the web I didn't found any solution... I will be grateful to you for your help.
Thanks so much,
Hippolyte

infinity on 1 Aug 2019
Hello,
You can refer this code
z = [0,9.99999953400000e-05,0.000219999991200000,0.000363999986200000,0.000536799977600000,0.000744159964100000,0.000992991944000000,0.00129159031500000,0.00164990835500000,0.00207988999700000,0.00259586796200000,0.00321504151600000,0.00395804977800000,0.00484965968200000,0.00591959154100000,0.00720350974200000,0.00874421157500000,0.0105930537600000,0.0120000000000000];
v1 = [0,0.00325397021600000,0.00701125750100000,0.0113216226500000,0.0162260455500000,0.0217406444100000,0.0278199545200000,0.0343086455900000,0.0409495839500000,0.0474558756800000,0.0535523554300000,0.0589996739700000,0.0636245292600000,0.0673341282000000,0.0701542241100000,0.0722371526700000,0.0737607573100000,0.0748035608600000,0.0750795450500000];
v2 = [0,0.00366377801700000,0.00791285352500000,0.0128013963500000,0.0183614737800000,0.0245786185200000,0.0313532627000000,0.0384678052200000,0.0456212766800000,0.0525096049400000,0.0588643367300000,0.0644851997700000,0.0692671336100000,0.0731802814600000,0.0762524020600000,0.0785669430900000,0.0802213325000000,0.0812515896100000,0.0814175550000000] ;
v3 =[0,0.00520340445400000,0.0111829468300000,0.0179754812400000,0.0255750928700000,0.0338966479500000,0.0427261762400000,0.0517216408400000,0.0604887040700000,0.0686496402800000,0.0758974471700000,0.0820465381000000,0.0870401606400000,0.0909297100700000,0.0938583875400000,0.0960021018900000,0.0974763949700000,0.0983292911300000,0.0984241062400000] ;
plot(v1,z,v2+0.2,z,v3+0.4,z)
legend('x = 0 mm','x = 50 mm','x = 100mm')
xticks([0 0.1 0.2 0.3 0.4 0.5])
xticklabels({'0','0.1','0','0.1','0','0.1'})
xlabel('velocity [m/s]')
ylabel('height [m]')
##### 2 CommentsShowHide 1 older comment
infinity on 2 Aug 2019
Hello,
I am still not get your point. But, maybe someone else can help you.
Good luck