[图文]CFD数值模拟原理-6-.doc

' ' aeue ? ? anb unb ? ? Pp' ? PE' ? Ae （修正方程） 邻近点上速度修正 压力修正（直接的力） ' ? a u 与点方程，直接解困难 近似取 nb nb ? 0 ? ? ? ' ' ' a u ? P ? P ? ?e e p E ? Ae ? Ae ? ' ' ' u ? P ? P ? p E? ? e a ? e 速度修正方程 ? ?v ' ? d ? P ' ? P ' ? d ? An n n p N n ? an ? ? 修正后的速度场能否满足连续性方程——检验 ? p’的 方程： u= u*+ u′,v= v* + v′代入连续性方程： ? n n ?1 ?p ? ?p ?x?y ? ? ?y ? ? ?x ? 0 ? ? u ?e ? ? ? u ? w ? ? ? v ?n ? ? ? v ? s ? ? ? ? ? ?t aP p ? aE p ? aW p ? aN p ? aS p ? b ' P ' E ' W ' N ' S ?aE ? ?e de ?y ? ?aN ? ?n d n ?y aw ? ? w d w ?y aS ? ? s d s ?y aP ? aE ? aW ? aN ? aS n n ?1 ?p ? ?p b? ? ?? ?u* ? ? ? ?u* ? ? ?y ? ?? ? v* ? ? ? ? v* ? ? ?x w e? s n? ? ? ?t 边界上p’＝0 质量守恒判据 §7－3 SIMPLE 算法 假定初始速度u0，v0 动量离散方程中各系数 假定一个压力初p* 求解动量方程 u*,v* 满足连方 输出 否 求压力修正方程 p‘ Semi Implicit Method for pressure Linked Equations 求解压力耦合方程的 半隐方法 p= p*+ p′, u= u*+ u′,v= v* + v′ 求解其它耦合问题 SIMPLE虽然开始是用于不可压问题，目前已成功应用于可 压缩流问题及有限元法。 N vn 例7－1 已知:pw=60 pE=40 ue=20 vn=7 uw ●w ●p ue ●E vs S uw=0.7(pw -pP vs =0.6(ps-pP 试求：pP uw vs 值（SIMPLE方法）