非稳态传热大作业.pdfVIP

传递过程原理大作业 宋劲 化13 班 2011011842 问题描述 如下图所示，一高温管式反应器由内径5cm，壁厚2.5cm 的陶瓷材料制成，为实现高温 条件下的反应，需要对反应器进行保温，保温层厚度为5cm，陶瓷和保温层的物性为： 陶瓷：ρ=2600 kg/m3 c =1150 J/(kg ·K ) k =3.0 W/(m ·K) 1 1 1 保温层：ρ=600 kg/m3 c =200 J/(kg ·K ) k =0.2 W/(m ·K) 2 2 2 反应器初始温度是300 K，内表面突然被暴露在1500K 的反应物流中，管壁和反应物流 间的对流换热系数为500 W/(m2 ·K) 。请给出保温层内表面温度随时间的变化规律，以及保 温层内温度分布，讨论保温层厚度的影响。（保温层外表面视为绝热） 图1 题图 题目分析 该问题是柱坐标系下的一维非稳态导热问题，比较特殊之处在于该体系分陶瓷层、保温 层两部分，共有内、中、外三个边界条件。 建立柱坐标系并传热微分方程进行化简，不难得到如下关系式： 2  1    (  2 ) t r r r T T 其中，无量纲温度=  T T 0  (T 1500K,T 300K )  0 根据题目给出的参数可校验发现该圆柱是厚壁物体，则必须解以上偏微分方程，考虑用 差分方程的方法给出数值解。 为保证结果的稳定性不受步长影响，选择用时间向后、空间一阶向后、空间二阶中心差 分方程。 n n1 n n1       i i O(t)  i i t i t t n n n n1       i i1 O(r)  i i r i r r 2 n n n n n n   2  2   i1 i1 i 2 i1     O(r)  i1 i r2 i (r)2 (r)2 以上上标n表示时间，下标i表示半径。 t、r分别为时间、半径的步长。 由此对于给定时间点，可以得到一个大规模方程组，自变量为各个r 位置的无量纲温度。 以上即可解决三个非边界点以外的所有差分点的列式，下面再给出边界点的列式。 1）r=0.0025m 处，陶瓷材料与热流界面处 在该位置可通过热平衡法得到其表达式：   C 2rr  hk p t 1 r 进而做差分近似处理得 n n-1 n n     1