Flujo de Poiseuille (Grupo 9B)

De MateWiki
Revisión del 17:21 7 dic 2022 de Jesus Berlanga Serrano (Discusión | contribuciones) (Campo de temperaturas)

Saltar a: navegación, buscar
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


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 off

3 Ecuación de Navier-Stokes para fluidos incompresibles

4 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 off

8.2 Curvas de nivel del campo

9 Gradiente de la temperatura

10 Caudal