Physic Labs

流体力学

纳维–斯托克斯方程

纳维–斯托克斯方程描述黏性流体的动量平衡,综合惯性、压强、体积力和黏性应力。

纳维–斯托克斯方程将牛顿第二定律应用于流体微元。方程形式取决于本构模型和假设;下式适用于均匀、不可压缩、黏度恒定的牛顿流体。

ρ(∂u∂t+(u⋅∇)u)=−∇p+μ∇2u+ρf,∇⋅u=0\rho\left(\frac{\partial\mathbf u}{\partial t}+(\mathbf u\cdot\nabla)\mathbf u\right)=-\nabla p+\mu\nabla^2\mathbf u+\rho\mathbf f,\qquad\nabla\cdot\mathbf u=0

定义: 动量方程项

ρ\rho 为密度,u\mathbf u 为速度,pp 为压强,μ\mu 为动力黏度,f\mathbf f 为单位质量体积力。左侧为局部与对流加速度;右侧包括压强梯度、黏性扩散和体积力。

研究平板间库埃特速度分布及黏度影响;此处是简单层流解,并非纳维–斯托克斯方程的一般解。

动量平衡

随体加速度包括局部变化 ∂tu\partial_t\mathbf u 与对流项 (u⋅∇)u(\mathbf u\cdot\nabla)\mathbf u。低 ReRe 时黏性会平滑速度梯度;高 ReRe 时对流和惯性更重要,流动通常更难预测。

基准解:库埃特流

ux(y)=UyH,0≤y≤Hu_x(y)=U\frac{y}{H},\qquad 0\le y\le H

两平行板相距 HH,下板静止、上板以速度 UU 移动;若沿板无压强梯度,牛顿流体的稳定层流速度分布为线性。无滑移边界条件使流体在壁面处速度与板相同。

例题: 间隙中的库埃特速度

两平板间距 H=2.0 mmH=2.0\,\mathrm{mm},上板以 U=0.10 m/sU=0.10\,\mathrm{m/s} 运动。求库埃特流的中点速度。

解答

u(H/2)=U(H/2)/H=U/2=0.050 m/su(H/2)=U(H/2)/H=U/2=0.050\,\mathrm{m/s}。

快速检测

各项都有明确物理意义,单位为单位体积的力:hoDmathbfu/Dt ho Dmathbf u/Dt 是惯性项,−ablap- abla p 是压强力,muabla2mathbfumu abla^2mathbf u 是黏性导致的动量扩散,homathbff homathbf f 是单位体积的体积力。不可压缩流体还须满足 ablacdotmathbfu=0 ablacdotmathbf u=0;若没有壁面无滑移等边界条件,解并不确定。

所示不可压缩方程基于哪种流体模型?

示例中平板间隙中点的库埃特速度是多少?

参考文献

  1. Frank M. White (2016). Fluid Mechanics
  2. G. K. Batchelor (1967). An Introduction to Fluid Dynamics