Couette flow between two plate

7 visualizaciones (últimos 30 días)
Garland Garland
Garland Garland el 15 de En. de 2022
Comentada: Garland Garland el 15 de En. de 2022
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 comentarios
Image Analyst
Image Analyst el 15 de En. de 2022
How about the plot function?
plot(velocity, 'b-');
Garland Garland
Garland Garland el 15 de En. de 2022
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

Iniciar sesión para comentar.

Respuesta aceptada

Torsten
Torsten el 15 de En. de 2022
Editada: Torsten el 15 de En. de 2022
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 comentarios
Torsten
Torsten el 15 de En. de 2022
Editada: Torsten el 15 de En. de 2022
Maybe making the plot this way is better:
plot(u,y)
xlabel('Velocity u')
ylabel('Plate height h')
Garland Garland
Garland Garland el 15 de En. de 2022
Great

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre Vector Fields en Help Center y File Exchange.

Etiquetas

Community Treasure Hunt

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

Start Hunting!

Translated by