扩散方程课件

,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,电子科技大学物理电子学院,周秀丽,第二篇 数学物理方程,本篇主要内容,：二阶线性偏微分方程的建立和求解,重点,：数学物理方程求解方法中的分离变量法和行波法.,特点,：加强物理模型和数学物理思想的介绍，以便充分了解模型的物理意义，有利于根据数学物理模型建立数学物理方程,第二篇 数学物理方程本篇主要内容：二阶线性偏微分方程的建立和,1,数学物理思想,数学物理方程（简称,数理方程,）是指从物理学及其它各门自然科学、技术科学中所导,出的函数方程，主要指偏微分方程和积分方程,数学物理方程所研究的内容和所涉及的领,域十分广泛，它深刻地描绘了,自然界中,的许多,物理现象,和,普遍规律,.,数学物理思想数学物理方程（简称数理方程）是指从物理学及其它各,2,声振动是研究声源与声波场之间的关系,热传导是研究热源与温度场之间的关系,泊松（S.D.Poisson 17811840，法国数学家）方程表示的是电势（或电场）和电荷分布之间的关系,定解问题,从物理规律角度来分析，数学物理定解问题表征的是场和产生这种场的源之间的关系,声振动是研究声源与声波场之间的关系热传导是研究热源与温度场之,3,根据分析问题的不同出发点，把数学物理问题分为,正向问题,和,逆向问题,.,不同出发点,正向问题，即为已知源求场,逆向问题，即为已知场求源.,前者是经典数学物理所讨论的主要内容.后者是高等数学物理（或称为现代数学物理）所讨论的主要内容,根据分析问题的不同出发点，把数学物理问题分为正向问题和逆向问,4,多数为二阶线性偏微分方程,振动与波（振动波，电磁波）传播满足,波动方程,热传导问题和扩散问题满足,热传导方程,静电场和引力势满足,拉普拉斯方程或泊松方程,数学物理方程的类型和所描述的物理规律,多数为二阶线性偏微分方程振动与波（振动波，电磁波）传播满足波,5,三类典型的数学物理方程,三类典型的数学物理方程,双曲型方程,波动方程为代表,抛物型方程,热传导方程为代表,椭圆型方程,泊松方程为代表,退化为拉普拉斯方程,三类典型的数学物理方程三类典型的数学物理方程双曲型方程波动方,6,分离变量法,偏微分方程,标准的常微分方程,标准解，即为各类特殊函数,三类数学物理方程的一种最常用解法,分离变量法偏微分方程标准的常微分方程标准解，即为各类特殊函数,7,第九章 数学建模-数学物理定解问题,9.1 数学建模-波动方程类型的建立,具有波动方程的数理方程的建立,弦的横振动,杆的纵振动,讨论,定解条件,传输线方程,第九章 数学建模-数学物理定解问题9.1 数学建模-,8,9.1.1波动方程的建立,1.弦的微小横振动,考察一根长为,且两端固定、水平拉紧的弦,讨论如何将这一物理问题转化为数学上的定解问题要确定弦的运动方程，需要明确：,确定弦的运动方程,（2）被研究的物理量遵循哪些物理定理？,牛顿第二定律.,（3）按物理定理写出数学物理方程（即建立泛定方程）,要研究的物理量是什么？,弦沿垂直方向的位移,9.1.1波动方程的建立1.弦的微小横振动考察一根长为且两,9,注意：,物理问题涉及的因素较多，往往还需要引入适当假设才能使方程简化,数学物理方程必须反映弦上任一位置上的垂直位移所遵循的普遍规律，所以考察点不能取在端点上，但可以取除端点之外的任何位置作为考察点,注意：,10,根据牛顿第二定律,方向运动的方程可以描述为,（,9.1.1,）,作用于小段,的纵向合力应该为零：,(9.1.2),仅考虑,微小,的横振动，,夹角,为很小的量，忽略,及其以上的高阶小量，则根据级数展开式有,根据牛顿第二定律方向运动的方程可以描述为（9.1.1）,11,注意到,:,故由图9.11得,这样，(9.1.1)和(9.1.2)简化为,注意到:故由图9.11得这样，(9.1.1)和(9.1.2,12,因此在微小横振动条件下，可得出,，弦中张力不随,而变，,可记为,故有,(9.1.5),变化量,可以取得很小，根据微分知识有下式成立,这样，,段的运动方程(9.1.5)就成为,(9.1.6),因此在微小横振动条件下，可得出，弦中张力不随而变，可记为,13,即为,（9.1.7）,上式即为弦作微小横振动的运动方程，简称为弦振动方程,其中,讨论：,（1）若设弦的重量远小于弦的张力，则上式(9.1.7)右端的重力加速度项可以忽略由此得到下列齐次偏微分方程：,（9.1.8）,称式（9.1.8）为弦的自由振动方程,即为 （9.1.7,14,(2)如果在弦的单位长度上还有横向外力,作用，则式(9.1.8)应该改写为,(9.1.9),式中,称为力密度,，为,时刻作用于,处单位质量上的横向外力,式（9.1.9）称为弦的受迫振动方程.,(2)如果在弦的单位长度上还有横向外力 作用，则式(9.1,15,2、均匀杆的纵振动,段的运动方程为,（9.1.10）,可得,（,9.1.11,）,这就是,杆的纵振动方程,2、均匀杆的纵振动段的运动方程为（9.1.10）可得,16,讨论,(1)对于均匀杆，,和,是常数，(9.11）可以改写成,（,9.,.12,）,其中,这与弦振动方程（9.8）具有完全相同的形式,（2）杆的受迫振动方程跟弦的受迫振动方程（9.9）,完全一样，,只是其中,应是杆的单位长度上单位,横截面积所受纵向外力,讨论(1)对于均匀杆，和是常数，(9.11）可以改写,17,3.传输线方程（电报方程）,（9.1.13）,同理可得：,(9.1.14),式（9.1.13）及（9.1.14）即为一般的传输线方程（或电报方程）,3.传输线方程（电报方程）（9.1.13）,18,(1)无失真线,(9.1.15),其中,(2)无损耗线,（9.1.16）,(9.1.17),具有与振动方程类似的数学形式，尽管它们的物理本质根本不同,(1)无失真线 (9.1.15)其中(2,19,(3)无漏导，无电感线,（9.1.18）,(9.1.19),它们具有与下节将讨论的一维热传导方程类似的数学形式，,尽管它们的物理本质根本不同,(3)无漏导，无电感线,20,9.1.2 波动方程的定解条件,定解条件：初始条件和边界条件,1.初始条件,波动方程含有对时间的二阶偏导数，它给出振动过程中每点,的加速度要确定振动状态，需知道开始时刻每点的位移和速度,波动方程的初始条件通常是,（,9.1.22,）,9.1.2 波动方程的定解条件定解条件：初始条件和边界条件1,21,例911,一根长为,的弦，两端固定于,和,在距离坐标原点为,的位置将弦沿着横向拉开距离,，如图9.5所示，然后放手任其振动，试写出初始条件。,x,u,o,b,l,h,图,9,.5,【,解,】初始时刻就是放手的那一瞬间，,按题意初始速度为零，即有,初始位移如图所示,例911 一根长为 的弦，两端固定于 和,在距离坐,22,2.边界条件,常见的线性边界条件分为三类：,第一类边界条件,直接规定了所研究的物理量在边界上的数值,第二类边界条件,规定了所研究的物理量在边界外法线方向上方向导数的数值,（,9.1.23,）,（,9.1.24,）,2.边界条件 常见的线性边界条件分为三类：第一类边界条件,23,第三类边界条件,规定了所研究的物理量及其外法向导数的线性组合在边界上的数值,（,9.1.25,）,其中,是时间,的已知函数，,为常系数,第三类边界条件 规定了所研究的物理量及其外法向导数的线性组合,24,9.2 数学建模热传导方程类型的建立,9.2.1数学物理方程热传导类型方程的建立,1.热传导方程,推导固体的热传导方程时，,需要利用能量守恒定律和关于热传导的傅里叶定律：,热传导的,傅里叶定律,：,时间内，通过面积元,流入小体积元的热量,与沿面积元外法线方向的温度变化率,成正比,也与,和,成正比，即：,（,9.2.1,）,式中,是导热系数,9.2 数学建模热传导方程类型的建立9.2.1数学物理方程,25,图9.8,取直角坐标系Oxyz,如图9.8,表示,t,时刻物体内任一点（x,y,z）处的温度,在,d,t,时间内通过,ABCD,面流入的热量为,同样，在,时间内沿y方向和,z,方向流入立方体的热量分别为,图9.8取直角坐标系Oxyz,如图9.8 表示t时刻物体内,26,在,t,到,时间内，小体积元的温度变化是,如果用,和,分别表示物体的,密度,和,比热,，,则根据能量守恒定律得热平衡方程,或写成,(9.2.2),在t到时间内，小体积元的温度变化是 如果用和分别表示物体的密,27,2.扩散方程,(9.2.3),其中,将一维推广到三维，即得到,(9.2.4),上述方程与一维热传导方程具有完全类似的形式,2.扩散方程,28,若外界有扩散源，且扩散源的强度为,这时，扩散方程应为,（,9.2.5,）,从上面的推导可知，热传导和扩散这两种不同的物理现象，,但可以用同一类方程来描述.,若外界有扩散源，且扩散源的强度为这时，扩散方程应为（,29,9.2.2 热传导（或扩散）方程的定解条件,1 初始条件,热传导方程的初始条件一般为,（,9.2.6,）,2 边界条件,第一类:,已知任意时刻,边界面,上的温度分布,(9.2.7),直接给出函数,u,在边界上的数值,所以是第一类边界条件.,9.2.2 热传导（或扩散）方程的定解条件 1 初始条,30,2.第二类,已知任意时刻,从外部通过边界流入物体内的热量。,设单位时间内通过边界上单位面积流入的热量为,.,考虑物体内以边界上面积元,为底的一个小圆柱体，,如图9.10所示.,图9.10,物体内部通过,流入小柱体的热量为,小柱体内温度升高,所需要的热量,随着柱高,趋于零而趋近于零,2.第二类 已知任意时刻从外部通过边界流入物体内的热量。,31,所以当,由热平衡方程给出:,考虑到,时,则得,(9.2.8),所以当由热平衡方程给出:考虑到时,则得,32,3.第三类,根据牛顿冷却定律:单位时间从周围介质传到边界上单位面积,的热量与表面和外界的温度差成正比,即,这里,是外界媒质的温度,.,为常数,与推导条件(9.2.11)相似,此时可得边界条件,(9.2.9),其中,3.第三类 根据牛顿冷却定律:单位时间从周围介质传到边界,33,9.3 数学建模稳定场方程类型的建立,9.3.1 数学建模稳定场方程类型的建立,1 静电场的电势方程,直角坐标系中,泊松方程,为,(9.3.1),若空间,中无电荷，即电荷密度,，上式成为,(9.3.2),称这个方程为,拉普拉斯方程,.,9.3 数学建模稳定场方程类型的建立9.3.1 数学建模,34,2.稳定温度分布,导热物体内的热源分布和边界条件不随时间变化,故热传导方程中对时间的偏微分项为零，从而热传导方程,（9.2.1)，(9.2.2)即为下列拉普拉斯方程和泊松方程.,(9.3.3),(9.3.4),2.稳定温度分布 导热物体内的热源分布和边界条件不随时间变,35,9.3.2 泊松方程和拉普拉斯方程的定解条件,泊松方程和拉普拉斯方程的定解条件不包含初始条件,而只有边界条件.,边界条件分为三类:,1、在边界上直接给定未知函数,即为第一类边界条件,2、在边界上给定未知函数导数的值,即为第二类边界条件,3、在边界上给定未知函数和它的导数的某种线性组合,即第三类边界条件.,9.3.2 泊松方程和拉普拉斯方程的定解条件 泊松方程和拉普,36,第一、二、三类边界条件可以统一地写成,(9.3.5),其中,是边界上的变点；,表示物理量,沿边界外法线方向的方向导数；,为常数，它们不同时为零,第一、二、三类边界条件可以统一地写成,37,9.4 数学物理定解理论,9.4.1 定解条件和定解问题的提法,边界条件的类型,除了前面我们介绍的,第一、第二、第三类边界条件,之外，还有其它边界条件，如,自然边界条件，衔接条件,周期性条件,和,无边界条件,9.4 数学物理定解理论 9.4.1 定解条件和定解问题的提,38,9.4.2 数学物理定解问题的适定性,(1)解的存在性,看所归结出来的定解问题是否有解；,(2)解的唯一性,看是否