Hi,
To visualize a doublet flow with a strength of 15 units at the origin in MATLAB, begin by setting the doublet strength to 15. Define the grid using “linspace” for x and y from -2 to 2 and create a 2D grid with “meshgrid”.
Calculate the velocity components U and V using the doublet equations and determine the stream function ‘psi’ for contour visualization. Use “quiver” for the vector plot and “contour” for the streamlines.
Ensure the plot has equal axis scaling and appropriate labels for clarity. Adjust grid resolution as needed for detail.
Please refer to the following example code of vector plot with contours showing doublet flow:
% Doublet strength
mu = 15;
% Grid setup
x = linspace(-2, 2, 100);
y = linspace(-2, 2, 100);
[X, Y] = meshgrid(x, y);
% Doublet flow components
U = -mu * (X.^2 - Y.^2) ./ (X.^2 + Y.^2).^2;
V = -2 * mu * X .* Y ./ (X.^2 + Y.^2).^2;
% Stream function for the doublet
psi = -mu * Y ./ (X.^2 + Y.^2);
% Plotting
figure;
hold on;
% Vector plot
quiver(X, Y, U, V, 'r');
% Contour plot
contour(X, Y, psi, 50, 'LineWidth', 1);
hold off;
You may refer to the following MathWorks documentation to know more about the functions mentioned above:
- linspace: https://www.mathworks.com/help/matlab/ref/double.linspace.html
- meshgrid: https://www.mathworks.com/help/matlab/ref/meshgrid.html
- quiver: https://www.mathworks.com/help/matlab/ref/quiver.html
- contour: https://www.mathworks.com/help/matlab/ref/contour.html
I hope the solution provided above is helpful.
