En dinámica de fluidos, las ecuaciones de Oseen o flujo de Oseen describen el flujo de un fluido viscoso e incompresible en números de Reynolds pequeños, tal como lo formuló Carl Wilhelm Oseen en 1910. El flujo de Oseen es una descripción mejorada de estos flujos, en comparación con el flujo de Stokes con la inclusión (parcial) de la aceleración convectiva.[1]
El trabajo de Oseen se basa en los experimentos de GG Stokes, que había estudiado la caída de una esfera a través de un fluido viscoso. Desarrolló un término de corrección, que incluía factores de inercia, para la velocidad del flujo utilizada en los cálculos de Stokes, para resolver el problema conocido como paradoja de Stokes. Su aproximación conduce a una mejora de los cálculos de Stokes.
Ecuaciones
Las ecuaciones de Oseen son, en el caso de un objeto que se mueve con una velocidad de flujo constante U a través del fluido, que está en reposo lejos del objeto, y en un marco de referencia unido al objeto:[1]
![{\displaystyle {\begin{aligned}-\rho \mathbf {U} \cdot \nabla \mathbf {u} &=-\nabla p\,+\,\mu \nabla ^{2}\mathbf {u} ,\\\nabla \cdot \mathbf {u} &=0,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e9c22bb7b32681ce1b1bfb0c36e58b9e2a9779a0)
donde
- u es la perturbación en la velocidad del flujo inducida por el objeto en movimiento, es decir , la velocidad total del flujo en el marco de referencia que se mueve con el objeto es − U + u ,
- p es la presión,
- ρ es la densidad del fluido,
- μ es la viscosidad dinámica,
- ∇ es el operador de gradiente, y
- ∇2 es el operador de Laplace.
Las condiciones de contorno para el flujo de Oseen alrededor de un objeto rígido son:
![{\displaystyle {\begin{aligned}\mathbf {u} &=\mathbf {U} &&{\text{en la superficie del objeto}},\\\mathbf {u} &\to 0&&{\text{y}}\quad p\to p_{\infty }\quad {\text{para}}\quad r\to \infty ,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b2ab329918cbf3fd47a83dace32dced74fe2a395)
con r la distancia desde el centro del objeto y p∞ la presión no perturbada lejos del objeto.
Ondas longitudinales y transversales[2]
Una propiedad fundamental de la ecuación de Oseen es que la solución general puede dividirse en ondas longitudinales y transversales.
Una solución
es una onda longitudinal si la velocidad es irrotacional y, por tanto, el término viscoso desaparece. Las ecuaciones se convierten en:
En consecuencia:
![{\displaystyle {\mathbf {u} _{\text{L}}}_{t}+U{\mathbf {u} _{\text{L}}}_{x}+{\frac {1}{\rho }}\nabla p=0,\quad \nabla \cdot \mathbf {u} _{\text{L}}=0,\quad \nabla \times \mathbf {u} _{\text{L}}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1c2ee0e6decf1aa5f8731de502f4f94d1de09830)
![{\displaystyle \mathbf {u} _{\text{L}}=\nabla \phi ,\quad \nabla ^{2}\phi =0,\quad p'=p-p_{\infty }=-\rho U\mathbf {u} _{\text{L}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f6cf2831ecc678885196f1f08c2d847cbb848f8d)
La velocidad se deriva de la teoría del potencial y la presión de las ecuaciones de Bernoulli linealizadas.
La solucuión
es una onda transversal si la presión
es idéntica a cero y el campo de velocidad es solenoidal. Las ecuaciones son:
![{\displaystyle {\mathbf {u} _{\text{T}}}_{t}+U{\mathbf {u} _{\text{T}}}_{x}=\nu \nabla ^{2}\mathbf {u} _{\text{T}},\quad \nabla \cdot \mathbf {u_{T}} =0.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/81e4d9ce4358381886be36cc45ac68096308984c)
Entonces la solución completa de Oseen viene dada por
![{\displaystyle \mathbf {u} =\mathbf {u} _{\text{L}}+\mathbf {u} _{\text{T}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/46eb888af344656cb2fb9ced842ac60bda66ad86)
un teorema de división debido a Horace Lamb.[3] La división es única si están especificadas las condiciones en el infinito (digamos que
)
Para ciertos flujos de Oseen, es posible dividir aún más la onda transversal en componente irrotacional y rotacional
Dejemos que
sea la función escalar que satisface
y desaparece en el infinito y, a la inversa, dejemos que
se dé de tal manera que
, entonces la onda transversal es:
![{\displaystyle \mathbf {u} _{\text{T}}=-{\frac {\nu }{U}}\nabla \chi +\chi \mathbf {i} ,\quad \mathbf {u} _{\text{1}}=-{\frac {\nu }{U}}\nabla \chi ,\quad \mathbf {u} _{\text{2}}=\chi \mathbf {i} .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0f9cd6771f5e4227c6f56c3fef25e082ccf8a5c8)
donde
viene determinada por
y
es el vector unitario. Sin embargo, tampoco
o
son transversales por sí mismas, pero
es transversal. Por lo tanto,
![{\displaystyle \mathbf {u} =\mathbf {u} _{\text{L}}+\mathbf {u} _{\text{T}}=\mathbf {u} _{\text{L}}+\mathbf {u} _{1}+\mathbf {u} _{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/00283c33f7856fc4affddcd6965b0cee269085b6)
El único componente rotacional es
.
Soluciones fundamentales[2]
La solución fundamental debido a una fuerza puntual singular incrustada en un flujo de Oseen es el Oseenlet. Las soluciones fundamentales de forma cerrada para los flujos inestables generalizados de Stokes y Oseen asociados con movimientos de traslación y rotación arbitrarios dependientes del tiempo se han derivado para los fluidos newtonianos[4] y fluidos micropolares.[5]
Usando la ecuación de Oseen, Horace Lamb pudo derivar expresiones mejoradas para el flujo viscoso alrededor de una esfera en 1911, mejorando la ley de Stokes hacia números de Reynolds algo más altos.[1] Además, Lamb derivó, por primera vez, una solución para el flujo viscoso alrededor de un cilindro circular.[1]
La solución a la respuesta de una fuerza singular
cuando no hay límites externos presentes se escribe como
![{\displaystyle U\mathbf {u} _{x}+{\frac {1}{\rho }}\nabla p-\nu \nabla ^{2}\mathbf {u} =\mathbf {f} ,\quad \nabla \cdot \mathbf {u} =0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/29861d97915dc9cb8c09eecf52cd2ac50cffdc28)
Si
donde
es la fuerza singular concentrada en el punto
y
es un punto arbitarrio y
es el vector dado, que da la dirección de la fuerza singular, entonces en ausencia de fronteras, la velocidad y la presión se deriva del tensor fundamental
y el vector fundamental
![{\displaystyle \mathbf {u} (q)=\Gamma (q,q_{o})\mathbf {a} ,\quad p'=p-p_{\infty }=\Pi (q,q_{o})\cdot \mathbf {a} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/68e517ec9e5cd657a86cbd2be46c6a8c321ab732)
Ahora bien, si
es una función arbitraria del espacio, la solución para un dominio no limitado es
![{\displaystyle \mathbf {u} (q)=\int \Gamma (q,q_{o})\mathbf {f} (q_{o})dq_{o},\quad p'(q)=\int \Pi (q,q_{o})\cdot \mathbf {f} (q_{o})dq_{o}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3ccbeeb97d28f52f4a3c842e8d045afbbc7adbe7)
donde
es el elemento de volumen/área infinitesimal alrededor del punto
.
Bidimensional
Sin pérdida de generalidad
tomado en el origen y
. Entonces el tensor fundamental y el vector son
![{\displaystyle \Gamma ={\begin{pmatrix}{\frac {\partial A}{\partial x}}&{\frac {\partial A}{\partial y}}\\{\frac {\partial A}{\partial y}}&-{\frac {\partial A}{\partial x}}\end{pmatrix}}+{\frac {1}{2\pi \nu }}e^{\lambda x}K_{o}(\lambda r){\begin{pmatrix}1&0\\0&1\end{pmatrix}},\quad \Pi ={\frac {\rho }{2\pi }}\nabla (\ln r)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/866f49988aafcc991b812bc586cc49e26aa3e996)
donde
![{\displaystyle \lambda ={\frac {U}{2\nu }},\quad r^{2}=x^{2}+y^{2},\quad A=-{\frac {1}{2\pi U}}\left[\ln r+e^{\lambda x}K_{o}(\lambda r)\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/63558e1e6128d9852c44500a37b1d6aecd68fe20)
donde
es la Función de Bessel# Funciones de Bessel modificadas de segunda especie de orden cero.
Tridimensional
Sin pérdida de generalidad
tomado en el origen y
. Entonces el tensor fundamental y el vector son
![{\displaystyle \Gamma ={\begin{pmatrix}{\frac {\partial A}{\partial x}}&{\frac {\partial B}{\partial x}}&{\frac {\partial C}{\partial x}}\\{\frac {\partial A}{\partial y}}&{\frac {\partial B}{\partial y}}&{\frac {\partial C}{\partial y}}\\{\frac {\partial A}{\partial z}}&{\frac {\partial B}{\partial z}}&{\frac {\partial C}{\partial z}}\end{pmatrix}}+{\frac {1}{4\pi \nu }}{\frac {e^{-\lambda (r-x)}}{r}}{\begin{pmatrix}1&0&0\\0&1&0\\0&0&1\end{pmatrix}},\quad \Pi =-{\frac {\rho }{4\pi }}\nabla \left({\frac {1}{r}}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9cc37338004adb58527ef4ca2763e1a62d872965)
donde
![{\displaystyle \lambda ={\frac {U}{2\nu }},\quad r^{2}=x^{2}+y^{2}+z^{2},\quad A={\frac {1}{4\pi U}}{\frac {1-e^{-\lambda (r-x)}}{r}},\quad B=-{\frac {1}{4\pi U}}{\frac {\left[1-e^{-\lambda (r-x)}\right]y}{r(r-x)}},\quad C=-{\frac {1}{4\pi U}}{\frac {\left[1-e^{-\lambda (r-x)}\right]z}{r(r-x)}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/74d65eb9a913571ada7747cbe4d0e30996fc1511)
Cálculos
Oseen consideró que la esfera estaba estacionaria y que el fluido fluía con una velocidad de flujo (
) a una distancia infinita de la esfera. Los términos de inercia se ignoraron en los cálculos de Stokes.[6] Es una solución límite cuando el número de Reynolds tiende a cero. Cuando el número de Reynolds es pequeño y finito, como 0,1, es necesario corregir el término inercial. Oseen sustituyó los siguientes valores de velocidad de flujo en las ecuaciones de Navier-Stokes.
![{\displaystyle u_{1}=u+u_{1}',\qquad u_{2}=u_{2}',\qquad u_{3}=u_{3}'.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/872760d9d0e24ab39fac1934133d6a497ea422e1)
Si se insertan en las ecuaciones de Navier-Stokes y se desprecian los términos cuadráticos de las cantidades primitivas, se obtiene la aproximación de Oseen:
![{\displaystyle u{\partial u_{1}' \over \partial x_{1}}=-{1 \over \rho }{\partial p \over \partial x_{1}}+\nu \nabla ^{2}u_{i}'\qquad \left({i=1,2,3}\right).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/faa5f84fbba41e3003036a9d2afcc4c5285a9759)
Como el movimiento es simétrico respecto al eje
y la divergencia del vector de vorticidad es siempre cero obtenemos:
![{\displaystyle \left(\nabla ^{2}-{U \over 2v}{\partial \over \partial x}\right)\chi =G(x)=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1460fb8b32145407423651ca50e683fa3ba84cb8)
la función
se puede eliminar añadiendo a una función adecuada en
, es la función de vorticidad, y la función anterior se puede escribir como:
![{\displaystyle {U \over v}{\partial u' \over \partial x}=\nabla ^{2}u'}](https://wikimedia.org/api/rest_v1/media/math/render/svg/92dc015f8baf9a04b8be202d41f25f0a68e044fb)
y por una cierta integración la solución para
es:
![{\displaystyle e^{-Ux \over 2v}\chi ={{Ce^{-U\operatorname {Re} \over 2v}} \over \operatorname {Re} }}](https://wikimedia.org/api/rest_v1/media/math/render/svg/089b40e621bc2b154cff5fcd79a7d042fb811d9c)
por lo que al dejar que
sea la "dirección preferente" produce:
![{\displaystyle \varphi ={A_{0} \over \operatorname {Re} }+A_{1}{\partial \over \partial x}{1 \over \operatorname {Re} }+A_{2}{\partial ^{2} \over \partial x^{2}}{1 \over \operatorname {Re} }+\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/40e970286d0f554891ee2952ffa4606f6f67e652)
entonces aplicando las tres condiciones de contorno se obtiene
![{\displaystyle C=-{3 \over 2}Ua,\ A_{0}=-{3 \over 2}va,\ A_{1}={1 \over 4}Ua^{3}\ {\text{, etc.}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/052769f20005769152852d4504908ffc93261a7c)
el nuevo coeficiente de arrastre mejorado pasa a ser:
![{\displaystyle C_{\text{d}}={12 \over \operatorname {Re} }\left(1+{3 \over 8}\operatorname {Re} \right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3df92f18f41cc5ef1cb65a1dd1ffc8875d296e01)
y finalmente, cuando se resolvió la solución de Stokes sobre la base de la aproximación de Oseen, se demostró que la resultante fuerza de arrastre viene dada por
![{\displaystyle F=6\pi \,\mu \,au\left(1+{3 \over 8}\operatorname {Re} \right),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a5362ed3fd04940f5a6c4635b62c60e3d06271fd)
donde
es el número de Reynolds basado en el radio de la esfera, ![{\displaystyle a}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ffd2487510aa438433a2579450ab2b3d557e5edc)
es lafuerza hidrodinámica
es la velocidad de flujo
es la viscosidad del fluido
La fuerza de la ecuación de Oseen difiere de la de Stokes en un factor de
![{\displaystyle 1+{3 \over 8}\operatorname {Re} .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/46e487ee42fc960f8151d578457d208ccb959468)
Error en la solución de Stokes
Las ecuaciones de Navier Stokes dicen:[7]
![{\displaystyle {\begin{aligned}\triangledown u'&~=0\\u\triangledown u'&~=-\triangledown p+\nu \triangledown ^{2}u',\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/803e8082eb9bcc501b0bcef5d480670f01966b4e)
pero cuando el campo de velocidad es
![{\displaystyle {\begin{aligned}u_{y}&=u\cos \theta \left({1+{a^{3} \over 2r^{3}}-{3a \over 2r}}\right)\\u_{z}&=-u\sin \theta \left({1-{a^{3} \over 4r^{3}}-{3a \over 4r}}\right).\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b0ff2d7acbad18d1ddcc6d63986287dfddeedff7)
En el campo lejano
≫ 1, la tensión viscosa está dominada por el último término. Es decir:
![{\displaystyle \triangledown ^{2}u'=O\left({a^{3} \over r^{3}}\right).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cc5c5d956d2973d876bfc417f05e63d8feb58b86)
El término de inercia está dominado por el término
![{\displaystyle u{\partial u' \over \partial z_{1}}\sim O\left({a^{2} \over r^{2}}\right).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a65ab38757892b6150f431ccd68a2fabe33251e0)
El error viene dado entonces por la relación:
![{\displaystyle u{{\partial u' \over \partial z_{1}} \over {\nu \triangledown ^{2}u'}}=O\left({r \over a}\right).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/092a2b4c00c970c22f8e729ca88c1163fb3e5b06)
Esto se convierte en ilimitado para
≫ 1, por lo tanto la inercia no puede ser ignorada en el campo lejano. Tomando el rizo, la ecuación de Stokes da
Dado que el cuerpo es una fuente de vorticidad,
se volvería logaritmo ilimitado para grandes
Esto es ciertamente antifísico y se conoce como paradoja de Stokes.
Solución para una esfera en movimiento en un fluido incompresible
Considérese el caso de una esfera sólida que se mueve en un líquido estacionario con una velocidad constante. El líquido se modela como un fluido incompresible, es decir, con densidad constante, y ser estacionario significa que su velocidad tiende a cero a medida que la distancia desde la esfera se acerca al infinito.
Para un cuerpo real habrá un efecto transitorio debido a su aceleración cuando comienza su movimiento; sin embargo, después de un tiempo suficiente, tenderá a cero, por lo que la velocidad del fluido en todas partes se acercará a la obtenida en el caso hipotético en el que el cuerpo ya se está moviendo por un tiempo infinito.
Por lo tanto, supongamos una esfera de radio a que se mueve a una velocidad constante
, en un fluido incompresible que está en reposo en el infinito. Trabajaremos en coordenadas
que se mueven junto con la esfera con el centro de coordenadas ubicado en el centro de la esfera. Se tiene:
![{\displaystyle {\begin{aligned}{\vec {u}}\left(\left\|{{\vec {x}}_{m}}\right\|=a\right)&={\vec {U}}\\{\vec {u}}\left(\left\|{{\vec {x}}_{m}}\right\|\rightarrow \infty \right)&\rightarrow 0\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/111f2610974948267e81674ba96ad26cdbc8f666)
Dado que estas condiciones de contorno, así como la ecuación de movimientos, son invariantes en el tiempo (es decir, no cambian al desplazar el tiempo
) cuando se expresan en las coordenadas
, la solución depende del tiempo sólo a través de estas coordenadas.
Las ecuaciones de movimiento son las ecuaciones de Navier-Stokes definidas en las coordenadas del marco de reposo
. Mientras que las derivadas espaciales son iguales en ambos sistemas de coordenadas, la derivada temporal que aparece en las ecuaciones satisface:
![{\displaystyle {\frac {\partial {\vec {u}}\left({\vec {x}},t\right)}{\partial t}}=\sum _{i}{{\frac {d{x_{m}}_{i}}{dt}}{\frac {\partial {\vec {u}}\left({\vec {x}}_{m}\right)}{\partial {x_{m}}_{i}}}}=-\left({\vec {U}}\cdot {\vec {\nabla }}_{m}\right){\vec {u}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bb87f8702a901c322ba1291cacd7382191fe6ee5)
donde la derivada
es con respecto a las coordenadas móviles
. En adelante omitiremos el subíndice m.
La aproximación de Oseen se resume en despreciar el término no lineal en
. Así, las ecuaciones de Navier-Stokes incompresibles se convierten en:
![{\displaystyle \left({\vec {U}}\cdot {\vec {\nabla }}\right){\vec {u}}+\nu \nabla ^{2}{\vec {u}}={\frac {1}{\rho }}{\vec {\nabla }}p}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fab1d835ddd673f981dab4964d9cb01773dc530a)
para un fluido que tiene densidad ρ y Viscosidad cinemática ν = μ/ρ (μ es la Viscosidad dinámica). p es la presión.
Debido a la Ecuación de continuidad para un fluido incompresible
, la solución se puede expresar mediante un potencial vectorial
. Este resulta estar dirigido en la dirección
y su magnitud es equivalente a la función de flujo utilizada en problemas bidimensionales. Resulta ser:
![{\displaystyle {\begin{aligned}\psi &=Ua^{2}\left(-{\frac {a}{4r^{2}}}\sin \theta +3{\frac {1-\cos \theta }{r\sin \theta }}{\frac {1-e^{-{\frac {Rr}{4a}}(1+\cos \theta )}}{R}}\right)\\{\vec {u}}&={\vec {\nabla }}\times (\psi {\hat {\varphi }})={\frac {1}{r\sin \theta }}{\frac {\partial }{\partial \theta }}\left(\psi \sin \theta \right){\hat {r}}-{\frac {1}{r}}{\frac {\partial }{\partial r}}\left(r\psi \right){\hat {\theta }}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b84bb8f23483c8a016b79584c254f72a869554a2)
donde
es el número de Reynolds para el flujo cercano a la esfera.
Nótese que en algunas notaciones
se sustituye por
para que la derivación de
a partir de
sea más similar a su derivación a partir de la función de flujo en el caso bidimensional (en coordenadas polares).
Elaboración
puede expresarse de la forma siguiente
![{\displaystyle \psi =\psi _{1}+\psi _{2}-\psi _{2}e^{-kr(1+\cos \theta )}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ca8c5c7a4c86871a2723b26fb47c19b591dc19d9)
donde:
![{\displaystyle {\begin{aligned}\psi _{1}&\equiv -{\frac {Ua^{3}}{4r^{2}}}\sin \theta \\\psi _{2}&\equiv {\frac {3Ua^{2}}{Rr}}{\frac {1-\cos \theta }{\sin \theta }}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d555f1591d98dcf574871cf2db447dbbcfc8ce85)
, so that
.
El vector Laplaciano de un vector del tipo
reads:
.
Por lo tanto, se puede calcular que:
![{\displaystyle {\begin{aligned}\nabla ^{2}\left(\psi _{1}{\hat {\varphi }}\right)&=0\\\nabla ^{2}\left(\psi _{2}{\hat {\varphi }}\right)&=0\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a98aa29874dc59ab7ba884d7c9418fe7e4373456)
Por lo tanto,
![{\displaystyle {\begin{aligned}\nabla ^{2}{\vec {\psi }}&=-\nabla ^{2}\left(\psi _{2}e^{-kr(1+\cos \theta )}{\hat {\varphi }}\right)\\&=-\left(\psi _{2}\nabla ^{2}e^{-kr(1+\cos \theta )}+2{\frac {\partial \psi _{2}}{\partial r}}{\frac {\partial }{\partial r}}e^{-kr(1+\cos \theta )}+{\frac {2}{r^{2}}}{\frac {\partial \psi _{2}}{\partial \theta }}{\frac {\partial }{\partial \theta }}e^{-kr(1+\cos \theta )}\right){\hat {\varphi }}\\&=-{\frac {6Ua^{2}}{R}}\sin \theta \left({\frac {k^{2}}{r}}+{\frac {k}{r^{2}}}\right)e^{-kr(1+\cos \theta )}{\hat {\varphi }}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f53c55ee5218bd22a874548fdaa64904525d35a0)
Así la vorticidad es:
![{\displaystyle {\vec {\omega }}\equiv {\vec {\nabla }}\times {\vec {u}}=-\nabla ^{2}{\vec {\psi }}={\frac {6Ua^{2}}{R}}\sin \theta \left({\frac {k^{2}}{r}}+{\frac {k}{r^{2}}}\right)e^{-kr(1+\cos \theta )}{\hat {\varphi }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a9fa469d795ff030205d23dd95b7e22e30226c85)
donde hemos utilizado la vanishing of the divergence de
para relacionar el vector laplaciano y un doble curl.
El lado izquierdo de la ecuación del movimiento es el rizo de lo siguiente:
![{\displaystyle \left({\vec {U}}\cdot {\vec {\nabla }}\right){\vec {\psi }}+\nu \nabla ^{2}{\vec {\psi }}=\left({\vec {U}}\cdot {\vec {\nabla }}\right){\vec {\psi }}-\nu {\vec {\omega }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a0691dd2defc27dbcd65255106e1eb50fcba86ff)
Calculamos la derivada por separado para cada término en
.
Obsérvese que:
![{\displaystyle {\vec {U}}=U\left(\cos \theta {\hat {r}}-\sin \theta {\hat {\theta }}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dda8c304ba711e2d052d3226f583fe7806db8550)
y también
![{\displaystyle {\begin{aligned}{\frac {\partial \psi _{2}}{\partial r}}&=-{\frac {1}{r}}\psi _{2}\\\sin \theta {\frac {\partial \psi _{2}}{\partial \theta }}&=\psi _{2}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/db976b21da78f95d396603c7c0ee2763c2c793b2)
Por lo tanto, tenemos:
![{\displaystyle {\begin{aligned}\left({\vec {U}}\cdot {\vec {\nabla }}\right)\left(\psi _{1}{\hat {\varphi }}\right)&=U\left(\cos \theta {\frac {\partial \psi _{1}}{\partial r}}-{\frac {1}{r}}\sin \theta {\frac {\partial \psi _{1}}{\partial \theta }}\right){\hat {\varphi }}={\frac {3U^{2}a^{3}}{4r^{3}}}\sin \theta \cos \theta {\hat {\varphi }}\\\left({\vec {U}}\cdot {\vec {\nabla }}\right)\left(\psi _{2}{\hat {\varphi }}\right)&=U\left(\cos \theta {\frac {\partial \psi _{2}}{\partial r}}-{\frac {1}{r}}\sin \theta {\frac {\partial \psi _{2}}{\partial \theta }}\right){\hat {\varphi }}=-U{\frac {1}{r}}(1+\cos \theta )\psi _{2}{\hat {\varphi }}=-{\frac {3U^{2}a^{2}}{Rr^{2}}}\sin \theta {\hat {\varphi }}\\\left({\vec {U}}\cdot {\vec {\nabla }}\right)\left(-\psi _{2}e^{-kr(1+\cos \theta )}{\hat {\varphi }}\right)&=-e^{-kr(1+\cos \theta )}\left(\left({\vec {U}}\cdot {\vec {\nabla }}\right)\left(\psi _{2}{\hat {\varphi }})+\psi _{2}\left({\vec {U}}\cdot {\vec {\nabla }}\right)\left(-kr(1+\cos \theta \right){\hat {\varphi }}\right)\right)\\&=U\psi _{2}e^{-kr(1+\cos \theta )}\left({\frac {1}{r}}(1+\cos \theta )+\cos \theta {\frac {\partial (kr(1+\cos \theta ))}{\partial r}}-{\frac {1}{r}}\sin \theta {\frac {\partial (kr(1+\cos \theta ))}{\partial \theta }}\right){\hat {\varphi }}\\&=U\psi _{2}(1+\cos \theta )\left({\frac {1}{r}}+k\right)e^{-kr(1+\cos \theta )}{\hat {\varphi }}={\frac {3U^{2}a^{2}}{R}}\sin \theta \left({\frac {1}{r^{2}}}+{\frac {k}{r}}\right)e^{-kr(1+\cos \theta )}{\hat {\varphi }}\\&={\frac {U}{2k}}{\vec {\omega }}=\nu {\vec {\omega }}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1d46e5df6aeb3bdf2409756495bf3345a7717b0e)
Combinando todos los términos que tenemos:
![{\displaystyle \left({\vec {U}}\cdot {\vec {\nabla }}\right){\vec {\psi }}+\nu \nabla ^{2}{\vec {\psi }}=\left({\frac {3U^{2}a^{3}}{4r^{3}}}\sin \theta \cos \theta -{\frac {3U^{2}a^{2}}{Rr^{2}}}\sin \theta \right){\hat {\varphi }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d624a51f6462a48307f5668a8aaaf36a6f7f81b2)
Tomando el rizo, encontramos una expresión que es igual a
por el gradiente de la siguiente función, que es la presión:
![{\displaystyle p=p_{0}-{\frac {3\mu Ua}{2r^{2}}}\cos \theta +{\frac {\rho U^{2}a^{3}}{4r^{3}}}\left(3\cos ^{2}\theta -1\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ebc27609f5acbe2fcd5c908e1a98df1b6524baff)
donde
es la presión en el infinito,
.es el ángulo polar originado desde el lado opuesto al punto de estancamiento frontal (
donde es el punto de estancamiento frontal).
Además, la velocidad se deriva tomando el rizo de
:
![{\displaystyle {\vec {u}}=U\left[-{\frac {a^{3}}{2r^{3}}}\cos \theta +{\frac {3a^{2}}{Rr^{2}}}-{\frac {3a^{2}}{R}}\left({\frac {1}{r^{2}}}+{\frac {k}{r}}[1-\cos \theta ]\right)e^{-kr(1+\cos \theta )}\right]{\hat {r}}-U\left[{\frac {a^{3}}{4r^{3}}}\sin \theta +{\frac {3a^{2}}{Rr}}k\sin \theta e^{-kr(1+\cos \theta )}\right]{\hat {\theta }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/92697fd8a92a753ba8bbffc2da167ca626d9d1b8)
Estos p y u satisfacen la ecuación de movimiento y, por tanto, constituyen la solución de la aproximación de Oseen.
Modificaciones a la aproximación de Ossen
Se puede cuestionar, sin embargo, si el término de corrección fue elegido por casualidad, porque en un marco de referencia que se mueve con la esfera, el fluido cerca de la esfera está casi en reposo, y en esa región la fuerza de inercia es despreciable y la ecuación de Stokes está bien justificado. [6] Lejos de la esfera, la velocidad del flujo se aproxima a u y la aproximación de Oseen es más precisa.[6] Pero la ecuación de Oseen se obtuvo aplicando la ecuación para todo el campo de flujo. Esta pregunta fue respondida por Proudman y Pearson en 1957,[8] quien resolvió las ecuaciones de Navier-Stokes y dio una solución de Stokes mejorada en la vecindad de la esfera y una solución de Oseen mejorada en el infinito, y emparejó las dos soluciones en una supuesta región común de su validez. Obtuvieron:
![{\displaystyle F=6\pi \,\mu \,aU\left(1+{3 \over 8}\operatorname {Re} +{9 \over 40}\operatorname {Re} ^{2}\ln \operatorname {Re} +{\mathcal {O}}\left(\operatorname {Re} ^{2}\right)\right).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e4f3b0b0bf1729959176c79e9ef8e8ff091ba6e9)
Aplicaciones
El método y la formulación para el análisis del flujo a un número de Reynolds muy bajo es importante. El movimiento lento de pequeñas partículas en un fluido es común en bioingeniería. La formulación de arrastre de Oseen puede utilizarse en relación con el flujo de fluidos en varias condiciones especiales, como: contener partículas, sedimentación de partículas, centrifugación o ultracentrifugación de suspensiones, coloides y sangre mediante el aislamiento de tumores y antígenos.[6] El fluido ni siquiera tiene que ser un líquido, y las partículas no tienen que ser sólidas. Puede utilizarse en varias aplicaciones, como la formación de smog y la atomización de líquidos.
El flujo sanguíneo en vasos pequeños, como los capilares, se caracteriza por números de Reynolds y Womersley pequeños. Un vaso de 10 µm de diámetro con un flujo de 1 milímetro/segundo, viscosidad de 0,02 poises para la sangre, densidad de 1 g/cm 3 y frecuencia cardiaca de 2 Hz, tendrá un número de Reynolds de 0,005 y un número de Womersley de 0.0126. A estos pequeños números de Reynolds y Womersley, los efectos viscosos del fluido se vuelven predominantes. Comprender el movimiento de estas partículas es esencial para la administración de fármacos y el estudio de los movimientos de metástasis de los cánceres.
Referencias
- ↑ a b c d Batchelor (2000), §4.10, pp. 240–246.
- ↑ a b Lagerstrom, Paco Axel. Laminar flow theory. Princeton University Press, 1996.
- ↑ Lamb, Horace. Hydrodynamics. Cambridge university press, 1932.
- ↑ Shu, Jian-Jun; Chwang, A.T. (2001). «Generalized fundamental solutions for unsteady viscous flows». Physical Review E 63 (5): 051201. Bibcode:2001PhRvE..63e1201S. PMID 11414893. S2CID 22258027. arXiv:1403.3247. doi:10.1103/PhysRevE.63.051201.
- ↑ Shu, Jian-Jun; Lee, J.S. (2008). «Fundamental solutions for micropolar fluids». Journal of Engineering Mathematics 61 (1): 69-79. Bibcode:2008JEnMa..61...69S. S2CID 3450011. arXiv:1402.5023. doi:10.1007/s10665-007-9160-8.
- ↑ a b c Fung (1997)
- ↑ Mei (2011)
- ↑ Proudman y Pearson (1957)
Bibliografía
- Oseen, Carl Wilhelm (1910), «Über die Stokes'sche formel, und über eine verwandte Aufgabe in der Hydrodynamik», Arkiv för matematik, astronomi och fysik, vi (29) .
- Batchelor, George (2000), An introduction to fluid dynamics, Cambridge Mathematical Library (second paperback edición), Cambridge University Press, ISBN 978-0-521-66396-0, MR 1744638 .
- Fung, Yuan-cheng (1997), Biomechanics: Circulation (2nd edición), New York, NY: Springer-Verlag .
- Mei, C.C. (4 de abril de 2011), «Oseen's improvement for slow flow past a body», Advanced Environmental Fluid Mechanics (Web.Mit.edu), consultado el 28 de febrero de 2013 .
- Proudman, I.; Pearson, J.R.A. (1957), «Expansions at small Reynolds numbers for the flow past a sphere and a circular cylinder», Journal of Fluid Mechanics 2 (3): 237-262, Bibcode:1957JFM.....2..237P, doi:10.1017/S0022112057000105 .
Datos: Q2033185