求高手分析下十分简单的变径圆管流动的控制方程

laoyangsky

我是刚学CFD的，不太清楚。
假设：没有热交换；层流；稳态；密度和μ为常数；水平流动；不考虑重力。
求解。
我的想法是：利用柱坐标系考虑。Vr、Vθ为0，Vz只与r，z有关
连续方程为：ρə（Vz）/əz=0
动量方程：r和θ分量为0，只有z分量。
ρ（Vr*ə（Vz）/ər+Vz*ə（Vz）/əz）=-ə（P）/əz+μ[(1/r)*ə（r*Vz）/ər+ə（Vz^2）/əz^2]
能量方程不用考虑。

真心麻烦高手帮我看看 我的分析有没有错。