Ver el código de «Flujo potencial alrededor de un cilindro circular»

[[File:Inviscid flow around a cylinder.gif|thumb|350px]]
En [[matemáticas]], '''flujo potencial alrededor de un cilindro circular''' es una solución clásica para el [[flujo de fluidos|flujo]] de un fluido [[Viscosidad| no viscoso]], [[Compresibilidad| incompresible]] alrededor de un cilindro que es transversal al flujo. Lejos del cilindro, el flujo es unidireccional y uniforme. El flujo no tiene [[vorticidad]] y, por tanto, el [[Velocidad de flujo| campo de velocidades]] es [[Campo vectorial conservativo|irrotacional]] y puede modelizarse como un [[flujo potencial]]. A diferencia de un fluido real, esta solución indica una [[Resistencia (fluidos)|resistencia]] neta nula sobre el cuerpo, resultado conocido como [[Paradoja de D'Alembert]].

== Solución matemática ==
[[Image:CylinderVelocityPressure.png|thumb|350px|Colores: campo de presión. {{color|rojo|'''Rojo'''}} es alta y {{color|azul|'''azul'''}} es baja. Vectores de velocidad]]
[[Image:CylinderVelocityPressureQuadrant.png|thumb|350px|Vista en primer plano de un cuadrante del flujo. Colores: campo de presión. {{color|rojo|'''Rojo'''}} es alta y {{color|azul|'''azul'''}} es baja. Vectores de velocidad]].
[[Image:CylinderStreamfunctionPotential.png|thumb|350px|Campo de presión (colores), función de flujo ('''negro''') con intervalo de contorno de {{math|0.2''Ur''}} de abajo a arriba, potencial de velocidad ('''blanco''') con intervalo de contorno {{math|0.2''Ur''}} de izquierda a derecha]].
Un cilindro (o disco) de radio {{mvar|R}} se coloca en un flujo bidimensional, incompresible y no viscoso. El objetivo es encontrar el vector velocidad constante {{math|'''V'''}} y la [[presión]] {{mvar|p}} en un plano, con la condición de que lejos del cilindro el vector velocidad (relativo a [[vector unitario]] {{math|'''i'''}} y {{math|'''j'''}}) sea:<ref>{{cite book|last=Batchelor|first=George Keith|title=An Introduction to Fluid Dynamics|publisher=Cambridge University Press|date=2000|isbn=9780521663960}}}</ref>

:<math>\mathbf{V}=U\mathbf{i}+0\mathbf{j} \,,</math>

donde {{mvar|U}} es una constante, y en el límite del cilindro

<math>\mathbf{V}\cdot\mathbf{hat n}=0 \,,</math>

donde {{math|'''n̂'''}} es el [[normal de superficie|vector normal]] a la superficie del cilindro. El flujo aguas arriba es uniforme y no tiene vorticidad. El flujo es no viscoso, incompresible y tiene [[densidad]] de masa constante {{mvar|ρ}}. Por lo tanto, el flujo permanece sin vorticidad y se dice que es ''irrotacional'', con {{math|∇ × '''V''' {{=}} 0}} en todas partes. Siendo irrotacional, debe existir un [[potencial de velocidad]] {{mvar|φ}}:

:<math>\mathbf{V}=\nabla\phi \,.</math>

Siendo incompresible, {{math|∇ - '''V''' {{=}} 0}}, por lo que {{mvar|φ}} debe satisfacer la [[ecuación de Laplace]]:

:<math> \nabla^2\phi=0 \,.</math>

La solución para {{mvar|φ}} se obtiene más fácilmente en [[coordenadas polares]] {{mvar|r}} y {{mvar|θ}}, relacionadas con las [[coordenadas cartesianas]] convencionales por {{math|''x'' {{=}} ''r'' cos ''θ''}} y {{math|''y'' {{=}} ''r'' sen ''θ''}}. En coordenadas polares, la ecuación de Laplace es:

:<math>\frac{1}{r}\frac{\partial}{\partial r}\left(r \frac{\partial \phi}{\partial r}\right) + \frac{1}{r^2}\frac{\partial^2 \phi}{\partial \theta^2} = 0 \,.</math>

La solución que satisface las [[condiciones de contorno]]s es<ref>{{cite book|last=Acheson|first=David J.|title=Elementary Fluid Dynamics|url=https://archive.org/details/elementaryfluidd0000ache|publisher=Oxford University Press|date=1990|isbn=9780198596790|page=130ff}}</ref>

:<math>\phi(r,\theta)=Ur\left(1+\frac{R^2}{r^2}\right)\cos\theta\,.</math>

Las componentes de la velocidad en coordenadas polares se obtienen a partir de las componentes de {{math|∇''φ''}} en coordenadas polares:

:<math>V_r=\frac{\partial \phi}{\partial r} = U\left(1-\frac{R^2}{r^2}\right)\cos\theta</math>

y

:<math>V_\theta=\frac{1}{r}\frac{\partial \phi}{\partial \theta} = - U\left(1+\frac{R^2}{r^2}\right)\sin\theta \,.</math>

Al ser no viscosa e irrotacional, la ecuación de Bernoulli permite obtener la solución para el campo de presiones directamente a partir del campo de velocidades:

:<math> p=\tfrac12\rho\left(U^2-V^2\right) + p_\infty,</math>

donde las constantes {{mvar|U}} y {{math|''p''<sub>∞</sub>}} aparecen para que {{math|''p'' → ''p''<sub>∞</sub>}} lejos del cilindro, donde {{math|''V'' {{=}} ''U''}}. Usando ''V''<sup>2</sup> {{=}} ''V''<sub>''r''</sub><sup>2</sup> + ''V''{{su|b=''θ''|p=2}},

:<math> p=\tfrac12\rho U^2\left(2\frac{R^2}{r^2}\cos(2\theta)-\frac{R^4}{r^4}\right) + p_\infty \,. </math>

En las figuras, el campo coloreado denominado "presión" es un gráfico de

:<math> 2 \frac{p - p_\infty}{\rho U^2} =2\frac{R^2}{r^2}\cos(2\theta)-\frac{R^4}{r^4,}.</math>

En la superficie del cilindro, es  decir, cuando {{math|''r'' {{=}} ''R''}}, la presión varía desde un máximo de 1 (mostrado en el diagrama en {{color|rojo|'''rojo'''}}) en los puntos de estancamiento en {{math|''θ'' {{=}} 0}} y {{math|''θ'' {{=}} π}} a un mínimo de -3 (mostrado en {{color|azul|'''azul'''}}) en los lados del cilindro, en {{math|''θ'' {{=}} {{sfrac|π|2}} y {{math|''θ'' {{=}} {{sfrac|3π|2}}. 

Del mismo modo, {{mvar|''V''}} varía de {{math|''V'' {{=}} 0}}}} en los puntos de estancamiento a {{math|''V'' {{=}} 2''U''}}}} en los laterales, en la baja presión.

=== Función de flujo ===
Siendo el flujo incompresible, se puede encontrar una [[función de corriente]] tal que

:<math>\mathbf{V}=\nabla\psi \times \mathbf{k} </math>

De esta definición se deduce, utilizando [[Identidad vectorial| identidades vectoriales]],

:<math>\mathbf{V}\cdot\nabla{\psi}=0 \,.</math>

Por tanto, un contorno de valor constante de {{mvar|ψ}} será también una línea de corriente, una línea tangente a {{math|'''V'''}}. Para el flujo más allá de un cilindro, encontramos:

:<math> \psi= U \left( r - \frac{R^2}{r} \right) \sin\theta,. </math>
== Interpretación física ==
La [[ecuación de Laplace]] es lineal, y es una de las [[ecuaciones diferenciales parciales]]s más elementales. Esta sencilla ecuación da la solución completa tanto para {{math|'''V'''}} como para {{mvar|p}} debido a la restricción de irrotacionalidad e incompresibilidad. Una vez obtenida la solución para {{math|'''V'''}} y {{mvar|p}}, se puede observar la consistencia del gradiente de presión con las aceleraciones.

La [[presión dinámica]] en el punto de estancamiento aguas arriba tiene un valor de {{math|{{sfrac|1|2}}''ρU''2}}. valor necesario para decelerar el flujo libre de velocidad {{mvar|U}}. Este mismo valor aparece en el punto de estancamiento aguas abajo, esta alta presión es de nuevo necesaria para desacelerar el flujo a velocidad cero. Esta simetría surge sólo porque el flujo es completamente sin fricción.

La baja presión en los lados del cilindro es necesaria para proporcionar la [[aceleración centrípeta]] del flujo:

:<math>\frac{\partial p}{\partial r}=\frac{\rho V^2}{L} \,,</math>

donde {{mvar|L}} es el radio de curvatura del flujo.{{citation needed|date=Enero 2011}} Pero {{math|''L'' ≈ ''R''}}, y {{math|''V'' ≈ ''U''}}. La integral de la ecuación de la aceleración centrípeta sobre una distancia {{math|Δ''r'' ≈ ''R''}} dará por tanto

:<math>p-p_\infty \approx - \rho U^2 \,.</math>
La solución exacta tiene, para la presión más baja,

:<math>p-p_\infty = - \tfrac32 \rho U^2 \,.</math>

La baja presión, que debe estar [[Presente (tiempo)|presente]] para proporcionar la aceleración centrípeta, también aumentará la velocidad del flujo a medida que el fluido se desplaza de valores más altos a valores más bajos de presión. Así encontramos la velocidad máxima en el flujo, {{math|''V'' {{=}} 2''U''}}, en la baja presión en los lados del cilindro.
== Comparación con el flujo de un fluido real a través de un cilindro ==
La simetría de esta solución ideal tiene un punto de estancamiento en la parte trasera del cilindro, así como en la parte [[Delantero|delantera]]. La distribución de presiones en las caras anterior y posterior son idénticas, lo que da lugar a la peculiar propiedad de que la [[Arrastre (física)|arrastre]] sobre el cilindro sea nula, propiedad conocida como [[Paradoja de D'Alembert]]. A diferencia de un fluido ideal no viscoso, un [[Viscosidad| flujo viscoso]] que pase por un cilindro, por pequeña que sea la viscosidad, adquirirá una fina [[capa límite]] adyacente a la superficie del cilindro. Se producirá la [[Separación de la capa límite]], y existirá una [[Estela (rastro)|estela]] de arrastre en el flujo detrás del cilindro. La presión en cada punto del lado de la estela del cilindro será menor que en el lado aguas arriba, lo que provocará una fuerza de arrastre en la dirección aguas abajo.
Un valor de {{math|''V'' > ''U''}} es consistente con la conservación del volumen del fluido. Con el cilindro bloqueando parte del flujo, {{mvar|V}} debe ser mayor que {{mvar|U}} en algún punto del plano que pasa por el centro del cilindro y es transversal al flujo.
==Expansión Janzen-Rayleigh==
El problema del flujo potencial compresible sobre cilindro circular fue estudiado por primera vez por O. Janzen en 1913<ref>O. JANZEN, Beitrag zu eincr Theorie der stationaren Stromung kompressibler Flussigkeiten. Phys. Zeits., 14 (1913)</ref> y por [[Lord Rayleigh]] en 1916<ref>Rayleigh, L. (1916). I. On the flow of compressible fluid past an obstacle. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 32(187), 1-6.</ref> con pequeños efectos compresibles. Aquí, el parámetro pequeño es el cuadrado del [[número de Mach]] <math>\mathrm{M}^2=U^2/c^2\ll 1</math>, donde {{mvar|c}} es la [[velocidad del sonido]]. Entonces la solución a la aproximación de primer orden en términos del potencial de velocidad es

:<math>\phi(r,\theta) = Ur\left(1+ \frac{a^2}{r^2}\right)\cos\theta - \mathrm{M}^2 \frac{Ur}{12} \left[\left( \frac{13 a^2}{r^2} - \frac{6 a^4}{r^4} + \frac{a^6}{r^6}\right) \cos\theta + \left(\frac{a^4}{r^4} - \frac{3a^2}{r^2} \right) \cos 3\theta\right]+ \mathrm{O}\left(\mathrm{M}^4\right) \,</math>

donde <math>a</math> es el radio del cilindro.
==Flujo potencial sobre un cilindro circular con ligeras variaciones==
El análisis de perturbaciones regulares para un flujo alrededor de un cilindro con ligeras perturbaciones en las configuraciones puede encontrarse en [[Milton Van Dyke]] (1975).<ref name=vandyke>{{cite book|last=Van Dyke|first=Milton|title=Perturbation Methods in Fluid Mechanics|publisher=Parabolic Press|date=1975}}{{page needed|date=December 2017}}</ref> En lo que sigue, {{mvar|ε}} representará un pequeño parámetro positivo y {{mvar|a}} es el radio del cilindro. Para análisis y discusiones más detalladas, se remite a los lectores al libro de [[Milton Van Dyke]] de 1975 ''Perturbation Methods in Fluid Mechanics''.<ref name=vandyke/>

===Cilindro ligeramente distorsionado===
Aquí el radio del cilindro no es {{math|''r'' {{=}} ''a''}}, sino una forma ligeramente distorsionada {{math|''r'' {{=}} ''a''(1 - ''ε'' sin<sup>2</sup> ''θ'')}}. Entonces la solución por aproximación de primer orden es

:<math>\psi(r,\theta) = Ur\left(1- \frac{a^2}{r^2}\right)\sin\theta + \varepsilon \frac{Ur}{2} \left( \frac{3a^2}{r^2}\sin \theta - \frac{a^4}{r^4} \sin 3 \theta \right) + \mathrm{O}\left(\varepsilon^2\right)</math>
===Círculo ligeramente pulsante===
Aquí el radio del cilindro varía con el tiempo ligeramente así que {{math|''r'' {{=}} ''a''(1 + ''ε''&thinsp;''f''(''t''))}}. Entonces la solución por aproximación de primer orden es

:<math>\psi(r,\theta,t) = Ur\left(1- \frac{a^2}{r^2}\right)\sin\theta + \varepsilon Ur\left( \frac{a^2}{Ur} \theta f'(t) - \frac{2 a^2}{r^2} f(t) \sin \theta\right) + \mathrm{O}\left(\varepsilon^2\right)</math>
===Flujo con ligera vorticidad===
En general, la velocidad de la corriente libre {{mvar|U}} es uniforme, en otras palabras {{math|''ψ''' {{=}} ''Uy''}}, pero aquí se impone una pequeña vorticidad en el flujo exterior.
====Cizallamiento lineal====
Aquí se introduce un cizallamiento lineal en la velocidad.

:<math>\begin{align}
    \psi &= U \left(y + \frac{1}{2} \varepsilon \frac{y^2}{a}\right)\,, \\[3pt]
  \omega &= - \nabla^2 \psi = - \varepsilon \frac{U}{a} \quad \text{as } x\rightarrow -\infty\,,
\end{align}</math>

donde {{mvar|ε}} es el parámetro pequeño. La ecuación de gobierno es

:<math>\nabla^2\psi = - \omega(\psi) \,.</math>

Entonces la solución por aproximación de primer orden es
:<math>\psi(r,\theta) = Ur\left(1- \frac{a^2}{r^2}\right)\sin\theta + \varepsilon \frac{Ur}{4} \left( \frac{r}{a}(1-\cos 2\theta) + \frac{a^3}{r^3} \cos 2 \theta - \frac{a}{r}\right) + \mathrm{O}\left(\varepsilon^2\right) \,.</math>
====Cizallamiento parabólico====
Aquí se introduce una cizalla parabólica en la velocidad exterior.

:<math>\begin{align}
\psi &= U \left(y + \tfrac16 \varepsilon \frac{y^3}{a^2}\right)\,, \\
\omega &= - \nabla^2 \psi = - \varepsilon U \frac{y}{a^2} \quad \text{as } x\rightarrow-\infty \,.
\end{align}</math>

Entonces la solución a la aproximación de primer orden es
:<math>\psi(r,\theta) = Ur\left(1- \frac{a^2}{r^2}\right)\sin\theta + \varepsilon \frac{Ur}{6} \left( \frac{r^2}{a^2}\sin^2\theta - 3 r \ln r \sin \theta + \chi \right) + \mathrm{O}\left(\varepsilon^2\right) \,,</math>

donde {{mvar|χ}} es la solución homogénea de la ecuación de Laplace que restablece las condiciones de contorno.
===Cilindro ligeramente poroso===

Sea {{math|''C''<sub>ps</sub>}} el coeficiente de presión superficial para un cilindro impermeable:

:<math>C_\mathrm{ps} = \frac{p_s-p_\infty}{\tfrac{1}{2}\rho U^2}=1-4 \sin^2\theta=2\cos 2\theta -1 \,,</math>

donde {{math|''p''<sub>s</sub>}} es la presión superficial del cilindro impermeable. Sea ahora {{math|''C''<sub>pi</sub>}} el coeficiente de presión interna dentro del cilindro, entonces una ligera velocidad normal debida a la ligera porosidad viene dada por

:<math>\frac{1}{r}\frac{\partial \psi}{\partial \theta} = \varepsilon U \left(C_\mathrm{pi} - C_\mathrm{ps}\right) = \varepsilon U \left(C_\mathrm{pi} +1 - 2\cos 2\theta\right) \quad \text{at } r=a \,,</math>

pero la condición de flujo neto cero

:<math>\int_0^{2\pi} \frac{1}{r}\frac{\partial \psi}{\partial \theta} \,\mathrm{d}\theta = 0</math>

requiere que {{math|''C''<sub>pi</sub> {{=}} -1}}. Por lo tanto,
:<math>\frac{\partial \psi}{\partial \theta} = - 2\varepsilon r U \cos 2 \theta \quad \text{at } r=a \,.</math>

Entonces la solución a la aproximación de primer orden es

:<math>\psi(r,\theta) = Ur\left(1- \frac{a^2}{r^2}\right)\sin\theta - \varepsilon U \frac{a^3}{r^2} \sin 2\theta+ \mathrm{O}\left(\varepsilon^2\right) \,.</math>
===Cuasi-cilindro corrugado===
Si el cilindro tiene radio variable en la dirección axial, el eje {{mvar|z}}, {{math|''r'' {{=}} ''a'' (1 + ''ε'' sen {{sfrac|''z''|''b''}})}, entonces la solución a la aproximación de primer orden en términos del potencial de velocidad tridimensional es

:<math>\phi(r,\theta,z) = Ur\left(1+ \frac{a^2}{r^2}\right)\cos\theta - 2\varepsilon U b \frac{K_1\left(\frac{r}{b}\right)}{K_1'\left(\frac{r}{b}\right)} \cos\theta \sin \frac{z}{b} + \mathrm{O}\left(\varepsilon^2\right) \,,</math>

donde {{math|''K''<sub>1</sub>({{sfrac|''r''|''b''}})<}} es la [[Función de Bessel| función de Bessel modificada de primer orden]].
== Véase también ==
* [[Transformada de Joukowsky]]
* [[Condición de Kutta]]
* [[Efecto Magnus]]

==Referencias ==
{{listaref}}

{{Control de autoridades}}
[[Categoría:Dinámica de fluidos]]