Analisis Diferencial - Placas Inclinadas

ECUACIONES DE NAVIER‐STOKES EN
COORDENADAS CARTESIANAS
DV
∇ · V =0 y ρ DT = -∇ P + ρ g + μ∇2 V 
Ecuación de Continuidad
∂u ∂v ∂w
∂x + ∂y + ∂z =0
Cantidad de Movimiento en dirección X
∂2 u ∂2 u ∂2 u
ρ( ∂u
∂t +u
∂u ∂u ∂u ∂P
∂x +v ∂y +w ∂z )=‐ ∂x +ρgx+μ( ∂x2 + ∂y2 + ∂z2 )
Cantidad de Movimiento en dirección Y

∂2 v ∂2 v ∂2 v
ρ( ∂v
∂t +u
∂v
∂x +v ∂v
∂y +w ∂v
∂z )=‐ ∂P
∂y +ρg y +μ( ∂x2 + ∂y2 + ∂z2 )
Cantidad de Movimiento en dirección Z

∂2 w ∂2 w ∂2 w
ρ( ∂w
∂t +u
∂w
∂x +v ∂w
∂y +w ∂w
∂z )=‐ ∂P
∂z +ρg z +μ( ∂x2 + ∂y2 + ∂z2 )
2) ENUNCIADO
Una película de agua de espesor constante (h) fluye en movimiento paralelo (laminar) entre dos placas inclinadas como lo
muestra la figura. La placa superior se mueve con velocidad Uo mientras que la inferior está quieta. Ignore la variación de la
presión en la dirección X.
a) Simplifique las ecuaciones de Navier Stokes, indicando claramente los supuestos hechos para encontrar las ecuaciones que
gobiernan el movimiento.
b) Encuentre el perfil de velocidades. Uo=2 m/s, h=0.04 m3 y ∅=0.1° μ = 10 -3 kg / (m * s), ρ = 1000 kg / m3
Printed by Wolfram Mathematica Student Edition

2 12. Analisis Diferencial E.C. y N.S. - Placas Inclinadas.nb
HIPÓTESIS
1. Placas infinitas en dirección X
∂
2.Flujo Permanente ∂t
=0
3. Flujo laminar (paralelo) unidireccional V=0, W=0.
4. Flujo Incompresible y Newtoniano μ=cte y ρ=cte.
∂P
5. No hay gradiente de presión en dirección X ∂x
=0
∂V
6. No hay variación de la velocidad en Z ∂z
=0
7. La gravedad posee componentes x, y
SIMPLIFICACIÓN DE LAS ECUACIONES

12. Analisis Diferencial E.C. y N.S. - Placas Inclinadas.nb 3
CONTINUDAD
3 3y 6
∂u ∂v ∂w ∂u
+ + =0 entonces tenemos que =0
∂x
∂x ∂y ∂z
CANTIDAD DE MOVIMIENTO EN DIRECCIÓN X

continuidad continuidad
2 3 3 5 3 6
∂u ∂u ∂u ∂u ∂P ∂2 u ∂2 u ∂2 u
ρ +u +v +w =- + ρ gx + μ + +
∂t ∂x ∂y ∂z ∂x ∂ x2 ∂ y2 ∂ z2
CANTIDAD DE MOVIMIENTO EN DIRECCIÓN Y
3 3 3 3 3 3
2
∂v ∂v ∂v ∂v ∂P ∂2 v ∂2 v ∂2 v
ρ +u +v +w =- + ρ gy + μ + +
∂t ∂x ∂y ∂z ∂y ∂ x2 ∂ y2 ∂ z2
CANTIDAD DE MOVIMIENTO EN DIRECCIÓN Z

3 3 3 6 6 3
2 3 3
∂w ∂w ∂w ∂w ∂P ∂2 w ∂2 w ∂2 w
ρ +u +v +w =- + ρ gz + μ + +
∂t ∂x ∂y ∂z ∂z ∂ x2 ∂ y2 ∂ z2
Solución de la ecuaciones
Cantidad de movimiento en X
∂2 u
0 = ρgx + μ
∂ y2
∂2 u
0 = ρ(gsen(θ)) + μ Integrando respecto a "y" y despejando
∂ y2
ρ(gsen(θ)) ∂u
- y + C1 = Integrando nuevamente con respecto a "y"
2μ ∂y
ρ(gsen(θ)) y2
- + C1 y + C2 = U
2μ 2
Cantidad de movimiento en Y
∂P
0=- + ρ * gy
∂y
∂P
0=- + ρ * g * cos(θ)
∂y

∂P
= ρ * g * cos(θ) + C1 Despejando
∂y
P = g y ρ cos(θ) + Po Integrando con respecto a "y"
Cantidad de movimiento en z
0=0
Determinación del perfil de velocidades

y=0 V=0 (ecuación 1)
y=h V=‐Uo (ecuación 2)
Teniendo en cuenta la ecuación (1)
ρ(g * sen(θ)) y2
- * + C1 * y + C2 = U
μ 2
0 + 0 + C2 = 0 entonces C2 = 0 Reemplazando
ρ(g * sen(θ)) y2
U=- * + C1 * y (3)
μ 2
Teniendo en cuenta la ecuación (1) y (2)
ρ(g * sen(θ)) h2
- * + C1 * h = -Uo
μ 2
ρ(g * sen(θ)) h Uo
C1 = * - Reemplazando
μ 2 h
ρ(g * sen(θ)) y2 ρ(g * sen(θ)) h Uo

- * + * - *y = U
μ 2 μ 2 h
ρ(g * sen(θ)) Uo
U= * yh - y2 - *y
2*μ h
Sabiendo además que : Uo = 2 m / s, h = 0.04 mm y ∅ = 0.1 ° μ = 10-3 kg  (m * s), ρ = 1000 kg  m3
1000 (9.81 * sen(0.1 °)) 2

U= * (0.04 y - y2) - *y
2 * 10-3 0.04

12. Analisis Diferencial E.C. y N.S. - Placas Inclinadas.nb 5
U = -8561 * y2 + 292.43 * y
Datos
h = 0.04; ϕ = 0.1 °; μ = 10-3; ρ = 1000; g = 9.81; Uo = 2;
ρ g Sin[ϕ] Uo
U= y h - y2 - y
2μ h
-50. y + 8560.84 0.04 y - y2
Expand[U]
292.433 y - 8560.84 y2
PlotU, y, 0, h, Frame → True, GridLines → Automatic
-1
-2
0.00 0.01 0.02 0.03 0.04
Determinación del Caudal
Teniendo en cuenta la definición del caudal Q = ∫ ∫ V ⅆA

h
Q=b * ∫o V ⅆy
h
Q= ∫o (-8561 * y2 + 292.43 * y) ⅆy Con b=1.0
y3 y2
Q=-8561 * 3 + 292.43 * 2 Reemplazando
0.043 0.042
Q=-8561 * 3 + 292.43 * 2 = 0.051 m3 / s
Caudal por unidad de ancho
h
 U ⅆy
0
0.0513156
Velocidad Media
Q ∫ ∫ V ⅆA Q 0.0510 m3 / s
Vmedia = = = = = 1.28 m / s
A A B*H 1.0 m * 0.04 m
Tensor Tasa Deformación

Por definición tenemos que:

∂U 1
∂x 2  ∂U
∂y
+ ∂V
∂x

ϵij = 1 ∂U ∂V ∂V
2  ∂y + ∂x  ∂x
∂U
∂x
=0
∂U
∂y
= -2 * 8561 * y + 292.43 = -2 * 8561 * 0.04 + 292.43 = -392.45
∂V
∂x
=0
∂V
∂y
=0
Por tanto tenemos que:
1
0 2 (-392.45) 0 -196.2
ϵij = 1 = 
2 (-392.45) 0 -196.2 0

Analisis Diferencial - Placas Inclinadas

Cargado por

Copyright:

Formatos disponibles

Analisis Diferencial - Placas Inclinadas

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Analisis Diferencial - Placas Inclinadas

Cargado por

Copyright:

Formatos disponibles

ECUACIONES DE NAVIER‐STOKES EN

Cantidad de Movimiento en dirección Y

Cantidad de Movimiento en dirección Z

b) Encuentre el perfil de velocidades. Uo=2 m/s, h=0.04 m3 y ∅=0.1° μ = 10 -3 kg / (m * s), ρ = 1000 kg / m3

Printed by Wolfram Mathematica Student Edition

3. Flujo laminar (paralelo) unidireccional V=0, W=0.

4. Flujo Incompresible y Newtoniano μ=cte y ρ=cte.

7. La gravedad posee componentes x, y

SIMPLIFICACIÓN DE LAS ECUACIONES

Printed by Wolfram Mathematica Student Edition

CANTIDAD DE MOVIMIENTO EN DIRECCIÓN X

CANTIDAD DE MOVIMIENTO EN DIRECCIÓN Y

CANTIDAD DE MOVIMIENTO EN DIRECCIÓN Z

Printed by Wolfram Mathematica Student Edition

Determinación del perfil de velocidades

Teniendo en cuenta la ecuación (1)

Teniendo en cuenta la ecuación (1) y (2)

ρ(g * sen(θ)) y2 ρ(g * sen(θ)) h Uo

Sabiendo además que : Uo = 2 m / s, h = 0.04 mm y ∅ = 0.1 ° μ = 10-3 kg  (m * s), ρ = 1000 kg  m3

1000 (9.81 * sen(0.1 °)) 2

Printed by Wolfram Mathematica Student Edition

Determinación del Caudal

Teniendo en cuenta la definición del caudal Q = ∫ ∫ V ⅆA

Tensor Tasa Deformación

Printed by Wolfram Mathematica Student Edition

Printed by Wolfram Mathematica Student Edition

También podría gustarte