Flujo de Couette entre dos tubos concéntricos (Grupo 35)
| Trabajo realizado por estudiantes | |
|---|---|
| Título | Flujo de Couette entre dos tubos concéntricos. Grupo 35 |
| Asignatura | Teoría de Campos |
| Curso | 2023-24 |
| Autores |
|
| Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura | |
Consideramos el flujo de un fluido incompresible a través de dos cilindros concéntricos de manera que el interior se mueve con una velocidad angular constante en sentido antihorario mientras que el exterior se encuentra fijo. Suponiendo que ambos cilindros tienen su eje en [math]OX_3[/math] y pintamos la sección transversal [math](x_3 = 0)[/math] el cilindro exterior queda proyectado sobre la la circunferencia ρ = 2 y el interior sobre la circunferencia ρ = 1. La velocidad
angular cilindro interior es ω > 0.
Contenido
1 Mallado de la sección transversal
Primero debemos representar el fluido definido en el enunciado acorde a los parámetros especificados en el enunciado. Las directrices requeridas son las siguientes: Tomar los ejes (comando axis) en el rectángulo [math](x, y) ∈ [-5;5] × [-5;5][/math] y para ver desde arriba view(2).
rho=1:0.1:3; % Intervalos de rho [1,3]
theta=0:0.1:2*pi; % Intervalo de theta [0,2*pi]
[RHO,THETA]=meshgrid(rho,theta); % Creación de matrices de los parámetros
figure(1)
x=RHO.*cos(THETA); % Parametrización de x
y=RHO.*sin(THETA); % Parametrización de y
mesh(x,y,0*x) % Dibujo de la figura con la coordenadas ya parametrizadas
axis([-5,5,-5,5]) % Región de vista del dibujo
view(2) % Perspectiva desde arriba
2.1 Definición del campo de velocidades
La velocidad de las partículas del fluido viene dada por [math]\vec u(\rho,\theta)=f(\rho)\vec e_\theta[/math] y su presión p (constante). Sabemos que [math](\vec u, p)[/math] satisface la ecuación de Navier-Stokes estacionaria:
y que se desprecia el primer termino, parte de la ecuación estacionaria, tenemos la siguiente igualdad: [math]\mu\triangle\vec u=\vec 0[/math] .
2.2 Cálculo del Laplaciano
En esta igualdad vamos a usar el Laplaciano [math]\triangle \vec u[/math] en cilíndricas ya que es así como nos lo da el enunciado, por lo tanto la formula se nos queda en: [math]\triangle\vec u = \triangledown(\triangledown\cdot\vec u)-\triangledown\times(\triangledown\times\vec u)[/math]
Usamos esta alternativa ya que no hemos visto el laplaciano de un campo de vectores cuando ese campo esta definido en la base cartesiana, tambien porque con esta ecuacion podemos usar la divergencia y el rotacional, los cuales son requeridos a lo largo del problema.
2.2.1 Cálculo de la divergencia
Calculamos la divergencia de [math]\vec u[/math]: [math]\triangledown\cdot\vec u=\frac{1}{\rho}[\frac{\partial}{\partial\rho}(0)+{\frac{\partial}{\partial\theta}(f(\rho))}+\frac{\partial}{\partial z}(0)=0][/math] Esto nos indica la condición de incompresibilidad, que sucede al ser la divergencia igual a 0.
2.2.2 Cálculo del rotacional
Ahora procedemos a calcular el rotacional de [math]\vec u[/math], y posteriormente el rotacional del rotacional:
[math]\ \triangledown \times \vec{u}= \frac{1}{\rho} \cdot \left|\begin{matrix} \vec e_\rho & \rho\cdot \vec e_\theta & \vec e_z \\ \frac{\partial}{\partial \rho } & \frac{\partial}{\partial \theta} & \frac{\partial}{\partial z} \\ 0 & \rho\cdot f(\rho ) & 0 \end{matrix}\right| = \left [ \frac{f\left ( \rho \right )}{\rho } + \frac{\partial \left ( f\left ( \rho \right ) \right )}{\partial \rho }\right ] \bar{e_{z}}= \left ( \frac{f\left ( \rho \right )}{\rho } + {f}'\left ( \rho \right )\right ) \bar{e_{z}} [/math] ;
[math]\ \triangledown \times \left ( \triangledown \times \vec{u} \right )=\frac{1}{\rho }\begin{vmatrix} \vec{e}_{\rho }&\rho \cdot \vec{e} _{\theta } & \vec{e}_{z}\\ \frac{\partial }{\partial \rho } & \frac{\partial }{\partial \theta } & \frac{\partial }{\partial z} \\ 0& 0 &\frac{f\left ( \rho \right )}{\rho }+{f}'\left ( \rho \right ) \end{vmatrix}= -\frac{1}{\rho }\left [ -\frac{1}{\rho ^{2}}f\left ( \rho \right ) +\frac{{f}'\left ( \rho \right )}{\rho }+{f}''\left ( \rho \right )\right ]\left ( \rho \right ) \vec{e}_{\theta } = \left [\frac{1}{\rho ^{2}}f\left ( \rho \right )-{\frac{1}{\rho}f'\left ( \rho \right )} - f''\left ( \rho \right )\right ] \vec{e}_{\theta } [/math]
2.2.3 Cálculo final
Con el rotacional y el rotacional del rotacional, tendríamos todos los cálculos necesarios para obtener el Laplaciano [math]\triangle \vec u[/math] , y por lo tanto la ecuación de Navier Stokes.
[math]\ \triangle \vec{u}=-\triangledown \times \triangledown \times \vec{u} =0-\left [\frac{1}{\rho ^{2}}f\left ( \rho \right )-{\frac{1}{\rho}f'\left ( \rho \right )} - f''\left ( \rho \right )\right ] \vec{e}_{\theta } = -\frac{1}{\rho }\left [ -\frac{1}{\rho ^{2}}f\left ( \rho \right ) +\frac{{f}'\left ( \rho \right )}{\rho }+{f}''\left ( \rho \right )\right ]\left ( \rho \right ) \vec{e}_{\theta } = -\frac{1}{\rho }\left [-\frac{1}{\rho}f\left ( \rho \right )+{f'\left ( \rho \right )} +\rho f''\left ( \rho \right )\right ] \vec{e}_{\theta } [/math]
[math] → \ \triangle \vec{u}= -\frac{1}{\rho }\left [ -\frac{1}{\rho ^{2}}f\left ( \rho \right ) +\frac{{f}'\left ( \rho \right )}{\rho }+{f}''\left ( \rho \right )\right ]\left ( \rho \right ) \vec{e}_{\theta } = \left [\frac{1}{\rho ^{2}}f\left ( \rho \right )-{\frac{1}{\rho}f'\left ( \rho \right )} - f''\left ( \rho \right )\right ] \vec{e}_{\theta } [/math]
La ecuación sería la siguiente:
[math]\mu \triangle \vec{u}= 0 \rightarrow \left [\frac{1}{\rho ^{2}}f\left ( \rho \right )-{\frac{1}{\rho}f'\left ( \rho \right )} - f''\left ( \rho \right )\right ]\vec{e}_{\theta }=\vec{0} [/math]
Después de calcular la Ec. de Navier-Stokes:
Tenemos la siguiente ecuación diferencial: [math]\ \frac{\partial }{\partial \rho }\left ( \rho \frac{\partial f\left ( \rho \right )}{\partial \rho } \right )= \frac{f(\rho)}{\rho} [/math]
2.3.1 Comprobante de que f(ρ) cumple la ecuación diferencial
[math]\ \frac{\partial }{\partial \rho }\left ( \rho \frac{\partial f\left ( \rho \right )}{\partial \rho } \right )= \frac{\partial }{\partial \rho } (\rho\cdot f'(\rho))=f'(\rho)+\rho(f''(\rho)) \rightarrow f'(\rho)+\rho(f''(\rho))=\frac{f(\rho)}{\rho} [/math]
Si pasamos el cociente [math]\frac{f(\rho)}{\rho}[/math] al otro lado de la igualdad, vemos que coincide con la ecuación calculada anteriormente.
2.3.2 Comprobante que la función es solución
Ahora vamos a comprobar si la función [math]f(\rho)=a\rho+\frac{b}{\rho}[/math] es solución de la ecuación diferencial mencionada anteriormente, en el caso de ser así, la velocidad de la frontera del fluido coincida con la frontera de los cilindros, tanto interior como exterior.
[math]\ \frac{\partial }{\partial \rho }\left ( \rho \frac{\partial f\left ( \rho \right )}{\partial \rho } \right )= \frac{\partial }{\partial \rho } (\rho\cdot (a-\frac{b}{\rho^2}))= \frac{\partial }{\partial \rho }(\rho\cdot a-\frac{b}{\rho }) = a + \frac{b}{\rho ^2} [/math]
[math]f'(\rho)= a-\frac{b}{\rho^2}[/math]
[math]f''(\rho)=\frac{2b}{\rho^3}[/math]
[math]\ \frac{1}{\rho}(a\cdot\rho +\frac{b}{\rho })=a-\frac{b}{\rho ^2}+\rho (\frac{2b}{\rho^3}) = a+\frac{b}{\rho^2} [/math]
Se comprueba que la función mencionada si que es solución, ya que coincide con la ecuación diferencial calculada:
[math]\ \frac{1}{\rho}(a\cdot\rho +\frac{b}{\rho })=a+\frac{b}{\rho ^2}\rightarrow \frac{f(\rho)}{\rho}=f'(\rho)+\rho f''(\rho) [/math]
2.3.3 Determinar los valores "a" y "b"
Recordamos que la formula de la velocidad viene dada por [math]\vec u(\rho,\theta)=f(\rho)\vec e_\theta[/math] y [math]f(\rho)=a\rho+\frac{b}{\rho}[/math] .
Para hallar estos valores, tenemos que adaptar la función a las características del problema, es decir, darle los valores de [math]\rho[/math] de las paredes de la tubería en las dos condiciones, los cuales son para [math]\rho =1[/math] y para [math]\rho =2[/math], e igualar la función a las velocidades que tienen los cilindros exterior(w) e interior(0):
Condición 1: Si [math]ρ=1 \Longrightarrow \vec{u}=\omega \vec{e_θ} \Longrightarrow f(1)\vec{e_θ}=\omega \vec{e_θ} \Longrightarrow (a+b)\vec{e_θ}=\omega \vec{e_θ}\Longrightarrow
a+b=\omega[/math]
Condición 2: Si [math]ρ=2 \Longrightarrow \vec{u}= \vec{0} \Longrightarrow f(2)\vec{e_θ}=\vec{0} \Longrightarrow (2a+\frac{b}{2})\vec{e_θ}=\vec{0}\Longrightarrow 2a+\frac{b}{2}= 0[/math]
Obtenemos como solucion de cada ecuación, un sistema de dos ecuaciones, iniciando con que [math]\ (a = -\frac{b}{4}\omega)[/math]
Con lo que tenemos finalmente obtenemos que: [math]\ a = -\frac{1}{3}\omega [/math] y [math]\ b = \frac{4}{3}\omega [/math]
Quedando la ecuación diferencial de Navier-Stokes de siguiente forma:
