库埃特流动和泊肃叶流动精选文档.pptVIP

2020/5/12 31 按 常微分方程 的记法并利用 μ = μ （ T ） d μ /d y= (d μ /dT)(dT/d y) d (k dT/ d y ） / dy+ τ (du/dT)(dT/d y) = 0 边界条件 y=0,T=Tw; y=D,T=Te 两点边界值问题 整理后方程组为 ? （ μ ? u/ ? y ） / ? y= ?τ / ? y=0 d (k dT/ d y ） / dy+ τ (du/dT)(dT/d y) = 0 τ = μ （ u e /D ） 边界条件 y=0,T=Tw; y=D,T=Te 2020/5/12 32 数值解法 —— 打靶法 1. 假设 τ , 如 τ = μ （ u e /D ）、 u(y) —— 不可压解为 初值 2. 边界条件 y=0,T=Tw; y=D,T=Te 变为 3. 边界条件 y=0,T=Tw;y=0, (dT/d y)w 已知 4. 通常为不可压流解 5. 求解 6. d (k dT/ d y ） / dy+ τ (du/dT)(dT/d y) = 0 3. 直到 y=D, 检查 T=Te?, 否则重假设 (dT/d y)w ， 回到 2 步计算，重复直到 T=Te 4. 得到了 T=T(y) 2020/5/12 33 4. T=T(y)_ —— 〉 μ = μ(y) 5. du/dy= τ / μ ，利用 1 步的 τ 和 6. 边界条件 y=0,u=0 求解上式 6. 直到 y=D, 检查 u=ue ？，否则，另假设 τ 进行 5 步，反复直到 u=ue 7. 用 6 步得到的 τ ，重新进行 2-7, 得到最终 收敛的 τ （大循环） 2020/5/12 34 时间相关有限差分法 二维非定常流 N-S 方程 2020/5/12 35 能量方程 2020/5/12 36 MacCormack 法 网格划分 流入边界条件， x=0, 用不可压流的解给 u,v,p,T 初始条件， t=0, u,v,p,T. 用不可压流解或 均匀值 2020/5/12 37 预测 - 校正法，以 x 轴动量方程为例 预测步 2020/5/12 38 校正步 用向后查分 2020/5/12 39 重 复 预 测 - 校 正 在 下 一 时 间 步 ， 得 到 结 果 2020/5/12 40 上述格式为 显示有限差分法，存在稳定性 限制 格式的稳定性 CFL 数，稳定性条件 Δ t x = Δ x/(u+a), Δ t y = Δ y/(v+a) 全场取最小值 2020/5/12 41 可压库埃特流结果 冷壁，速度和温度型 速度型影响小，温度型影响大 2020/5/12 42 绝热壁 速度型影响和温度型影响比冷壁大 2020/5/12 43 解析分析 一些量随温度的变化 Pr 近似常数 2020/5/12 1 第 16 章 一些特殊流动 N-S 方程精确解 库埃特流动和泊肃叶流动 2020/5/12 2 16.1 引言 平行流：流线 是直的，且互相平行 本章内容： N-S 方程精确解，包括 库埃特和泊肃叶流 基本概念：表面摩擦 (skin fraction) 热传导 (heat transfer) 恢复因子 (recovery factor) 雷诺比拟 (Renolds analogy) 2020/5/12 3 库埃特流动 2020/5/12 4 16.2 COUETTE FLOW( 库埃特流 ): GENERRAL DISCUSSION 边界条件 At y=D: u=u e , T=T e 。 上界流体和与动平板间的摩擦剪力 ? e 热传导 q e 平行流线（ Paralel streamline ） At y=0: u=0 , T=Tw 。 下界流体和与动平板间的摩擦剪力 ? w 热传导 q w 2020/5/12 5 y 方向的热流量 q y =-k dT/dy (15.2) 流动方向从温度高的壁面流向温度低的 冷壁：热量传输从流体到壁面 热壁：热量传输从壁面到流体 温度场： 1. 平板上下温度一般不同 , 产生 温度梯度 2. 动能由摩擦消耗变成内能 , 内能的变化由 温度升高显示出来（粘性耗散） 2020/5/12 6 无穷长平行流动特点 任何特性沿 x 方向不变 ,( 任何量如果变化 , 就会变到无穷大 或者无穷小 ) v=w=0 ? u/ ? x= ? T/ ? x= ? p/ ? x=0 考虑定常流，应用到粘性流动方程组： 2020/5/12 7 得到： x-momentum equation: ? / ? y( ?? u/ ? y)=0