数值传热学部分习题答案.doc

? ? ( ρuw ? ( ρvw ? ? ?w ? ? ? ?w ? ? + = ?η ? + ?η ? ? ?y ?x ? ?x ? ?y ? ? ?x ? ?y ? ? 2 2 ?? ψ + ? ψ ? w = 0 ? ?x 2 ?y 2 ? 由假设的流函数ψ 求解涡量控制方程， 然后由求解获得的涡量 w 通过流函数控制方程更新 0 得到ψ ，重复直到收敛。因为ψ , w 有相应控制方程且互相耦合，所以上述步骤可以得到合 * 理的ψ ，最终可以得到速度 u , v 。在涡量－流函数法中，不存在速度及压力关系隐含在连 续方程的问题，所以在获得收敛的流函数后，便可以利用 Poisson 方程： ?? ? 2ψ ?2 p ?2 p ?? ρ 2 + = ? ?x 2 ?x 2 ?y 2 ? ?? 求解 p 分布； ?? ? 2ψ ? ? 2 ?? ?? ?y ? ? ? 2ψ ? ? ? ? ? ??? ? ? ?x?y ? 2 ? ? ? ? 8－3 坐标右图所示，将ψ i , 2 ,ψ i , 3 分别在 (i , n 点 Taylor 展开： ψ i , 2 = ψ i ,1 + ψ i ,3 ?ψ ?y δy + i ,1 ? 2ψ ?y 2 (δy 2 ? 3ψ + 3 2! ?y i ,1 (δy 2 + O (δy 3 3 ! i ,1 ( 2δy 3 + O (δy 4 3! i ,1 （1） ?ψ = ψ i ,1 + ?y ? 2ψ (2δy + 2 ?y i ,1 ? 2ψ ?y 2 ( 2δy 2 ? 3ψ + 3 2! ?y i ,1 ?u ?y （2） 因为 ?ψ ?y = u i ,1 = 0 ， i ,1 = i ,1 = wi ,1 ，所以 8 × (1 ? ( 2 得： i ,1 8ψ i , 2 ? ψ i ,3 = 7ψ i , 2 + 2(δy 2 wi ,1 + O(δy 4 整理即可得： wi ,1 = ? 7ψ i ,1 + 8ψ i , 2 ? ψ i ,3 2δy 2 + O(δy 2 8－4 将流函数ψ i , 2 在 (i,1 点展开，有： ψ i , 2 = ψ i ,1 + ?ψ ?y δy i ,1 2 + ? 2ψ ?y 2 (δy / 2 2 + O (δy 3 2 ! i ,1 26 因为 ?ψ ?y = u i ,1 = 0 ， i ,1 ? 2ψ ?y 2 = i ,1 ?u ?y = wi ,1 ，所以上式整理为： i ,1 ψ i , 2 = ψ i ,1 + 进而可得 B 节点法中 Thom 公式的表达： (δy 2 wi ,1 + O(δy 3 8 wi ,1 = 8(ψ i , 2 ? ψ i ,1 δy 2 + O(δy 将流函数ψ i , 2 在 (i,1 点展开成具有三阶截差的 Taylor 公式可得： ψ i,2 且有： ?ψ = ψ i ,1 + ?y δy i ,1 ? 2ψ + 2 ?y 2 (δy / 2 2 ? 3ψ + 3 2! ?y i ,1 (δy / 2 3 + O (δy 4 3! i ,1 ? 3ψ ?y 3 = i ,1 ? ? ? 2ψ ? 2 ?y ? ? ?y ? ?w wi , 2 ? wi ,1 ? ? = ?y = (δy / 2 + O (δy ? 将流函数的各阶导数表达式代入 Taylor 公式得： ψ i , 2 = ψ i ,1 + ? wi , 2 ? wi ,1 ? (δy 3 (δy 2 4 ? wi ,1 + ? O ( δ y + ? (δy / 2 ? 48 + O(δy 8 ? ? 12(ψ i , 2 ? ψ i ,1 (δy 2 整理可得 B 节点法中的 Woods 公式： wi ,1 = ? 1 wi , 2 + O(δy 2 2 8－5 对直角坐标系中，对流换热的有效压力为： p eff = p + ρ c gy sin θ + ρ c gx cos θ 对柱坐标系中，对流换热的有效压力为： p eff = p + ρ c gx 对极坐标系中，对流换热的有效压力为： p eff = p + ρ c gr cos θ 27 9－1 由势流的伯努利方程有： pe = p + 1 2 ρu 2 由工作条件可得空气密度 ρ = 1.205 ，所以可得： p = pe ? 由 u 2 1 2 1 ρu = 100000 ? × 1.205 × 50 2 = 98493.75 Pa 2 2 / u = 5% ，可得： u 2 = (5% × u 2 = (0.05 * 50 2 = 6.25 所以： 有效压力： p t = ρu ; 2 = 1.205 * 6.25 = 7.