[工学]高等计算流体力学讲义8.pdfVIP

高等计算流体力学讲义（8） 第三章 不可压缩流动的数值方法 §1 基本方程及其性质 一、基本方程 考虑不可压缩 NS 方程： ∇iu 0 （1） ∂ρu +∇i uu f =−∇ +∇i (ρ ) ρ p τ （2 ） ∂t 其中粘性应力为， τ 2μS （3 ） S 1 (=∇u +∇uT ) 2 如果粘性系数为常数， ∇iτ μ=Δu （4 ） 常粘性系数不可压缩 NS 方程可以写为： ∇iu 0 ∂u ， +∇i(uu) f =−∇p +υΔu ∂t 其中υ μ/ ρ 为运动粘性系数。NS 方程也可以写为无量纲化形式 ∇iu 0 ∂u 1 +∇i(uu) f =−∇p + Δu ∂t Re 其中 ρ 已经吸收到p 中（p 代表p / ρ ）。不可压缩方程的边界条件为： 固体壁面：u uwall ， 进口条件： u uin ， 1 ∂u 出口条件： 0 。 ∂n 不可压缩方程中的压力场可以相差任一常数而对速度场无影响，所以压力场只是在相差任意 常数的条件下是确定的。为了确定全场压力值，还应指定流场中某一点的压力。 二、不可压 N －S 方程的特点： （1） 方程为二阶偏微分方程，二阶项中包含参数μ （粘性系数）。 边界层、分离、湍流… （2 ） 方程是非线性的，表现为对流项∇i(uu) 。 ∂u 对一维问题，非线性项为u 。假定u 的波数为k 的Fourier 分量为 ∂x u u(t)sin kx （5 ） ∂u 1 2 1 2 则：u u sin 2kx 。即振幅由u → u ；波数由k →2k 。也就是说，振幅呈现非线 x ∂ 2 2 性变化，且可以产生高频成分。粘性的作用，使得解的结构进一步复杂化，考虑模型方程 2 ∂u 1 ∂ u 2