Couette flow between two plate
显示 更早的评论
Hi, how can i plot the velocity profile for steady couette flow between two plates? d^2u/d^2y=(1/2*miu)*dp/dx for dp/dx=[_2,_1,0,1,2] , miu=1
2 个评论
Image Analyst
2022-1-15
How about the plot function?
plot(velocity, 'b-');
Garland Garland
2022-1-15
采纳的回答
Underlying equation is
d^2u/dy^2 = -G/mu with G = -dp/dx = const
function main
dpdx = [-2 -1 0 1 2];
mu = 1.0;
U = 0.01;
H = 0.05;
dy = 0.001;
np = numel(dpdx);
y = (0:dy:H).';
ny = numel(y);
u = zeros(ny,np);
for ip=1:np
u(:,ip) = -dpdx(ip)/(2*mu)*y.*(H-y) + U*y/H;
end
plot(y,u)
end
3 个评论
Garland Garland
2022-1-15
Maybe making the plot this way is better:
plot(u,y)
xlabel('Velocity u')
ylabel('Plate height h')
Garland Garland
2022-1-15
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Fluid Network Interfaces Library 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
