Visualización de campos en fluidos (Grupo C12)

1 Introducción. Ecuación de Navier-Stokes

Vamos a considerar el flujo de un fluido incompresible a través de un canal con paredes rectas. Trabajaremos en el plano con coordenadas cartesianas. Para empezar definimos un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1] El codigo matlab para realizarlo es el siguiente

x=0:0.1:4;    %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);     %Hacemos el mallado
mesh(xx,yy,xx*0) 
axis([0,4,-1,2])  %Definimos los ejes

Mallado de puntos

%Dibujo del campo
x=0:0.1:4;    %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);   
%SÛlido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%solido despuÈs de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])

.......

%Dibujo del campo
x=0:0.1:4;    %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);   
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
figure(3)
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy))
uy=0*xx;
quiver(xx,yy,ux,uy)
axis([0,4,-1,2])

Campo de presiones

Campo de presiones proyectado sobre z=0

Campo de velocidades

2 Campos de velocidades, presiones y temperatura

Este artículo está en versión beta. El autor de este artículo no lo ha terminado todavía, por favor no lo edites hasta que elimine este mensaje.