Fiche NavierStokes2
Fiche NavierStokes2
Fiche NavierStokes2
∂v
∂ρ/∂t + div(ρv) = 0 (continuity equation) ρ + ρ(v.−→
grad)v = ρg − −→ p + η∆v
grad (Navier-Stokes equation)
∂t
In cylindrical coordinates z
continuity equation:
P
∂ρ 1 ∂ 1 ∂(ρuϕ ) ∂(ρuz ) ρ z y
+ (ρrur ) + + =0
∂t r ∂r r ∂ϕ ∂z
x φ illustration: Wikimedia/Cylindrical_Coordinates
∂ur u2ϕ
" ! #
∂ur ∂ur uϕ ∂ur 1 ∂p 1 ∂ ∂ur 1 ∂ 2 ur ∂ 2 ur ur 2 ∂uϕ
+ ur + + uz − =− + gr + ν r + 2 2
+ 2
− 2 − 2
∂t ∂r r ∂ϕ ∂z r ρ ∂r r ∂r ∂r r ∂ϕ ∂z r r ∂ϕ
" ! #
2 2
∂uϕ ∂uϕ uϕ ∂uϕ ∂uϕ ur uϕ 1 ∂p 1 ∂ ∂uϕ 1 ∂ uϕ ∂ uϕ 2 ∂ur uϕ
+ ur + + uz + =− + gϕ + ν r + 2 + + 2 − 2
∂t ∂r r ∂ϕ ∂z r ρr ∂ϕ r ∂r ∂r r ∂ϕ2 ∂z 2 r ∂ϕ r
" ! #
∂uz ∂uz uϕ ∂uz ∂uz 1 ∂p 1 ∂ ∂uz 1 ∂ 2 uz ∂ 2 uz
+ ur + + uz =− + gz + ν r + 2 +
∂t ∂r r ∂ϕ ∂z ρ ∂z r ∂r ∂r r ∂ϕ2 ∂z 2
stress tensor:
" #
∂ur ∂ uϕ 1 ∂ur
τrr = −p + 2η τrϕ =η r +
∂r ∂r r r ∂ϕ
! !
1 ∂uϕ ur ∂uϕ 1 ∂uz
τϕϕ = −p + 2η + τϕz =η +
r ∂ϕ r ∂z r ∂ϕ
!
∂uz ∂ur ∂uz
τzz = −p + 2η τϕz =η +
∂z ∂z ∂r
In spherical coordinates z
θ
P
ρ
continuity equation:
y
∂ρ 1 ∂ 2 1 ∂(ρuϕ ) 1 ∂
+ 2 ρr ur + + (sin(θ)ρuθ ) = 0 φ
∂t r ∂r r sin(θ) ∂ϕ r sin(θ) ∂θ
x illustration: Wikimedia/Spherical_Coordinates