Flujo de Poiseuille Grupo 30

Trabajo realizado por estudiantes
Título	Deformaciones de una placa plana. Grupo 30-O
Asignatura	Teoría de Campos
Curso	2024-25
Autores	Ivan Ortega Perez Natalia Esteban Tezanos Ana España Franco Abdallah Attar Altarazi Guillermo Rodriguez Navadijos
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura

1 Introducción

La ley de Poiseuille, también conocida como ley de Hagen-Poiseuille, describe el flujo laminar estacionario de un líquido incompresible. En este caso, analizaremos el flujo de un líquido incompresible a través de una tubería cilíndrica con un radio de 3, lo que implica una sección transversal circular constante. Este flujo depende del gradiente de presión y del radio de la tubería.

Para desarrollar este análisis, hemos utilizado el software Matlab, que nos ha permitido representar gráficamente los resultados, como secciones transversales y gradientes, de manera visual. Esto facilita al lector una mejor comprensión de la Ley de Poiseuille, ayudándole a interpretar y entender sus implicaciones de forma clara y didáctica.

2 Dibujar un mallado 2D de la sección longitudinal

Se trata de graficar la sección longitudinal de la tubería en coordenadas [math] x_{1} = 0 [/math], [math] \left ( \rho,z \right )\epsilon \left [ 0,4 \right ]\times \left [ 0,10 \right ]. [/math]

x=0:0.05:2;  %Creamos Vectores
y=0:0.2:10;
[XX,YY]=meshgrid(x,y);  %Creamos Malla
mesh(XX,YY,0*XX);  %Representamos la sección
axis([0,4,0,10]);  %Rango de los ejes
xlabel('ρ') ;
ylabel('z') ;
view(2);
title ('Malla de la Sección Longitudinal');

4 Ecuación de Navier-Stokes

La velocidad de las partículas de nuestro fluido viene dada por el campo [math]\vec{u}(\rho,\theta,z)= f\left(\rho\right)\vec{e_{z}}[/math], y la presión por [math]p\left(x,y\right)=p_{1}+\left (p_{2}-p_{1}\right)(z-1)/4 [/math], donde [math] p_{1} [/math] es la presión en los puntos [math] z=1 [/math], [math] p_{2} [/math] la presión en los puntos [math] z=5 [/math].

Ambas magnitudes, [math] \left ( \vec{u},\rho \right ) [/math], cumplen la ecuación estacionaria de Navier-Stokes, independiente del tiempo, donde [math] \mu [/math] es el coeficiente de viscosidad de fluido: [math] \left ( \vec{u}\cdot \triangledown \right )\vec{u}+\triangledown p=\mu\Delta\vec{u}, [/math] Comprobamos que [math]f\left ( \rho \right ) [/math] satisface la siguiente ecuación diferencial, despreciendo la parte convectiva, que corresponde con el primer término. [math] \frac{1}{\rho} \frac{\partial }{\partial \rho }\left ( \rho \frac{\partial f\left ( \rho \right )}{\partial \rho } \right )=\frac{p_{2}-p_{1}}{\mu}. [/math] Resolvemos multiplicando por [math] \rho [/math] e integrando 2 veces:

1) Multiplicamos por [math] \rho [/math]

[math] \frac{\partial }{\partial \rho}\left ( \rho \frac{\partial f\left ( \rho \right )}{\partial \rho } \right )= \rho\frac{\left ( p_{2}-p_{1} \right )}{\mu} [/math]

2) Integramos

[math] \int \frac{\partial }{\partial \rho}\left ( \rho \frac{\partial f\left ( \rho \right )}{\partial \rho } \right ) d\rho = \frac{\left ( p_{2}-p_{1} \right )}{\mu}\int\rho d\rho [/math]

[math] \rho \frac{\partial f\left ( \rho \right )}{\partial \rho } = \frac{\left ( p_{2}-p_{1} \right )}{\mu} \cdot \frac{\rho^{2}}{2} [/math]

[math]\frac{\partial f\left ( \rho \right )}{\partial \rho } =\frac{p_{2}-p_{1}}{2\mu}\cdot \rho[/math]

3)Integramos por segunda vez

[math] \int \left ( \frac{\partial f\left ( \rho \right )}{\partial \rho } \right ) d\rho = \frac{\left ( p_{2}-p_{1} \right )}{2\mu}\int\rho d\rho [/math]

[math] f\left(\rho \right )=\frac{\left ( p_{2}-p_{1} \right )}{2\mu} \cdot \left (\frac{\rho^{2}}{2} + c \right ) [/math]

Para darle valor a la constante, usamos el dato de que en [math] \rho = 3 [/math] la velocidad es cero, por tanto, como la velocidad es [math] f\left ( \rho \right )\vec{e_{z}}, [/math] entonces [math] f\left ( 3 \right ) [/math] debe ser cero. La [math] f\left ( \rho \right ) [/math] nos queda: [math] f\left(\rho \right )=\frac{\left ( p_{2}-p_{1} \right )}{4\mu} \cdot \left (\frac{\rho^{2}}{1} - 9 \right ) [/math]

5 Verificación de la condición de incompresibilidad

Dado que el agua es un líquido incompresible, su volumen debe permanecer constante, lo que implica que su densidad no varía. Para garantizar esta propiedad, se verifica que la divergencia del campo de velocidades sea cero. Esto se debe a que, en un fluido, la divergencia del campo de velocidades en un punto refleja la variación de la densidad del fluido en ese punto. [math] \triangledown \cdot \vec{u}=0 [/math] [math]\vec{u}(\rho,\theta,z)= f\left(\rho\right)\vec{e_{z}}[/math]

[math]\triangledown\cdot\vec{u} \left(\rho,\theta,z\right)=\frac{1}{\rho} \left ( \frac{\partial }{\partial \rho}\left ( \rho \cdot u_{\rho } \right )+\frac{\partial }{\partial \theta}\left ( u_{\theta } \right ) + \frac{\partial }{\partial z}\left (\rho \cdot u_{z}\right)\right) \triangledown\cdot\vec{u} \left(\rho,\theta,z\right)=\frac{1}{\rho}\left ( \frac{\partial }{\partial z}\left(\rho \cdot f\left(\rho\right)\right)\right)[/math]

\triangledown\cdot\vec{u} \left(\rho,\theta,z\right)= \frac{1}{\rho} \frac{\partial}{\partial \rho} \left[ \frac{(p_2 - p_1)}{4\mu} \rho \left( \rho^2 - 9 \right) \right] = 0.</math> (para cualquier valor de [math]\rho[/math]).