第三章第三节 一维非稳态导热的分析解
传热学传热学--第三章 第三节 一维非稳态导热问题
传热学--第三章第三节一维非稳态导热问题 §3 — 3 一维非稳态导热的分析解 本节介绍第三类边界条件下：无限大平板、无限长圆柱、球的分析解及应用。如何理解无限大物体，如：当一块平板的长度、宽度>> 厚度时，平板的长度和宽度的边缘向四周的散热对平板内的温度分布影响很少，以至于可以把平板内各点的温度看作仅是厚度的函数时，该平板就是一块“无限大”平板。若平板的长度、宽度、厚度相差较小，但平板四周绝热良好，则热量交换仅发生在平板两侧面，从传热的角度分析，可简化成一维导热问题。 一、无限大平板的分析解 已知：厚度的无限大平板，初温t0，初始瞬间将其放于温度为的流体中，而且> t0，流体与板面间的表面传热系数为一常数。 试确定在非稳态过程中板内的温度分布。 解：如图3-5 所示，平板两面对称受热，所以其内温度分布以其中心截面为对称面。对 于x 0 的半块平板，其导热微分方程：(0 （边界条件） （边界条件） 对偏微分方程分离变量求解得： （3-10 ） 其中离散值是下列超越方程的根，称为特征值。 其中Bi 是以特征长度为的毕渥数。 由此可见：平板中的无量纲过余温度与三个无量纲数有关：以平板厚度一半为特 征长度的傅立叶数、毕渥数及即：（3-12) 二、非稳态导热的正规状况阶段 1 、平板中任一点的过余温度与平板中心的过余温度的关系 前述得到的分析解是一个无穷级数，计算工作量大，但对比计算表明，当Fo>0.2 时，采用该级数的第一项与采用完整的级数计算平板中心温度的误差小于1% ，因此，当Fo>0.2 时，采用以下简化结果：（3-13 ） 其中特征值之值与Bi 有关。 由上式（3-13 ）可知：Fo>0.2 以后平板中任一点的过余温度(x ，τ) 与平板中心的过余温度(0 ，τ)=（τ ）之比为：（3-14 ） 此式反映了非稳态导热过程中一种很重要的物理现象：即当Fo>0.2 以后，虽然(x ，τ) 与（τ ）各自均与τ 有关，但其比值则与τ 无关，而仅取决于几何位置（x/ ）及边界条件（Bi ）。也就是说，初始条件的影响已经消失，无论初始条件分布如何，只要 一维非稳态导热问题的数 值解 Prepared on 22 November 2020 计算传热学程序报告 题目：一维非稳态导热问题的数值解 姓名： 学号： 学院：能源与动力工程学院 专业：工程热物理 日期：2014年5月25日 一维非稳态导热问题数值解 求解下列热传导问题： 1.方程离散化 对方程进行控制体积分得到： 非稳态项：选取T 随x 阶梯式变化，有 扩散项：选取一阶导数随时间做显示变化，有 进一步取T 随x 呈分段线性变化，有 e P E e x T T x T )()(δ-=??，w W P w x T T x T )()(δ-=?? 整理可以得到总的离散方程为： 2.计算空间和时间步长 取空间步长为： h=L/N 网格Fourier 数为： 2 20x t x t F ??=??=α（小于时稳定） 时间步长为： 3.建立温度矩阵与边界条件 T=ones(N+1,M+1) T(:,1)=Ti(初始条件温度都为0） T(1,:)=To(边界条件x=0处温度为1） T(N+1,:)=Te(边界条件x=L 处温度为0） 4.差分法求解温度 由离散方程可得到： 转化为相应的温度矩阵形式： 5.输入界面 考虑到方程的变量，采用inputdlg函数设置5个输入变量，对这5个变量设置了默认值，如图1所示。在计算中可以改变不同的数值，得到不同的结果，特别注意稳定条件的临界值是。根据设置的默认值，得到的计算结果如图2所示。 图1matlab变量输入界面 图2默认值的计算结果 6.结果分析 根据上面的分析，给出了程序的输入界面，以及默认值状态下的数值解。可以通过改变不同的输入值，得到需要的分析结果，总结出了下面4点结论：（1）取F0=，得到一维非稳态导热结果如下图所示 图2F0=时一维非稳态导热 从图中可以看出，对于长度L=1的细杆，初始时刻t=0时温度为0，边界条件x=0时，T=1,边界条件x=1时，T=0。随着时间的增加，温度从x=0通过导热的形式传递到x=1，不同时刻不同位置杆的温度都不同，并且随着时间的增加，杆的温度也逐渐增加。 （2）取F0=，可以得到不同位置的温度响应曲线，如下图所示 图3F0=时不同x位置处的温度响应 图中红色曲线代表x=位置的温度瞬态响应，黑色曲线代表x=位置的温度瞬态响应，蓝色曲线代表x=位置的温度瞬态响应。从图中可以看出，随着x的增加，曲线与x轴的交点值越大，温度开始传递到该位置的所需的时间越长。 一维非稳态导热的数值计算 一、实验名称 一维非稳态导热的数值计算 二、实验内容 一块无限大平板（如图3所示），其一半厚度为L=0.1m ，初始温度T 0=1000℃，突然将其插入温度T ∞=20℃的流体介质中。平板的导热系数λ=34.89W/m ℃，密度ρ=7800 kg/m 3，比热c=0.712310 J/kg ℃，平板与介质的对流换热系数为h=233W/m 2.℃，求平板内各点的温度分布。 三、实验编程 #include for(k=0; k<=1000; k++) {for(i=1; i<=L-1; i++) for(j=1; j<=L-1; j++) {T[i][j] = T[i][j] + (a/4)*(T[i+1][j] + T[i][j+1] + T[i-1][j] + T[i][j-1] - 4*T[i][j]); } } printf(" a = %lf\n", a); printf("T[x][y] = ...\n"); for(i=0; i<=L; i++) for(j=0; j<=L; j++) {printf("%.1lf\t", T[i][j]); if(j == L) putchar(10); } return 0; } 四、运行结果 第三章非稳态导热分析解法 本章主要要求： １、重点内容： ① 非稳态导热的基本概念及特点； ② 集总参数法的基本原理及应用； ③ 一维及二维非稳态导热问题。 2 、掌握内容： ① 确定瞬时温度场的方法； ② 确定在一时间间隔内物体所传导热量的计算方法。 3 、了解内容：无限大物体非稳态导热的基本特点。 许多工程问题需要确定：物体内部温度场随时间的变化，或确定其内部温度达某一极限值所需的时间。如：机器启动、变动工况时，急剧的温度变化会使部件因热应力而破坏。因此，应确定其内部的瞬时温度场。钢制工件的热处理是一个典型的非稳态导热过程，掌握工件中温度变化的速率是控制工件热处理质量的重要因素；金属在加热炉内加热时，要确定它在炉内停留的时间，以保证达到规定的中心温度。 §3—1 非稳态导热的基本概念 一、非稳态导热 1 、定义：物体的温度随时间而变化的导热过程称非稳态导热。 2 、分类：根据物体内温度随时间而变化的特征不同分： 1 2 ）物体的温度随时间而作周期性变化 如图 3-1 所示，设一平壁，初值温度 t 0 ，令其左侧的表面温 度突然升高到 并保持不变，而右侧仍与温度为 的空气接触，试分 析物体的温度场的变化过程。 首先，物体与高温表面靠近部分的温度很快上升，而其余部分仍 保持原来的 t 0 。 如图中曲线 HBD ，随时间的推移，由于物体导热温度变化波及范 围扩大，到某一时间后，右侧表面温度也逐渐升高，如图中曲线 HCD 、 HE 、 HF 。 最后，当时间达到一定值后，温度分布保持恒定，如图中曲线 HG （若 λ=const ，则 HG 是直线）。 由此可见，上述非稳态导热过程中，存在着右侧面参与换热与不参 与换热的两个不同阶段。 （ 1 ）第一阶段（右侧面不参与换热） 温度分布显现出部分为非稳态导热规律控制区和部分为初始温度区的混合分布，即：在此阶段物体温度分布受 t 分布的影响较大，此阶段称非正规状况阶段。 （ 2 ）第二阶段，（右侧面参与换热） 当右侧面参与换热以后，物体中的温度分布不受 to 影响，主要取决于边界条件及物性，此时，非稳态导热过程进入到正规状况阶段。正规状况阶段的温度变化规律是本章讨论的重点。 2 ）二类非稳态导热的区别：前者存在着有区别的两个不同阶段，而后者不存在。 3 、特点； 非稳态导热过程中，在与热流量方向相垂直的不同截面上热流量不相等，这是非稳态导热区别于稳态导热的一个特点。 计算传热学程序报告 题目：一维非稳态导热问题的数值解 姓名： 学号：学院：能源与动力工程学院专业：工程热物理日期：2014年5月25 日 一维非稳态导热问题数值解 求解下列热传导问题: 1. 方程离散化 对方程进行控制体积分得到: 非稳态项：选取T随x阶梯式变化，有扩散项：选取一阶导数随时间做显示变化，有进一步取T随x呈分段线性变化，有 T ()e x T E T P (T \ T P T W ，( )w (x)x ( X) 整理可以得到总的离散方程为: 2. 计算空间和时间步长 取空间步长为: h=L/N 网格Fourier数为： F。一^ —V (小于时稳定) x x 时间步长为： 3. 建立温度矩阵与边界条件 T=o nes(N+1,M+1) T(:,1)=Ti(初始条件温度都为0) T(1,:)=To(边界条件x=0处温度为1) T(N+1,:)=Te(边界条件x=L处温度为0) 4. 差分法求解温度 由离散方程可得到： 转化为相应的温度矩阵形式： 5. 输入界面 考虑到方程的变量，采用inputdlg函数设置5个输入变量，对这5个变量设置了默认值，如图1所示。在计算中可以改变不同的数值，得到不同的结果，特别注意稳定条件的临界值是。根据设置的默认值，得到的计算结果如图2所示。图1matlab变量输入界面 图2 默认值的计算结果 6. 结果分析根据上面的分析，给出了程序的输入界面，以及默认值状态下的数值解。可以通过改变不同的输入值，得到需要的分析结果，总结出了下面4 点结论： （1）取F o=，得到一维非稳态导热结果如下图所示 图2F。二时一维非稳态导热 从图中可以看出，对于长度L=1 的细杆，初始时刻t=0 时温度为0，边界条件 x=0时，T=1,边界条件x=1时，T=0。随着时间的增加，温度从x=0通过导热的形式传递到x=1，不同时刻不同位置杆的温度都不同，并且随着时间的增加，杆的温度也逐渐增加。 （2）取F o=，可以得到不同位置的温度响应曲线，如下图所示 图3F o=时不同X位置处的温度响应 图中红色曲线代表x=位置的温度瞬态响应，黑色曲线代表x=位置的温度瞬态响应，蓝色曲线代表X=位置的温度瞬态响应。从图中可以看出，随着X的增加，曲线与X 轴的交点值越大，温度开始传递到该位置的所需的时间越长。随着x 的增加，温度响应曲线的变化速率越慢，最终的达到的温度也越低。 （3）取F o=，得到不同位置的温度响应曲线如下图所示 图4F o=时不同X位置处的温度响应 图中三条曲线分别是X=，X=，x=位置的温度瞬态响应。与图3的F o=进行对比，两种情况下的F o值不同，F o值越大表明热扩散系数的值越大。从图中可以 看出热扩散系数对于导热的影响，尸0=时，与F o=相比较，各位置开始响应时所需的时间较长，而且各位置响应曲线的变化速率较小，最终的达到的温度也较低，说明了热扩散系数越小，热传导越慢，传递效率越低。 （4）取F o=，得到非稳定的数值解如图所示 图5F o二时一维非稳态导热 图6F o=时不同X位置处的温度响应 从图中可以看出，对于显示格式的离散方程，并不是所有的F o值都能得到有意义的解，必须要求F o＜时才能得到稳定的数值解，当F o＞时，会出现物理上不真实的解。 传热学C 程序源 二维稳态导热的数值计算 2.1物理问题 一矩形区域，其边长L=W=1，假设区域内无内热源，导热系数为常数，三个边温度为T1=0，一个边温度为T2=1，求该矩形区域内的温度分布。 2.2 数学描述 对上述问题的微分方程及其边界条件为：2222T T 0x y ??+=?? x=0，T=T 1=0 x=1，T=T 1=0 y=0，T=T 1=0 y=1，T=T 2=1 该问题的解析解：112121(1)sin n n n sh y T T n L x n T T n L sh W L ππππ∞=??? ?---????=? ?-????? ??? ∑ 2.3数值离散 2.3.1区域离散 区域离散x 方向总节点数为N ，y 方向总节点数为M ，区域内任一节点用I,j 表示。 2.3.2方程的离散 对于图中所有的内部节点方程可写为：2222,,0i j i j t t x y ??????+= ? ??????? 用I,j 节点的二阶中心差分代替上式中的二阶导数，得： +1,,-1,,+1,,-1222+2+0i j i j i j i j i j i j T T T T T T x y --+= 上式整理成迭代形式：()()22 ,1,-1,,1,-12222+2() 2()i j i j i j i j i j y x T T T T T x y x y ++=++++ (i=2,3……,N-1)，(j=2,3……,M-1) 补充四个边界上的第一类边界条件得：1,1j T T = (j=1,2,3……,M) ,1N j T T = (j=1,2,3……,M) ,1i j T T = (i=1,2,3……,N) !*************************************************************** ! 一维非稳态导热通用程序(不变部分) ! This is a general purpose program to solve 1-D diffusion ! problem in the form of: ! ρcdt/dz=1/a(x)d/dx(a(x)Γdt/dx)+s !******************2003.7 revised******************************** !...................Define Variables......................... MODULE VARIABLES INTEGER,PARAMETER::L1=130 REAL,DIMENSION(L1):: X,XF,XM,XP,R,RF,AP REAL,DIMENSION(L1)::AE,AW,CN,T,TA REAL,DIMENSION(L1)::TG,GM,RC INTEGER:: K=1,KM=1,KP=1,OM=1 INTEGER:: JB,JE,KE,KI,KF,KN,KR,KT,LS,MD,M1,M2,NF REAL:: AEC,AI,BE,BI,DF,DS,DT,EP,EX,PW,TU,TM,XE,XI REAL:: A1,A2,T2,TC,SC,SP,RO,TE,DN,LM END !...............................................Main Program........................................... PROGRAM MAIN USE V ARIABLES IMPLICIT NONE INTEGER I OPEN(1,FILE="q.dat") OPEN(2,FILE="temp.dat") NF=1！（求解变量指标） KN=1！（非稳态问题输出局部变量，输出一次加1） TU=0！（当前时间） 150 KT=1！（非线性问题迭代次数） CALL Speci !First to specify the problem CALL Grid!Set up grid points 200 CALL Difsor!Specify the diff-coeff and source term 220 CALL InterOutput!Output intermediate results CALL Coeff!Set up coefficients of discretization equation CALL TDMA!Solve the algebraic equation by TDMA IF(LS.EQ.2.OR.LS.EQ.4) THEN ！（对非线性问题） IF(DF.GT.EP) THEN！（如果最大偏差大于允许值） DO I=1,M1 TA(I)=TA(I)+OM*(T(I)-TA(I))！（采用亚松弛方式将当前解付给上一次迭代值） END DO DF=0！（最大偏差置零） 一维非稳态热传导热源反问题研究 摘要 本文是关于热传导的正反问题的研究，即利用偏微分方程中典型热传导方程 t时刻温度分布与热源位置。 求解含有内热源的金属细杆 本文从解偏微分方程出发，由已知条件最终得出温度分布函数及热源位置函数并建立了两个数学模型。 模型一：利用偏微分方程及初始温度分布函数建立了一段时间后的温度分布与热源强度、位置之间的数学模型，最终解出一段时间后长杆上的温度分布。 模型二：通过一类抛物型偏微分方程模型，解决已知初始温度分布函数、一段时候后的温度分布函数及热源强度的确定热源位置和中间任意时刻的温度分布函数。 u x t,即t时刻的温度根据模型一建立偏微分方程组，用分离变量法求解(,) 分布函数，并通过Matlab中的PDE(偏微分方程)工具箱求解偏微分方程组，且使解可视化。 u x T,结合抛物型方程，运用根据模型二依然建立偏微分方程组，通过测得(,) 离散正则法，确定热源位置，并通过论证说明问题的唯一性和确定性，给出反问题的数值解法。最后再简单介绍差分法解决热传导在非稳态导热问题中的应用。 最后是结论部分，主要总结本文的结果并提出一些尚待进一步研究的问题，以及研究该反问题的应用前景。 相同t不同x的温度变化曲线相同x不同t的温度变化曲线 一维非稳态热传导热源反问题研究 一、问题的提出 在金属细秆的传热过程中，温度差是导致其发生必要条件，有无热源决定传导效率的高低。从一维非稳态传导问题的数学模型和初始条件出发，经过对有内热源问题的进一步分析，在初始温度分布已知的情况下，对分布函数的处理显得很关键。对热源反问题的处理中，我们的问题是如何寻找某种合理的附件条件，通过已知方程来解决方程右端的热源的具体位置并使其具有唯一性。本文利用微分方程并建立了满足温度分布的数学物理模型，从理论上导出了温度分布函数和热源位置的求解，并借助计算机软件画出了温度分布图。 二、问题的分析 对于热传导问题，为了使函数解决起来更容易，对于细秆的初始温度分布() g x我们可以设它在区间[0,L]连续，那么() g x可以展成正弦或余弦级数，对于有内热源的处理，由于细秆边界条件是齐次的，我们采用叠加原理把一根金属细秆的导热问题分解为有热源的具有其次边界条件的稳态导热问题和一个非稳态 其次问题，则原问题的解为 (,)1(,)2() u x t u x t u x =+。 对于源反问题的解决有如下3个问题: 1、反问题的唯一性:附加条件给得是否合理，也就是说，这个附加条件是否可以唯一确定热源的具体位置。 2、反问题的稳定性:反演所得到的热源的具体位置，该热源是否是连续地依赖于测量数据() h t? 3、反问题的数值解法:如何用可行的数值方法反演该热源的具体位置。用离散正则法将温度分布离散化，由已知初始温度分布再利用计算机软件得出热源位置 三、模型假设 1、金属细杆边界与外界无热量交换，即与外界绝缘 第三章非稳态导热分析解法 本章主要要求： 1、重点内容: ①非稳态导热的基本概念及特点； ②集总参数法的基本原理及应用; ③一维及二维非稳态导热问题。 2 、掌握内容：①确定瞬时温度场的方法； ②确定在一时间间隔内物体所传导热量的计算方法。 ３、了解内容：无限大物体非稳态导热的基本特点。 许多工程问题需要确定:物体内部温度场随时间的变化,或确定其内部温度达某一极限值所需的时间。如：机器启动、变动工况时,急剧的温度变化会使部件因热应力而破坏。因此,应确定其内部的瞬时温度场。钢制工件的热处理是一个典型的非稳态导热过程，掌握工件中温度变化的速率是控制工件热处理质量的重要因素；金属在加热炉内加热时，要确定它在炉内停留的时间，以保证达到规定的中心温度。 §3—１非稳态导热的基本概念 一、非稳态导热 1 、定义:物体的温度随时间而变化的导热过程称非稳态导热。 ２、分类：根据物体内温度随时间而变化的特征不同分: 1 )物体的温度随时间的推移逐渐趋于恒定值，即： 2 )物体的温度随时间而作周期性变化 如图3－1 所示，设一平壁，初值温度ｔ 0 ,令其左侧的表面温 度突然升高到并保持不变,而右侧仍与温度为的空气接触，试分 析物体的温度场的变化过程。 首先，物体与高温表面靠近部分的温度很快上升，而其余部分仍 保持原来的 t ０。 如图中曲线 HＢD ,随时间的推移，由于物体导热温度变化波及范 围扩大,到某一时间后，右侧表面温度也逐渐升高,如图中曲线HCD 、 HE 、ＨＦ。 最后，当时间达到一定值后,温度分布保持恒定,如图中曲线HG （若λ=const ，则ＨG 是直线）。 由此可见,上述非稳态导热过程中,存在着右侧面参与换热与不参 与换热的两个不同阶段。 （１ )第一阶段(右侧面不参与换热) 温度分布显现出部分为非稳态导热规律控制区和部分为初始温度区的混合分布,即：在此阶段物体温度分布受 t 分布的影响较大，此阶段称非正规状况阶段。 ( ２）第二阶段,（右侧面参与换热) 当右侧面参与换热以后,物体中的温度分布不受to 影响,主要取决于边界条件及物性,此时,非稳态导热过程进入到正规状况阶段。正规状况阶段的温度变化规律是本章讨论的重点。 2 )二类非稳态导热的区别:前者存在着有区别的两个不同阶段,而后者不存在。 3 、特点; 非稳态导热过程中,在与热流量方向相垂直的不同截面上热流量不相等，这是非稳态导热区别于稳态导热的一个特点。 原因：由于在热量传递的路径上，物体各处温度的变化要积聚或消耗能量,所以，在热流量传递的方向上。 二、非稳态导热的数学模型 1 、数学模型 非稳态导热问题的求解规定的｛初始条件，边界条件｝下，求解导热微分方程。 2 、讨论物体处于恒温介质中的第三类边界条件问题 在第三类边界条件下,确定非稳态导热物体中的温度变化特征与边界条件参数的关系。 已知:平板厚 2 、初温 to 、表面传热系数ｈ、平板导热系数，将其突然置于温度为的流体中冷却。 试分析在以下三种情况: <<１／h 、>>1/h 、=1/h 时，平板中温度场的变化。 1 ） 1/ｈ<< 因为 1/h 可忽略,当平板突然被冷却时,其表面温度就被冷却到，随着时间 非稳态导热㈠ 分离变量法在一维非稳态导热问题中的应用 一、无限大平壁一维非稳态导热（参考文献[1]PP45-50） 大平壁在等温介质中的冷却：常物性、无内热源、第三类边界条件（见图1）。 图1 令：f t t -=θ，导热微分方程为：0,022 τδθτθx x a ??=?? 初始条件为：()x f x =≤≤=θδτ,0,0；边界条件： δ δθ θ λτδθ τ====??-==??=x x x h x x x x |,0,0|, 0,00 假设解的形式为()()()ττθΓ?=x X x , 则两个常微分方程： 方程(3—2—3)的解是 由下表知，特征方程（3-2-4）的为()()x A x X m m m εcos =， 特征值为方程()λδεεh m m =tan 的正根，范数()()[] λ λεδλεεh h h N m m m +++=22 22 21 则根据()[]() m L m N dx x f A β??= 特征函数 ，得待定常数为： ，其中，εδβ= 二、半无限大物体一维非稳态导热（参考文献[2]PP40-45） 常物性、无内热源、第二、三类边界条件 典型问题：一半无限大物体，∞≤≤x 0，初始温度为F (x)，当时间τ＞0时，x ＝0的边界上以对流方式向温度为零度的介质传输热量，如图2所示。 该问题的数学描述为： 三、多维的齐次问题（参考文献[2]PP49-57） 典型问题：矩形截面的柱体，为二维非稳态导热，材料为常物性，物体内没有内热源，边界条件如图3所示。 τ??=??+??t a y t x t 1222 2， 0 ,00, 0,0042=+??====+??==??=t H y t b y t y t H x t a x x t x 时，时，时，时， ， ()y x F t ,0==时，τ 图3矩形截面柱体的二维非稳态导热 令()()()()τθΓ=y Y x X y x t ,,，分离方程如下： ,0, 0022=+====+''X H dx dX a x A dx dX x X X ） （β， 0,0,00 42=+====+''Y H dy dY b x B Y y Y Y ）（γ， ()()()C e τ γβατ22+-=Γ 上述问题的完全解为：()()()( )τ γβαγβτ2 2,,,,11 n m e y Y x X C y x t n m m mn n +-∞=∞ =∑∑ = 求解待定系数C mn 后，得： ()()()()()()()()()y d x d y x F y Y x X e y Y x X N N y x t a b n m n m m n n m n m ''''''=??∑∑+-∞=∞ =0011,,,,,1,,22γβγβγβττγβα 上式中出现的本征函数、本征值及范数可从表1-2中直接查得，即： ()()x x X m m ββcos ,=， ()( ) 2 22 22222 1H H a H N m m m +++=βββ，且m β为方程()2tan H a m m =ββ的正根： ()()y y X n n γγsin ,=， () () n n n n n n H H b H N +++=222 22 1 γγγ，且n γ为方程()4cot H b n n -=γγ的正根： ()( ) ()()()()()()()y d x d y x F y x y x H H b H H H a H e y x t a x b y n m n m n n n n n m m m n n m ' '''''?? ++++++=? ? ∑∑ ='='+-∞ =∞ =00 2 2222 22 222211 ,sin cos sin cos 4,,2 2γβγβγ γββττ γβα 四、某些非齐次或非线性问题的处理思路 1、线性、齐次多维非稳态热传导问题（参考文献[2]PP57-61） 对线性、齐次多维非稳态热传导问题，可以象一维问题那样，用分离变量法求解，其结果必定是二重 或三重级数，不便于计算和应用。在一定条件下，规则物体中齐次多维非稳态热传导问题的解可以简单地 计算传热学程序报告题目：一维非稳态导热问题的数值解姓名： 学号： 学院：能源与动力工程学院 专业：工程热物理 日期：2014年5月25日 一维非稳态导热问题数值解 求解下列热传导问题： 1.方程离散化 对方程进行控制体积分得到： 非稳态项：选取T 随x 阶梯式变化，有 扩散项：选取一阶导数随时间做显示变化，有 进一步取T 随x 呈分段线性变化，有 e P E e x T T x T )()(δ-=?? ， w W P w x T T x T )()(δ-=?? 整理可以得到总的离散方程为： 2.计算空间和时间步长 取空间步长为： h=L/N 网格Fourier 数为： 220x t x t F ??=??= α（小于0.5时稳定） 时间步长为： 3.建立温度矩阵与边界条件 T=ones(N+1,M+1) T(:,1)=Ti (初始条件温度都为0） T(1,:)=To (边界条件x=0处温度为1） T(N+1,:)=Te (边界条件x=L 处温度为0） 4.差分法求解温度 由离散方程可得到： 转化为相应的温度矩阵形式： 5.输入界面 考虑到方程的变量，采用inputdlg 函数设置5个输入变量，对这5个变量设置了默认值，如图1所示。在计算中可以改变不同的数值，得到不同的结果，特别注意稳定条件的临界值是0.5。根据设置的默认值，得到的计算结果如图2所示。 图1 matlab 变量输入界面 图2 默认值的计算结果 6.结果分析 根据上面的分析，给出了程序的输入界面，以及默认值状态下的数值解。可以通过改变不同的输入值，得到需要的分析结果，总结出了下面4点结论： （1）取F =0.48，得到一维非稳态导热结果如下图所示 图2 F0=0.48时一维非稳态导热从图中可以看出，对于长度L=1的细杆，初始时刻t=0时温度为0，边界条件x=0时，T=1,边界条件x=1时，T=0。随着时间的增加，温度从x=0通过导热的形式传递到x=1，不同时刻不同位置杆的温度都不同，并且随着时间的增加，杆的温度也逐渐增加。 （2）取F =0.48，可以得到不同位置的温度响应曲线，如下图所示 图3 F0=0.48时不同x位置处的温度响应图中红色曲线代表x=0.1位置的温度瞬态响应，黑色曲线代表x=0.2位置的温度瞬态响应，蓝色曲线代表x=0.4位置的温度瞬态响应。从图中可以看出，随着x的增加，曲线与x轴的交点值越大，温度开始传递到该位置的所需的时间越长。随着x的增加，温度响应曲线的变化速率越慢，最终的达到的温度也越低。 （3）取F =0.25，得到不同位置的温度响应曲线如下图所示 图4 F0=0.25时不同x位置处的温度响应图中三条曲线分别是x=0.1，x=0.2，x=0.4位置的温度瞬态响应。与图3 的F 0=0.48进行对比，两种情况下的F 值不同，F 值越大表明热扩散系数 的值 越大。从图中可以看出热扩散系数对于导热的影响，F 0=0.25时，与F =0.48相 比较，各位置开始响应时所需的时间较长，而且各位置响应曲线的变化速率较小， 最终的达到的温度也较低，说明了热扩散系数越小，热传导越慢，传递效率越低。 （4）取F 0= 0.51，得到非稳定的数值解如图所示 图5 F0=0.51时一维非稳态导热 图6 F0=0.51时不同x位置处的温度响应 从图中可以看出，对于显示格式的离散方程，并不是所有的F 值都能得到 有意义的解，必须要求F 0<0.5时才能得到稳定的数值解，当F >0.5时，会出现 二维稳态导热的数值计算 2.1物理问题 一矩形区域，其边长L=W=1，假设区域内无内热源，导热系数为常数，三个边温度为T1=0，一个边温度为T2=1，求该矩形区域内的温度分布。 2.2 数学描述 对上述问题的微分方程及其边界条件为：2222T T 0x y ??+=?? x=0，T=T 1=0 x=1，T=T 1=0 y=0，T=T 1=0 y=1，T=T 2=1 该问题的解析解：112121(1)sin n n n sh y T T n L x n T T n L sh W L ππππ∞=??? ?---??? ?=? ?-?? ??? ??? ∑ 2.3数值离散 2.3.1区域离散 区域离散x 方向总节点数为N ，y 方向总节点数为M ，区域内任一节点用I,j 表示。 2.3.2方程的离散 对于图中所有的内部节点方程可写为：2222,,0i j i j t t x y ???? ??+= ? ??????? 用 I,j 节点的二阶中心差分代替上式中的二阶导数，得： +1,,-1,,+1,,-1 2 2 2+2+0i j i j i j i j i j i j T T T T T T x y --+ = 上式整理成迭代形式：()()22 ,1,-1,,1 ,-12222+2()2() i j i j i j i j i j y x T T T T T x y x y ++=++++ (i=2,3……,N-1)，(j=2,3……,M-1) 补充四个边界上的第一类边界条件得：1,1j T T = (j=1,2,3……,M) ,1N j T T = (j=1,2,3……,M) ,1i j T T = (i=1,2,3……,N) ,2i M T T = (i=1,2,3……,N)一维非稳态导热问题的数值解
一维非稳态导热的数值计算
第三章非稳态导热分析解法
一维非稳态导热问题的数值解
一维非稳态导热的数值计算
一维非稳态导热通用程序
一维非稳态热传导热源反问题研究
第三章非稳态导热分析解法
第03次课J-分离变量法在一维非稳态导热问题中的应用
一维非稳态导热问题的数值解
传热学上机C程序源答案之一维非稳态导热的数值计算