Flujo elemental

Flujo elemental es una colección de flujos básicos a partir de los cuales es posible construir flujos más complejos por superposición. Algunos de los flujos reflejan casos específicos y restricciones como incompresible o irrotacional, o ambos, como en el caso del flujo potencial.[1]

Flujo uniforme bidimensional

Uniforme
Flujo potencial. Líneas de corriente para un flujo uniforme ideal

Dada una velocidad uniforme de un fluido en cualquier posición del espacio:

V 0 = v 0 cos ( θ 0 ) e x + v 0 sin ( θ 0 ) e y {\displaystyle \mathbf {V_{0}} =v_{0}\cos(\theta _{0})\mathbf {e} _{x}+v_{0}\sin(\theta _{0})\mathbf {e} _{y}}

Este flujo es incompresible porque la velocidad es constante, las primeras derivadas de las componentes de la velocidad son cero y la divergencia total es cero: v = 0 {\displaystyle \nabla \cdot \mathbf {v} =0}

Dado que la circulación es siempre cero el flujo es también irrotacional, podemos deducirlo del teorema de la circulación de Kelvin y del cálculo explícito de la vorticidad:

ω z = v x y v y x = 0 {\displaystyle \omega _{z}={\frac {\partial v_{x}}{\partial y}}-{\frac {\partial v_{y}}{\partial x}}=0}

Al ser incompresible y bidimensional, este flujo se construye a partir de una función de flujo:

v x = ϕ x {\displaystyle v_{x}=-{\frac {\partial \phi }{\partial x}}}
v y = ϕ y {\displaystyle v_{y}=-{\frac {\partial \phi }{\partial y}}}

de donde

ψ = v 0 sin ( θ 0 ) x + v 0 cos ( θ 0 ) y {\displaystyle \psi =-v_{0}\sin(\theta _{0})x+v_{0}\cos(\theta _{0})y}

y en coordenadas cilíndricas

v r = ϕ r {\displaystyle v_{r}={\frac {\partial \phi }{\partial r}}}
v θ = 1 r ϕ θ {\displaystyle v_{\theta }={\frac {1}{r}}{\frac {\partial \phi }{\partial \theta }}}

de donde

ψ = v 0 r sin ( θ θ 0 ) {\displaystyle \psi =-v_{0}r\sin(\theta -\theta _{0})}

Como es habitual la función de flujo está definida hasta un valor constante que aquí tomamos como cero. También podemos confirmar que el flujo es irrotacional a partir de:

2 ψ = 0 {\displaystyle \nabla ^{2}\psi =0}

Siendo irrotacional, la función potencial es en cambio:

v x = ϕ x {\displaystyle v_{x}=-{\frac {\partial \phi }{\partial x}}}
v y = ϕ y {\displaystyle v_{y}=-{\frac {\partial \phi }{\partial y}}}

y por lo tanto

ϕ = v 0 cos ( θ 0 ) x v 0 sin ( θ 0 ) y {\displaystyle \phi =-v_{0}\cos(\theta _{0})x-v_{0}\sin(\theta _{0})y}

y en coordenadas cilíndricas

v r = ϕ r {\displaystyle v_{r}={\frac {\partial \phi }{\partial r}}}
v θ = 1 r ϕ θ {\displaystyle v_{\theta }={\frac {1}{r}}{\frac {\partial \phi }{\partial \theta }}}
ϕ = v 0 r cos ( θ θ 0 ) {\displaystyle \phi =-v_{0}r\cos(\theta -\theta _{0})}

Fuente bidimensional de líneas de flujo

Fuente puntual
Flujo potencial. Líneas de corriente para una fuente lineal ideal

El caso de una línea vertical que emite a velocidad fija una cantidad constante de fluido Q por unidad de longitud es una fuente lineal. El problema tiene una simetría cilíndrica y puede tratarse en dos dimensiones en el plano ortogonal.

Las fuentes lineales y los sumideros lineales (abajo) son flujos elementales importantes porque desempeñan el papel de monopolo(s) para fluidos incompresibles (que también pueden considerarse ejemplos de campos solenoidales es decir, campos libres de divergencia). Los patrones de flujo genéricos también pueden descomponerse en términos de expansión multipolar, del mismo modo que para los campos eléctrico y magnético, donde el monopolo es esencialmente el primer término no trivial (por ejemplo, constante) de la expansión.

Este patrón de flujo también es irrotacional e incompresible.

Se caracteriza por una simetría cilíndrica:|

v = v r ( r ) e r {\displaystyle \mathbf {v} =v_{r}(r)\mathbf {e} _{r}}

Donde el flujo total saliente es constante

S v d S = 0 2 π ( v r ( r ) e r ) ( e r r d θ ) = 2 π r v r ( r ) = Q {\displaystyle \int _{S}\mathbf {v} \cdot d\mathbf {S} =\int _{0}^{2\pi }(v_{r}(r)\,\mathbf {e} _{r})\cdot (\mathbf {e} _{r}\,r\,d\theta )=\!2\pi \,r\,v_{r}(r)=Q}

Por lo tanto,

v r = Q 2 π r {\displaystyle v_{r}={\frac {Q}{2\pi r}}}

Se deriva de una función de flujoion

ψ ( r , θ ) = Q 2 π θ {\displaystyle \psi (r,\theta )=-{\frac {Q}{2\pi }}\theta }

o de una función potencial

ϕ ( r , θ ) = Q 2 π ln r {\displaystyle \phi (r,\theta )=-{\frac {Q}{2\pi }}\ln r}

Sumidero de línea bidimensional

El caso de una línea vertical que absorbe a un ritmo fijo una cantidad constante de fluido Q por unidad de longitud es un sumidero lineal. Todo es igual al caso de una línea fuente a parte del signo negativo. v r = Q 2 π r {\displaystyle v_{r}=-{\frac {Q}{2\pi r}}} Esto se deriva de una función de flujo

ψ ( r , θ ) = Q 2 π θ {\displaystyle \psi (r,\theta )={\frac {Q}{2\pi }}\theta }

o a partir de una función potencial

ϕ ( r , θ ) = Q 2 π ln r {\displaystyle \phi (r,\theta )={\frac {Q}{2\pi }}\ln r}

Dado que los dos resultados son iguales a una parte de un signo menos, podemos tratar de forma transparente tanto las fuentes como los sumideros lineales con la misma corriente y funciones potenciales permitiendo que Q asuma tanto valores positivos como negativos y quitando el signo menos en la definición de Q.

Doblete bidimensional o fuente de línea dipolar

Líneas de corriente de flujo potencial para una línea ideal de doblete o dipolo

Si consideramos una línea fuente y una línea sumidero a una distancia d podemos reutilizar los resultados anteriores y la función de flujo será

ψ ( r ) = ψ Q ( r d / 2 ) ψ Q ( r + d / 2 )   d ψ Q ( r ) {\displaystyle \psi (\mathbf {r} )=\psi _{Q}(\mathbf {r} -\mathbf {d} /2)-\psi _{Q}(\mathbf {r} +\mathbf {d} /2)\ \simeq \mathbf {d} \cdot \nabla \psi _{Q}(\mathbf {r} )}

La última aproximación es de primer orden en d.

Dado d = d [ cos ( θ 0 ) e x + sin ( θ 0 ) e y ] = d [ cos ( θ θ 0 ) e r + sin ( θ θ 0 ) e θ ] . {\displaystyle \mathbf {d} =d[\cos(\theta _{0})\mathbf {e} _{x}+\sin(\theta _{0})\mathbf {e} _{y}]=d[\cos(\theta -\theta _{0})\mathbf {e} _{r}+\sin(\theta -\theta _{0})\mathbf {e} _{\theta }].}

Queda ψ ( r , θ ) = Q d 2 π sin ( θ θ 0 ) r {\displaystyle \psi (r,\theta )=-{\frac {Qd}{2\pi }}{\frac {\sin(\theta -\theta _{0})}{r}}}

La velocidad es entonces

v r ( r , θ ) = Q d 2 π cos ( θ θ 0 ) r 2 {\displaystyle v_{r}(r,\theta )={\frac {Qd}{2\pi }}{\frac {\cos(\theta -\theta _{0})}{r^{2}}}}

v θ ( r , θ ) = Q d 2 π sin ( θ θ 0 ) r 2 {\displaystyle v_{\theta }(r,\theta )={\frac {Qd}{2\pi }}{\frac {\sin(\theta -\theta _{0})}{r^{2}}}}

Y el potencial en cambio

ϕ ( r , θ ) = Q d 2 π cos ( θ θ 0 ) r {\displaystyle \phi (r,\theta )={\frac {Qd}{2\pi }}{\frac {\cos(\theta -\theta _{0})}{r}}}

Línea de vórtice bidimensional

Flujo potencial. Líneas de corriente para una línea de vórtice ideal

Este es el caso de un filamento de vórtice que gira a velocidad constante, existe una simetría cilíndrica y el problema puede resolverse en el plano ortogonal.

Al igual que en el caso anterior de las fuentes lineales, las líneas de vórtice desempeñan el papel de monopolos para flujo irrotacional.

También en este caso el flujo es a la vez irrotacional e incompresible y por lo tanto un caso de flujo potencial.

Se caracteriza por una simetría cilíndrica:

v = v θ ( r ) e θ {\displaystyle \mathbf {v} =v_{\theta }(r)\,\mathbf {e} _{\theta }}

Donde la circulación total es constante para cada línea cerrada alrededor del vórtice central

v d s = 0 2 π ( v θ ( r ) e θ ) ( e θ r d θ ) = 2 π r v θ ( r ) = Γ {\displaystyle \oint \mathbf {v} \cdot d\mathbf {s} =\int _{0}^{2\pi }(v_{\theta }(r)\,\mathbf {e} _{\theta })\cdot (\mathbf {e} _{\theta }\,r\,d\theta )=\!2\pi \,r\,v_{\theta }(r)=\Gamma }

y es cero para cualquier línea que no incluya el vórtice.

Por lo tanto, v θ = Γ 2 π r {\displaystyle v_{\theta }={\frac {\Gamma }{2\pi r}}}

Esto se deriva de una función de corriente

ψ ( r , θ ) = Γ 2 π ln r {\displaystyle \psi (r,\theta )={\frac {\Gamma }{2\pi }}\ln r}

o de una función potencial

ϕ ( r , θ ) = Γ 2 π θ {\displaystyle \phi (r,\theta )=-{\frac {\Gamma }{2\pi }}\theta }

Que es dual al caso anterior de una fuente lineal

Flujo potencial bidimensional genérico

Dado un flujo bidimensional incompresible que además es irrotacional tenemos:

2 ψ = 0 {\displaystyle \nabla ^{2}\psi =0}

Que es en coordenadas cilíndricas[2]

1 r r ( r ψ r ) + 1 r 2 2 ψ θ 2 = 0 {\displaystyle {\frac {1}{r}}{\frac {\partial }{\partial r}}\left(r{\frac {\partial \psi }{\partial r}}\right)+{\frac {1}{r^{2}}}{\frac {\partial ^{2}\psi }{\partial \theta ^{2}}}=0}

Buscamos una solución con variables separadas:

ψ ( r , θ ) = R ( r ) Θ ( θ ) {\displaystyle \psi (r,\theta )=R(r)\Theta (\theta )}

lo que da

r R ( r ) d d r ( r d R ( r ) d r ) = 1 Θ ( θ ) d 2 Θ ( θ ) d θ 2 {\displaystyle {\frac {r}{R(r)}}{\frac {d}{dr}}\left(r{\frac {dR(r)}{dr}}\right)=-{\frac {1}{\Theta (\theta )}}{\frac {d^{2}\Theta (\theta )}{d\theta ^{2}}}}

Dado que la parte izquierda sólo depende de r y la parte derecha sólo depende de θ {\displaystyle \theta } , las dos partes deben ser iguales a una constante independiente de r y θ {\displaystyle \theta } . La constante será positiva. Por tanto, r d d r ( r d d r R ( r ) ) = m 2 R ( r ) {\displaystyle r{\frac {d}{dr}}\left(r{\frac {d}{dr}}R(r)\right)=m^{2}R(r)}

d 2 Θ ( θ ) d θ 2 = m 2 Θ ( θ ) {\displaystyle {\frac {d^{2}\Theta (\theta )}{d\theta ^{2}}}=-m^{2}\Theta (\theta )}

La solución de la segunda ecuación es una combinación lineal de e i m θ {\displaystyle e^{im\theta }} y e i m θ {\displaystyle e^{-im\theta }} . Para tener una velocidad de un solo valor (y también una función de corriente de un solo valor) m debe ser un número entero positivo.

Por lo tanto la solución más genérica viene dada por

ψ = α 0 + β 0 ln r + m > 0 ( α m r m + β m r m ) sin [ m ( θ θ m ) ] {\displaystyle \psi =\alpha _{0}+\beta _{0}\ln r+\sum _{m>0}{\left(\alpha _{m}r^{m}+\beta _{m}r^{-m}\right)\sin {[m(\theta -\theta _{m})]}}}

En cambio, el potencial viene dado por

ϕ = α 0 β 0 θ + m > 0 ( α m r m β m r m ) cos [ m ( θ θ m ) ] {\displaystyle \phi =\alpha _{0}-\beta _{0}\theta +\sum _{m\mathop {>} 0}{(\alpha _{m}r^{m}-\beta _{m}r^{-m})\cos {[m(\theta -\theta _{m})]}}}

Referencias

  1. Oliver, David (14 de marzo de 2013). El corcel peludo de la física: Mathematical Beauty in the Physical World. Springer Science & Business Media. ISBN 978-1-4757-4347-0. 
  2. operador de Laplace

Bibliografía

  • Fitzpatrick, Richard (2017), Theoretical fluid dynamics, IOP science, ISBN 978-0-7503-1554-8 .
  • Faber, T.E. (1995), Fluid Dynamics for Physicists, Cambridge university press, ISBN 9780511806735 .

Bibliografía adicional

  • Batchelor, G.K. (1973), An introduction to fluid dynamics, Cambridge University Press, ISBN 978-0-521-09817-5 .
  • Chanson, H. (2009), Applied Hydrodynamics: An Introduction to Ideal and Real Fluid Flows, CRC Press, Taylor & Francis Group, Leiden, The Netherlands, 478 pages, ISBN 978-0-415-49271-3 .
  • Lamb, H. (1994) [1932], Hydrodynamics (6th edición), Cambridge University Press, ISBN 978-0-521-45868-9 .
  • Milne-Thomson, L.M. (1996) [1968], Theoretical hydrodynamics (5th edición), Dover, ISBN 978-0-486-68970-8 .

Enlaces externos

  • Richard Fitzpatrick University of Texas, Austin (2017). «Fluid Mechanics». University of Texas, Austin. Consultado el 7 de febrero de 2018. 
Control de autoridades
  • Proyectos Wikimedia
  • Wd Datos: Q48817327
  • Wd Datos: Q48817327