一维非稳态热传导热源反问题研究摘要本文是关于热传导的正反问题的研究，即利用偏微分方程中典型热传导方程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、金属细杆边界与外界无热量交换，即与外界绝缘2、热源强度在整个时间段里始终保持常量。

3、在求解源反问题时，热源分布相对细杆长度来说，可假设为点热源。

四、符号说明(,)u x t :温度分布函数，即x 处在t 时刻的温度(,)f x t ：热源强度a 2：热扩散系数，单位为m 2/sL ：细杆长度()g x :初始温度分布函数，即t=0时的杆在x 处温度其他运算过程中使用符号在步骤中说明，再不赘述。

五、模型的建立与求解5.1 建立热传导微分方程，并求出温度分布函数(,)u x t由题意可的模型如下：222()u u a f x t x∂∂=+∂∂ 0<x<L t>0 （1） (0,)(,)0u t u L t == 0x ≥ （2） (,0)()u x g x = 0x L ≤≤ （3） 根据可以把有内热源的非稳态导热问题分解为有热源的具有齐次边界条件的稳态导热问题和一个非稳态齐次问题，即原方程的解为：(,)1(,)2()u x t u x t u x =+ （4）式中，1(,)u x t 是如下非稳态的解：(方程组A)22211u u a t x ∂∂=∂∂ （A-1）1(0,)1(,)0u t u l t == （A-2）1(,0)()2()()u x g x u x p x =-= （A-3） 2()u x 是如下稳态问题的解：（方程组B ）2222()0u a f x x∂+=∂ 0<x<L 2(0)2()0u u l ==对方程组A:利用分离变量法，假定偏微分方程的解是两个独立变数的乘积，即设1(,)()()u x t X x Y t =•代入微分方程（A-1）中，可得'2''XY a YX = 或'''2X Y X a Y==-λ(λ)为常数 （A-4）（若为+λ，则推导后所得出的解，其结果将对λ取任何值都不能满足边界条件），再由边值条件，有u1(0,t)=X(0)Y(t)=0 u1(L,t)=X(L)Y(t)=0必有(0)()0x x L ==。

由（A-4），有''0X X +λ=；X(0)=X(L)=0 （A-5）；'20Y a Y +λ= （A-6）；至此，通过分离变量，我们把微分方程的边值问题转化成为常微分方程的边值问题。

要是(A-5)有非零解，只要取λ=λn =222n L π（n=1，2，3···）对应的非零解为sin n x Lπ，类似方程(A-6)的解 2()()exp(())n n a Y t Y t t Lπ==- 于是得到方程组A 的一组线性无关的解： 2(,)sin exp(())n n x n a u x t t L L ππ=⋅- （n=1,2,3···）通过适当取值，设p(x)是函数sin n x L π的有限线性组合 即1()sin Nn n x p x c L π=∑最终得到21(,)sin exp(())N n n x n a u x t c t L L ππ=-∑ 其中，02()sin L n n x c p x dx L L π=⎰2121(,)()sin sin exp(())N n x n x n a u x t p x dx t L L L L πππ⎡⎤=-⎢⎥⎣⎦∑ 显然，(,)u x t 满足边值条件(,0)()u x p x =对方程组B ：即考虑有内热源的情况下，利用傅里叶变换求解2222()0u a f x x∂+=∂ 0<x<L 2(0)2()0u u l ==的解，得出()2()()F L u x F x x L -=+ 式中，()F x 是21()f x a-的二次积分。

最后，利用叠加定理将方程组A 与B 合并，即为所建模型温度分布函数的解： 21(,)1(,)2()2()()sin sin exp(())()N u x t u x t u x n x n x n a F L p x dx t F x x L L L L L πππ=+-⎡⎤=-++⎢⎥⎣⎦∑5.2 已知T 时刻的温度分布和初始温度分布，求解热源位置。

对于求源问题，由已知初始温度分布和一段时间后细秆的温度分布，那即为(,)()u x T f x =来求热源位置。

a 、离散正则化一维问题的一般提法为：(,)(,)(,)0,(0,)(,0)()0(0,)()(0,)(0,)(,)0(0,)t xx x u x t u x t f x x t T u x g x x u t t t T u t u l t t T =+→>∈⎧⎪=→>⎪⎨=ψ→∈⎪⎪==→∈⎩ （A ） 式中，边界条件(0,)()(0,)x u t t t T =ψ→∈当(,)f x t ，()g x ，()t ψ已知的情况下，原方程变为一维热传导方程的正问题，问题5.1就是属于一类一维热传导问题，对此借助参考文献[1]得出另一解为： ()000(,)2(,)()()(()())t t u x t f d k x t g k x t k x t d ∞=ττ--τψτ+ξ-ξ, ++ξ,ξ⎰⎰⎰ （B ）其中，21(,)exp()4x k x t t =- 对于已知(0,)(,)0u t u l t ==，将其换为附加条件(0,)()0u t t t T =φ→<< （C ）此时根据已知条件，原问题就构成一维热传导方程的寻源反问题，把C 式代入B 式中可得到：()000()2(0,)()()(()())t t f d t k t d g k t k t d ∞ττ=φ+-τψττ-ξ-ξ, +ξ,ξ⎰⎰⎰(D) 以下为使用离散正则法解决源反问题：式（D ）可写成第一类积分算子方程的形式：,,Af F f Z F U =∈∈ （E ）其中，A 为第一类积分算子，Z 为解空间，U 为数据空间。

由于原问题是不适定的，故将积分算子离散化后得到的线性方程组常常是病态的，而且随着维数的增加其病态状况更加严重。

这时，直接求解离散方程显然是不行的，必须将其正则化后求解才能得到稳定的近似解。

离散正则化的过程可归结如下：1)将积分方程离散化得,,h h h h h h h A f F f W F L =∈∈其中， h A 为积分算子的离散矩阵；h f 为方程的离散解； h F 为近似数据；h为离散步长；2)形成相应的法方程T h T h h h h A A f A F = 3)法方程的正则化'()()T h h h h a a A A C f CF +α=-其中，C 为单位矩阵；4)利用偏差原理决定正则参数()a α=δ，使得((h a φδ))=0，其中δ为实际数据与测量数据之间的误差，((h h h ha A f F 2φδ))=||-||-δ； 5)求得离散正则解1()()h T T h a h h h f A A C A F -δ=+αb 、对离散正则模型进行数值实验（注：本文实验模拟特殊函数借助文献<2>） 设T 时刻细秆的温度分布：T=1，2()0.00001(1),01T g x x x x =-<<3414(exp(1))),01324t t t t t t φ()=-++-<<4.5(),01t t T ψ=<< 则根据拟最优准则<4> 得 345()325T t t t F t =-+由式E 得22(1)T f t t =-。

在实际应用中，得到的只是一些离散的数据，并且存在着一定的误差故给()T F t 一个扰动δ作为式E 的右端项()f F t ，应用离散正则化方法来求解()f F t 作为右端项的第一类积分算子方程，求解结果为 ()f f t 。