热传导方程扩散方程高教知识.ppt

《热传导方程扩散方程高教知识.ppt》由会员分享，可在线阅读，更多相关《热传导方程扩散方程高教知识.ppt（45页珍藏版）》请在三一文库上搜索。

1、从不同的物理模型出发，建立数学物理中三类从不同的物理模型出发，建立数学物理中三类 典型方程典型方程根据系统边界所处的物理条件和初始状态列出根据系统边界所处的物理条件和初始状态列出 定解条件定解条件提出相应的定解问题提出相应的定解问题第一章第一章 数学建模和基本原理介绍数学建模和基本原理介绍1全面分析1.1 1.1 数学模型的建立数学模型的建立数学模型建立的一般方法：数学模型建立的一般方法：确定所研究的物理量；确定所研究的物理量；建立适当的坐标系；建立适当的坐标系；划划出出研研究究小小单单元元，根根据据物物理理定定律律和和实实验验资资料料写写出出 该该单单元元与与邻邻近近单单元元的的相相互互作作

2、用用，分分析析这这种种相相互互 作作用用在在一一个个短短时时间间内内对对所所研研究究物物理理量量的的影影响响，表达为数学式表达为数学式;简化整理，得到方程。简化整理，得到方程。2全面分析2 热传导动方程热传导动方程第一节第一节 热传导方程的导出和定解条件热传导方程的导出和定解条件一、热传导方程的导出：一、热传导方程的导出：给定一空间内物体给定一空间内物体 ，设其上的点，设其上的点 在时刻在时刻 的温度为的温度为 。模型：模型：问题：问题：研究温度研究温度 的运动规律。的运动规律。3全面分析分析：（两个物理定律和一个公式）分析：（两个物理定律和一个公式）1 1、热量守恒定律、热量守恒定律:2 2

3、傅里叶、傅里叶（Fourier）热传导定律热传导定律:温度变温度变化吸收化吸收的热量的热量通过边通过边界流入界流入的热量的热量热源放热源放出的热出的热量量为热传导系数。为热传导系数。3 3、热量公式、热量公式:4全面分析任取物体任取物体 内一个由光滑闭曲面内一个由光滑闭曲面 所围成的区域所围成的区域 ，研究物体在该区域，研究物体在该区域 内热量变化规律。内热量变化规律。热传导方程的推导：热传导方程的推导：热量热量守恒守恒定律定律区域区域 内各点的温度从时刻内各点的温度从时刻 的温度的温度 改变为时刻改变为时刻 的温度的温度 所吸收（或所吸收（或放出）的热量，应放出）的热量，应等于等于从时刻从

4、时刻 到时刻到时刻 这段这段时间内通过曲面时间内通过曲面 流入（或流出）流入（或流出）内的热内的热量和热源提供（或吸收）的热量之和。即量和热源提供（或吸收）的热量之和。即 内温度变化所需要的热量内温度变化所需要的热量 =通过曲面通过曲面 流入流入 内的内的热量热量 +热源提供的热量热源提供的热量 下面分别计算这些热量下面分别计算这些热量5全面分析（1）内温度变化所需要的能量内温度变化所需要的能量那么包含点那么包含点 的体积微元的体积微元 的温度从的温度从 变为变为 所需要的热量为所需要的热量为 设物体设物体的比热（单位质量的物体温度改变的比热（单位质量的物体温度改变所需要的热量为所需要的热量为

5、密度为密度为 整个整个 内温度变化所需要的能量内温度变化所需要的能量6全面分析（2）通过曲面）通过曲面 进入进入 内的热量内的热量由傅里叶热传导定律，从由傅里叶热传导定律，从 到到 这段时间内通过这段时间内通过 进入进入 内的热量为内的热量为由高斯公式由高斯公式知知7全面分析（3）热源提供的热量）热源提供的热量用用 表示热源强度，即单位时间内从单位表示热源强度，即单位时间内从单位体积内放出的热量，则从体积内放出的热量，则从 到到 这段时间内这段时间内 内热内热源所提供的热量为源所提供的热量为由热量守恒定律得：由热量守恒定律得：由由 及及 的任意性知的任意性知8全面分析三维无热源热传导方程：三维

6、无热源热传导方程：三维有热源的热传导方程：三维有热源的热传导方程：（均匀且各向同性物均匀且各向同性物体，即体，即 都为常数的物体）都为常数的物体）其中其中称为非齐次项（自由项）。称为非齐次项（自由项）。通常称（通常称（1.5）为）为非齐次的热传导方程非齐次的热传导方程，而称（，而称（1.6）为为齐次热传导方程齐次热传导方程。9全面分析二、定解条件（初始条件和边界条件）二、定解条件（初始条件和边界条件）初始条件：初始条件：边界条件：边界条件：1 1、第一边界条件、第一边界条件（Dirichlet 边界条件）边界条件）特别地：特别地：时，物体表面保持恒温。时，物体表面保持恒温。10全面分析2 2、

7、第二边界条件、第二边界条件（Neumann 边界条件）边界条件）特别地：特别地：时，表示物体绝热。时，表示物体绝热。3 3、第三边界条件、第三边界条件(D-N 混合边界条件混合边界条件)其中：其中：表示表示 沿边界沿边界 上的单位外法线方向上的单位外法线方向 的方向导数的方向导数注：注：11全面分析注意第三边界条件的推导：注意第三边界条件的推导：研究物体与周围介质在物体表面上的热交换问题研究物体与周围介质在物体表面上的热交换问题 把一个温度变化规律为把一个温度变化规律为 的物体放入的物体放入 空空气介质中，已知与物体表面接触处的空气介质温度为气介质中，已知与物体表面接触处的空气介质温度为 ，它

8、与物体表面的温度，它与物体表面的温度 并不相同。这给出了并不相同。这给出了第三边界条件的提法。第三边界条件的提法。热传导热传导试验定试验定律或牛律或牛顿定律顿定律从物体流到介质中的热量和两者的温差成正比从物体流到介质中的热量和两者的温差成正比：其中比例常数其中比例常数 称为称为热交换系数热交换系数12全面分析流过物体表面流过物体表面 的流量可以从物质内部（傅里叶的流量可以从物质内部（傅里叶定律）和外部介质（牛顿定律）两个方面来确定：定律）和外部介质（牛顿定律）两个方面来确定：或或即得到（即得到（1.10）：）：13全面分析例例 长为长为l 的均匀杆，两端有恒定的均匀杆，两端有恒定热热流流进进入

9、其入，其强强度度为为，写出，写出这这个个热传导问题热传导问题的的边边界条件。界条件。在边界上有：在边界上有：若端点是绝热的，则解：x=l处：xq0q0nnx=0处：14全面分析三、定解问题三、定解问题定义定义1 在区域在区域上，由偏微分方程、初上，由偏微分方程、初始条件和边界条件中的其中之一组成的定解问题称为始条件和边界条件中的其中之一组成的定解问题称为初边值问题或混合问题初边值问题或混合问题。例如三维热传导方程的第一初边值问题为：例如三维热传导方程的第一初边值问题为：15全面分析始条件组成的定解问题称为始条件组成的定解问题称为初值问题或柯西问题初值问题或柯西问题。例如三维热传导方程的初值问

10、题为：例如三维热传导方程的初值问题为：定义定义2 在区域在区域上，由偏微分方程和初上，由偏微分方程和初16全面分析2 2、上述边界条件形式上与波动方程的边界条件上述边界条件形式上与波动方程的边界条件一样，但表示的物理意义不一样；一样，但表示的物理意义不一样；3 3、热传导方程的初始条件只有一个，而波动方热传导方程的初始条件只有一个，而波动方程有两个初始条件。程有两个初始条件。1 1、热传导、热传导方程不仅仅描述热传导现象，也可以方程不仅仅描述热传导现象，也可以刻画分子、气体的扩散等，也称扩散方程；刻画分子、气体的扩散等，也称扩散方程；注注4 4、除了三维热传导方程外，物理上，除了三维热传导方程

11、外，物理上，温度的分温度的分布在同一个界面上是相同的布在同一个界面上是相同的，可得，可得一维热传导方一维热传导方程：程：而对于薄片的热传导，而对于薄片的热传导，可得可得二维热传导方程：二维热传导方程：17全面分析3 3 拉普拉斯方程拉普拉斯方程当我们研究物理中的各类现象，如振动、热传导、当我们研究物理中的各类现象，如振动、热传导、扩散等的扩散等的稳定稳定过程时，由于表达该物理过程的物过程时，由于表达该物理过程的物理量理量 不随时间变化而变化，因此不随时间变化而变化，因此 .如果我们考虑的是一个稳定的热场，则可以得到如果我们考虑的是一个稳定的热场，则可以得到不随时间变化而变化的温度不随时间变化而

12、变化的温度 所满足的方所满足的方程：程：方程方程(*)(*)称为三维称为三维拉普拉斯拉普拉斯(Laplace)(Laplace)方程方程或者或者调和方程调和方程，它通常表示成为，它通常表示成为 或者或者 的形式。的形式。18全面分析拉普拉斯方程和泊松方程不仅描述稳定状态下温拉普拉斯方程和泊松方程不仅描述稳定状态下温度的分布规律，而且也描述稳定的浓度分布及静度的分布规律，而且也描述稳定的浓度分布及静电场的电位分布等物理现象。电场的电位分布等物理现象。其中其中如果我们考虑有源的稳定热场，则可以得到方程：如果我们考虑有源的稳定热场，则可以得到方程：非齐次方程非齐次方程 通常叫做通常叫做泊松泊松(Po

13、isson)(Poisson)方程方程，记作，记作或者或者19全面分析1 1、DilichletDilichlet问题。问题。2 2、NeumannNeumann问题。问题。2 2、NeumannNeumann问题。问题。3 3、第三边值问题。第三边值问题。20全面分析波动方程（双曲型）波动方程（双曲型）声波、电磁波、杆的振声波、电磁波、杆的振 动；动；热热传传导导方方程程（抛抛物物型型）热热传传导导，物物质质扩扩散散时时 的的浓浓度度变变化化规规律律,土壤力学土壤力学 中的渗透方程中的渗透方程;LaplaceLaplace方程方程 （椭圆型）（椭圆型）稳定的浓度分布稳定的浓度分布,静静电电场

14、场的的电电位位,流体的势。流体的势。总总 结：结：21全面分析1.3 1.3 定解问题的提法定解问题的提法初始条件和边界条件通称为初始条件和边界条件通称为定解条件定解条件。定解问题定解问题是指泛定方程和相应定解条件的结合体。是指泛定方程和相应定解条件的结合体。泛定方程和相应初始条件构成的定解问题称为泛定方程和相应初始条件构成的定解问题称为初值初值问题问题或者或者柯西柯西(Cauchy)(Cauchy)问题问题。22全面分析波方程的Cauchy问题由泛定方程和相应边界条件构成的定解问题称为由泛定方程和相应边界条件构成的定解问题称为边值问题。边值问题。Laplace方程的边值问题23全面分析由偏微

15、分方程和相应的初始条件及边界条件构成由偏微分方程和相应的初始条件及边界条件构成的定解问题称为的定解问题称为混合问题混合问题。热传导方程的混合问题热传导方程的混合问题24全面分析例例 设弦的两端固定于设弦的两端固定于x=0 和和x=l，弦的初始位移，弦的初始位移如下图，初速度为零，求弦满足的定解问题。如下图，初速度为零，求弦满足的定解问题。解：解：25全面分析一个定解问题的一个定解问题的适定性适定性(Well-posedness)(Well-posedness)包含以包含以下几个方面：下几个方面：1 1）解的）解的存在性存在性，即所提的定解问题是否有解；，即所提的定解问题是否有解；3 3）解的）

16、解的稳定性稳定性，即看定解问题的解是否连续依赖，即看定解问题的解是否连续依赖定解条件。也就是说，当定解条件有微小变动时，定解条件。也就是说，当定解条件有微小变动时，引起解的变动是否足够小。若是，则称解是稳定的，引起解的变动是否足够小。若是，则称解是稳定的，否则称解是不稳定的。否则称解是不稳定的。2 2）解的）解的唯一性唯一性，即所提的定解问题是否有唯一的，即所提的定解问题是否有唯一的解；解；26全面分析数理方程的一些基本概念数理方程的一些基本概念(1)偏微分方程偏微分方程 含有未知多元函数及其偏导数的方程，如含有未知多元函数及其偏导数的方程，如其中其中是未知多元函数，是未知多元函数，而而 是未

17、知是未知变变量；量；为为的偏的偏导导数数.有有时为时为了了书书写方便，通常记写方便，通常记27全面分析(2)方程的阶方程的阶 偏微分方程中未知函数偏导数的最高阶数称偏微分方程中未知函数偏导数的最高阶数称为方程的阶为方程的阶(3)方程的次数方程的次数 偏微分方程中最高阶偏导数的幂次数称为偏微分方程中最高阶偏导数的幂次数称为偏微分方程的次数偏微分方程的次数(4)线性方程线性方程 一个偏微分方程对未知函数和未知函数的所一个偏微分方程对未知函数和未知函数的所有（组合）偏导数的有（组合）偏导数的幂次数幂次数都是一次的，就称为线性方程，都是一次的，就称为线性方程，高于一次以上的方程称为非线性方程高于一次以

18、上的方程称为非线性方程28全面分析(5)准线性方程准线性方程 一个偏微分方程，如果仅对方程中所有最一个偏微分方程，如果仅对方程中所有最高阶偏导数是线性的，则称方程为准线性方程高阶偏导数是线性的，则称方程为准线性方程(6)自由项自由项 在偏微分方程中，不含有未知函数及其偏导数的在偏微分方程中，不含有未知函数及其偏导数的项称为自由项项称为自由项29全面分析5 5、微分方程的解、微分方程的解 古典解：如果将某个函数 u 代入偏微分方程中，能使方程成为恒等式，且方程中出现的偏导数都连续，则这个连续函数就是该偏微分方程的古典解。通解：解中含有相互独立的和偏微分方程阶数相同的任意常数的解。特解：通过定解条

19、件确定了解中的任意常数后得到的解。形式解：未经过严格数学理论验证的解为形式解。6 6、求解方法、求解方法分离变量法、特征线法、格林函数法30全面分析例例2.1 设 在直线R上具有二阶连续导数，验证 在 平面上都是 的古典解.解解 直接计算可得 代，到方程中即得结论成立.类似可证也是方程的古典解.31全面分析32全面分析(4)按未知函数及其导数的系数是否变化分为常系数按未知函数及其导数的系数是否变化分为常系数和变系数微分方程和变系数微分方程；(5)按自由项是否为零分为按自由项是否为零分为齐次方程齐次方程和和非齐次方程非齐次方程3 3、微分方程一般分类、微分方程一般分类 (1)按自变量的个数，分为

20、二元和多元方程按自变量的个数，分为二元和多元方程；(2)按未知函数及其导数的幂次，分为线性微分方程和按未知函数及其导数的幂次，分为线性微分方程和 非线性微分方程非线性微分方程；(3)按方程中未知函数导数的最高阶数，分为一阶、二阶按方程中未知函数导数的最高阶数，分为一阶、二阶 和高阶微分方程和高阶微分方程；33全面分析判断下列方程的类型思考34全面分析35全面分析一般二阶线性偏微分方程（n个自变量）两个自变量二阶线性偏微分方程的一般形式 线性方程的叠加原理线性方程的叠加原理36全面分析称形如的符号为微分算子。37全面分析二阶偏微分方程可简写为定解条件定解条件可简写为可简写为38全面分析 几种不同

21、的原因的综合所产生的效果等于这些不同原因单独产生的效果的累加。(物理上)39全面分析线性方程的解具有叠加特性 4 4、叠加原理、叠加原理 叠加原理的叠加原理的物理意义物理意义：几种不：几种不同的原因的综合所产生的效果同的原因的综合所产生的效果等于这些不同原因单独产生的等于这些不同原因单独产生的效果的累加。效果的累加。（以热传导方程为例）（以热传导方程为例）叠加原理叠加原理I设设是下面方程的解：是下面方程的解：在在G内收敛并且对内收敛并且对t可逐项求导一次，对可逐项求导一次，对x可逐项求可逐项求导两次，则和函数在导两次，则和函数在G内仍然是（内仍然是（1 1）的解）的解.若级数若级数 也就是说，如果也就是说，如果 是（是（1 1）的解，则其无限线性组合也是解。）的解，则其无限线性组合也是解。40全面分析叠加原理II41全面分析叠加原理III 设设是下面方程的解：是下面方程的解：若若在积分号下对在积分号下对 t 求导一次，对求导一次，对 x 可求导两次，则可求导两次，则在在G上是以下方程的解：上是以下方程的解：42全面分析叠加原理IV43全面分析44全面分析例 非齐次波动方程的Cauchy问题的解等于问题(I)和问题(II)的解之和45全面分析