有限容积法在冰水相变过程中的数值计算.pdfVIP

有限容积法在冰水相变过程中的数值计算 1 1 2 许彬 ，张敏 ，John C.Chai 1 南京理工大学动力工程学院，江苏南京（210094 ） 2 南洋理工大学机械与航天学院，新加坡（639798 ） 摘 要：用有限容积法，在结构化网格中求解冰水相变过程的温度场并寻找界面运动规律， 求解算例包括一维单区域问题和二维第一类边界问题。在一维问题中，将带有误差函数的精 确解同数值解进行了比较，并得到了满意的一致性。这些结果充分证明了该方法的正确性， 并且为工程应用提供了有益的参考。 关键词：相变过程，两相界面，有限容积法 1. 引言 相变过程是自然界和工程技术领域中的一种常见现象，冰层的融化大地的冰冻，铸件的 凝固食品的冷冻，乃至地球自身的形成等都是一些典型的实例。近年来，航天技术、能量贮 存技术, 以及生物工程技术的发展，促进了相变问题的研究，相变传热技术也因此而得到越 [1-6] 来越广泛的应用 。 对于相变过程的研究，P. Rath[7]等人对工业界中电子仪器精刻过程中的化学蚀刻过程进 行了类似的研究。本文采用有限容积法，在结构化网格中，对冰水相变过程进行数值计算和 研究分析，通过两个算例，寻找相变过程的基本规律，并得到有参考价值的结果。 2. 基本控制方程和离散方法 [6] 非稳态的导热控制方程为 ： ∂T ρcp div(kgradT ) +S (2.1) ∂t 对于二维问题，方程（2.1 ）可以简化为如下形式： ⎛ ⎞ ∂T ∂ ∂⎛ T ⎞ ∂ ∂T ⎜ ⎟ ρcp ⎜k ⎟+ k ＋S (2.2) ⎜ ⎟ t x x ∂y ∂y ∂ ∂ ∂⎝ ⎠ ⎝ ⎠ 方程（2.2 ）在全隐式时间积分方案下的离散方程为： a T a T +a T +a T +a T +b (2.3) p p E E W W N N S S 相邻节点包括 W、E 、S 和 N ，其中: ke ∆y kw ∆y kn ∆x ks ∆x a a a a E ( ) , W ( ) , N ( ) , S ( ) (2.4) δx δx δx δx e w n s a a +a +a +a +a0 −S ∆x∆y (2.5) p E