钢球的非稳态冷却过程的数值模拟.pdf

钢球的非稳态冷却过程的数值模拟 魏晓 （山东豪迈化工技术） 摘要：非稳态导热是当设备处于变动工作条件下时，其内部温度场随时间变动， 处于不稳定状态。本文对一钢球的非稳态冷却过程，用近似分析法和数值解法进 行求解对比，发现数值计算结果更精确，能更真实地反应物理过程的实际情况。 1.导热问题概述 导热是3种热量传递方式(导热、对流、辐射)中的一种，工程中的导热问题 包括两种情况：稳态导热和非稳态导热。稳态导热是整个过程中物体各点的温度 不随时间变化，始终保持一个温度；非稳态导热是当设备处于变动工作条件下时， 其内部温度场随时间变动，处于不稳定状态。自然界和工程上许多导热过程为非 稳态t f()，例如：冶金、热处理与热加工中工件被加热或冷却，锅炉、内燃机 等装置起动、停机、变工况，自然环境温度，供暖或停暖过程中墙体与室内空气 温度等。 傅里叶定律揭示导热问题的基本规律：即在导热现象中，单位时间内通过给 定截面的热量，正比例于垂直该截面方向上的温度变化率和截面面积，而热量传 递的方向与温度升高的方向相反。由傅里叶定律并结合能量守恒定律，建立了导 热微分方程。等号左侧的项为非稳态项，右侧中间3项为导热项，最右侧的一项 为源项。 t t  t t  t t  t  c          x  x y  y z  z 传统方法求解导热问题实际上就是对导热微分方程在定解条件下的积分求 解，这种方法获得的解称为分析解。但工程技术中遇到的许多几何形状或者边界 条件复杂的导热问题，由于数学上的困难而无法得到其分析解。近几十年来，随 着计算机技术的迅速发展，对物理问题进行离散求解的数值方法发展得十分迅 速，并得到广泛应用。 数值解法的基本思想是：把原来在时间、空间坐标系中连续的物理量的场， 用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关于这些 值的代数方程来获得离散点上被求物理量的值，这些离散点上被求物理量值的集 合称为该物理量的数值解。 2 问题描述 一个直径为5cm 的钢球，放入加热设备中加热至723K （450℃），假设此时 整个钢球温度均匀，然后突然被置于温度为303K （30℃）的空气中。 2 已知钢球表面与周围环境的表面传热系数为24W/(m .K)；钢球的比热容为 3 480J/(kg.K)，密度为7753kg/m ，导热系数为33W/(m.K)。求钢球表面和中心冷 却到573K （300℃）所需要的时间。 3 求解 对于非稳态传热问题，通常用分析解法、近似分析法、数值解法来求解。其 中分析解法包括分离变量法、积分变换、拉普拉斯变换；近似分析法包括集总参 数法和积分法；数值解法包括有限差分法、蒙特卡洛法、有限元法、分子动力学 模拟。以下分别采用集总参数法和有限元两种方法对该问题进行求解并对比。 3.1 近似分析法——集总参数法 集总参数法是当导热体内部热阻忽略不计时，即Bi→0时研究非稳态导热的 Bi 一种方法。毕渥数 数表示物体内部导热热阻和外部对流热阻的比值，其表达 h 式为Bi 。其实用判别条件是Bi＜0.1，这一判别式产生的依据是使整个导  热体中温度的不均匀性在5％以内。此时，温度仅为时间的函数，而与空间坐标 没有关系。 对于该钢球非稳态冷却问题，首先判断能否用集总参数法求解： h(V/A) 毕渥数为B 0.00606 0.033，可以采用集总参数法。因  2 2 hA 24W/