用拉普拉斯变换求解非稳态导热.doc

用拉普拉斯变换求解非稳态导热 06数学系 骆逸弘 P在一个学期的偏微分方程课程的学习之后，我们可以粗浅的认识到，如果一个微分方程中出现多元函数的偏导数，或者说如果未知函数和几个变量有关，而且方程中出现未知函数对几个变量的导数，那么这种微分方程就是偏微分方程。 在科学技术日新月异的发展过程中，人们研究的许多问题用一个自变量的函数来描述已经显得不够了，不少问题有多个变量的函数来描述。比如，从物理角度来说，物理量有不同的性质，温度、密度等是用数值来描述的叫做纯量；速度、电场的引力等，不仅在数值上有不同，而且还具有方向，这些量叫做向量；物体在一点上的张力状态的描述出的量叫做张量，等等。这些量不仅和时间有关系，而且和空间坐标也有联系，这就要用多个变量的函数来表示。应该指出，对于所有可能的物理现象用某些多个变量的函数表示，只能是理想化的，如介质的密度，实际上“在一点”的密度是不存在的。而我们把在一点的密度看作是物质的质量和体积的比当体积无限缩小的时候的极限，这就是理想化的。介质的温度也是这样。这样就产生了研究某些物理现象的理想了的多个变量的函数方程，这种方程就是偏微分方程。也许这样的定义很宽泛而且显得并不专业，但是不能否认偏微分方程在各个 今天我主要想谈的问题是拉普拉斯变换法在非稳态导热中的应用。 首先，简单的来看一下拉普拉斯变换的概念。 设函数f(t)当t≥0时有定义，而且积分∫f(t)e-stdt（s时一个复参量）在s的某一域内收敛，则由此积分所确定的函数称为函数f(t)的拉普拉斯变换.记为F(s) =L[ f (t)]=∫f(t)e-stdt。 拉普拉斯变换有许多非常好的性质，如线性性质、微分性质、积分性质、位移性质、延迟性质、初值定理和终值定理、卷积定理等。这些性质在解题时非常重要。在利用拉普拉斯变换求解导热问题时，关键的一步是把变换后的函数从复变量s区域变回到时间变量t区域的逆变换。而许多逆变换都可直接或利用性质转化之后通过查拉普拉斯变换表得到。 利用拉普拉斯变换求解非稳态导热问题（偏微分方程）的一般步骤：（1）根据问题建立数学模型（偏微分方程）； （2）将温度看作时间t的函数，对方程及定解条件关于t取拉普拉斯变换，把偏微分方程和定解条件化为象函数的常微分方程的定解问题； （3）解常微分方程，求出象函数U(x,s)； （4）取拉普拉斯逆变换，求出温度函数u(x,t)。 在具体的来看一看边界热流为常数的非稳态导热问题。 一个半无限大物体(x ≥0)的初始温度为零，当时间t0时，在x=0的边界上有恒定热流 qw的作用，试求t0时物体中的温度分布。 【分析】设u 表示物体中的温度，x表示坐标，t表示时间，λ表示导热系数，则温度是时间及坐标的函数，即u=u(x,t)。据题意，该问题的数学模型为 ?u∕?t=a?2u∕?x2 0x+∞,t0 u=0 x≥0,t=0 ?u∕?x=-qw∕λ x=0,t0 u=0 x→+∞,t0 对上述定解问题关于t取Laplace变换，并利用微分性质和初始条件可得 L[u(x,t)] =U(x,s) L[?u∕?t]= sU(x,s)-u(x,0)=sU(x,s) L[?u∕?x]= dU(x,s)∕dx L[?2u∕?x2]= d2U(x,s)∕dx2 于是，上述方程组转化为 d2U(x,s)∕dx2-sU/a=0 0x+∞ dU/dx=-qw∕λs x=0 U(x,s)=0 x→+∞ 这是一个二阶常系数线性齐次微分方程的边值问题。由特征值法易求得通解为 U(x,s)=c1e√sx/a +c2e-√sx/a 由边界条件可知 c1=0, c2= qw√(a/s)∕λs 对上式取拉普拉斯逆变换并查表得 U(x,t)= qw√a∕λ[2√(t/∏)e-x2/4at-xerfc(x/2√at)/√a] 当然，以上的认识是很粗浅的，而且偏微分方程在其他的领域中也在发挥着自己的作用，随着数学的不断发展，偏微分方程也会有更广阔的应用空间。 （由于我的word没有安装公式编辑器，所以文中的公式都是手工打的，看起来有点乱，请老师见谅！）