Flujo incompresible


Flujo incompresible

Flujo incompresible

Un flujo se clasifica en compresible e incompresible, dependiendo del nivel de variación de la densidad del fluido durante ese flujo. La incompresibilidad es una aproximación y se dice que el flujo es incompresible si la densidad permanece aproximadamente constante a lo largo de todo el flujo. Por lo tanto, el volumen de todas las porciones del fluido permanece inalterado sobre el curso de su movimiento cuando el flujo o el fluido es incompresible. En esencia, las densidades de los líquidos son constantes y así el flujo de ellos es típicamente incompresible. Por lo tanto, se suele decir que los líquidos son sustancias incompresibles. Ejemplo: una presión de 210 atm hace que la densidad del agua liquida a 1 atm cambie en sólo 1 por ciento. Cuando se analizan flujos de gas a velocidades altas, la velocidad del flujo a menudo se expresa en términos del número adimensional de Mach que se define como


En donde c es la velocidad del sonido cuyo valor es de 346 m/s en el aire a temperatura ambiente al nivel del mar. Se dice que un flujo es sónico cuando Ma=1, subsónico cuando Ma<1, supersónico cuando Ma>1, e hipersónico cuando Ma>>1. Los flujos de líquidos son incompresibles hasta un nivel alto de exactitud, pero el nivel de variación de la densidad en los flujos de gases y el nivel consecuente de aproximación que se hace cuando se modelan estos flujos como incompresibles depende del número de Mach. Con frecuencia, los flujos de gases se pueden aproximar como incompresibles si los cambios en al densidad se encuentran por debajo de alrededor de 100 m/s. Así el flujo de un gas no es necesariamente compresible.

Contenido

La ecuación de Bernoulli y un criterio para el flujo incompresible

Una de las ecuaciones más utilizadas en mecánica de fluidos es la ecuación de Bernoulli:

(0) \frac{P}{\rho}+\frac{V^2}{2}+gz = \mbox{cte.}

Se demostrara que en el limite de números de Mach muy pequeños, la ecuación isoenergética e isoentrópica para la presión se vuelve idéntica a la ecuación Bernoulli. Creando un criterio para decidir si el flujo de un gas se puede tratar como incompresible. Considerando un flujo estacionario sin esfuerzo cortante, trabajo en el eje o transferencia de calor. A estas condiciones, la presión de estancamiento es constante. Se supondrá que los cambios en elevación son despreciables. Si el fluido es incompresible, la presión en cualquier lugar se puede calcular a partir de la ecuación de Bernoulli en la forma de presión (Flujo incompresible):

(1) P = P_0 -\frac{\rho V^2}{2}

Si el fluido es compresible y un gas ideal, las presiones estática y de estancamiento están relacionadas por medio de (Flujo compresible):

(2) P = P_0 \left(1+\frac{k-1}{2}M^2\right)^\frac{k}{k-1}

Si la consideración se restringe a números de Mach menores que 1, se puede expandir el término del número de Mach es una serie infinita empleando el teorema binomial de Newton:

(3) P = P_0 \left(1+\frac{k-1}{2}M^2\right)^\frac{k}{k-1} \approx 
P_0 \left( 1+ \frac{k}{2}M^2 + \frac{k}{8}M^4 + O(M^6) \right)

De la ecuación:

(4) M^2 = \frac{V^2}{kP}

Se tiene:

(5) P_0 \approx P \left[1- \frac{V^2}{2}\left( 1+ \frac{M^2}{4} \right)\right]

Si el número de Mach es pequeño, entonces M2/4 es pequeño comparado con 1 y se puede escribir que:

(6) P_0 \approx P \left[1- \frac{V^2}{2}\right]

En consecuencia, la ecuación Bernoulli es una aproximación a la relación de presión del flujo isoenergetico e isoentropico para números de Mach pequeños. Lo preciso de esta aproximación depende de lo pequeño del número de Mach. La ecuación (5) muestra que a bajos números de Mach el error es proporcional a M2/4. si se deseara limitar el error al emplear la ecuación Bernoulli para el calculo de la presión a no más del 2 por ciento, entonces:

M < \sqrt{4 \cdot 0.02} \approx 0.283

No hay nada especial en el error del 2 por ciento. Para estimaciones gruesas, un error del 5 por ciento podría ser aceptable, en cuyo caso el número de Mach debe ser menor que 0.45. El criterio más ampliamente utilizado para el límite entre el flujo compresible y el incompresible coloca el umbral del número de Mach en 0.3: En general se puede suponer que un flujo con M < 0.3 sea incompresible.

La ecuación de Navier-StoKes para flujo isotérmico incompresible

Por definición el tensor de esfuerzo es linealmente proporcional al tensor de razón de formación. Para flujo incompresible (ρ = constante), también se supone flujo aproximadamente isotérmico sabiendo que los cambios locales en temperatura son pequeños o inexistentes; esto elimina la necesidad de una ecuación diferencial de conservación de energía. Una consecuencia de la última suposición es que las propiedades del fluido, como viscosidad dinámica μ y la viscosidad cinemática v, también son constantes. Con dichas suposiciones se puede demostrar que el tensor de esfuerzo viscoso se reduce a:

Tensor de esfuerzo viscoso para un fluido newtoniano incompresible con propiedades constantes:

(7) Tij = 2μεij

donde εij es el tensor de razón de deformación. La ecuación (7) muestra que el esfuerzo es linealmente proporcional a la deformación. En coordenadas cartesianas, se mencionan las nueve componentes del tensor de esfuerzo viscoso, seis de las cuales son independientes debido a su simetría:

(8) 
T_{ij} =
\begin{pmatrix}
  T_{xx} & T_{xy} & T_{xz} \\
  T_{yx} & T_{yy} & T_{yz} \\
  T_{zx} & T_{zy} & T_{zz}
\end{pmatrix}
=
\begin{pmatrix}
  2 \mu \frac{\partial u}{\partial x} &
  \mu \left( \frac{\partial u}{\partial y} + \frac{\partial v}{\partial x} \right) &
  \mu \left( \frac{\partial u}{\partial z} + \frac{\partial w}{\partial x} \right) \\
  \mu \left( \frac{\partial v}{\partial x} + \frac{\partial u}{\partial y} \right) &
  2 \mu \frac{\partial v}{\partial y} &
  \mu \left( \frac{\partial v}{\partial z} + \frac{\partial w}{\partial y} \right) \\
  \mu \left( \frac{\partial w}{\partial x} + \frac{\partial u}{\partial z} \right) &
  \mu \left( \frac{\partial w}{\partial y} + \frac{\partial v}{\partial z} \right) &
  2 \mu \frac{\partial w}{\partial z}
\end{pmatrix}

En coordenadas cartesianas, el tensor de esfuerzo de la ecuación de fluidos en movimiento se convierte por lo tanto en:

(9) 
\sigma_{ij} =
\begin{pmatrix}
  -P & 0 & 0 \\
  0 & -P & 0 \\
  0 & 0 & -P
\end{pmatrix}
=
\begin{pmatrix}
  2 \mu \frac{\partial u}{\partial x} &
  \mu \left( \frac{\partial u}{\partial y} + \frac{\partial v}{\partial x} \right) &
  \mu \left( \frac{\partial u}{\partial z} + \frac{\partial w}{\partial x} \right) \\
  \mu \left( \frac{\partial v}{\partial x} + \frac{\partial u}{\partial y} \right) &
  2 \mu \frac{\partial v}{\partial y} &
  \mu \left( \frac{\partial v}{\partial z} + \frac{\partial w}{\partial y} \right) \\
  \mu \left( \frac{\partial w}{\partial x} + \frac{\partial u}{\partial z} \right) &
  \mu \left( \frac{\partial w}{\partial y} + \frac{\partial v}{\partial z} \right) &
  2 \mu \frac{\partial w}{\partial z}
\end{pmatrix}

Ahora se sustituye la ecuación (8) en las tres componentes cartesianas de la ecuación de Cauchy. Considere primero la componente x, se convierte en:

(10) 
\rho \frac{Du}{Dt}
 = -\frac{\partial P}{\partial x}
 + \rho g_x
 + 2 \mu \frac{\partial^2 u}{\partial x^2}
 + \mu \frac{\partial}{\partial y} \left( \frac{\partial v}{\partial x} + \frac{\partial u}{\partial y} \right)
 + \mu \frac{\partial}{\partial z} \left( \frac{\partial w}{\partial x} + \frac{\partial u}{\partial z} \right)

Dado que la presión consiste sólo de un esfuerzo normal, únicamente aporta un término a la ecuación (10). Sin embargo, ya que el tensor de esfuerzo viscoso consiste tanto de esfuerzos normal como de corte, aporta tres términos. También en tanto las componentes de velocidad sean funciones suaves de x, y y z, el orden de diferenciación es irrelevante. Par ejemplo, la primera parte del último término en la ecuación (10) se puede reescribir como:


\mu \frac{\partial}{\partial z} \left( \frac{\partial w}{\partial x} \right) =
\mu \frac{\partial}{\partial x} \left( \frac{\partial w}{\partial z} \right)

Después de cierto reordenamiento inteligente de los términos viscosos en la ecuación (10):


\begin{align}
\rho \frac{Du}{Dt}
&= -\frac{\partial P}{\partial x} + \rho g_x + \mu \left[
  \frac{\partial^2 u}{\partial x^2} +
  \frac{\partial}{\partial x} \frac{\partial u}{\partial x} +
  \frac{\partial}{\partial x} \frac{\partial v}{\partial y} +
  \frac{\partial^2 u}{\partial y^2} +
  \frac{\partial}{\partial x} \frac{\partial w}{\partial z} +
  \frac{\partial^2 u}{\partial z^2}
\right] \\
&= -\frac{\partial P}{\partial x} + \rho g_x + \mu \left[
  \frac{\partial}{\partial x} \left(
    \frac{\partial u}{\partial x} +
    \frac{\partial v}{\partial y} +
    \frac{\partial w}{\partial z}
  \right) +
  \frac{\partial^2 u}{\partial x^2} +
  \frac{\partial^2 u}{\partial y^2} +
  \frac{\partial^2 u}{\partial z^2}
\right]
\end{align}

El término entre paréntesis es cero debido a la ecuación de continuidad para flujo incompresible.

\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} + \frac{\partial w}{\partial z} = 0

También se reconocen los últimos tres términos como el Laplaciano de la componente de velocidad u en coordenadas cartesianas. Por lo tanto, la componente x de la ecuación de cantidad de movimiento se escribe como:

(11) \rho \frac{Du}{Dt} = -\frac{\partial P}{\partial x} + \rho g_x + \mu \nabla^2 u

De manera similar se escriben las componentes y y z de la ecuación de cantidad de movimiento como:

(12) \rho \frac{Dv}{Dt} = -\frac{\partial P}{\partial y} + \rho g_y + \mu \nabla^2 v

Y

(13) \rho \frac{Dw}{Dt} = -\frac{\partial P}{\partial z} + \rho g_z + \mu \nabla^2 w

respectivamente. Por último, se combinan las tres componentes en una ecuación vectorial; el resultado es la ecuación de Navier-Stokes para flujo incompresible con viscosidad constante.

Ecuación de Navier-Stokes:

(14) \rho \frac{D \vec V}{Dt} = -\vec \nabla P + \rho \vec g + \mu \nabla^2 \vec V

Aunque los componentes de la ecuación Navier-stokes se dedujeron en coordenadas cartesianas, la forma vectorial de la ecuación es válida en cualquier sistema coordenado ortogonal. Esta famosa ecuación recibe su nombre en honor al ingeniero francés Louis Marie Henri Navier (1785-1836) y al matemático inglés Sir George Gabriel Stokes (1819-1903), quienes desarrollaron los términos viscosos, aunque de manera independiente. La ecuación de Navier-Stokes es la base de la mecánica de fluidos. Puede parecer suficientemente inocua, pero es una ecuación diferencial parcial de segundo orden, no-lineal e inestable. Si fuera posible resolver esta ecuación para flujos de cualquier geometría, seria más sencillo. Por desgracia, las soluciones analíticas no se obtienen excepto para campos de flujo muy simples. La ecuación tiene cuatro incógnitas (tres componentes de velocidad y la presión), aunque sólo representa tres ecuaciones (tres componentes puesto que es una ecuación vectorial). Obvio, es necesaria otra ecuación para solucionar el problema. La cuarta ecuación es la ecuación de continuidad para flujo incompresible:

\vec \nabla \cdot \vec V = 0

Antes de intentar resolver ese conjunto de ecuaciones diferenciales, es necesario elegir un sistema coordenado y expandir las ecuaciones en dicho sistema coordenado.


Ecuaciones de continuidad y de Navier-Stokes en coordenadas cartesianas

La ecuación de continuidad y la ecuación de Navier-Stokes se expanden en coordenadas cartesianas (x, y, z) y (u, v, w): Ecuación de continuidad de flujo incompresible:

\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} + \frac{\partial w}{\partial z} = 0

Componente x de la ecuación de Navier-Stokes de flujo incompresible:


\rho \left( \frac{\partial u}{\partial t} + u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y} + w \frac{\partial u}{\partial z} \right)
 = -\frac{\partial P}{\partial x} + \rho g_x + \mu \left( \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2} \right)

Componente y de la ecuación de Navier-Stokes de flujo incompresible:


\rho \left( \frac{\partial v}{\partial t} + u \frac{\partial v}{\partial x} + v \frac{\partial v}{\partial y} + w \frac{\partial v}{\partial z} \right)
 = -\frac{\partial P}{\partial y} + \rho g_y + \mu \left( \frac{\partial^2 v}{\partial x^2} + \frac{\partial^2 v}{\partial y^2} + \frac{\partial^2 v}{\partial z^2} \right)

Componente Z de la ecuación de Navier-Stokes de flujo incompresible:


\rho \left( \frac{\partial w}{\partial t} + u \frac{\partial w}{\partial x} + v \frac{\partial w}{\partial y} + w \frac{\partial w}{\partial z} \right)
 = -\frac{\partial P}{\partial z} + \rho g_z + \mu \left( \frac{\partial^2 w}{\partial x^2} + \frac{\partial^2 w}{\partial y^2} + \frac{\partial^2 w}{\partial z^2} \right)

Ecuaciones de continuidad y de Navier-Stokes en coordenadas cilíndricas

La ecuación de continuidad y la ecuación de Navier-Stokes se expanden en coordenadas cilíndricas (r, θ, z) y (ur, uθ, uz):

Ecuación de continuidad de flujo incompresible:

\frac{1}{r} \frac{\partial (ru_r)}{\partial r} + \frac{1}{r} \frac{\partial (u_\theta)}{\partial \theta} + \frac{\partial (u_z)}{\partial z} = 0

Componente r de la ecuación de Navier-Stokes de flujo incompresible:


\rho \left(
   \frac{\partial u_r}{\partial t}
 + u_r \frac{\partial u_r}{\partial r}
 + \frac{u_\theta}{r} \frac{\partial u_r}{\partial \theta}
 + \frac{u^2_\theta}{r}
 + u_z \frac{\partial u_r}{\partial z}
 \right)
 = -\frac{\partial P}{\partial r}
 + \rho g_r
 + \mu \left[
   \frac{1}{r} \frac{\partial}{\partial r} \left( r \frac{\partial u_r}{\partial r} \right)
 + \frac{u_r}{r^2}
 + \frac{1}{r^2} \frac{\partial^2 u_r}{\partial \theta^2}
 + \frac{2}{r^2} \frac{\partial u_r}{\partial \theta}
 + \frac{\partial^2 u_r}{\partial z^2}
 \right]

Componente θ de la ecuación de Navier-Stokes de flujo incompresible:


\rho \left(
   \frac{\partial u_\theta}{\partial t}
 + u_r \frac{\partial u_\theta}{\partial r}
 + \frac{u_\theta}{r} \frac{\partial u_\theta}{\partial \theta}
 + \frac{u_r u_\theta}{r}
 + u_z \frac{\partial u_\theta}{\partial z}
 \right)
 = -\frac{1}{r} \frac{\partial P}{\partial \theta}
 + \rho g_\theta
 + \mu \left[
   \frac{1}{r} \frac{\partial}{\partial r} \left( r \frac{\partial u_\theta}{\partial r} \right)
 + \frac{u_\theta}{r^2}
 + \frac{1}{r^2} \frac{\partial^2 u_\theta}{\partial \theta^2}
 + \frac{2}{r^2} \frac{\partial u_\theta}{\partial \theta}
 + \frac{\partial^2 u_\theta}{\partial z^2}
 \right]

Componente Z de la ecuación de Navier-Stokes de flujo incompresible:


\rho \left(
   \frac{\partial u_z}{\partial t}
 + u_r \frac{\partial u_z}{\partial r}
 + \frac{u_\theta}{r} \frac{\partial u_z}{\partial \theta}
 + u_z \frac{\partial u_z}{\partial z}
 \right)
 = -\frac{\partial P}{\partial z}
 + \rho g_z
 + \mu \left[
   \frac{1}{r} \frac{\partial}{\partial r} \left( r \frac{\partial u_z}{\partial r} \right)
 + \frac{1}{r^2} \frac{\partial^2 u_z}{\partial \theta^2}
 + \frac{\partial^2 u_z}{\partial z^2}
 \right]

Los términos adicionales en ambos lados de las componentes r y θ de la ecuación de Navier-Stokes surgen debido a la naturaleza especial de las coordenadas cilíndricas. De esta manera, conforme se mueve en la dirección θ, el vector unitario er , también cambia de dirección; por lo tanto, las componentes r y θ se acoplan.

A continuación, citaremos las seis componentes independientes del tensor de esfuerzo viscoso en coordenadas cilíndricas:


\tau_{ij} =
\begin{pmatrix}
  \tau_{rr} & \tau_{r\theta} & \tau_{rz} \\
  \tau_{\theta r} & \tau_{\theta\theta} & \tau_{\theta z} \\
  \tau_{zr} & \tau_{z\theta} & \tau_{zz}
\end{pmatrix} =
\begin{pmatrix}
  2 \mu \frac{\partial u_r}{\partial r} &
 \mu \left[ r \frac{\partial}{\partial r} \left( \frac{u_\theta}{r} \right) \frac{\partial u_r}{\partial \theta} \right] &
 \mu \left( \frac{\partial u_r}{\partial z} + \frac{\partial u_r}{\partial r} \right) \\
 \mu \left[ r \frac{\partial}{\partial r} \left( \frac{u}{r} \right) + \frac{1}{r} \frac{\partial u_r}{\partial \theta} \right] &
 2 \mu \left( \frac{1}{r} \frac{\partial u_r}{\partial \theta} + \frac{\partial u_r}{\partial r} \right) &
 \mu \left( \frac{\partial u_\theta}{\partial z} + \frac{1}{r} \frac{\partial u_z}{\partial \theta} \right) \\
 \mu \left( \frac{\partial u_r}{\partial z} + \frac{\partial u_r}{\partial r} \right) &
 \mu \left( \frac{\partial u_\theta}{\partial z} + \frac{1}{r} \frac{\partial u_z}{\partial \theta} \right) &
  2 \mu \frac{\partial u_z}{\partial z}
\end{pmatrix}

La aplicación de las ecuaciones diferenciales de movimiento tanto en coordenadas cartesianas como en cilíndricas. Existen dos tipos de problemas para los que son útiles las ecuaciones diferenciales (de continuidad y de Navier-Stokes):

  • Cálculo de campo de presión para un campo de velocidad conocido.
  • Cálculo de campos de velocidad y presión para un flujo de geometría conocida y condiciones de frontera conocidas.

Por simplicidad, sólo se considera flujo incompresible, cuando se eliminan el cálculo de ρ como una variable. Además, la forma de la ecuación de Navier-Stokes sólo es válida para fluidos newtonianos con propiedades constantes (viscosidad, conductividad térmica, entre otras). Para finalizar, se suponen variaciones de temperatura despreciables, de modo que T no es una variable. Quedan cuatro variables o incógnitas (presión más tres componentes de velocidad) y se tienen cuatro ecuaciones diferenciales.

Bibliografía

  • Yunus A. Cengel, John M. Cimbala, "MECANICA DE FLUIDOS- FUNDAMENTOS Y APLICACIONES". MCGRAWHILL, MARZO 2006.
  • P.GERHART, R. GROSS, J. HOCHSTEIN,"FUNDAMENTOS DE MECANICA DE FLUIDOS",2DA.EDICION, ED. ADDISON-WESLEY IBEROAMERICANA. EUA.1992
Obtenido de "Flujo incompresible"

Wikimedia foundation. 2010.

Mira otros diccionarios:

  • Flujo sanguíneo — Saltar a navegación, búsqueda Contenido 1 Definición. 2 Introducción 2.1 Valores normales en el humano 2.2 Índice c …   Wikipedia Español

  • Flujo en tubería — Saltar a navegación, búsqueda Uno de los aspectos de la dinámica de fluidos es el comportamiento de los flujos de fluidos, es decir, el movimiento de estos últimos. Contenido 1 Algunas características generales del flujo de fluidos 2 La ecuación… …   Wikipedia Español

  • Flujo compresible — Saltar a navegación, búsqueda Ondas de choque. Todos los fluidos son compresibles, incluyendo los líquidos. Cuando estos cambios de volumen son demasiado grandes se opta por considerar el flujo como compresible (que muestran una variación… …   Wikipedia Español

  • Fluido incompresible — Saltar a navegación, búsqueda Un fluido incompresible es cualquier fluido cuya densidad siempre permanece constante con el tiempo, y tiene la capacidad de oponerse a la compresión del mismo bajo cualquier condición. De hecho, todos los fluidos… …   Wikipedia Español

  • Coeficiente de presión — El coeficiente de presión es un número adimensional que describe la presión relativa a través de un campo de flujo en dinámica de fluidos. El coeficiente de presión es usado en aerodinámica e hidrodinámica. Cualquier punto inmerso en el flujo de… …   Wikipedia Español

  • Principio de Bernoulli — Para el teorema matemático enunciado por Jakob Bernoulli, véase Teorema de Bernoulli. Esquema del Principio de Bernoulli. El principio de Bernoulli, también denominado ecuación de Bernoulli o Trinomio de Bernoulli, describe el comportamiento de… …   Wikipedia Español

  • Mecánica de fluidos — La mecánica de fluidos es la rama de la mecánica de medios continuos (que a su vez es una rama de la física) que estudia el movimiento de los fluidos (gases y líquidos) así como las fuerzas que los provocan.[1] La característica fundamental que… …   Wikipedia Español

  • Mecánica — La Mecánica comprende el estudio de las máquinas (Polea simple fija). Para otros usos de este término, véase Mecánica (desambiguación). La mecánica (Griego Μηχανική y de latín mechanìca o arte de construir una máquina) es la rama de la física que …   Wikipedia Español

  • Helicidad (mecánica de fluidos) — Para otros usos de este término, véase helicidad. En mecánica de fluidos, la helicidad es una magnitud física asociada al flujo que mide en cada punto del mismo la proyección de la vorticidad sobre la velocidad, dado que la combinación de un giro …   Wikipedia Español

  • Campo solenoidal — Saltar a navegación, búsqueda Un campo solenoidal es un campo vectorial v cuya divergencia es cero: Esta condición se satisface si v es derivable de un potencial vector, por ejemplo A, ya que: Ya que entonces se cumple automáticamente que …   Wikipedia Español


We are using cookies for the best presentation of our site. Continuing to use this site, you agree with this.