(1.6)--§7 热传导方程与扩散方程模型建模方法.doc

§7热传导方程与扩散方程模型建模方法

（案例一）传导方程模型建模方法

一问题与分析：建立空间中某个物体的热传导过程中其每个位置处温度变化过程的数学模型。

二模型假设：1、假设空间物体占据位置为G；

2、建立空间直角坐标系；

3、在时刻，在物体G的位置处物体的温度是；

三模型建立与计算：

1、通过微曲面温度变化与热量传递的关系模型

物理定律模型：传热学中傅里叶实验定律表明，一个物体在很小的时间段内沿着内法向方向流过一个微小曲面的微热量与物体的温度沿着曲面的法线方向的方向导数成正比，即有：

这实际上就是表明：在一个小的时间段内，由于不同位置的温度不同，导致发生热量的流动，分析的对象是选取任意的过点的小的曲面，那么流过曲面的热量值，等于温度沿法线方向的导数值，乘以面积、时间长度、热传导系数等。由于时间很短，实际上方向导数只是考虑了开始时刻的导数值。

其中是在点位置处的热传导系数，就是一个比值；由于热量总是从温度高的一侧流向低的一侧，因此和异号。

2、通过封闭曲面热流量平衡模型

2.1无热源的平衡模型：假设物体G中任取一个封闭的曲面，其所包围的空间区域为，则在间隔内流进闭曲面内部的热量为：；这里表示的是温度分布沿单位外法线方向的方向导数。这个关系成立的基础是在这个物体G的表面上产生了热量，就是说表面上的温度在发生变化，并且这个变化能够传导到物质的内部，这个温度应该是已知的，至于如何产生的，可以有不同的机理。

根据温度数据的记录，时间间隔中，区域G内每一点的温度产生了变化值，这是因为吸收热量产生的，在这一点的附近单位体积的小区域内吸收的热量为，这里是比热，是密度；则根据三重积分的定义可知：整个区域G吸收的热量是：

注意这个量是根据物质的各个部分的温度值决定的，只要整个物质每个时刻的温度是已知的，就能计算出来这个热量值，至于这些热量来自于何处不需要考虑。

因此，根据上述的计算分析得到：

再根据高斯公式

得到：

交换积分顺序，进一步得到下面的积分结果，就是：

实际上这个物质G可以选取任何一个小的区域和任何一个时间段上来计算，因此上面的等式使得任何空间小区域和时间段上都是成立的，因此就有：

最终得到局部任何一点处的数量关系：

这就是温度分布模型。只要有温度差异，就会产生热量的流动传递。

进一步，如果物体是各向同性的、均匀的，则都是常数，则化为方程

得到标准的热传导方程。

模型计算：可以利用数学计算软件里的现成模块进行计算，也可以利用差分格式建立每个点位置和相应时刻相应的差分方程进行计算，可以获得离散近似解的形式。此处具体的差分格式形式省略。

四建模方法点评：利用微小元素分析计算建模方法，选定分析计算对象的适当的微小形式部分，通过利用物理原理，利用产生现象的来源数据进行的计算，建立相关数值的计算模式，然后在利用相关数值的记录结果进行的相应计算，使得二者相等，就可以得到微分方程模型。

（案例二）带有热源的温度分布模型

一问题与分析:建立数学模型，表达与分析带有热源的系统中，温度的分布状况规律，即每个点处的温度值满足的变化关系与规律。

二模型假设：1、在物体G内部单位时间单位体积中产生的热量为；

2、记：，

3、系统内每点的温度值为；

三模型建立与计算：

1、带有热源的新的平衡关系模型为：这个关系可以通过

可以得到：

这就是带有热源的热传导方程模型。

模型计算：可以利用数学计算软件里的现成模块进行计算，也可以利用差分格式建立每个点位置和相应时刻相应的差分方程进行计算，可以获得离散近似解的形式。此处具体的差分格式形式省略。

四建模方法点评：利用已经建立起来的热传导模型。增加辅助的热源因素，建立起新的反映温度平衡关系的数学模型。

（案例三）带有约束条件的热传导温度分布模型

一问题与分析:建立数学模型，表达与分析带有约束条件的系统中，温度的分布状况规律，即每个点处的温度值满足的变化关系与规律。

二模型假设：1、系统内每点的温度值为；

2、给出区域系统的边界上的每个点处的温度数值；

3、给出在区域边界；每个点处温度沿着特定射线方向的变化率，包括方向导数；

4、给出初始时刻各个位置处的温度值分布；

三模型建立与计算：

1、第一类边界模型：假设满足物体开始各点的温度分布：；物体在边界上的各个时刻的温度分布：，就是物体与外界接触部分的温度，就会源源不断地产生里面的温度分布，因此提出如下的边界条件模型：

，

这种形式的条件称之为第一类边界条件，又称之为狄利克雷(Dirichlet)边界条件。

模型计算：可以利用数学计算软件里的现成模块进行计算，也可以