Diferencia entre revisiones de «Visualización de campos en fluidos (Grupo C12)»
De MateWiki
(→Introducción. Ecuación de Navier-Stokes) |
(→Introducción. Ecuación de Navier-Stokes) |
||
| Línea 35: | Línea 35: | ||
[[Archivo:C-12-21.jpg|200px|thumb|right|.......]] | [[Archivo:C-12-21.jpg|200px|thumb|right|.......]] | ||
| − | {{matlab|codigo=}} | + | |
| + | |||
| + | {{matlab|codigo= | ||
| + | %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]) | ||
| + | }} | ||
| + | |||
| + | [[Archivo:C-12-311.jpg|200px|thumb|right|.......]] | ||
| + | [[Archivo:C-12-321.jpg|200px|thumb|right|.......]] | ||
| + | [[Archivo:C-12-332.jpg|200px|thumb|right|.......]] | ||
== Campos de velocidades, presiones y temperatura == | == Campos de velocidades, presiones y temperatura == | ||
Revisión del 18:25 6 dic 2013
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
%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])
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.
|
