《数学物理方程——4 热传导方程的建立》由会员分享，可在线阅读，更多相关《数学物理方程——4 热传导方程的建立（13页珍藏版）》请在金锄头文库上搜索。
1、数学物理方法数学物理方法抛物型方程的建立抛物型方程的建立2.3 抛物型方程的建立抛物型方程的建立推导固体的热传导方程时， 需要利用能量守恒定律和关于热传导的Fourier定律：推导固体的热传导方程时， 需要利用能量守恒定律和关于热传导的Fourier定律：热传导的Fourier定律定律：dt时间内，通过面积元时间内，通过面积元dS流入小体积元的热量流入小体积元的热量dQ与沿面积元外法线方向的温度变化率与沿面积元外法线方向的温度变化率u n 成正比成正比dSdt也与和成正比，即：也与和成正比，即：式中式中k是导热系数是导热系数1.热传导方程1.热传导方程( , , )udQk x y zdsdt
2、ku dSdtn= = ?数学物理方法数学物理方法抛物型方程的建立抛物型方程的建立其中其中 n 为曲面为曲面 ds 的外法向向量的外法向向量，k为热传导系数为热传导系数。故从故从211.ttSuQkds dtn= 12tt到到这段时刻流入曲面内部的热量为这段时刻流入曲面内部的热量为又，区域又，区域内温度升高吸收的热量为内温度升高吸收的热量为221( , , , )( , , , ).Qcu x y z tu x y z tdxdydz=数学物理方法数学物理方法抛物型方程的建立抛物型方程的建立21( , , , )ttucx y z t dt dxdydzt=21( , , , )ttucx y
3、 z t dxdydz dtt= 12.QQ=其中其中c为比热为比热,为密度。由能量守恒定律为密度。由能量守恒定律，有有221( , , , )( , , , ).Qcu x y z tu x y z tdxdydz=数学物理方法数学物理方法抛物型方程的建立抛物型方程的建立()SSukdSk u dSk u dxdydzn= =?.k udxdydz=21ttucdxdydz dtt 21.ttk udxdydz dt= （流入的热量流入的热量）又由又由Gauss公式有公式有故有故有（吸收的热量吸收的热量）数学物理方法数学物理方法抛物型方程的建立抛物型方程的建立,uck ut= 222 22
4、222,uuuuaautxyz=+=222222,uuuuxyz =+2.kac=因此有即其中因此有即其中数学物理方法数学物理方法抛物型方程的建立抛物型方程的建立( , , , ),F x y z t2,uauft= +.Ffc=若物体内部有热源若物体内部有热源，设单位时间内，设单位时间内，单位体积内所产生的热量为单位体积内所产生的热量为则易得相应的热传导方程为则易得相应的热传导方程为其中其中数学物理方法数学物理方法抛物型方程的建立抛物型方程的建立0,u t=0u =uf =如果我们考虑的是稳恒的温度场，如果我们考虑的是稳恒的温度场，即即u与时间与时间t 无关无关，温度分布达到某种动态平衡状态
5、温度分布达到某种动态平衡状态，则有则有这时上述两个方程分别变为这时上述两个方程分别变为Laplace方程方程和和Poisson方程方程数学物理方法数学物理方法抛物型方程的建立抛物型方程的建立解长为 l 的柱形管，一端封闭一端开放，管外空 气中含有某种浓度为 u0的气体向管内扩散，试写出该扩散问题的定解问题。xyoxldxx + +设管长方向为 x 轴，浓度函数为 u(x, t)，x = 0 端封闭，x = l 端开放，管的横截面积为 S 。由扩散定律知：xuDq = =( () )( () )SdxuuSdtqqtdttxdxx = = + + +dt 时间内流入微元 dv = Sdx 内的气
6、体分子满足的方程为：()()SdxuuSdtDxu xu tdttxdxx= +两边同除以Sdxdt ，整理得2. 扩散方程2. 扩散方程022 = = xuDtu数学物理方法数学物理方法抛物型方程的建立抛物型方程的建立若外界有扩散源，且扩散源的强度为，若外界有扩散源，且扩散源的强度为，( , , , )f x y z t这时，扩散方程应为这时，扩散方程应为222 2 222( , , , )uuuuaf x y z ttxyz+=从上面的推导可知，热传导和扩散这两种不同的物理现象， 但可以用同一类方程来描述. 从上面的推导可知，热传导和扩散这两种不同的物理现象， 但可以用同一类方程来描述.
7、其中，其中，2.aD=数学物理方法数学物理方法抛物型方程的建立抛物型方程的建立2.3热传导(或扩散)方程的定解条件2.3热传导(或扩散)方程的定解条件( , , ),x y z0.tu=( , , , ),f x y z t.Suf=1. 初值条件 若开始时刻物体的温度分布为则初值条件 若开始时刻物体的温度分布为则 u 应满足初值条件应满足初值条件2. 边值条件分三种情况： （边值条件分三种情况： （1）给出任何时刻物体表面的温度值： 若传热过程中边界上的温度为已知的函数则有）给出任何时刻物体表面的温度值： 若传热过程中边界上的温度为已知的函数则有数学物理方法数学物理方法抛物型方程的建立抛物型
8、方程的建立0.Su n=1u11()dQk uu dsdt=.uqknS= （2）给出任何时刻物体边界曲面的热流密度函数：）给出任何时刻物体边界曲面的热流密度函数：若传热过程中，物体若传热过程中，物体G与周围介质绝热，与周围介质绝热，即即S上的热量流量为零，则有上的热量流量为零，则有（3）若传热过程中，物体）若传热过程中，物体G与周围介质之间 有热量交换，以表示邻接处介质的温度，与周围介质之间 有热量交换，以表示邻接处介质的温度，由由Newton冷却定律知，冷却定律知，（流出（流出G的热量）的热量）数学物理方法数学物理方法抛物型方程的建立抛物型方程的建立而单位时间内通过而单位时间内通过 G 的
9、任一部分边界的任一部分边界S111().SQk uu ds =S流入周围介质的热量为的热量为流入周围介质的热量为的热量为单位时间流过单位时间流过2.SuQkdsn =12QQ=S由由及及的任意性知的任意性知11()0Sukk uun+=若外部介质 温度为零，则若外部介质 温度为零，则10.Sukk un+=数学物理方法数学物理方法抛物型方程的建立抛物型方程的建立例题例题解长为 l 的均匀杆的导热问题 （1）杆的两端温度保持零度 （2）杆的两端均绝热 （3）杆的一端恒温零度，另一端绝热试写出三种情况下的边界条件。0, 00= = =lxxuu（1）xukq = =（2）杆长方向的热量流动由傅里叶定律知，热流密度 两端绝热，既无热量流动，所以设 u(x, t) 为杆的温度函数0, 00= =lxxxu xu（3）或0, 00= = = = lxxxuu0, 00= = = =lx xuxu以上均为齐次边界条件。
