Tarea 6 Bulk Viscosity
Tarea 6 Bulk Viscosity
Tarea 6 Bulk Viscosity
ρ ( ∂∂tu +u ∙ ∇ u)=−∇ P+ ρ g
La ecuación se puede resolver con la ayuda del tensor de viscosidades:
T =−P I + τ
Sustituyendo en la ecuación de Navier y en términos de la derivada total:
Dv
ρ = ρ g−∇ ∙ ( P I ) +∇ ∙ τ
Dt
En dónde para obtener la divergencia del tensor (PI)
∂ ∂ ∂
i [
∂x
( P I xx ) +
∂y
( ∂z ]
P I yx ) + ( P I zx )❑
∇ ∙ ( P I ) =+ j ∂ ( P I xy ) + ∂ ( P I yy ) + ∂ ( P I zy )
∂x [ ∂y ∂z ]
∂ ∂ ∂
+k
∂x [( P I xz ) + ( P I yz ) + ( P I zz )
∂y ∂z ]
Debido a que el tensor partió de la matriz identidad, los elementos fuera de la diagonal
principal son cero, los de la diagonal son uno:
∂P ∂P ∂P
∇ ∙ ( P I ) =i +j +k =∇ P
∂x ∂y ∂z
ρ ( ∂∂tu +u ∙ ∇ u)=−∇ P+ ρ g+ ∇ ∙ τ
El tensor τ se puede definir como:
2 ∂y
[(
τ ij =2 μ d ij + k− μ
3 ∂x )( )]
d kk δ ij
En términos de gradiente
∂ vi ∂ v j 2
τ ij =μ ( +
∂ x j ∂ xi ) [( ) ]
+ κ− μ δ ij
3
2 ∂ vk ∂ vi ∂ v j
(
τ ij = k − μ
3 ∂ xk )( )
δ ij +2 μ +
∂ x j ∂ xi [ ]
∂ ∂ 2 ∂ vk ∂ vi ∂ v j
∂ xk
(τ ij )=
∂ xk
k− μ
3 ((
∂ xk
+μ +
∂ x j ∂ xi )( ) [ ])
∂ ∂ 2 ∂ vi ∂ ∂vj
∂ xk
(τ ¿ ¿ij)=
∂ xj
k− μ+ μ
3 [(
∂ xi
+μ
∂ xi ∂ x i
¿
[ ]
) ]
∂ ∂ 1 ∂ vi ∂ ∂ vi
∂ xk
(τ ¿ ¿ij)=
∂ xj
k+ μ
3 ∂ xi (
+μ
∂ xi ∂ x i
¿) [ ]
∂ 1 ∂ ∂ vi ∂ ∂v j
∂ xk
(τ ¿ ¿ij)= k + μ (
3 ∂ x j ∂ xi
+μ )
∂ xi ∂ xi
¿
[ ] [ ]
1
( )
∇ ∙ τ= k + μ ∇ [ ∇ ∙ v ] + μ ∇ ∙ [ ∇ v ]
3
1
[( ) ]
∇ ∙ τ=∇ ∙ k + μ ∇ v + μ ∇ 2 v
3
1
( )
∇ ∙ T ii =−∇ P+ k + μ grad ∙÷ ( v ) + μ ∇2 v
3
∂u 1
ρ ( ∂t ) (
+u ∙ ∇ u =−∇ P+ ρ g+ k + μ grad ∙÷( v ) + μ ∇ 2 v
3 )
Como conclusión, el término de la viscosidad del bulk presentada en el artículo no es correcta. El
término es 1/3.