Ecuaciones de Cantidad de Movimiento
Ecuaciones de Cantidad de Movimiento
Ecuaciones de Cantidad de Movimiento
TF 1221
FENÓMENOS DE TRANSPORTE I
GUIA DE
ECUACIONES DE CONTINUIDAD
Y CANTIDAD DE MOVIMIENTO
INTRODUCCIÓN
OBJETIVOS
Terminal ¾ Proporcionar al estudiante los conceptos básicos de
Mecánica de Fluidos y su aplicación a problemas de
interés práctico.
ÍNDICE
Página
Introducción 2
Objetivos 3
Índice 4
Bibliografía 18
Ecuaciones de continuidad y cantidad de movimiento. 5
LA ECUACIÓN DE CONTINUIDAD
∆y ∆x
∆x∆y∆z
∂ρ
∂t
[ ] ( ) y − (ρv y ) y + ∆y ⎤⎥⎦
= ∆y∆z ( ρv x ) x − ( ρv x ) x + ∆x + ∆x∆z ⎡ ρv y
⎢⎣
[
+ ∆x∆y ( ρv z ) z − ( ρv z ) z + ∆z ] (2)
Ecuaciones de continuidad y cantidad de movimiento. 6
∂ρ ⎛∂ ∂ ∂ ⎞
∂t
= −⎜⎜ ( ρv x ) + ( )
ρv y + ( ρv z )⎟⎟ (3)
⎝ ∂x ∂y ∂z ⎠
Ésta es la ecuación de continuidad y puede escribirse en una forma más
conveniente utilizando notación generalizada:
∂ρ
= −(∇ • ρ v ) (4)
∂t
El término (∇ • ρ v ) se denomina divergencia de ρv, en donde la
∂ρ ∂
+ (ρv x ) + ∂ ρv y + ∂ (ρv z ) = 0
( )
∂t ∂x ∂y ∂z
∂ρ 1 ∂
+ (ρrv r ) + 1 ∂ (ρv θ ) + ∂ (ρv z ) = 0
∂t r ∂r r ∂θ ∂z
∂ρ 1 ∂
+ 2
∂t r ∂r
(
ρr 2 v r +
1
) ∂
r sen θ ∂θ
(ρv θ sen θ) + 1 ∂ (ρv φ ) = 0
r sen θ ∂φ
Ecuaciones de continuidad y cantidad de movimiento. 8
LA ECUACIÓN DE MOVIMIENTO
Para un elemento diferencial de volumen ∆x∆y∆z, se puede escribir el
siguiente balance de cantidad de movimiento.
(τzx)]z+∆z
z
(τyx)]y+∆y
(τxx)]x (τxx)]x+∆x
(τyx)]y
y
(τzx)]z
[ ]
∆y∆z ( ρv x v x ) x − ( ρv x v x ) x + ∆x + ∆x∆z ⎡ ρv y v y
⎢⎣
( ) y − (ρv y v y ) y + ∆y ⎤⎥⎦
(5)
[
+ ∆x∆y (ρv z v z ) z − ( ρv z v z ) z + ∆z ]
De igual forma, la velocidad con que el componente x de la cantidad de
movimiento entra por transporte molecular por la cara situada en x es τxx]x∆y∆z,
y con la que sale por x+ ∆x es τxx]x+∆x∆y∆z; para las demás caras se pueden
escribir expresiones similares. Hay que tomar en cuenta que τij es la densidad de
flujo de cantidad de movimiento j a través de una cara perpendicular al eje i.
Sumando estas seis contribuciones, se obtiene:
[
∆y∆z τ xx x
− τ xx x + ∆x
] + ∆x∆z ⎡⎢⎣τ yx y − τ yx ⎤ + ∆x∆y τ
y + ∆y ⎥
⎦
[ zx z − τ zx z + ∆z
] (6)
ecuación por ∆x∆y∆z, y tomando el limite cuando ∆x, ∆y, ∆z tiende a cero, se
tiene
⎛ ∂v x ∂v ∂v ∂v ⎞ ∂P ⎛ ∂τ xx ∂τ yx ∂τ zx ⎞
ρ ⎜⎜ + v x x + v y x + v z x ⎟⎟ = − −⎜ + + ⎟ + ρg x (8)
⎝ ∂t ∂x ∂y ∂z ⎠ ∂x ⎜⎝ ∂x ∂y ∂z ⎟
⎠
Los componentes y y z, se pueden obtener de forma análoga. Ésta
ecuación se puede escribir en notación vectorial:
Ecuaciones de continuidad y cantidad de movimiento. 10
∂
∂t
( )
(ρ v ) = −(∇ • ρ vv ) − (∇ • P ) − ∇ • τ + ρ g (9)
donde:
∂
(ρ v ) es la velocidad de acumulación de cantidad de movimiento
∂t
(∇ • ρ vv ) es la velocidad de ganancia de cantidad de movimiento por
convección
(∇ • P ) es la fuerza de presión que actúa sobre el elemento
(∇ • τ ) es la velocidad de ganancia de cantidad de movimiento por
transporte viscoso
ρg es la fuerza de gravedad que actúa sobre el elemento
LA ECUACION DE MOVIMIENTO EN
COORDENADAS RECTANGULARES (x, y, z):
En función de τ:
Componente x:
⎛ ∂v ∂v x ∂v x ∂v x ⎞ ∂P ⎛ ∂τ xx ∂τ yx ∂τ zx ⎞
ρ⎜⎜ x + v x + vy + vz ⎟⎟ = − −⎜ + + ⎟ + ρg x
⎝ ∂t ∂x ∂y ∂z ⎠ ∂x ⎜⎝ ∂x ∂y ∂z ⎟
⎠
Componente y:
⎛ ∂v y ∂v y ∂v y ∂v y ⎞ ∂P ⎛ ∂τ xy ∂τ yy ∂τ zy ⎞
ρ⎜⎜ + vx + vy + vz ⎟=−
⎟ − ⎜⎜ + + ⎟ + ρg y
⎟
⎝ ∂t ∂x ∂y ∂z ⎠ ∂y ⎝ ∂x ∂y ∂z ⎠
Componente z:
⎛ ∂v ∂v z ∂v z ∂v ⎞ ∂P ⎛ ∂τ xz ∂τ yz ∂τ zz ⎞
ρ⎜⎜ z + v x + vy + vz z ⎟⎟ = − −⎜ + + ⎟ + ρg z
⎝ ∂t ∂x ∂y ∂z ⎠ ∂z ⎜⎝ ∂x ∂y ∂z ⎟
⎠
Componente x:
⎛ ∂v ∂v x ∂v x ∂v x ⎞ ∂p ⎛⎜ ∂ 2 v x ∂ 2 v x ∂ 2 v x ⎞
⎟ + ρg x
ρ⎜⎜ x + v x + vy + vz ⎟⎟ = − + µ + +
⎝ ∂t ∂x ∂y ∂z ⎠ ∂x ⎜⎝ ∂x 2 ∂y 2 ∂z 2 ⎟
⎠
Componente y:
⎛ ∂v y ∂v y ∂v y ∂v y ⎞ ∂p ⎛ ∂ 2v y ∂ 2v y ∂ 2v y ⎞
ρ⎜⎜ + vx + vy + vz ⎟=− + µ ⎜ + + ⎟ + ρg
∂t ∂x ∂y ∂z ⎟ ∂y ⎜ ∂ 2
∂ 2
∂z 2 ⎟ y
⎝ ⎠ ⎝ x y ⎠
Componente z:
⎛ ∂v ∂v z ∂v z ∂v ⎞ ∂P ⎛ ∂ 2v z ∂ 2v z ∂ 2v z ⎞
ρ⎜⎜ z + v x + vy + v z z ⎟⎟ = − + µ⎜ + + ⎟ + ρg z
⎜ ∂x 2 ⎟
⎝ ∂t ∂x ∂y ∂z ⎠ ∂z ⎝ ∂y 2 ∂z 2 ⎠
Ecuaciones de continuidad y cantidad de movimiento. 12
En función de τ:
Componente r:
⎛ ∂v r ∂v r v θ ∂v r v θ2 ∂v ⎞ ∂P ⎛ 1 ∂ 1 ∂τ rθ τ θθ ∂τ rz ⎞
ρ⎜ + vr + − + vz r ⎟ = − − ⎜⎜ (rτ rr ) + − + ⎟⎟ + ρgr
⎜ ∂t ∂r r ∂θ r ⎟
∂z ⎠ ∂r ⎝ r ∂r r ∂θ r ∂z ⎠
⎝
Componente θ:
⎛ ∂v ∂v θ v θ ∂v θ v r v θ ∂v θ ⎞ 1 ∂P ⎛ 1 ∂ 2 1 ∂τ θθ ∂τ zθ ⎞
ρ⎜⎜ θ + Vr + + + vz v ⎟⎟ = − − ⎜⎜ 2 (r τ rθ ) + + ⎟⎟ + ρg θ
⎝ ∂t ∂r r ∂θ r ∂z ⎠ r ∂θ ⎝ r ∂r r ∂θ ∂z ⎠
Componente z:
⎛ ∂v ∂v z v θ ∂v z ∂v z ⎞ ∂P ⎛ 1 ∂ 1 ∂τ θz ∂τ zz ⎞
ρ⎜⎜ z + v r + + vz ⎟⎟ = − − ⎜⎜ (rτ rz ) + + ⎟⎟ + ρg z
⎝ ∂t ∂r r ∂θ ∂z ⎠ ∂z ⎝ r ∂r r ∂θ ∂z ⎠
Componente θ:
⎛ ∂v θ ∂v v ∂v θ v r v θ ∂v ⎞ 1 ∂P
ρ ⎜⎜ + vr θ + θ + + vz θ ⎟⎟ = −
⎝ ∂t ∂r r ∂θ r ∂z ⎠ r ∂θ
⎡ ∂ ⎛1 ∂ ⎞ 1 ∂ vθ
2
2 ∂v r ∂ v θ ⎤
2
+ µ⎢ ⎜ (rv θ )⎟ + 2 2 + 2 + ⎥ + ρgθ
⎣⎢ ∂r ⎝ r ∂r ⎠ r ∂θ r ∂θ ∂z 2 ⎥⎦
Componente z:
⎛ ∂v ∂v z v θ ∂v z ∂v z ⎞ ∂P ⎡1 ∂ ⎛ ∂v z ⎞ 1 ∂ vz ∂ vz ⎤
2 2
ρ⎜⎜ z + v r + + vz ⎟⎟ = − + µ⎢ ⎜⎜ r ⎟⎟ + 2 + ⎥ + ρg z
⎝ ∂t ∂r r ∂θ ∂z ⎠ ∂z ⎣⎢ r ∂r ⎝ ∂r ⎠ r ∂θ ∂z 2 ⎦⎥
2
Ecuaciones de continuidad y cantidad de movimiento. 13
LA ECUACIÓN DE MOVIMIENTO EN
COORDENADAS ESFÉRICA (r, θ, Φ)
En función de los esfuerzos (τ):
Componente r:
⎛ ∂v ∂v r v θ ∂v r v φ ∂v r v θ + v φ
2 2 ⎞
ρ⎜ r
+ vr + + − ⎟
⎜ ∂t ∂r r ∂θ r sen θ ∂φ r ⎟
⎝ ⎠
∂τ rφ τ θθ + τ φφ ⎞
=−
∂P ⎛ 1 ∂ 2
− ⎜⎜ 2
∂r ⎝ r ∂r
r τ rr +( 1 ∂
r sen θ ∂θ
)
(τ rθ Senθ) + 1
r sen θ ∂φ
−
r
⎟ + ρgr
⎟
⎠
Componente θ:
⎛ ∂v ∂v θ v θ ∂v θ v φ ∂v θ v r v θ v φ ⋅ cot θ ⎞⎟
2
ρ⎜ θ
+ vr + + + −
⎜ ∂t ∂r r ∂θ r sen θ ∂φ r r ⎟
⎝ ⎠
∂τ θφ τ rθ cot θ
=−
1 ∂P ⎛ 1 ∂ 2
− ⎜⎜ 2 ⋅
r ∂θ ⎝ r ∂r
(
r τ rθ +
1 ∂
r sen θ ∂θ
)(τ θθ sen θ) + 1
r sen θ ∂φ
+
r
−
r
⎞
τ φφ ⎟⎟ + ρg θ
⎠
Componente φ:
⎛ ∂v φ ∂v φ v θ ∂v φ v φ ∂v φ v r v φ v φ v θ cot θ ⎞
ρ⎜⎜ + vr + + + + ⎟
⎟
⎝ ∂t ∂r r ∂θ r sen θ ∂φ r r ⎠
1 ∂τ θφ 1 ∂τ φφ τ rφ 2 cot θ
=−
1 ∂P ⎛ 1 ∂ 2
− ⎜⎜ 2
r sen θ ∂φ ⎝ r ∂r
r τ rφ + (
r ∂θ
+ )
r sen θ ∂φ
+
r
+
r
⎞
τ θφ ⎟⎟ + ρg φ
⎠
∂P ⎛ 2 2 ∂v θ 2 2 ∂v φ ⎞
=− + µ⎜⎜ ∇ 2 v r − 2 v r − 2 − 2 v θ cot θ − 2 ⎟ + ρg r
⎟
∂r ⎝ r r ∂θ r r sen θ ∂φ ⎠
Ecuaciones de continuidad y cantidad de movimiento. 14
Componente θ
⎛ ∂v ∂v θ v θ ∂v θ v φ ∂v θ v r v θ v φ cot θ ⎞⎟
2
ρ⎜ θ
+ vr + + + −
⎜ ∂t ∂r r ∂θ r sen θ ∂φ r r ⎟
⎝ ⎠
1 ∂P ⎛ 2 ∂v r v 2 cos θ ∂v φ ⎞
=− + µ⎜⎜ ∇ 2 v θ + 2 − 2 θ2 − 2 ⎟ + ρg θ
r ∂θ ⎝ r ∂θ r sen θ r sen 2 θ ∂φ ⎟⎠
Componente φ
⎛ ∂v φ ∂v φ v θ ∂v φ v φ ∂v φ v r v φ v φ v θ cot θ ⎞
ρ⎜⎜ + vr + + + − ⎟
⎟
⎝ ∂t ∂r r ∂θ r sen θ ∂φ r r ⎠
1 ∂P ⎛ v 2 ∂v r 2 cos θ ∂v θ ⎞
=− + µ⎜⎜ ∇ 2 v φ − 2 θ 2 + 2 + 2 ⎟⎟ + ρg φ
r sen θ ∂φ ⎝ r sen θ r sen θ ∂φ r sen 2 θ ∂φ ⎠
donde
1 ∂ ⎛ 2 ∂a ⎞ 1 ∂ ⎛ ∂a ⎞ 1 ∂ ⎛ ∂a 2 ⎞
∇a 2 = ⎜ r ⎟ + ⎜ sen θ ⎟ + ⋅ ⎜ ⎟
r 2 ∂r ⎝ ∂r ⎠ r 2 sen θ ∂θ ⎝ ∂θ ⎠ r 2 sen 2 θ ∂θ ⎜⎝ ∂φ 2 ⎟
⎠
Ecuaciones de continuidad y cantidad de movimiento. 15
⎡ ∂v ⎤
τ xx = −µ ⎢2 x − (∇ • v )⎥
2
⎣ ∂x 3 ⎦
⎡ ∂v y 2 ⎤
τ yy = −µ ⎢2 − (∇ • v )⎥
⎣ ∂y 3 ⎦
⎡ ∂v ⎤
τ zz = −µ ⎢2 z − (∇ • v )⎥
2
⎣ ∂z 3 ⎦
⎡ ∂v ∂v y ⎤
τ xy = τ yx = −µ ⎢ x + ⎥
⎣ ∂y ∂x ⎦
⎡ ∂v ∂v y ⎤
τ zy = τ yz = −µ ⎢ z + ⎥
⎣ ∂y ∂z ⎦
⎡ ∂v ∂v ⎤
τ zx = τ xz = −µ ⎢ x + z ⎥
⎣ ∂z ∂x ⎦
donde
∂v y
(∇ • v ) = ∂v x + +
∂v z
∂x ∂y ∂z
Ecuaciones de continuidad y cantidad de movimiento. 16
⎡ ∂v ⎤
τ rr = −µ ⎢2 r − (∇ • v )⎥
2
⎣ ∂r 3 ⎦
⎡ ⎛ 1 ∂v θ v r ⎞ 2 ⎤
τ θθ = −µ ⎢2⎜⎜ + ⎟⎟ − (∇ • v )⎥
⎣ ⎝ r ∂θ r ⎠ 3 ⎦
⎡ ∂v ⎤
τ zz = −µ ⎢2 z − (∇ • v )⎥
2
⎣ ∂z 3 ⎦
⎡ ∂ ⎛v ⎞ 1 ∂v r ⎤
τ rθ = τ θr = −µ ⎢r ⎜⎜ θ ⎟⎟ + ⎥
⎣ ∂r ⎝ r ⎠ r ∂θ ⎦
⎡ ∂v 1 ∂v z ⎤
τ zθ = τ θz = −µ ⎢ θ + ⎥
⎣ ∂z r ∂θ ⎦
⎡ ∂v ∂v ⎤
τ zr = τ rz = −µ ⎢ r + z ⎥
⎣ ∂z ∂r ⎦
donde
∂v θ ∂v z
(∇ • v ) = 1 ∂ (rv r ) + 1 +
r ∂r r ∂θ ∂z
Ecuaciones de continuidad y cantidad de movimiento. 17
⎡ ∂v ⎤
τ rr = −µ ⎢2 r − (∇ • v )⎥
2
⎣ ∂r 3 ⎦
⎡ ⎛ 1 ∂v θ v r ⎞ 2 ⎤
τ θθ = −µ ⎢2⎜⎜ + ⎟⎟ − (∇ • v )⎥
⎣ ⎝ r ∂θ r ⎠ 3 ⎦
⎡ ⎛ 1 ∂v φ v r v θ cot θ ⎞ 2 ⎤
τ φφ = −µ ⎢2⎜⎜ + + ⎟ − (∇ • v )⎥
⎟ 3
⎣⎢ ⎝ r sen θ ∂φ r r ⎠ ⎦⎥
⎡ ∂ ⎛v ⎞ 1 ∂v r ⎤
τ rθ = τ θr = −µ ⎢r ⎜⎜ θ ⎟⎟ + r ∂θ ⎥
⎣ ∂r ⎝ r ⎠ ⎦
⎡ sen θ ∂ ⎛ v φ ⎞ 1 ∂v θ ⎤
τ φθ = τ θφ = −µ ⎢ ⎜ ⎟+ ⎥
⎜ ⎟
⎢⎣ r ∂θ ⎝ sen θ ⎠ r sen θ ∂φ ⎥⎦
⎡ 1 ∂v r ∂ ⎛ vφ ⎞⎤
τ φr = τ rφ = −µ ⎢ + r ⎜⎜ ⎟⎥
⎟
⎣⎢ r sen θ ∂φ ∂r ⎝ r ⎠⎦⎥
donde
∂v φ
(∇ • v ) = 1 ∂ 2
(
r vr + )1 ∂
(v θ sen θ) + 1
r ∂r
2
r sen θ ∂θ r sen θ ∂φ
Ecuaciones de continuidad y cantidad de movimiento. 18
BIBLIOGRAFÍA