Diferencia entre revisiones de «Tubos Concéntricos»

De MateWiki
Saltar a: navegación, buscar
Línea 60: Línea 60:
 
<center><math> f(1) = 0 \Longrightarrow a + b = 0 </math></center>
 
<center><math> f(1) = 0 \Longrightarrow a + b = 0 </math></center>
 
<center><math> f(2) = w\cdot\rho \Longrightarrow 2a + \frac{b}{2} = 2w </math></center>
 
<center><math> f(2) = w\cdot\rho \Longrightarrow 2a + \frac{b}{2} = 2w </math></center>
 +
Resolvemos el sistema, resultando:
 +
<center><math> a = \frac{4w}{3} </math></center>
 +
<center><math> b = -\frac{4w}{3} </math></center>

Revisión del 02:44 3 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 OX3 y pintamos la sección transversal (x3 = 0) el cilindro exterior queda proyectado sobre la circunferencia [math] \rho = 0 [/math] y el interior sobre la circunferencia [math] \rho = 1 [/math]. La velocidad angular del cilindro exterior es [math] \omega \gt 0 [/math].

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 general de Navier-Stokes que satisface nuestro fluido es la siguiente:

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

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


Ahora bien, conocemos lo siguiente:

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


Teniendo todo esto en cuenta, nuestra fórmula 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 en coordenadas cartesianas es relativamente sencillo y se describe de la siguiente manera:

[math] \triangle\vec{u} = \triangle(u[/math]1[math]\vec{\imath} + u[/math]2[math]\vec{\jmath} + u[/math]3[math]\vec{k}) = \triangle u[/math]1[math]\vec{\imath} + \triangle u[/math]2[math]\vec{\jmath} + \triangle u[/math]3[math]\vec{k}[/math]

Sin embargo, nos encontramos en coordenadas cilíndricas, lo cuál complica los cálculos. Para calcular el laplaciano vectorial podemos utilizar la siguiente igualdad:

[math] \Delta\vec{u} = \nabla(\nabla\cdot\vec{u}) - \nabla\times(\nabla\times\vec{u}) [/math]

Sin embargo, calculísticamente es más sencillo mediante la fórmula del laplaciano vectorial en coordenadas cilíndricas, la cuál es la siguiente:

[math] \Delta\vec{u} = \vec{e}_\rho(\Delta u_\rho - \frac{u_\rho}{\rho^2} - \frac{2}{\rho^2}\frac{\partial u_\theta}{\partial\theta}) + \vec{e}_\theta(\Delta u_\theta - \frac{u_\theta}{\rho^2} + \frac{2}{\rho^2}\frac{\partial u_\rho}{\partial\theta}) + \vec{e}_z \Delta u_z [/math]

A priori, esta puede parecer más tediosa de calcular; sin embargo, una vez despreciados los términos que son igual a cero, queda de la siguiente manera:

[math] \Delta\vec{u} = \vec{e}_\theta (\Delta u_\theta - \frac{u_\theta}{\rho^2}) [/math]

Procedemos a calcular [math] \Delta u_\theta [/math]:

[math] \Delta u_\theta = \frac{1}{\rho}\frac{\partial}{\partial\rho}(\rho\frac{\partial u_\theta}{\partial\rho}) + \frac{1}{\rho^2}\frac{\partial^2 u_\theta}{\partial\theta^2} + \frac{\partial^2 u_\theta}{\partial z^2} \Longrightarrow \Delta u_\theta = \frac{f'(\rho)}{\rho} + f''(\rho) [/math]

Por lo que:

[math] \Delta\vec{u} = [-\frac{f(\rho)}{\rho^2} + \frac{f'(\rho)}{\rho} + f''(\rho)]\cdot\vec{e}_\theta = \vec{0} [/math]

2.3 Obtención de la ecuación diferencial

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

[math] -\frac{f(\rho)}{\rho^2} + \frac{f'(\rho)}{\rho} + f''(\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}(\rho\frac{\partial f(\rho)}{\partial\rho}) = \frac{f(\rho)}{\rho} [/math]

Desarrollamos la ecuación:

[math] -\frac{f(\rho)}{\rho} + f'(\rho) + \rho f''(\rho) = 0 [/math]

Nos damos cuenta de que esta ecuación es igual a la que hemos obtenido nosotros multiplicada por [math] \rho [/math], 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]

Calculamos [math] f'(\rho) [/math] y [math] f''(\rho) [/math]:

[math] f'(\rho) = a - \frac{b}{\rho^2} [/math]
[math] f''(\rho) = \frac{2b}{\rho^3} [/math]

Sustituimos en nuestra ecuación diferencial:

[math] -\frac{f(\rho)}{\rho^2} + \frac{f'(\rho)}{\rho} + f''(\rho) = 0 \Longrightarrow -\frac{a}{\rho} - \frac{b}{\rho^3} + \frac{a}{\rho} - \frac{b}{\rho^3} + \frac{2b}{\rho^3} = 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) = w\cdot\rho \Longrightarrow 2a + \frac{b}{2} = 2w [/math]

Resolvemos el sistema, resultando:

[math] a = \frac{4w}{3} [/math]
[math] b = -\frac{4w}{3} [/math]
Obtenido de «https://mat.caminos.upm.es/w/index.php?title=Tubos_Concéntricos&oldid=49436»