INSTITUTO POLITECNICO NACIONAL
ESCUELA SUPERIOR DE INGENIERIA QUIMICA E INDUSTRIAS
EXTRACTIVAS
FLUJO DE FLUIDOS
Trabajo de investigación: Ecuaciones de Navier Stokes
Grupo:2IV41
Fecha de entrega: 09/11/2020
Introducción
Las ecuaciones de Navier-Stokes son ecuaciones que describen el movimiento
tridimensional del fluido viscoso. Las ecuaciones de Navier-Stokes proporcionan
una muy buena aproximación de la realidad.
En dinámica de fluidos , las ecuaciones de Navier-Stokes son ecuaciones que
describen el movimiento tridimensional de sustancias fluidas viscosas . Estas
ecuaciones llevan el nombre de Claude-Louis Navier (1785-1836) y George
Gabriel Stokes (1819-1903). En situaciones en las que no hay fuertes gradientes
de temperatura en el fluido, estas ecuaciones proporcionan una muy buena
aproximación de la realidad.
Las ecuaciones de Navier-
Stokes consisten en
una ecuación de continuidad dependiente del
tiempo para la conservación de la masa ,
tres ecuaciones de conservación del momento dependientes del tiempo y
una ecuación de conservación de la energía dependiente del tiempo . Hay cuatro
variables independientes en el problema, las coordenadas espaciales x, y y z de
algún dominio, y el tiempo t.
Las ecuaciones de Navier-Stokes también son de gran interés en un sentido
puramente matemático. Desafortunadamente, el carácter altamente intermitente e
irregular del flujo turbulento complica todos los análisis . Todavía no se ha
demostrado que en tres dimensiones siempre existan soluciones, o que, si existen,
entonces son fluidas. De hecho, a menudo se dice que la solución general de las
ecuaciones de Navier-Stokes con turbulencias es considerado último problema no
resuelto de la físico matemática clásica
Origen de las ecuaciones
Dado que la ecuación de Cauchy no es muy útil como se encuentra por que el
tensor esfuerzo contiene nueve componentes de los cuales seis son
independientes. Debido a esto, si contamos la densidad del fluido y las tres
componentes de la velocidad juntamos diez incógnitas y como sabemos para
poder resolver el sistema el número de incógnitas debe ser el mismo que el
número de ecuaciones. Es por eso que ecuaciones como la de continuidad que
solo tiene una y como la de Cauchy que solo tiene tres no nos permiten solucionar
el sistema.
Para obtener las ecuaciones faltantes, denominadas ecuaciones constitutivas, se
debe realizar un análisis de esfuerzos de presión y esfuerzos viscosos.
Dado que el único esfuerzo que actúa sobre un fluido en reposo es de la presión
hidrostática, tensor esfuerzo se expresa como:
( )( )
σxx σxy σxz −P 0 0
σij= σyx σyy σyz = 0 −P 0
σzx σzy σzz 0 0 −P
Donde la presión hidrostática es la presión
termodinámica obtenida mediante una ecuación de
estado, la cual esta relacionada con la temperatura y
la densidad del fluido. Cabe resaltar que podemos
denostar es P como presión hidrostática cuando
hablamos de fluido compresible pero cuando
hablamos de fluidos incompresibles debemos
cambiar esa P como presión mecánica.
Por otro lado, cuando el fluido está en movimiento, se ejerce sobre el un esfuerzo
viscoso o desviatorio Este esfuerzo nos permite establecer una relación con la
velocidad y la viscosidad. De tal forma que :
( )( )( )
σxx σxy σxz −P 0 0 τxx τxy τxz
σij= σyx σyy σyz = 0 −P 0 + τyx τyy τyz
σzx σzy σzz 0 0 −P τzx τzy τzz
Deducción de la ecuación de Navier Stokes para flujo isotérmico incompresible
El estudio de flujo con sestas ecuaciones se limita a fluido newtonianos, donde
podemos deducir que el tensor de esfuerzo ( τ ) es directamente proporcional al
tensor de deformación(ε). También cuando hablamos de un flujo isotérmico
incompresible, deducimos que ρ=cte., T=cte., µ=cte. y por lo tanto v=cte. Bajo
estas suposiciones podemos demostrar que:
τij=2 μεij
Y por lo tanto
( )
2μ ( ∂∂ ux )
( ∂ y ∂x)
∂u ∂v
+ μ μ( +
∂z ∂ x )
∂u ∂w
( ) ∂∂wx ∂∂yu ∂ w∂ y∂ v ∂ z ∂ w∂ y )
τxx τxy τxz
τyz = μ ( + ) 2 μ( ) μ( +
∂ v ∂u ∂v ∂v ∂w
τij= τyx τyy
τzx τzy τzz
μ(
∂x ∂ z )
μ(
∂ y ∂z )
2 μ(
∂z )
+ +
Y por lo tanto en el tensor de esfuerzo
( )
2μ ( ∂u∂x )
( + ) μ( +
∂u ∂ v
∂ y ∂x
μ
∂z ∂x)
∂u ∂ w
( ) ∂∂wx ∂∂uy ∂ w∂ y∂ v ∂ z ∂ w∂ y )
−P 0 0
0 + μ( + ) 2 μ( ) μ( +
∂v ∂u ∂v ∂v ∂w
σij= 0 −P
0 0 −P
μ(
∂x ∂ z ) ( ∂ y ∂z )
2 μ(
∂z )
+ μ +
Y ahora esta expresión se sustituye en la ecuación de Cauchy. Considerando
primero la componente x
( )
2
Du −∂ P ∂ y ∂ ∂ v ∂u ∂ ∂w ∂u
ρ = + ρgx+2 μ 2 + μ + +μ ( + )
Dt ∂x ∂x ∂x ∂x ∂ y ∂ z ∂x ∂ z
Nótese que dado que la presión es un esfuerzo normal solo aparece en un térmico
de la ecuación mientras que el esfuerzo viscoso tiene esfuerzos normales y de
corte aporta tres términos a la ecuación
Después reordenamos la ecuación, de tal forma que :
ρ
Du −∂ P
Dt
=
∂x
+ ρgx+ μ
[ ( + + + +)
∂ ∂ u ∂ v ∂ w ∂ 2 u ∂ 2 u ∂2 u
+
∂ x ∂ x ∂ y ∂ z ∂ x 2 ∂ y2 ∂ z 2 ]
Y como el término entre paréntesis es cero debido a la ecuación de continuidad y
los últimos tres dígitos de conocen como el laplaciano de la velocidad u en
coordenadas cartesianas podemos decir que:
Du −∂ P 2
ρ = + ρgx+ μ ∇ u
Dt ∂x
Y para Y y Z
Dv −∂ P 2
ρ = + ρgy + μ ∇ v
Dt ∂y
Dw −∂ P 2
ρ = + ρgz+ μ ∇ w
Dt ∂z
Y al combinar las tres ecuaciones en una sola obtenemos la ecuación de Navier
Stokes para un flujo incompresibles con viscosidad constante
D⃗
V
=− ⃗
∇ P+ ρ ⃗g + μ ∇ ⃗
2
ρ V
Dt
Coordenadas de la ecuación de Navier Stokes
Podemos expresar tanto la ecuación de continuidad en coordenadas cartesianas,
coordenadas cilíndricas y coordenadas esfericas según sea la necesidad de
nuestro problema.
A continuación:
Ecuación de continuidad y de Navier Stokes en coordenadas cartesianas
∂u ∂ v ∂w
+ + =0
∂x ∂ y ∂z
ρ ( ∂u
∂t
+u
∂u
∂x
+v
∂u
∂y
+w
∂z )
∂u −∂ P
=
∂x
+ ρgx+ μ (
∂ 2 u ∂2 u ∂2 u
+ +
∂ x 2 ∂ y 2 ∂ z2 )
( ) ( )
2 2 2
∂v ∂v ∂v ∂ v −∂ P ∂u ∂ u ∂ u
ρ +u + v +w = + ρgy+ μ 2
+ 2+ 2
∂t ∂x ∂y ∂z ∂y ∂x ∂ y ∂ z
ρ ( ∂w
∂t
+u
∂w
∂x
+v
∂w
∂y
+w
∂ w −∂ P
∂z
=
∂z )+ ρgz + μ 2(
∂2 u ∂ 2 u ∂ 2 u
+ 2+ 2
∂ x ∂ y ∂z )
Ecuación de continuidad y de Navier Stokes en coordenadas cilíndricas
1 ∂ ( rur ) 1 ∂ ( uθ ) 1 ∂ ( uz )
+ + =0
r ∂r r ∂θ r ∂ z
( ) ( ( ) )
2 2 2
∂ ur ∂ ur uθ ∂ ur u θ ∂ ur −∂ P 1 ∂ ∂ ur ur 1 ∂ u 2 ∂u ∂ ur
ρ + ur + − +uz = + ρgθ + μ r − + − +
∂t ∂ r r ∂θ r ∂z ∂r r ∂r ∂r r r ∂ θ2 r 2 ∂θ ∂ z 2
ρ ( ∂u
∂t
+ur
∂ uθ uθ ∂u uruθ
+
∂ r r ∂θ
−
r
+uz
∂ uθ −∂ P
∂z
=
∂θ )
+ ρgθ + μ
1 ∂
r
r ∂r ∂r ( ( )
∂u uθ 1 ∂2 u 2 ∂ur ∂2 u
− 2 + 2 2− 2 +
r r ∂ θ r ∂ θ ∂ z2 )
( ) ( ( ) )
2 2
∂ uz ∂uz uθ ∂uz ∂ uz −∂ P 1 ∂ ∂ uz 1 ∂ uz ∂ uz
ρ +ur + +uz = + ρgz + μ r + 2 +
∂t ∂r r ∂ θ ∂z ∂z r ∂r ∂r r ∂θ 2 ∂ z 2
Conclusión
Como se puede ver, las ecuaciones de Navier-
Stokes son ecuaciones diferenciales parciales no
lineales de segundo orden, sus soluciones se han
encontrado para una variedad de problemas
interesantes de flujo viscoso. Se pueden usar
para modelar el clima, las corrientes oceánicas, el flujo de aire alrededor de un
perfil y el flujo de agua en una tubería o en un reactor. Las ecuaciones de Navier-
Stokes en sus formas completas y simplificadas ayudan con el diseño de aviones
y automóviles, el estudio del flujo sanguíneo, el diseño de reactores nucleares y
muchas otras cosas.
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.
Aunque las ecuaciones de Navier-Stokes tienen solo un número limitado de
soluciones analíticas conocidas, son susceptibles de modelado computarizado de
cuadrícula fina. La herramienta principal disponible para su análisis es el análisis
CFD . CFD es una rama de la mecánica de fluidos que utiliza análisis numéricos y
algoritmos para resolver y analizar problemas que involucran flujos de fluidos
turbulentos. Es ampliamente aceptado que las ecuaciones de Navier-
Stokes (o las ecuaciones de Navier-Stokes simplificadas promediadas por
Reynolds ) son capaces de exhibir soluciones turbulentas, y estas ecuaciones son
la base de esencialmente todos los códigos CFD. Ahora es posible lograr
resultados CFD aproximados, pero realistas, para una amplia variedad de flujos
viscosos complejos de dos y tres dimensiones.
Referencias
(Yunus Cengel, 2006, págs. 426-430)Mecánica de Fluidos, Fundamentos y aplicaciones