Diferencia entre revisiones de «Usuario:Bernardo Rodríguez»
De MateWiki
(→Campo de presiones y velocidades) |
(→Campo de presiones y velocidades) |
||
| Línea 50: | Línea 50: | ||
[[Archivo: pvijk.png ]] | [[Archivo: pvijk.png ]] | ||
| − | Se | + | Se utilizó el siguiente código para el campo de presiones |
{{Matlab|codigo= | {{Matlab|codigo= | ||
| Línea 71: | Línea 71: | ||
[[Archivo: pumpumpum.png ]] | [[Archivo: pumpumpum.png ]] | ||
| + | |||
| + | Por otro lado, seutilizó el siguiente código para el campo de velocidades | ||
| + | |||
| + | {{Matlab|codigo= | ||
| + | |||
| + | x=0:0.1:8; | ||
| + | y=-1:0.1:1; | ||
| + | [xx,yy] = meshgrid(x,y); | ||
| + | figure(1) | ||
| + | fx=-0.5*yy.^2+0.5; | ||
| + | fy=0; | ||
| + | quiver(xx,yy,fx,fy) | ||
| + | axis([0,4,-2,2]) | ||
| + | view(2) | ||
| + | |||
| + | |||
| + | |||
| + | }} | ||
| + | |||
| + | El cual nos da la siguiente gráfica: | ||
| + | |||
| + | |||
| + | [[Archivo: bienrapido.png ]] | ||
Revisión del 17:43 29 nov 2019
| Trabajo realizado por estudiantes | |
|---|---|
| Título | Trabajo 7: Visualización de campos escalares y vectoriales en fluidos (Grupo C-11) |
| Asignatura | Teoría de Campos |
| Curso | 2019-20 |
| Autores | Bernardo Rodríguez, Eduardo Martín Velásquez, Joaquín Ramón López |
| Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura | |
Contenido
1 Enunciado
Visualización de campos escalares y vectoriales en fluidos. 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.
2 Mallado
Dibujar un mallado que represente los puntos interiores del rectángulo [0, 8] × [−1, 1] ocupado por un fluido. Fijar los ejes en la regi´on [0, 4] × [−2, 2]. Para ello se utilizó el siguiente código:
x=0:0.1:8;
y=-1:0.1:1;
[xx,yy] = meshgrid(x,y);
figure(1)
mesh(xx,yy,0*xx)
axis([0,4,-2,2])
view(2)
Del cual se obtiene la siguiente gráfica.
4 Campo de presiones y velocidades
Suponiendo valores para:
p1=2 p2=1 μ=1
Se utilizó el siguiente código para el campo de presiones
x=0:0.1:8;
y=-1:0.1:1;
[xx,yy] = meshgrid(x,y);
figure(1)
f=3-xx;
surf(xx,yy,f)
axis([0,4,-2,2])
view(2)
colorbar
El cual nos da la siguiente gráfica:
Por otro lado, seutilizó el siguiente código para el campo de velocidades
x=0:0.1:8;
y=-1:0.1:1;
[xx,yy] = meshgrid(x,y);
figure(1)
fx=-0.5*yy.^2+0.5;
fy=0;
quiver(xx,yy,fx,fy)
axis([0,4,-2,2])
view(2)
El cual nos da la siguiente gráfica:
