非齐次边界条件下物体导热导温性能分析

1、非齐次边界条件下物体导热导温性能分析王钊摘要：对有不同热物性的物体，在首端恒定热流密度的第二类边界条件下，进行一维分析，物体的末端定义为第三类边界条件，即规定了边界上的物体与周围流体间的表面换热系数以及周围流体的温度。采 用Dofort-Frankel差分格式进行matlab编程求解，分析首端温度变化速度和最后时间几何体的温度分布，以及稳态下的几何体温度分布，比较不同的热物性对物体导热和导温性能的影响。【关键词】：微分方程；matlab；导热；导温.引言自然科学与工程技术中种种运动发展过程与 平衡现象各自遵守一定的规律。这些规律的定量表 述一般地呈现为关于含有未知函数及其导数的方 程，并且大多

2、是偏微分方程。同时大多没有解析解， 但是可以用用数值方法解出有意义的解。数值解法 需要对连续性的问题加以离散化，从而求出离散节 点的数值解。matlab是目前最流行、应用最广泛的科学和工 程计算软件，可以求出偏微分方程的数值解，来对 问题进行定量的分析。本文采用Dofort-Frankel差分 格式，进行matlab编程求解在第二类和第三类边界 条件共同作用下的物体内的温度分布，比较不同的 热物性对物体导热速度的影响。.一维导热微分方程与差分方程模型在直角坐标系中，对一维热传导现象的研究可 以得到二阶抛物型偏微分方程，本文分析这种类型 方程的混合问题，假设 入、P、c为常数后，方程可 以理解为

3、在区域RNoExMOSET内求函数u（x,t）满足方程:：2U2 二0,-二：二 x :二二,t 0 E%11)及初始条件和边界条件：Uy = ： x ,0 _x _1u 0,t t ,u 1,t =g2 t ,0 tT建立差分格式对x-t平面进行二维网格划分，取水平方向的h为x方向的空间步长，取竖直方向的 t为时间步x： = jh t = k长，使用二组平行直线：j和k ，可以将矩形区域 b,TL0,L分割成矩形网格。在上网格内点外，对 FE分别采用向前、向后及中心差金商公式，对 &采用二阶中心差商公式，一维热传 导方程 Dofort-Frankel格式表示为：.kt .k 1. k . k

4、 1- k 1. kjnakT ：一 2，h227h个 f (2)由Taylor公式可以知道，方程与一维热传导2 h2方程相容，截断误差分别为,1下面对定解条件进行离散化，对于初始条件及 第一类边界条件，可以直接得到:t0 =3初始条件为：tj tJ,边界条件为:U,j -u0, j 1,ju0, j - f1,j hUN,j -UN _1,j；2, jUN,j =f2,jh其中,1,j = 1 j . , 2,j = ,2 j . , f1 j .= f1,j , f2 j .= f2,j差分格式的求解令r =ah2为网格比，则公式(2)可以结合初始条件、边界条件得到差分方程组：k 1k 1

5、 kk(1+2r)tj + =(12r)tj +2rtj /+ 2rtj + 3 7式(3)是3层格式，需要事先知道前面2个时间层上的解，第1层上的解可以通过对初始条件 离散得到，第2层上的解可以使用的任意差分格式 得到，其余层就可以用这个格式求。.几何模型与物理参数物体的几何模型定义为长方体，端面为0.01mX0.01m,长度为0.1m。左0.01mX0.01m端面为第二类边界条件，恒定热流，定义为0.05W;右0.01mX0.01m端面为第三类边界条件，指定端面与流体 间的表面传热系数为5斗刀2 K流体的温度指定为与物体的初始温度一致为 0C。几何体的模型始下图1所示：与流体对流四个面更热

6、恒定热流 r=诙1 / ,1u X. 工 x J F标匕/本文选取四种不同的材料比较物体内温度分布，分别为纯铜、纯铝、纯铁和聚氯乙烯，假设四 种材料在实验的过程中，物性参数保持不变，物性如下表1所示:表1四种材料的物性参数C / (m k1P kg/m3cJ/(kg k )12 fs/m 纯铜39889303868660纯铝236271090210358纯铁80787045544760PC0.21300130084500004.纯铜和PC对比实验编写matlab主程序，格式如下所示:function U x t=temp_change_DF(uX,uT,M,N,C)函 数会调用定义的材料物性：

7、导热系数、比热和密度。uX=0.1是指一维空间长度，M是指长度方向的分 区间数，本文取10区间，共11点；uT是指计算时 间的长度，本文取100秒，N指时间方向的等分区间数，本文100区间；即matlab程序的计算次数为：1111 次计算纯铜的温度分布如下图 2所示，图中显示随时间的变化，几何体的一维温度分布:铜在一边固定热流一边自然对流下的图像U化 变 度 温 的 铜 维0.250.20.150.10.0501000.1图2纯铜材料随时间变化的温度分布图图1实验物体的相同几何模型分析数值解的数据，发现从第 11秒开始，几何体的末端才有温度变化。对导热系数最低的PC为00.04m。对比图2和图

8、3,可以看到纯铜和 PC材料进行同样的分析，温度分布如下图3所示:PC在一边固定热流一边自然对流下的图像时间变量t00空间变量x图3PC材料随时间变化的温度分布图U化变度温的CP维一在100秒时的温度分别为 0.2 C和28 C,在同一个图上面，差别过大使得纯铜的数据看不明显，所以后面分析就不放入PC的数据。5.纯铜、纯铝和纯铁对比实验下面分析纯铜、纯铝和纯铁在首端点的温升，如下图6所示:分析数值解的数据，温度传递的非常缓慢。查端点随时间变化的温度变化图如图4所示：2827.52726.52625.5250102030405060708090100时间变化t图4PC材料首端随时间变化的温度变化

9、图PC材料苜厢在恒M热流0.05W的作用下，温度随时间变化图图6各材料首端随时间变化的温度变化图从图中可以看到，由于材料的导热系数太小,所以首端点的温度相对要比较高才能传递恒定的热流。查看100秒时，几何体的温度分布图: 302520尸T化15变 度 温10PC材料100秒时温度分布图在0.04m 处，温度为0.006 C000.010.020.030.040.050.060.070.080.090.1长度变化x图5PC材料100秒时的温度分布图从图中并结合数值解的数据可以看到，铜从第 10秒开始，末端点的温度开始变化，同时首端点温 度的开始线性升高；铝也是从第10秒开始，末端点的温度开始变化

10、，同时首端点温度开始线性升高； 铁是从第25秒开始，末端点的温度开始变化，同 时首端点温度开始近线性升高。从图6中还可以看到铜温度线性升高部分的斜率最小，铝其次，铁最大；这说明铜材在首端点 也就是热流导入点的温度上升最缓慢，表明铜的导 热能力最好，铝其次，铁的导热能力最差。下面分析纯铜、纯铝和纯铁在第100秒的温度分布，如下图所示：分析图5的数值解数据，可以看到，由于 PC材料的导温系数过小，在 100秒时，温度传递长度炯材科温度分布 铝材料温度分布 铁材料温度分布0.010.020.030.040.050.060.070.080.090.100.010.020.030.040.050.060

11、.070.080.090.1长度变化x铜材料温度分布铝材料温度分布铁材料温度分布PC材料温度分布0.40.35 0.3 T 0.25 化 变 度 0.2 温0.15 0.1 0.05长度变化x图7各材料在100秒时随时间变化的温度分布图从图7中可以看到，铜材沿长度方向，温度分 布最均衡，分析数据温度在(0.1424 C0.2053C), 说明铜的导温能力最好；铝材沿长度方向温度分布 较陡，分析数据温度在(0.1956C0.3017C)即铝 材的导温能力比铜差；铁材沿长度方向温度分布最 陡，分析数据温度在(0.0698 C0.3676 C)即铁的 导温能力在这三者里最差。.数据验证计算与稳态温度

12、分布下面通过简略的数值验算数值解，分析纯铜、 纯铝和纯铁在50 c时的温度分布，如下图：0.35 00.010.020.030.040.050.060.070.080.090.1长度变化x图8各材料在50秒时随时间变化的温度分布图以铜材料为例，50秒内，首端传递的总热量 为：0.05WX50秒=2.5焦耳，假设热量没有通过自 然对流散热末端，则这些热量可以使铜的温度升高 0.0725 C,与图中铜的数据进行验证，可以认为数当稳态时，可以计算出物体的温度分布与首端 温度。稳态时，导入的热量等于对流换热量，所以 末端的温度为100 C,各种材料在稳态时的温度分 布如下图所示：100.8100.71

13、00.6T 100.5 布 分 度 100.4 温 I 态 100.3 稳100.2100.1100图9铜、铝和铁在稳态的温度分布图 350 300T 250 布 分 度 温 时200 态 稳15010000.010.020.030.040.050.060.070.080.090.1长度变化x图10PC在稳态的温度分布图即在稳态时，纯铜，纯铝和纯铁的首端温度均在100 c左右；而PC的首端温度为350 o即对导热系数从80398这个范围，最终稳态时的首末端温差相差不大，即导热能力相当。.总结由表1中热物性表，比较各材料的导热系数，我们可以看到，铜的导热能力最好，铝其次，铁的导热能力最差。所以铜材在首端点也就是热流导入点的温度上升最缓慢。PC的导热能力最差，所以PC首端点的温升最快，温度也最高。通过分析首端点值解正确的温度上升速度，可以相对比较不同材料的导热系由表1中热物性表，比较各材料的导温系数, 我们可以看到，铜的导温系数最高，铝其次，铁的导温系数最低。所以在100秒的时间之后，铜材内 的温度分布最平滑，铝其次，铁材内的温度分布最 陡。对于稳态计算的结果，通过图 9和图10可以看 到，最终影响稳态温度的因素是导热系数，密度和 比热参数对稳态温度分布没有影响，导热系数越小， 首末二端的温差越大，PC的温差达到250 C。对于散热器来说，选择导温系数大的材料，意 味着达到热平衡所