|
|
| (No se muestran 67 ediciones intermedias de 4 usuarios) |
| Línea 1: |
Línea 1: |
| − | {{ TrabajoED | Flujo de Poiseuille (Grupo 10) | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC23/24|2023-24]] | Lucía Domínguez Álvarez; Eduardo Juarranz del Valle; Adrián Díaz Gadea; Pablo Amado Silva; Carmen Pardos Martínez }}
| |
| | | | |
| − | ===Introducción===
| |
| − | La ley de Poiseuille o también llamada ley de Hagen-Poiseuille es una ley que permite determinar el flujo laminar estacionario de un líquido incomprensible, en nuestro caso consideraremos el flujo de un líquido incomprensible en una tubería cilíndrica de radio 2, teniendo así una sección circular constante. Por lo tanto, dicho flujo vendrá determinado por el gradiente de presión y el radio.
| |
| − |
| |
| − | Para la realización y creación de este artículo, hemos hecho uso del programa informático Matlab para poder representar gráficamente todos los resultados, como por ejemplo gráficos de secciones, gradientes... de una manera visual para así poder ayudar al usuario a comprender la Ley de Poiseuille ayudándole a comprender las interpretaciones de dicha ley.
| |
| − |
| |
| − | [[Archivo:Sección_longitudinal_tubería.png|350px|miniaturadeimagen|right|''Sección longitudinal'']]
| |
| − |
| |
| − | ===Sección Longitudinal de la Tubería===
| |
| − | Mallado de la represenación de la sección longitudinal de la tubería <math> x_{1} = 0 </math>, <math> \left ( \rho,z \right )\epsilon \left [ 0,3 \right ]\times \left [ 0,10 \right ]. </math>
| |
| − | {{matlab|codigo=
| |
| − | 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,3,0,10]); %Rango de los ejes
| |
| − | xlabel('ρ') ;
| |
| − | ylabel('z') ;
| |
| − | view(2);
| |
| − | title ('Malla de la Sección Longitudinal');
| |
| − | }}
| |
| − | ===Ecuación de Navier-Stokes===
| |
| − |
| |
| − | ===Campo de presiones y campo de velocidades.===
| |
| − | A continuación calcularemos el campo de presiones (campo escalar) y el campo de velocidades (campo vectorial).
| |
| − | Suponiendo que <math> p_{1}=4, p_{2}=1 </math> y <math>\mu=1. </math>
| |
| − | Donde <math> p_{1} </math> es la presión en los puntos <math> z=4 </math>, <math> p_{2} </math> la presión en los puntos <math> z=1 </math> y <math> \mu </math> el coeficiente de viscosidad del fluido.
| |
| − |
| |
| − | =====Campo de presiones.=====
| |
| − |
| |
| − | Para calcular el ''' Campo de presiones''' hacemos uso de la ecuación de presión <math>p\left ( x,y \right ) </math>.
| |
| − |
| |
| − | <math>p\left ( x,y \right )=p_{1}+\left ( p_{2}-p_{1} \right )\left ( z-1 \right ){/2}=4+\left ( 1-4 \right )\left ( z-1 \right ){/2}=\frac{-3z}{2}+\frac{11}{2} </math>
| |
| − |
| |
| − | ======Representación del campo de presiones.======
| |
| − | Para poder representar el campo de presiones debemos estudiar su comportamiento de la presión frente a la altura. Como podemos ver hay una relación lineal entre ambas. Cuanto más aumenta la profundidad, más aumenta la presión y de igual forma, si disminuye una, también lo hace la otra.
| |
| − |
| |
| − |
| |
| − | {{matlab|codigo=
| |
| − | clc;
| |
| − | clear all;
| |
| − | z=0:0.1:10;
| |
| − | f=(-3)*z/2+(11/2); %Falta comprobar
| |
| − | plot(z,f)
| |
| − | xlabel('Variación de altura');
| |
| − | ylabel('Variación de presión');
| |
| − | title(' Gráfica del campo de presiones');
| |
| − | }}
| |
| − |
| |
| − | =====Campo de velocidades.=====
| |
| − | Para calcular el ''' Campo de velocidades''' haremos uso de la ecuación que representa la velocidad de las partículas del fluido. Es decir por <math>\vec{u}(\rho,\theta,z) </math>
| |
| − |
| |
| − | <math>\vec{u}(\rho,\theta,z)= f\left(\rho\right)\vec{e_{z}}= </math>
| |
| − |
| |
| − | ======Representación del campo de velocidades.======
| |
| − | Para proceder con la representación del campo de velocidades debemos tener en cuenta que este es de carácter vectorial. Por lo que para representarlo usando Matlab o Octave deberemos pasar primero de coordenadas cilíndricas (ya que se trata de una tubería cilíndrica) a coordenadas cartesianas.
| |
| − |
| |
| − | {{matlab|codigo=
| |
| − | clear all;
| |
| − | clc;
| |
| − | x=0:0.1:3;
| |
| − | y=0:0.1:10;
| |
| − | [xx,yy]=meshgrid(x,y);
| |
| − | ux=.; %Cambio a Coord. Cartesianas
| |
| − | uy=.; %Cambio a Coord. Cartesianas
| |
| − | hold on % Debo rellenar la parte de arriba del código
| |
| − | quiver(xx,yy,ux,uy)
| |
| − | axis([0,5,0,12])
| |
| − | hold off
| |
| − | colorbar;
| |
| − | view(2)
| |
| − | title(' GRÁFICA DEL CAMPO DE VELOCIDADES')
| |
| − | }}
| |
| − |
| |
| − | ===Líneas de corriente del campo.===
| |
| − |
| |
| − |
| |
| − | ===Puntos con módulo de velocidad máxima.===
| |
| − | La ley de Poiseuille establece que el caudal de un fluido incompresible en un tubo cilíndrico de radio constante es proporcional a la diferencia de presión entre los extremos del tubo, al radio del tubo y a la inversa de la viscosidad del fluido. La velocidad del fluido en un punto del tubo se puede calcular a partir del caudal y del radio del tubo.
| |
| − | El módulo de la velocidad del fluido es:
| |
| − | <center>''añadir ecuación''</center>
| |
| − |
| |
| − | Para encontrar los puntos donde el módulo de la velocidad del fluido es máximo, debemos derivar la expresión anterior con respecto a ρ. Obtenemos:
| |
| − | <center><math> \frac{\partial \vec{u}} {\partial \rho}=\frac{-\rho}{2}\vec{e_z} </math></center>
| |
| − | <center><math> \frac{\partial \vec{u} }{\partial \theta}=0 </math></center>
| |
| − | <center><math> \frac{\partial \vec{u} }{\partial z}=0 </math></center>
| |
| − |
| |
| − | La derivada es igual a cero cuando ρ = 0. Por lo tanto, el módulo de la velocidad del fluido es máximo en el punto ρ = 0, que corresponde al eje del tubo.
| |
| − | Es decir, el módulo de la velocidad del fluido es máximo en el eje del tubo y su valor es proporcional al radio del tubo y a la inversa de la viscosidad del fluido.
| |
| − |
| |
| − | ====Comportamiento del módulo de la velocidad.====
| |
| − | La siguiente gráfica muestra el comportamiento del módulo de la velocidad del fluido en función de ρ. Se muestra la curva de la variación de la velocidad del fluido en función del radio de la sección de la tubería cilíndrica. Podemos ver como a medida que aumenta el radio, disminuye la velocidad. Por lo que concluimos que la velocidad máxima se situa en el eje de la tubería.
| |
| − | {{matlab|codigo=
| |
| − | rho=0:0.1:2;
| |
| − | f=abs((-(rho.^2/4)+1); % Me falta cambiar ec.
| |
| − | plot(rho,f);
| |
| − | title('Comportamiento del módulo del campo de velocidades');
| |
| − | xlabel('Radio de la sección');
| |
| − | ylabel('Variación de la velocidad');
| |
| − | axis([0,2,0,3.5])
| |
| − | }}
| |
| − |
| |
| − |
| |
| − | ===Rotacional===
| |
| − |
| |
| − | ===Campo de temperaturas y curvas de nivel===
| |
| − |
| |
| − | ===Gradiente de temperatura===
| |
| − |
| |
| − | ===Caudal===
| |
| − |
| |
| − |
| |
| − | [[Categoría:Teoría de Campos]]
| |
| − | [[Categoría:TC23/24]]
| |
