用c语言编程数值法求解传热,251论述相变传热问题一种数值解法.doc

251论述相变传热问题一种数值解法.doc

论述相变传热问题的一种数值解法

张丽蓉 解国珍

(北京建筑工程学院，环境与能源工程学院，北京，100044)

摘要：具有相变的传热过程，是一个物体三相间相互转换的过程，相变传热存在着一个相间移动的边界面,因而相变传热的计算比单纯热传导更为复杂。本文以水的凝固问题为例，采用编程的方法对水的凝固问题进行了数值求解，得出了在第一类边界条件下相变界面位置变化与界面移动速度变化以及两相中的温度变化等结果，并对相变过程进行分析，在其实际运用中具有重义。

关键词：相变 热传导 数值计算

Discussion the phase-change heat transfer problems in a numerical solution

Zhang Lirong，Xie Guozhen

(School of Environment and Energy Engineering，BUCEA,Beijing 100044)

Abstract: With a phase-change heat transfer process is an object between three phases conversion process,there is a phase-change heat transfer and phase moving boundary surface, and thus the calculation of phase-change heat transfer is more complicated than the simple heat conduction. In this paper,the solidification of water as an example,using programming methods of solidification of water conducted a numerical solution obtained in the first position of the boundary conditions of Phase Transformation interface changes and interfaces as well as the moving speed two-phase temperature and other results,and analysis of phase-change process in its practical application has a important meaning.

Key words: phase-change heat transfer numerical solution

前言

具有可移动边界的传热问题,出现在许多有实际意义的领域中,如金属的熔化与凝固,半导体材料的扩散,火箭重返大气层时表面保护材料的烧蚀,冷冻手术，天体物理,气象学,生态学和医学等,而在化学工业中，如汽化与冷凝,升华与凝华,蒸发,溶解以及某些具有相变的化学反应等[1]。正因为如此，文献[2]中指出“在这些领域中进行深入的研究仍然是十分必要的。”在国外这类问题对许多科学家和工程师都具有吸引力[3]-[4]。

对于相变问题，求解方法可以分为解析法(包括精确分析和近似分析)、数值方法(包括有限差分法、有限元法和有限容积法)。这种问题是1891年J.Stefan在研究地极冰块厚度时首先提出的[5]-[6],故又称Stefan问题，或自由边界问题,可动边界问题,其基本的和共同的特点是有一个分界面，把两个不同物性的区域分开,而这个分界面以时间的某个函数而运动,在这个分界面处有热量和质量的传递,有的还伴有化学反应。这类问题在数学上是用抛物型的偏微分方程描述。求解这类问题必须确定分界面运动的方式和速度,这样提出的问题是非线性的,尽管具有明显的实用价值，求解是相当困难的。现有的精确解很少，最基本的精确解法是由Franz Neumann于1860年研究一维半无限大物体相变问题时得到的，将Neumann法稍加推广可解决以下相变问题:(1)固、液相密度不同的相变问题;(2)定热流边界相变问题;(3)接触相变问题;(4)存在相变区的相变问题;(5)多相变问题，除了这些少数精确解之外，所有的课题必须用计算机数值求解或计算机模拟机求解[6]，本文给出了采用迭代方法通过计算机语言编程求解这类问题的具体解法和步骤，并以水的凝固为例，在计算机上得到了冷却速度和相界面的位置以及介质中的温度分布规律，因此对于研究边界条件和凝固时间对水的凝固的影响，找到凝固的理想边界条件和时间条件具有一定意义。

具有相变传热过程的数学描述

关于相变传热过程(即Stefan)的均匀温度，