Diferencia entre revisiones de «Tubos Concéntricos»

Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/mat/public_html/w/includes/diff/DairikiDiff.php on line 434

Revisión del 20:22 6 dic 2022

Flujo de Couette entre dos tubos concéntricos. Vamos a considerar el flujo de un fluido incompresible a través de dos cilindros concéntricos de manera que el exterior se mueve con velocidad angular constante en sentido antihorario mientras que el interior está fijo. Si suponemos que ambos cilindros tienen su eje OX₃ y pintamos la sección transversal (x₃ = 0) el cilindro exterior queda proyectado sobre la circunferencia [math] \rho = 2 [/math] y el interior sobre la circunferencia [math] \rho = 1 [/math]. La velocidad angular del cilindro exterior es [math] \omega \gt 0 [/math].

Contenido

1 Representación de la sección transversal
2 Cálculo de Velocidades
3 Representación del campo de velocidades
4 Representación de las líneas de corriente del campo
5 Cálculo de la velocidad de fluido máxima
6 Cálculo y dibujo del rotacional [math] \vec{u} [/math]
7 Representación de la temperatura

1 Representación de la sección transversal

2 Cálculo 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 es constante. Sabemos que [math] (\vec{u},p) [/math] satisface la ecuación de Navier-Stokes estacionaria.

2.1 Interpretación física del problema

La ecuación de Navier-Stokes que satisface nuestro fluido es la siguiente:

[math] \frac{\partial\vec{u}}{\partial t} + (\vec{u}\cdot\nabla)\vec{u} + \nabla p = \mu\Delta\vec{u} [/math]

Donde p es la presión y [math] \mu [/math] es la viscosidad del fluido.

Ahora bien, conocemos lo siguiente:

Que es estacionario, por lo cual el diferencial de densidad en función del tiempo es igual a cero.
Que su presión es constante, por lo que el campo de presiones, [math] \nabla p [/math], es cero.
Que debemos despreciar el segundo término, [math] (\vec{u}\cdot\nabla)\vec{u} [/math], que es la parte convectiva, por lo que podemos suponer que la viscosidad es mucho mayor.

Teniendo todo esto en cuenta, la ecuación de Navier-Stokes queda de la siguiente manera:

[math] \mu\Delta\vec{u} = \vec{0} \Longrightarrow \Delta\vec{u} = \vec{0} [/math]

2.2 Cálculo del laplaciano vectorial del campo de velocidades

El cálculo del laplaciano vectorial en coordenadas cartesianas es relativamente sencillo y se describe de la siguiente manera:

[math] \Delta\vec{u} = \Delta (u_1\vec{\imath} + u_2\vec{\jmath} + u_3\vec{k}) = \Delta u_1\vec{\imath} + \Delta u_2\vec{\jmath} + \Delta u_3\vec{k} [/math]

Sin embargo, en este caso nos encontramos en coordenadas cilíndricas, lo cual complica los cálculos. Para calcular el laplaciano vectorial podemos escribir el campo en la base cartesiana pero las componentes se pueden quedar en coordenadas cilíndricas y realizar el laplaciano en dichas coordenadas. Sabemos que:

[math] \vec{e}_\theta = -sen\theta\vec{\imath} + cos\theta\vec{\jmath} [/math]

Por lo que el campo queda de la siguiente manera:

[math] \vec{u} = f(\rho)\cdot(-sen\theta\vec{\imath} + cos\theta\vec{\jmath}) \Longrightarrow \vec{u} = -f(\rho)sen\theta\vec{\imath} + f(\rho)cos\theta\vec{\jmath} [/math]

Aplicamos el laplaciano según la fórmula descrita anteriormente. Para ello debemos calcular el laplaciano de u₁ y u₂ en coordenadas cilíndricas:

[math] \Delta u = \frac{1}{\rho}\frac{\partial}{\partial\rho}\left ( \rho\frac{\partial u}{\partial\rho} \right ) + \frac{1}{\rho^2}\frac{\partial ^2 u}{\partial\theta^2} + \frac{\partial ^2 u}{\partial z^2} [/math] [math] \Delta u_1 = -\frac{sen\theta}{\rho}\frac{\partial}{\partial\rho} \left ( \rho\frac{\partial f(\rho)}{\partial\rho}\right ) + \frac{sen\theta}{\rho^2}f(\rho) [/math] [math] \Delta u_2 = \frac{cos\theta}{\rho}\frac{\partial}{\partial\rho} \left ( \rho\frac{\partial f(\rho)}{\partial\rho}\right ) -\frac{cos\theta}{\rho^2}f(\rho) [/math]

Por lo que:

[math] \Delta\vec{u} = \left [ -\frac{sen\theta}{\rho}\frac{\partial}{\partial\rho} \left ( \rho\frac{\partial f(\rho)}{\partial\rho}\right ) + \frac{sen\theta}{\rho^2}f(\rho) \right ] \vec{\imath} + \left [ \frac{cos\theta}{\rho}\frac{\partial}{\partial\rho} \left ( \rho\frac{\partial f(\rho)}{\partial\rho}\right ) -\frac{cos\theta}{\rho^2}f(\rho) \right ] \vec{\jmath} = \vec{0} [/math]

Cambiamos las componentes a coordenadas cilíndricas:

[math] \Delta\vec{u} = \frac{1}{\rho} \left [ \frac{\partial}{\partial\rho} \left ( \rho\frac{\partial f(\rho)}{\partial\rho} \right ) -\frac{f(\rho)}{\rho^2}\right ] (-sen\theta\vec{\imath} + cos\theta\vec{\jmath}) \Longrightarrow \Delta\vec{u} = \frac{1}{\rho} \left [ \frac{\partial}{\partial\rho} \left ( \rho\frac{\partial f(\rho)}{\partial\rho} \right ) -\frac{f(\rho)}{\rho^2}\right ] \vec{e}_\theta [/math]

2.3 Obtención de la ecuación diferencial

Igualamos cada componente del vector a cero, y simplificamos, obteniendo así una ecuación diferencial:

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

2.4 Comprobaciones

El enunciado exige que comprobemos si la ecuación diferencial satisface ciertas condiciones:

Comprobar que [math] f(\rho) [/math] satisface la siguiente ecuación diferencial:

[math] \frac{\partial}{\partial\rho}\left ( \rho\frac{\partial f(\rho)}{\partial\rho}\right ) = \frac{f(\rho)}{\rho} [/math]

Nos damos cuenta de que esta ecuación es igual a la que hemos obtenido nosotros, por lo que podemos afirmar que satisface la ecuación diferencial.

Comprobar que [math] f(\rho) [/math] es solución, siendo [math] f(\rho) [/math]:

[math] f(\rho) = a\rho + \frac{b}{\rho} [/math]

Sustituimos [math] f(\rho) [/math] en la ecuación diferencial:

[math] \frac{\partial}{\partial\rho} \left ( \rho\frac{\partial \left ( a\rho + \frac{b}{\rho} \right ) }{\partial\rho} \right ) -\frac{a\rho + \frac{b}{\rho}}{\rho} = 0 \Longrightarrow \frac{\partial}{\partial\rho}\left ( a\rho - \frac{b}{\rho}\right ) - a - \frac{b}{\rho^2} = 0 \Longrightarrow a + \frac{b}{\rho^2} - a -\frac{b}{\rho} = 0 [/math]

Simplificando resulta [math] 0 = 0 [/math], por lo que [math] f(\rho) [/math] es solución.

2.5 Determinación de a y b en función de la velocidad

Se nos pide determinar a y b de manera que la velocidad en la frontera del fluido coincida con la de los cilindros interior y exterior. De esta manera, podemos determinar que:

[math] f(1) = 0 \Longrightarrow a + b = 0 [/math] [math] f(2) = \omega\cdot\rho \Longrightarrow 2a + \frac{b}{2} = 2\omega [/math]

Resolvemos el sistema, resultando:

[math] a = \frac{4\omega}{3} [/math] [math] b = -\frac{4\omega}{3} [/math]

2.6 Comprobación de la condición de incompresibilidad

Se nos pide que comprobemos la condición de incompresibilidad, es decir, que el agua siempre ocupa el mismo volumen, la cual es:

[math] \nabla\cdot\vec{u} = 0 [/math]

La divergencia de un campo vectorial en coordenadas cilíndricas viene dada por la siguiente expresión:

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

Operamos:

[math] \nabla\cdot\vec{u} = \frac{1}{\rho}\left\{\frac{\partial}{\partial\rho}(0) + \frac{\partial}{\partial\theta}(f(\rho)) + \frac{\partial}{\partial z}(0)\right\} \Longrightarrow \nabla\cdot\vec{u} = \frac{1}{\rho}(0 + 0 + 0) \Longrightarrow \nabla\cdot\vec{u} = 0 [/math]

Se cumple la condición de incompresibilidad.

3 Representación del campo de velocidades

Suponiendo que, [math] \omega = 1 [/math] y [math] \mu = 1 [/math]:

[math] a = \frac{4}{3} [/math] y [math] b = -\frac{4}{3} [/math]

Sustituimos en la función:

[math] f(\rho) = a\rho + \frac{b}{\rho} \Longrightarrow f(\rho) = \frac{4}{3}\left ( \rho - \frac{1}{\rho} \right ) [/math]

De manera que nuestro campo es:

[math] \vec{u} = f(\rho)\vec{e}_\theta \Longrightarrow \vec{u} = \frac{4}{3}\left ( \rho - \frac{1}{\rho} \right )\cdot\vec{e}_\theta [/math]

Representamos:

4 Representación de las líneas de corriente del campo

Para ello, necesitamos calcular el campo [math] \vec{v} [/math] que en cada punto es ortogonal a [math] \vec{u} [/math].

[math] \vec{v} = \vec{k}\times\vec{u} \Longrightarrow \vec{v} = \frac{4}{3}\left ( \rho - \frac{1}{\rho} \right ) (\vec{e}_z\times\vec{e}_\theta) \Longrightarrow \vec{v} = \frac{4}{3}\left ( \rho - \frac{1}{\rho} \right ) (-\vec{e}_\rho) \Longrightarrow \vec{v} = -\frac{4}{3}\left ( \rho - \frac{1}{\rho} \right )\cdot\vec{e}_\rho [/math]

4.1 Comprobación de irrotacionalidad de [math] \vec{v} [/math]

Comprobamos que, como dice el enunciado, el campo [math] \vec{v} [/math] es irrotacional debido a que la divergencia de [math] \vec{u} [/math] es nula, cosa que ya hemos calculado con anterioridad.

[math] \nabla\times\vec{v} = \begin{vmatrix} \vec{e}_\rho & \rho\vec{e}_\theta & \vec{e}_z \\ \frac{\partial}{\partial\rho} & \frac{\partial}{\partial\theta} & \frac{\partial}{\partial z} \\ -\frac{4}{3}\left ( \rho - \frac{1}{\rho} \right ) & 0 & 0 \end{vmatrix} \Longrightarrow \nabla\times\vec{v} = \vec{0} [/math]

Se comprueba la irrotacionalidad de [math] \vec{v} [/math].

4.2 Cálculo de [math] \psi [/math]

Conocemos que:

[math] \vec{v} = \nabla\psi \Longrightarrow -\frac{4}{3}\left ( \rho - \frac{1}{\rho} \right )\cdot\vec{e}_\rho = \frac{\partial\psi}{\partial\rho}\vec{e}_\rho + \frac{1}{\rho}\frac{\partial\psi}{\partial\theta}\vec{e}_\theta + \frac{\partial\psi}{\partial z}\vec{e}_z [/math]

Despejamos [math] \psi [/math]:

[math] \frac{\partial\psi}{\partial\rho} = -\frac{4}{3}\left ( \rho - \frac{1}{\rho} \right ) \Longrightarrow \psi = \int -\frac{4}{3}\left ( \rho - \frac{1}{\rho} \right ) \partial\rho [/math]

Calculamos la integral:

[math] \psi = \frac{4}{3} \left ( ln(\rho) - \frac{\rho^2}{2} \right ) [/math]

4.3 Representación de las líneas de corriente de [math] \vec{u} [/math]

5 Cálculo de la velocidad de fluido máxima

Para hallar cuándo el módulo de la velocidad del fluido es máximo debemos maximizar el módulo del campo de velocidades. El módulo del campo de velocidades es:

[math] |\vec{u}(\rho)| = \left | \left ( \frac{4}{3}\rho - \frac{4}{3\rho} \right ) \vec{e}_\theta \right | = \frac{4}{3} \left ( \rho - \frac{1}{ρ} \right ) [/math]

Derivamos [math] |\vec{u}(\rho)| [/math] e igualamos a 0:

[math] \frac{\partial |\vec{u}(\rho)|}{\partial\rho} = \frac{4}{3} \left ( 1 + \frac{1}{\rho^2} \right ) = 0 [/math]

Observamos que esta función no posee máximos ni mínimos, de manera que procedemos a analizar en el intervalo en el cual se encuentran nuestros radios [1,2]. Observamos además, que es continua en dicho intervalo, por lo que analizamos en los extremos de este:

Gráfico módulo campo de velocidades

a=4/3;                                   % Definicion del parámetro a en función de la velocidad angular
 b=-a;                                    % Definicion del parámetro b en función de la velocidad angular
 u=1:0.05:2;                              % Definicioón del parametro u (rho)
 v=0:0.05:2*pi;                           % Definición del parametro v (theta)
 [U,V] = meshgrid(u,v);                   % Matriz de los parámetros
 x = U.*cos(V);                           % Coordenada x (cartesiana)
 y = U.*sin(V);                           % Coordenada y (cartesiana)
 f = a*U+b./U;                            % Módulo de la velocidad es función del radio
 surf(x,y,f);                             % Dibujo del campo escalar
 axis([-3,3,-3,3])                        % Definición de los ejes
 view(2);                                 % Vista en 2D

6 Cálculo y dibujo del rotacional [math] \vec{u} [/math]

El cálculo del rotacional [math] \vec{u} [/math] ya ha sido calculado previamente en el ejercicio 2 para hallar el laplaciano. Por lo tanto, suponiendo que ω=1, se tiene que:

[math](\nabla\times\vec{u}) = [\frac{\frac{4}{3} ρ -\frac{4}{3ρ}}{ρ}+\frac{\partial(\frac{4}{3} ρ -\frac{4}{3ρ})}{\partial ρ}]\vec{e_z} = \frac{8}{3}\vec{e_z}[/math]

Para saber en qué puntos el rotacional será mayor se calculará la norma del rotacional.

[math] |\nabla\times\vec{u}| = |\frac{8}{3}\vec{e_z}| = \frac{8}{3}\vec{e_z} [/math]

Debido a que el campo del rotacional es constante, en todos los puntos el rotacional es igual. Esto queda demostrado en la siguiente gráfica.

7 Representación de la temperatura

Sabiendo que la temperatura del fluido viene dada por la expresión [math]T(ρ,θ) = 1+ρ^2 sin^2θ e^{−(ρ−3/2)2} [/math] se dibuja el campo de temperaturas y las curvas de nivel:

@@ Línea 103: / Línea 103: @@
 Derivamos <math> |\vec{u}(\rho)| </math> e igualamos a 0:
-<center><math> |\vec{u}(\rho)| = \frac{4}{3} \left ( 1 + \frac{1}{\rho^2} \right ) = 0 </math></center>
+<center><math> \frac{\partial |\vec{u}(\rho)|}{\partial\rho} = \frac{4}{3} \left ( 1 + \frac{1}{\rho^2} \right ) = 0 </math></center>
 Observamos que esta función no posee máximos ni mínimos, de manera que procedemos a analizar en el intervalo en el cual se encuentran nuestros radios [1,2]. Observamos además, que es continua en dicho intervalo, por lo que analizamos en los extremos de este: