Ecuaciones de Oseen

concepto de mecánica de fluidos

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

editar

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]

 

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:

 

con r la distancia desde el centro del objeto y p la presión no perturbada lejos del objeto.

Ondas longitudinales y transversales[2]

editar

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:

 
 

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:

 

Entonces la solución completa de Oseen viene dada por

 

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:

 

donde   viene determinada por  y   es el vector unitario. Sin embargo, tampoco   o   son transversales por sí mismas, pero   es transversal. Por lo tanto,

 

El único componente rotacional es  .

Soluciones fundamentales[2]

editar

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

 

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  

 

Ahora bien, si   es una función arbitraria del espacio, la solución para un dominio no limitado es

 

donde   es el elemento de volumen/área infinitesimal alrededor del punto  .

Bidimensional

editar

Sin pérdida de generalidad   tomado en el origen y  . Entonces el tensor fundamental y el vector son

 

donde

 

donde   es la Función de Bessel# Funciones de Bessel modificadas de segunda especie de orden cero.

Tridimensional

editar

Sin pérdida de generalidad   tomado en el origen y  . Entonces el tensor fundamental y el vector son

 

donde

 

Cálculos

editar

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.

 

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:

 

Como el movimiento es simétrico respecto al eje   y la divergencia del vector de vorticidad es siempre cero obtenemos:

 

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:

 

y por una cierta integración la solución para   es:

 

por lo que al dejar que   sea la "dirección preferente" produce:


 

entonces aplicando las tres condiciones de contorno se obtiene

 

el nuevo coeficiente de arrastre mejorado pasa a ser:

 

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

 

donde

  es el número de Reynolds basado en el radio de la esfera,  
  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

 

Error en la solución de Stokes

editar

Las ecuaciones de Navier Stokes dicen:[7]

 

pero cuando el campo de velocidad es

 

En el campo lejano   ≫ 1, la tensión viscosa está dominada por el último término. Es decir:

 

El término de inercia está dominado por el término

 

El error viene dado entonces por la relación:

 

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

editar

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:

 

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:

 

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:

 

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:

 

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

editar

  puede expresarse de la forma siguiente

 

donde:

 
 , so that  .

El vector Laplaciano de un vector del tipo   reads:

 .

Por lo tanto, se puede calcular que:

 

Por lo tanto,

 

Así la vorticidad es:

 

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:

 

Calculamos la derivada por separado para cada término en  .

Obsérvese que:

 

y también

 

Por lo tanto, tenemos:

 

Combinando todos los términos que tenemos:

 

Tomando el rizo, encontramos una expresión que es igual a   por el gradiente de la siguiente función, que es la presión:

 

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  :

 

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

editar

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:

 

Aplicaciones

editar

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

editar
  1. a b c d Batchelor (2000), §4.10, pp. 240–246.
  2. a b Lagerstrom, Paco Axel. Laminar flow theory. Princeton University Press, 1996.
  3. Lamb, Horace. Hydrodynamics. Cambridge university press, 1932.
  4. 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. 
  5. 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. 
  6. a b c Fung (1997)
  7. Mei (2011)
  8. Proudman y Pearson (1957)

Bibliografía

editar