Diferencia entre revisiones de «Flujo de Poiseuille (Grupo 9B)»
(→Curvas de nivel del campo) |
(→Gradiente de la temperatura) |
||
| Línea 93: | Línea 93: | ||
=== Gradiente de la temperatura === | === Gradiente de la temperatura === | ||
| + | {{matlab|codigo= | ||
| + | x=0:0.2:2; | ||
| + | y=0:0.2:10; | ||
| + | [X,Y]=meshgrid(x,y); | ||
| + | figure(1); | ||
| + | p=1+((X-(1/2)).^2).*(exp(-(Y-1).^2)); | ||
| + | [pX,pY]=gradient(p); | ||
| + | hold on | ||
| + | quiver(X,Y,pX,pY) | ||
| + | axis([0,3,0,10]); | ||
| + | title('Gradiente de temperatura'); | ||
| + | shading flat | ||
| + | grid on | ||
| + | hold off | ||
| + | }} | ||
| + | |||
=== Caudal === | === Caudal === | ||
Revisión del 17:34 7 dic 2022
| Trabajo realizado por estudiantes | |
|---|---|
| Título | Flujo de Poiseuille . Grupo 9-B |
| Asignatura | Teoría de Campos |
| Curso | 2022-23 |
| Autores | Jesus Berlanga Serrano,Iñigo Castells Gómez, Javier Azañedo Guisado, Adriana Fernández Rivas |
| Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura | |
Contenido
1 Introducción
La Ley de Poiseuille (o de Hagen-Poiseuille) es una ecuación hemodinámica fundamental. Debido a que la longitud de las tuberías y la viscosidad son relativamente constantes, el flujo viene determinado básicamente por el gradiente de presión y por el radio. Dicha ecuación está formulada para flujos laminares de fluidos homogéneos con viscosidad constante, si la velocidad del flujo es alta o si el gradiente de presión es elevado, se pueden generar remolinos o turbulencias que modifican el patrón del flujo.
Para la creación de este artículo, se ha hecho uso de programa informático Matlab para la representación gráfica de mallados, gradientes, etc.
2 Mallado de la sección
Este mallado de dimensión 2 muestra la mitad de la sección longitudinal de la tubería, centrada en el eje OZ.
x=0:0.2:2; %Vector X.
y=0:0.2:10; %Vector Y.
[xx,yy]=meshgrid(x,y); %Mallado XY.
hold on
mesh(xx,yy,0.*xx); %Representamos la tubería
axis([0,3,0,10]); %Rango de los ejes
xlabel('ρ') ;
ylabel('z') ;
view(2);
title ('Mallado de la sección');
hold off4 Campo de presiones y el campo de velocidades
5 Líneas de corriente
6 Comportamiento del módulo de la velocidad
7 Rotacional
x=0:0.2:2;
y=0:0.2:10;
[xx,yy]=meshgrid(x,y);
rot=abs(-(xx.^2/4)-1);
hold on
surf(xx,yy,rot);
colorbar;
view(2);
axis([0,3,0,10]);
title('Rotacional del campo');
hold off
8 Campo de temperaturas
8.1 Campo
x=0:0.2:2;
y=0:0.2:10;
[X,Y]=meshgrid(x,y);
rho=sqrt(X.^2+Y.^2);
hold on
T=1+((rho-(1/2)).^2).*(exp(-(Y-1).^2));
surf(X,Y,T);
view
axis([0,3,0,10]);
colorbar
title('Campo de temperaturas')
hold off8.2 Curvas de nivel del campo
x=0:0.2:2;
y=0:0.2:10;
[X,Y]=meshgrid(x,y);
rho=sqrt(X.^2+Y.^2);
hold on
T=1+((rho-(1/2)).^2).*(exp(-(Y-1).^2));
contour(X,Y,T,10);
axis([0,3,0,10]);
view(2);
title('Curvas de nivel de temperatura')
colorbar
hold off
9 Gradiente de la temperatura
x=0:0.2:2;
y=0:0.2:10;
[X,Y]=meshgrid(x,y);
figure(1);
p=1+((X-(1/2)).^2).*(exp(-(Y-1).^2));
[pX,pY]=gradient(p);
hold on
quiver(X,Y,pX,pY)
axis([0,3,0,10]);
title('Gradiente de temperatura');
shading flat
grid on
hold off