相变传热.ppt

相变传热 焦冬生 热科学与能源工程系 内容 引言 精确解（诺曼（Neumann）解） 分析方法—移动热源法 近似方法 数值方法 一、引言 相变过程其实就是传热传质过程； 这类传热现象的基本特征 由边界（固液交界面）移动引起的非线性（nonlinearity）化，使得此类问题变得更复杂，并且每一个问题均有其独特性。 引起数学处理较为困难的其他因素： 在相变前、进行中和后，其物理性质依赖于温度，而温度分布是三维的和瞬时的变化。 有时，溶解和凝结时发生的复杂又令人困惑的现象，使得传统的分析方法无法解决。 热传导方程 静止的均匀物体内含有热源的各向同性物体的热传导方程 直角坐标系 柱坐标和球坐标 边界条件和初始条件 第一类边界条件：边界上温度给定 第二类边界条件：边界边上温度的法向导数给定 第三类边界条件：边界边上温度和法向导数的线性组合给定 解析解——分离变量法 平板 解析解——分离变量法 此方程中，左边只是空间变量x的函数，右边只是时间变量t函数,要使等式成立，只有两边都等于同一个常数 时间变量方程 空间变量函数满足微分方程 特征方程的根为一对共轭复数根 通解为 温度的完全解可由上述分离方程的基本解按线性迭加原理构成，其形式为 这个解既满足热传导问题的微分方程，又满足边界条件，但是它并不一定满足初始条件。因此，将初始条件应用于上式可得 未知系数可根根据下述特征函数的正交性来确定： 我们用算子 对F(x)的两边进行运算，再根据正交性，可得 相变传热现象 物理现象（连续介质） 有单一相变温度和明确界面。 相变有一个温度范围，存在两相区。 相变传热模型 温度法 以温度为唯一的因变量，分别在固相和液相区建立能量守恒方程。 焓法 焓和温度共同作为因变量，无需分区建立控制方程 温度法 控制方程 界面耦合条件 以下情况不考虑速度场： 忽略密度变化的影响，液相内只有导热； 密度不同，但液相一直处于相变温度。 焓法对材料的密度和相变特性没作特殊假设。 积分形式的控制方程 把原来在两个活动区域及固液界面成立的方程组转换为在一个固定区域内成立的方程，无需跟踪界面。便于数值计算。 界面移动规律 半无限大物体在某个时间t ，固体层厚度为x(t)，暴露的表面保持在温度Ts，Ts比相变温度低。在交界面处释放出来的热量必须由热传导经过固体导出，假定液相中没有温度梯度。 被释放出来的热量经过固体层导出 解出增长层的表达式 二、诺曼（Neumann）解:半无限大物体相变问题 斯蒂芬（Stefan）问题 考虑一个分布在正x区域处于均匀温度Ti的液体， Ti高于其固态的熔解温度Tf。 在t=0时，处于x=0的液态表面突然降到Tw（ Tf ），并保持下去。 这样，凝固从x=0的表面开始，固液交界面则沿x正方向移动。 固液两相密度相同 这类问题通常称为斯蒂芬问题。 斯蒂芬（Stefan）问题的数学描述固相内数学描述液相内数学描述 假定固相和液相内温度分布的解分别为 解（13）满足微分方程、边界条件和初始条件（10-11） 。 带入界面条件（12）可得 为常数，得到系数A，B的表达式 通过固定边界位置，诺曼得到固相区和液相区的温度分布如下（具体方法参见M.N.奥齐西克《热传导》）： 表面热流 总传热量 误差函数 诺曼解的应用 获得了一些应用； 但求解问题的范围很小。 固液相密度不同的材料，固定边界将导致液相的整体运动，自由界面将导致两个界面的移动。 凝固过程（ ）液相在x方向的整体运动速度 液相的能量方程 固相的解同前 液相解 定热流边界相变 接触相变——界面温度恒定，溶解过程中溶解物被完全排除 存在相变区 假定两相区固相成分与距离成线性关系 相变热处理成内热源 多相变问题 半无限大液体初始温度为T，t=0后表面瞬间冷却并保持较低温度Tw，此时液体发生凝固。假设介质有两个相变温度Tm1，Tm2，且TiTm1Tm2Tw，相变热分别为hm1，hm2。固液密度相同 构造温度函数 带入边界温度条件 相界面的位置 利用相界面的能量守恒条件 三、分析方法—移动热源法 （Moving heat source method） 凝固（或熔解）过程中热量的释放（或吸收）可视为在固—液交界面处有一移动的平面热源（或汇）。从这一认识出发，可以把瞬态的相变问题在形式上考虑为具有移动平面热源的瞬态热传导问题。 用这一方法求解瞬态相变问题最早的是Lightfoot。 N.M.H.Lightfoot, The Effect of Lat