Usuario discusión:Camino del Blanco
Bienvenido a MateWiki! Esperamos que contribuyas mucho y bien. Probablemente desearás leer las páginas de ayuda. Nuevamente, bienvenido y diviértete! Herraiz (discusión) 21:58 30 nov 2013 (CET)
1 Visualización de campos escalares y vectoriales en fluidos. (Grupo 23-C)
 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.
|
2 Introducción
Vamos a estudiar el flujo de un fluido incompresible alrededor de un obstáculo con forma circular de radio igual a 1, ósea un círculo unidad. Para la realización de este trabajo usaremos tanto coordendas cartesianas como con coordenadas Cilíndricas (polares).
3 Mallado
La región ocupada por el fluido será el exterior del círculo unidad tal y como indica a continuación el mallado que hemos dibujado. Hemos tomado un mallado del anillo comprendido entre los radio 1 y 5 y centro el origen, con unos ejes dentro del intervalo [-4,4]x[-4,4].
El código matlab con el que lo hemos realizado es el siguiente:
u=linspace(1,5,100); % Creamos el vector u en el intervalo[1,5]
v=linspace(0,2*pi,100); % Creamos el vector v en el intervalo [0,2*pi]
[uu,vv]=meshgrid(u,v); % matrices de las coordenadas u y v
xx=uu.*cos(vv); % Parametrizamos en coordendas cartesianas y así definimos el anillo
yy=uu.*sin(vv);
mesh(xx,yy,0*xx) % Dibujamos el mallado
hold on
x=cos(v);
y=sin(v);
plot(x,y,'k','Linewidth',1) % Dibujamos el obstáculo: círculo unidad
axis([-4,4,-4,4]) % Definimos los ejes en [-4,4]x[-4,4]
view(2) % Visualizamos lo dibujado en el plano x e y
