It is usual in the field of mechanics of fluids to plot interest variables as function of the thickness and along the length, such as velocity profiles, as:
which shows the velocity profile across the boundary layer in an open flow, or:
which shows the velocity profile of internal flow. This type of plot is extremely helpful to visualize the distribution of the variable across the two dimensions. This is not restricted to velocity components, as temperature may also be shown similarly:
It is clearly in the above that the tempreature distribution is constant at the entrance and parabolic the fully developed region. This type of plot, for instance, is used to show the temperature distribution of lubricant flow in radial bearings, as in:
We can see in this plot the exact distribution of temperature at selected locations of the x axis.
How can one plot such distribution in MATLAB?
Suppose the simplest case:
subplot(121); quiver(X,Y,U,0*U);
subplot(122); surf(X,Y,U); view(2);
The quiver and surf functions let us visualize the variable distribution on space, but this distribution is not entirely clear either from the quiver or surf plots. We can enhance the quality of the figure in a different way:
figure; contourf(X,Y,U,100,'LineColor','none'); hold on; quiver(X,Y,U,0*U,'k');
And, although the colorscale and the arrows lenght give us an idea of how the distribution of U as function of y is at each x location, I would like to actually plot this distribution. I can plot it at several different tilesets:
tiledlayout(1,6,"TileSpacing","none");
nexttile; plot(U(:,ii),Y(:,ii)); xticks([]); yticks([]);
But this solution is not good. All plots look the same (we can certainly change the xlim at each plot to scale all plots properly), but one fundamental information is lot: the location of each plot along the x-axis. I even tried to add another axis to the background of this plot, to include the x distance, but this solution was not good, because it is extremely dificutl to align both axes programatically, as the variable may have a negative value, so it is not easy to align the 0 with the desired x position on the axis below. The second difficult with this approach is that the y variable may change in scale, so I don't know how this can be accomplished.
The quiver and surf solutions really are simple. But if the expressions are obtained analitically, how can we plot this distributions as desired?
I though on creating the quiver plot with a fine mesh on the y distribution, without arrow heads, and connect the arrow heads with a line to give me the graph of the variable as function of y. Is this possible?
I am using MATLAB R2021b, but have access to newer versions.
