Couette flow between two plate

12 次查看（过去 30 天）
Garland Garland
Garland Garland 2022-1-15
评论： Garland Garland 2022-1-15
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
Image Analyst 2022-1-15
How about the plot function?
plot(velocity, 'b-');
Garland Garland
Garland Garland 2022-1-15
U(y)=1/2*miu*(dp/dx)*y^2+(U_1/2*miu*(dp/dx))y The upper plate moves by U=0.01m/s

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Torsten
Torsten 2022-1-15
编辑：Torsten 2022-1-15
Underlying equation is
d^2u/dy^2 = -G/mu with G = -dp/dx = const
https://en.wikipedia.org/wiki/Couette_flow
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 个评论
显示 1更早的评论
Torsten
Torsten 2022-1-15
编辑：Torsten 2022-1-15
Maybe making the plot this way is better:
plot(u,y)
xlabel('Velocity u')
ylabel('Plate height h')

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Vector Fields 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by