Instituto Tecnológico de Saltillo

Departamento Metal-Mecánica
Proyecto final de unidad
“Aplicaciones de Navier Stokes”
DR. LUCIANO ELIEZER RAMÍREZ VIDAURRI
Integrantes:
Kevin Josué Cárdenas Hernández
Yessica Montserrat Banda Flores
Mirna Judith Esquivel Álvarez
Melisa Moncerrat Rodríguez Vargas
OBJETIVO

El objetivo de este trabajo es desarrollar y resolver el problema , así como aplicar las
herramientas, fórmulas y aprendizajes vistos y adquiridos en la clase, emplear las
Matemáticas. Emplear los métodos y procedimientos estadísticos y probabilísticos para
obtener conclusiones a partir del desarrollo del tema de NAVIER STOCKE.

INTRODUCCIÓN

Uno de los campos de la física más complicados de estudiar son los fluidos,
el comportamiento de gases y líquidos en movimiento. Comprender, por ejemplo, los
flujos de aire turbulento o los remolinos que se forman cuando el agua discurre por una
tubería o la sangre por una arteria son de suma importancia, tanto para la ingeniería
como para la medicina.

Las ecuaciones que rigen la dinámica de fluidos son las conocidas como ecuaciones de
Navier-Stokes, producto del francés constructor de puentes Claude-Louis Navier y del
matemático irlandés George Stokes. El primero en obtener estas ecuaciones fue el
francés en una época (1822)

Las ecuaciones de Navier-Stokes describen el movimiento de los fluidos “reales“ como

la velocidad, presión, temperatura y transporte del medio de flujo, mas mo los detalles
de las fluctuaciones de la tubería.
MARCO TEÓRICO

La resolución de las Ecuaciones de Navier-Stokes es uno de los problemas del milenio

recompensando por 1 millón de euros. Estas ecuaciones no son más que la aplicación
de los principios de conversación de masa, movimiento y energía a un fluido. Son tan
importantes porque su resolución implica un conocimiento absoluto de cómo un fluido se
va a mover partiendo unas condiciones iniciales y un gran uso que se le podría dar es el
de por ejemplo dar las previsiones meteorológicas exactas.

Sin embargo, el movimiento de fluido es tan caótico que las ecuaciones adquieren una
gran complejidad en su resolución y solo pueden resolverse numéricamente o
analíticamente por medio de muchas simplificaciones e idealizaciones del fluido.

DESCRIPCION LAGRANGIANA Y EULERIANA

En primer lugar, es importante saber que hay dos modos de describir el movimiento de
un fluido. Mediante la descripción Lagrangiana, las partículas del fluido son seguidas una
a una y de manera independiente a medida que estas se mueven. Como si fueran
canicas muy pequeñas.

Por otra parte, bajo la descripción Euleriana, se fija un volumen de control (fluido
contenido en un tramo de tubería, por ejemplo) y se analiza en conjunto las
características del fluido que fluye en su interior.

Es necesario combinar ambos puntos de vista para poder aplicar las leyes de newton al
fluido en el que debemos de considerar el movimiento independiente de las partículas
(lagrangiana) mientras que la observación y medición del comportamiento macroscópico
del fluido se realiza tomando volúmenes de control (euleriana).

DESCRIPCION DE LAGRANGIANA

Seguimos una partícula de fluido en su movimiento. Partiendo de unas condiciones

iniciales {𝑟⃗0, 𝑡0} y las fuerzas externas, empleando las ecuaciones de la dinámica
podríamos obtener la evolución de posición de la partícula con una formulación tal
que 𝑟⃗ = 𝑟⃗(𝑡; 𝑟⃗0, 𝑡0).

Una vez conocida la posición de cada una de las partículas, podemos definir el valor de
la magnitud en un punto de fluido como 𝐹 = 𝐹(𝑟⃗(𝑡; 𝑟⃗0, 𝑡0), 𝑡) que será la magnitud
asociada a la partícula estudiada.

DESCRIPCION EULERIANA

Ahora nos centramos en describir las propiedades de un fluido alrededor de una zona
centrada en 𝑥⃗ en un tiempo t. De esta forma, el valor de una magnitud expresada
según esta descripción como campo viene dada por 𝐹 = 𝐹(𝑥⃗, 𝑡)

UNION DE AMBAS DESCRIPCION

Para que ambas descripciones o modelos sean válidos, debe de haber una
equivalencia en los valores 𝐹 𝑐𝑢𝑎𝑛𝑑𝑜 𝑥⃗ = 𝑟⃗ de en el mismo instante de tiempo t, es
decir, se estudie el mismo punto del fluido:

𝐹(𝑟⃗(𝑡; 𝑟⃗0, 𝑡0), 𝑡) = 𝐹(𝑥⃗, 𝑡)𝑐𝑢𝑎𝑛𝑑𝑜 𝑥⃗ = 𝑟⃗(𝑡; 𝑟⃗0, 𝑡0)

Si aplicamos la derivada de dicha función con respecto al tiempo:

𝑑 𝜕𝐹 𝜕𝐹 𝜕𝑟1 𝜕𝐹 𝜕𝑟2 𝜕𝐹 𝜕𝑟3 𝑑

[𝐹(𝑟⃗(𝑡; 𝑟⃗0, 𝑡0), 𝑡)] = + ° + ° + ° = [𝐹(𝑥⃗, 𝑡)]
𝑑𝑡 𝜕𝑡 𝜕𝑟1 𝜕𝑡 𝜕𝑟2 𝜕𝑡 𝜕𝑟3 𝜕𝑡 𝑑𝑡
𝑑𝑟1
Teniendo en cuenta que 𝑥⃗ = 𝑟⃗ y sabiendo que = 𝑢1 (la derivada temporal de la
𝑑𝑡
posición en una coordenada es igual a la velocidad instantánea en dicha coordenada),
definimos la siguiente relación:

𝐷 𝜕𝐹 𝜕𝐹 𝜕𝐹 𝜕𝐹 𝜕𝐹
𝐹(𝑥⃗, 𝑡) = + °𝑢1 + °𝑢2 + °𝑢3 = +𝑢 ⃗⃗ 𝐹
⃗⃗°∇
𝐷𝑡 𝜕𝑡 𝜕𝑟1 𝜕𝑟2 𝜕𝑟3 𝜕𝑡
𝐷
Siendo 𝐷𝑡 definida como la derivada temporal material de la descripción Euleriana de un
𝑑
campo fluido, equivalente a la derivada temporal total 𝑑𝑡
de la descripción Lagrangiana
que expresa la variación temporal total de una magnitud que se desplaza en el tiempo.
 𝑑
Es decir, que la derivada temporal total se emplea cuando hablamos de partículas, y
𝑑𝑡
𝐷
la derivada temporal material 𝐷𝑡 cuando hablamos de volúmenes.

ECUACIONES DE NAVIER-STOKES

Las ecuaciones de conservación de masa y cantidad de movimiento definidas hasta

ahora sin hacer ninguna simplificación del fluido son:

MASA

𝜕𝜌 𝜕
+ (𝜌 ∙ 𝑢𝑖) = 0
𝜕𝑡 𝜕𝑥𝑖

CANTIDAD DE MOVIMIENTO

𝜕𝑢𝑖 𝜕𝑢𝑖 𝜕𝑃 𝜕 𝜕𝑢𝑖 𝜕𝑢𝑗 𝜕 𝜕𝑢𝑗

𝜌∙ + 𝜌 ∙ 𝑢𝑗 ∙ =− +𝜇 ( + ) +⋋ ( ) + 𝜌 ∙ 𝑔𝑖
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑖 𝜕𝑥𝑗

Si consideramos un flujo a densidad constante o incompresible, entonces la ecuación de

𝜕𝑢𝑖
conversación de masa se reduce a 𝜕𝑥𝑖 = 0 de forma que de la ecuación de conservación
de cantidad de movimiento desaparecen los siguientes términos:

𝜕 𝜕𝑢𝑗 𝜕 𝜕𝑢𝑗
∎𝜇 ( )=𝜇 ( )=0
𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑖 𝜕𝑥𝑗

𝜕 𝜕𝑢𝑗
∎⋋ ( )=0
𝜕𝑥𝑗 𝜕𝑥𝑗

De esta forma, para un flujo a densidad constante las ecuaciones de Navier-Stokes se

reducen a:

𝜕𝑢𝑖
=0
𝜕𝑥𝑖

𝜕𝑢𝑖 𝜕𝑢𝑖 𝜕𝑃 𝜕 2 𝑢𝑖
𝜌∙ + 𝜌 ∙ 𝑢𝑗 ∙ =− + 𝜇 2 + 𝜌 ∙ 𝑔𝑖
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥 𝑗
𝜇
Que, empleando vectores y definiendo la viscosidad cinemática como v: 𝜌:

𝜕𝑢 𝜕𝓋 𝜕𝑤
+ + =0
𝜕𝑥 𝜕𝒴 𝜕𝑧
 𝜕𝑢 𝜕𝑢 𝜕𝑢 𝜕𝑢 1 𝜕𝑃 𝜕 2𝑢 𝜕 2𝑢 𝜕 2
+𝑢∙ +𝑣∙ +𝓌∙ =− + 𝑣 ( 2 + 2 + 2 ) + 𝑔𝑥
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜌 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑧

𝜕𝑣 𝜕𝑣 𝜕𝑣 𝜕𝑣 1 𝜕𝑃 𝜕 2𝑣 𝜕 2𝑣 𝜕 2𝑣
+𝑢∙ +𝑣∙ +𝑤∙ =− + 𝑣 ( 2 + 2 + 2 ) + 𝑔𝓎
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜌 𝜕𝑦 𝜕𝑥 𝜕𝑦 𝜕𝑧

𝜕𝑤 𝜕𝑤 𝜕𝑤 𝜕𝑤 1 𝜕𝑃 𝜕 2𝑤 𝜕 2𝑤 𝜕 2𝑤
+𝑢∙ +𝑣∙ +𝑤∙ =− +𝑣( 2 + + ) + 𝑔𝑧
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜌 𝜕𝑧 𝜕𝑥 𝜕𝒴 2 𝜕𝑧 2

Si nos fijamos en los términos del sistema de ecuaciones diferenciales, como es lógico
son 4 ecuaciones (físicamente independientes entre si) para 4 incógnitas, siendo”””
posible””” matemáticamente obtener una solución para las 4 incógnitas:

𝑢(𝑥⃗, 𝑡), 𝑣(𝑥⃗, 𝑡), 𝑤(𝑥⃗, 𝑡), 𝑃(𝑥⃗, 𝑡)

Flujo de Momentum:

La ecuación de movimiento, expresa el principio de conservación de cantidad de

movimiento. Al sistema pude entrar cantidad de movimiento por transporte, de acuerdo
con la expresión newtoniana to no-newtonianal, de densidad de flujo de cantidad de
movimiento. También puede entrar cantidad de movimiento debido al movimiento global
del fluido, Las fuerzas que nos interesan son las fuerzas de presión (actuando sobre
superficie) y las fuerzas de gravedad (que actúan sobre todo el volumen).

Para un elemento diferencial de volumen se puede escribir el siguiente balance de

cantidad de movimiento, la ecuación es para sistemas en estado no estacionario:

(Velocidad de (Velocidad de
acumulación de entradas de
cantidad de cantidad de
movimiento) movimiento) (Suma de las fuerzas
(Velocidad de que actúan sobre el
salida de cantidad sistema)
del movimiento)

Enlizando la ecuación en la dirección x la cantidad de movimiento que entra y sale del

elemento de volumen indicado como se muestra en la figura anterior, se produce por dos
mecanismos como la convección y transporte molecular, es decir, debido al flujo global
del fluido y a causa de los gradientes de velocidad.

La velocidad con que entra debido a la cantidad de movimiento por convección en el eje
x de la cantidad de movimiento por la cara situada en x en las
coordenadas y la velocidad de salida sumándole las fuerzas que actuaran sobre el
sistema del cual se esté estudiando.

Para las demás caras se pueden escribir expresiones similares. Hay que tomar en cuenta
que las densidades de flujo para la cantidad de movimiento en el eje j a través de una
cara perpendicular al eje sumando las seis contribuciones se obtiene que:

Estas densidades de flujo de la cantidad de movimiento se consideran esfuerzos.

cortantes por lo tanto el esfuerzo normal actuara en una sola cara y será un esfuerzo
tangencial que actuara sobre la cara en dirección xy que resultara de las fuerzas
viscosas.

Ecuacion de Continuidad:

En el flujo la continuidad es consecuencia de la conservación de la materia, esto se

refiere a la constancia dei flujo a is largo del camino recorrido, su enunciado es. El flujo
de un fluido en movimiento es el mismo en dos puntos diferentes del camino recorrido
por el fluido.

Esta ecuación también recibe el nombre de ecuación de continuidad del fujo. Expresa
que la cantidad de masa por unidad de tiempo que ingresa por un purito debe ser igual
a la cantidad de masa por unidad de tiempo que sale por un punto recorrido del fluido

Si un fluido es un liquido no viscoso e incomprensibio su densidad permanece constante

durante el flujo, es por ello que se puede eliminar la densidad en ambos miembros de la
ecuación del fluja, per lo que la ecuación de continuidad del flujo se reduce a la ecuación
de continuidad del caudal del liquido.

La ecuación de continuidad se otitiene por la aplicación de un balance de materia a un

elemento que posee un volumen y se puede representar por medio de una diferencial a
través de la cual circula el fluido (figura 10).
La ecuación de continuidad nos permitirá calcular la velocidad de un fluido en un sistema
de conductos cerrados, esto se puede determinar por la velocidad de figo y masa como
se menciona en la formula general de continuidad:

(Velocidad de (Velocidad de entrada (Velocidad de salida de

acumulación de de materia) la materia)
materia)

La ecuacion de continuidad resultante es:

Para un fluido de densidad constante: (Fluido incomplensible)

Flujo Laminar de Fluidos newtonianos y no newtonianos:

Mediante balances de cantidad de movimiento se pueden resolver varios problemas de
flujo viscoso en estado estacionario. Los métodos más generales permiten describir
sistemas de flujo más complicados.

Es importante que primeramente revises y entiendas los distintos regimenes de flujo:

laminar, de transición y turbulento. El análisis adimensional es un método para verificar
ecuaciones y planificar experimentos sistemáticos, a partir del análisis dimensional se
obtiene una serie de grupos adimensionales que va a permitir obtener resultados
expenmentales en condiciones limitadas, a situaciones en que se tengan diferentes.
dimensiones geométricas, cinematicas y dinámicas, y muchas veces en casos en que
las propiedades del fluido y del flujo son distintas de las que se tuvieron durante los
experimentos

La importancia del analisis adimensional viene

dada por la dificultad del establecimiento de
ecuaciones en determinados fujos, ademas de la
dificultad de su resolución, siendo imposible
obtener retaciones empiricas.

El número de Reynolds representa la relación de

fuerzas inerciales y fuerzas viscosas. Permite
caracterizar la naturaleza del flujo, es decir, si se
trata de un flujo laminar o turbulento. Es un
número adimensional utilizado en mecánica de
fluidos, un fluido.

Flujo de estudio:

Un ejemplo de un flujo no laminar es un río suave: cada parte del río se mueve en la
misma dirección a la misma velocidad. Un fluido laminar es la fracturación de ese río, de
modo que diferentes partes del flujo se mueven en diferentes direcciones a diferentes
velocidades. Los físicos describen la formación de turbulencias como, primero, un
remolino en un flujo suave, y luego la formación de remolinos dentro de ese remolino, y
remolinos más finos dentro de esos remolinos: se desplazan hacia abajo, de modo que
el fluido se rompe en partes discretas, todos interactuando, cada uno moviéndose a su
manera.
Probelama 8

Un fluido viscoso e incompresible fluye entre las dos placas infinitas, verticales y
paralelas que se muestran en la figura. Determine, mediante el uso de las ecuaciones de
Navier-Stokes, una expresión para el gradiente de presión en la dirección del flujo.

Exprese su respuesta en términos de la velocidad media. Suponga que el flujo es laminar,

estable y completamente desarrollado.

1) Se tiene un sistema con coordenadas rectangulares.

2) Se tiene 2 placas fijas con fluido con una viscosidad con densidad rho.
3) Su velocidad en los extremos de las placas tendrá el comportamiento que se
muestra en la Fig.1.

Fig.1. Se tiene este perfil de velocidades ya que tiene paredes fijas y en ellas se
tiene una v=0.
4) Después se tiene la siguiente Fig.2, la cual muestra una distancia h entonces, a
la mitad hacia una de las placas tendría una distancia h/2, este punto lo
tomaremos como y=0

Y=0
Fig.2.
ℎ
2
Campo de velocidad o perfil de velocidad.
-Las placas están fijas, lo cual estas no generan movimiento. Entonces, quien es el agente de
el movimiento del fluido se genera debido a una diferencia de presiones.
-El fluido siempre a querido ir a un lugar donde hay menos presión entonces quedaría
representado como en la Fig.3.

P1

P2

Fig.3.
-El fluido tendrá esa dirección siempre y cuando P2 sea mayor que P1, o mejor dicho si el
gradiente de presiones en el eje x es menor a cero (0).
𝜕𝑃
𝑃2 > 𝑃1 Ó <0
𝜕𝑋
Es lo que genera movimiento de frio y a el gradiente de presiones.
ENCONTAR CAMPO DE VELOCIDADES
-Se tiene supuestos:
*Incompresible
*Newtoniano
*Estacionario
*Es un problema Bidimensional
𝜕𝑉𝑥
Vx=? 𝜕𝑦
=?

𝜕𝑉𝑥
Vy=0 𝜕𝑥
=?
𝜕𝑉𝑥
Vz=0 𝜕𝑧
=0

-ECUACIÓN DE CONTINUIDAD-
Variación de densidad con respecto al tiempo mas la divergencia de rho por campo de
velocidades esto es igual a cero(0).
 𝜕𝜌
+ 𝐷(𝑟 ∙ 𝑣⃗) = 0
𝜕𝑡
𝜕𝜌
Se elimina términos: es estacionario e incompresible no hay variación de densidad se elimina 𝜕𝑡
quedaría: divergencia de rho por 0 = 𝐷(𝑟 ∙ 𝑣⃗) al igual de que es incompresible y es una constante,
se saca rho de la divergencia y quedaría: la divergencia por 𝑣⃗ y es =0.
∴ 𝑟(𝐷 ∙ 𝑣⃗) = 0

Se elimina rho y queda: 𝐷 ∙ 𝑣⃗ = 0

Entonces esta nos dice que la variación con respecto a U con respecto a X más la derivada de
V con respecto a Y más la derivada de W con respecto a Z es igual a cero (0).
𝜕𝑢 𝜕𝑣 𝜕𝑤
+ + =0
𝜕𝑥 𝜕𝑦 𝜕𝑧
Se tiene movimiento del fluido en solo una dirección que es X, en la dirección Y no hay
movimiento del fluido, el movimiento del fluido en esta dirección no se mueve.

V=0
𝜕𝑢 𝜕𝑣 𝜕𝑤
Entonces= 𝜕𝑥
+ 𝜕𝑦 + 𝜕𝑧
=0

-Como es un problema bidimensional y el efecto del eje z es nulo por lo tanto no hay velocidad
en el eje z no hay variación, por lo tanto, se anula:
V=0 w=0

𝜕𝑢 𝜕𝑣 𝜕𝑤
+ + =0
𝜕𝑥 𝜕𝑦 𝜕𝑧
Por lo tanto, de la ecuación de continuidad nos queda que la derivada de U respecto a X es igual
a cero:
𝜕𝑢
=0
𝜕𝑥
Esto implica implica que U es una función de Y
𝜕𝑢
= 0 → 𝑢 = 𝑢(𝑦) Ecuación de continuidad.
𝜕𝑥
El perfil de velocidad va cambiando a medida que aumenta el eje Y al igual va cambiando su
forma.
-ECUACIÓN DE MOMENTUM-
𝜕 2 𝑉𝑥 𝜕𝑃
0=𝜇∙ − + 𝜌𝑔𝑥
𝜕𝑦 2 𝜕𝑥
𝜕 2 𝑉𝑥 1 𝜕𝑃
= ∙ − 𝜌𝑔𝑥
𝜕𝑦 2 𝜇 𝜕𝑥
*Se deriva para encontrar Vx respecto a Y
 𝜕 2 𝑉𝑥 1 𝜕𝑃
= ∙ − 𝜌𝑔𝑥
𝜕𝑦 2 𝜇 𝜕𝑥
𝜕𝑉𝑥 1 𝜕𝑃
= 𝜇 ∙ 𝜕𝑥 − 𝜌𝑔𝑥 (y) + C1.
𝜕𝑦

1 𝜕𝑃 𝑦2
𝑣𝑥 = ∙ − 𝜌𝑔𝑥 + 𝐶1(𝑌) + 𝐶2.
𝜇 𝜕𝑥 2
-CONDICIONES A LA FRONTERA-
ℎ
1era condición→y= 2 Vx=0

*La velocidad en la pared, el cual se tomo como v=0 como se muestra en la Fig.4., entonces Y
ℎ
es igual a 2, esta es llamada velocidad de no deslizamiento.

Vmáx.

Y=0
0
Fig.4.
-La condición de guardia que se debe utilizar es la misma condición a la frontera que se hace en
los flujos a través de tubería tiene el perfil de temperatura de la Fig.1.
*Donde en las paredes la velocidad es cero, pero en el centro en el y=0 se tiene la velocidad
máxima, entonces, decimos que:
La derivada de U con respecto a Y es igual a cero:
𝜕𝑢
𝜕𝑦
=0 𝑦=0 2da Condición a la frontera.

Nota: Estamos diciendo que la pendiente en este punto sería 0 y que al final la velocidad en esta
zona seria máxima.
o Sustituir:
−𝜌𝑔𝑥 𝜕𝑃
2 0= (0) + 𝐶1(0) + 𝐶2.
𝜇 𝜕𝑥
𝐶2 = 𝑂

1 2
1 −𝜌𝑔𝑥 𝜕𝑃 2ℎ 1
0= 𝜇 𝜕𝑥 2
( )+ 𝐶1(2 ℎ) + 𝐶2.

−𝜌𝑔𝑥 𝜕𝑃 ℎ 2 1
0= ( )+ 𝐶1( ℎ) + 𝐶2.
𝜇 𝜕𝑥 8 2

1 𝜌𝑔𝑥 𝜕𝑃 2 ℎ
(ℎ ) = 𝐶1( ) + 𝐶2.
8 𝜇 𝜕𝑥 2
 1 𝜌𝑔𝑥 𝜕𝑃 2
(ℎ ) = 𝐶1 + 𝐶2.
4ℎ 𝜇 𝜕𝑥
0
1 𝜌𝑔𝑥 𝜕𝑃 2
(ℎ ) + 𝐶2 = 𝐶1
4ℎ 𝜇 𝜕𝑥

PERFIL DE VELOCIDADES
Remplazamos C1 en U y tendremos nuestro perfil de velocidades
1 𝜌𝑔𝑥 𝜕𝑃 1 𝜌𝑔𝑥 𝜕𝑃
U(y)= 4𝑦 [ 𝜇 𝜕𝑥
(𝑦 2 )] − 4ℎ [ 𝜇 𝜕𝑥 (ℎ2 )]

Factorizamos:
1 𝜌𝑔𝑥 𝜕𝑃 1 1
U(y)= [ (𝑦 2 ) − (ℎ2 )] Perfil de velocidad
4 𝜇 𝜕𝑥 𝑦 ℎ

-ESFUERZO CORTANTE- en la pared.

Se trabaja en el siguiente plano.

Fig.2.
y

*2 placas fijas, con perfil de velocidad, tenemos

ℎ
el eje en el centro y=0 con distancia
2

*El esfuerzo cortante que vamos a calcularen la pared es (𝜏𝑥𝑦).

𝜏𝑥𝑦

ℎ
Estamos buscando el esfuerzo cortante en el punto 2 porque ahí es donde esta la pared, es cuál
es miu por la variación ósea la derivada de U con respecto a Y más la derivada de V con respecto
a X.
 ℎ 𝜕𝑢 𝜕𝑣
𝜏𝑥𝑦( ) = 𝜇( + )
2 𝜕𝑦 𝜕𝑥

~Pero la derivada de V no existía por lo tanto es igual a cero.

0
ℎ 𝜕𝑢 𝜕𝑣
𝜏𝑥𝑦( ) = 𝜇( + )
2 𝜕𝑦 𝜕𝑥

Por lo tanto, nuestro esfuerzo en la pared solo será:

ℎ 𝜕𝑢 ℎ
𝜏𝑥𝑦( ) = 𝜇 𝑒𝑛 𝑒𝑙 𝑝𝑢𝑛𝑡𝑜 ( ).
2 𝜕𝑦 2
𝜇 → 𝑉𝑖𝑠𝑐𝑜𝑠𝑖𝑑𝑎𝑑 𝑑𝑖𝑛𝑎𝑚𝑖𝑐𝑎.

-Ahora tocaría sacar la derivada de la variación de velocidad en Y.

𝜕𝑢 1 𝜕𝑃 1
= 𝜌𝑔𝑥 (− 2 ) (2𝑦)
𝜕𝑦 4 𝜕𝑥 𝑦
ℎ
𝜕𝑢 1 𝜕𝑃 −2 2
(𝑦 = ℎ) = ∙ 𝜌𝑔𝑥 ∙
𝜕𝑦 4𝜇 𝜕𝑥 (ℎ)2
2
𝜕𝑢 1 𝜕𝑃 1
(𝑦 = ℎ) = ∙ 𝜌𝑔𝑥 ∙ ∙ (−4 ∙ )
𝜕𝑦 4𝜇 𝜕𝑥 ℎ
1 𝜕𝑃 1
𝜏𝑥𝑦 = 𝜇 ∙ (−1 ∙ 𝜌𝑔𝑥 ∙ ∙ )
𝜇 𝜕𝑥 ℎ

𝜕𝑃 1
𝜏𝑥𝑦 = −1 ∙ 𝜌𝑔𝑥 ∙ ∙
𝜕𝑥 ℎ

•VELOCIDAD MEDIA•
-Se define como caudal partionario en un área transversal del fluido:
𝑄
𝑉̅ 𝑥 =
𝐴
-El caudal se define teniendo un campo de velocidad entonces queda definido como:
*La integral del campo de velocidades por un diferencial de área:
 Q=∫ 𝑢 ∙ 𝑑𝐴

∴ la velocidad media es igual a uno sobre el área por la integra del campo de velocidades por el
diferencial de área.
1
𝑉̅ 𝑥 = ∙ ∫ 𝑢 ∙ 𝑑𝐴
𝐴

*El área es aquella por el cual pasa un caudal (el área transversal).
-Lo primero que debemos hacer es darnos una coordenada en x, ósea dar el espesor de la placa,
lo que se quiere decir es que al tener una tubería de este estilo el cual pasa por aquí un fluido,
tiene una distancia de "'h" pero no se tiene distancia perpendicular la cual podemos dar un
espesor (espesor variable) "b” ya que cuando se tenga uno sea más fácil de sustituir y encontrar
un valor numérico.

b Figura 5.

Cuando nos demos el espesor, el área será h por el espesor por b.

Nota: como es rectangular el espesor se eliminará, lo contrario que en el circular que tiene que
ser dos figuras de área del perfil circular.
* Ahora integramos el perfil de velocidad y el diferencial de área.
1
𝑉̅ 𝑥 = ∫ 𝑢(𝑦) ∙ 𝑑𝐴
(ℎ ∙ 𝑏)

*Lo que pasa aquí es que el perfil de velocidad depende de Y que va a tener a integrar con
respecto a Y, por lo tanto, va a quedar nuestro diferencial como se muestra a continuación:
1
𝑉̅ 𝑥 = ∫ 𝑢(𝑦) ∙ (𝑏)𝑑𝑦
(ℎ ∙ 𝑏)

Espeso Variación de Y
r

*Ahora ya vamos a poder integrar el perfil de velocidad con respecto a su variable independiente.
 1
𝑉̅ 𝑥 = ∫ 𝑢(𝑦) ∙ (𝑏)𝑑𝑦
(ℎ ∙ 𝑏)

1
𝑉̅ 𝑥 = ∫ 𝑢(𝑦) ∙ 𝑑𝑦
ℎ

Pondremos los límites de integración los cuales van a ser de donde están definidos este dominio
ℎ ℎ
que es de − 2 𝑎 2 .
ℎ
1 2 1 𝜌𝑔𝑥 𝜕𝑃 1 2 1
𝑉̅ 𝑥 = ∫ [ (𝑦 ) − (ℎ2 )]𝑑𝑦
ℎ −ℎ 4 𝜇 𝜕𝑥 𝑦 ℎ
2

Integrar:
ℎ
1 1 𝜌𝑔𝑥 𝜕𝑃 2 1 2 1
𝑉̅ 𝑥 = ∙ ∫ [ (𝑦 ) − (ℎ2 )]𝑑𝑦
ℎ 4 𝜇 𝜕𝑥 −ℎ 𝑦 ℎ
2

ℎ ℎ
1 𝜌𝑔𝑥 𝜕𝑃 𝑦 2 2
𝑉̅ 𝑥 = ∙( | ) − (ℎ ∙ 𝑦)| 2 )
4ℎ 𝜇 𝜕𝑥 2 ℎ ℎ
− −
2 2

ℎ 2 ℎ 2
1 𝜌𝑔𝑥 𝜕𝑃 ( ) (− ) ℎ ℎ
𝑉̅ 𝑥 = ∙ [( 2 ) − ( 2 )] − [ℎ ∙ ( − (− ))]
4ℎ 𝜇 𝜕𝑥 2 2 2 2

0
ℎ2 ℎ2
1 𝜌𝑔𝑥 𝜕𝑃 ( 4 ) ( 4)
𝑉̅ 𝑥 = ∙[ − ] − [ℎ ∙ ℎ]
4ℎ 𝜇 𝜕𝑥 2 2

1 𝜌𝑔𝑥 𝜕𝑃
𝑉̅ 𝑥 = ∙ (−ℎ2 )
4ℎ 𝜇 𝜕𝑥

𝒉𝟐 𝝆𝒈𝒙 𝝏𝑷
̅𝒙 = −
𝑽
𝟒 𝝁 𝝏𝒙
𝝏𝒑 𝟒𝝁𝑽̅𝒙
= + 𝝆𝒈𝒙
𝝏𝒙 𝒉𝟐
Grafica

Conclusión
Las ecuaciones Navier-Stokes son una herramienta muy útil a la hora de analizar el
comportamiento de los fluidos (aunque con ciertos condicionamientos), ya que como
sabemos no siempre podremos predecir en su totalidad el comportamiento de éstos, pero
podremos tener una aproximación muy real y acertada de cómo sería.
Bibliografías:

Unidad 1. Transporte de cantidad de movimiento (momentum). (2017, 17 febrero).

Issuu.
https://issuu.com/jazinaruizhernandez/docs/unidad1.transportedecantidaddemovi
m

Luisyep, P. (2021, diciembre 4). Las Ecuaciones de Navier-Stokes. Ingenieria Basica.

https://ingenieriabasica.es/las-ecuaciones-de-navier-stokes/

MC Larrosa, J., & Perfil, V. T. mi. (s/f). Escuela de Economía. Blogspot.com.

Recuperado el 10 de diciembre de 2023, de https://eebb-
uns.blogspot.com/2018/02/matematica-la-complejidad-de-navier.html?m=1