[拼音]:paowuxing pianweifen fangcheng

[外文]:partial differential equation of parabolic type

简称抛物型方程,一类重要的偏微分方程。热传导方程是最简单的一种抛物型方程。

热传导方程

研究热传导过程的一个简单数学模型。根据热量守恒定律和傅里叶热传导实验定律导致热传导方程

   (1)

式中u是温度;Δ是拉普拉斯算符;α2是导温系数;

;k是热传导系数;с、ρ分别是比热和密度;

;F是外加热源密度。自然界还有很多现象同样可以用方程(1)来描述,例如分子在介质中的扩散过程等,因此方程(1)通常亦称为扩散方程。

定解问题

为了确定一个具体的热传导过程,除了列出方程(1)以外,还必须知道物体Ω的初始温度(初始条件)和在它的边界嬠Ω上所受到的外界的影响(边界条件)。

初始条件:

   (2)

边界条件,最通常的形式有三类。

第一边界条件(或称狄利克雷条件):

   (3)

即表面温度为已知函数。

第二边界条件(或称诺伊曼条件):

   (4)

式中n是嬠Ω的外法向,即通过表面的热量已知。

第三边界条件(或称罗宾条件):

(5)式中α≥0;即物体表面给定热交换条件。

除了以上三类边界条件外还可以在边界嬠Ω上给定其他形式的边界条件,如斜微商条件、混合边界条件等。

方程(1)连同初始条件(2)以及边界条件(3)、(4)、(5)中的任意一个一起构成了一个定解问题,根据边界条件的不同形式,分别称为第一、二、三边值问题,统称为热传导方程的初边值问题或混合问题。若Ω呏R3,则由方程(1)和初始条件(2)构成的定解问题称为热传导方程的初值问题或柯西问题。

基本解与格林函数

基本解是点热源的影响函数。如果在t=0时刻在(ξ,η,ζ)处给定单位点热源,即u0(x,y,z,0)=δ(ξ,η,ζ)(δ是狄喇克函数),则当t>0时由它引起的在全空间 R3的温度分布(即热传导方程(1)的解)称为热传导方程的基本解。通过傅里叶变换可以得到它的表达式。当t>0时

热传导方程初值问题(1)、(2)的解可通过叠加的步骤由基本解生成

对于一个有界区域Ω,若边界温度为零,在初始时刻在(ξ,η,ζ)处给定一个单位点热源u(x,y,z,0)=δ(ξ,η,ζ),当t>0时由它引起在Ω内的温度分布(即热传导方程的解)称为热传导方程第一边值问题的格林函数,记作G(x-ξ,y-η,z-ζ,t)。根据格林公式

,式中l*是l的共轭算子,

任意第一边值问题(1)、(2)、(3)的解都可通过格林函数表为

;格林函数可以通过基本解来表示:

这里

时是一个定义在捙×[0,∞)上的充分光滑函数。对于一维问题或Ω为立方体等特殊区域,格林函数可以通过分离变量法或镜像法去求得。

极值原理

一个内部有热源的热传导过程(即在方程(1)中ƒ≥0),它的最低温度一定在边界上或初始时刻达到,这就是所谓的极值原理。事实上,还可以有更强的结论:

(1)如果在t=T时在Ω内部某一点达到了最低温度,那么在这个时刻T以前(即t

(2)如果这个最低温度只在t=T时刻的某一边界点P达到,那么在这一点上

(n是嬠Ω的外法向),此即所谓的边界点引理。

极值原理与边界点引理在热传导方程的研究中有很多应用,它的一个最直接的推论就是导出了热传导方程初边值问题解的惟一性和稳定性。

至于初值问题(1)、(2)的解的惟一性,它与解在无穷远点的性态有关。如果对于初值问题(1)、(2),附加上无穷远点增长阶的限制

,这里A,M是任意给定正常数,那么由极值原理可以证明初值问题(1)、(2)的解必惟一。

解的正则性(光滑性)

若ƒ呏0,则由初值问题解的表达式可看出,若u0(x,y,z)有界连续,则初值问题(1)、(2)的解u(x,y,z,t)当t>0时都是无穷次连续可微的,而且关于空间变量x,y,z是解析的,关于时间变量t属于谢弗莱二类函数,即在|x|

当ƒ扝0时,热传导方程解的可微性质与ƒ的性质有关,例如为了得到热传导方程的古典解,除了需要假定ƒ(x,y,z,t)连续以外,还要求对x,y,z或对t是赫尔德连续的。

解的渐近性

如果边界上的温度以及热源密度与时间无关(即

),则热传导过程将趋于稳定状态,也就是当t→∞时,不管什么初始条件,物体内部温度总趋于同一个极限(稳定态的温度分布u(x,y,z)),它是椭圆边值问题:

的解。

解的半群性质

热传导是一个单向的不可逆过程,热总是由高温流向低温。如果边界温度为零,S(t)表示由初始时刻的温度场映到t时刻的温度场的线性解算子,即

,由于热传导的不可逆性质,因此算子族

具有半群性质:

(1)S(0)=I(I为恒同算子);

(2)S(t+τ)=S(t)S(τ)t,τ≥0;

(3)

。由泛函分析中的希尔-吉田定理,存在一个相应的无穷小生成子A,S(t)=e-tA,使得具有齐次边条件的第一边值问题(1)、(2)、(3)的解具有明显的表达式

,式中

线性和拟线性抛物型方程

。二阶线性偏微分方程

  (6)在区域Q内称为是抛物型的,如果存在常数α >0,使得对于任意ξ∈Rn,(x1,x2,…,xn,t)∈Q 有

。如果αij连续可微,那么(6)可改写为

   (7)的形式。(7)称为具有散度形式的抛物型方程,(6)称为非散度形式的抛物型方程。当

时,(6)与(7)是有区别的,不能互推。如果方程(6)、(7)中的系数和右端还依赖于u,墷u,则(6)和(7)称为拟线性抛物型方程。

抛物型方程和椭圆型方程的研究有相似的地方,它们互相影响、互为借鉴。椭圆型方程理论很多结果在抛物型方程中都有相应的定理,例如先验估计、极值原理等。

拟线性蜕化抛物型方程

考虑在绝热过程中气体通过多孔介质的流动,这个过程可由下述方程来刻画:

,式中m>1,u是气体密度,通常研究它的非负解。由于当u=0时方程蜕化,因此它是一个拟线性蜕化抛物型方程。对于这个问题的系统理论研究是从 1957年开始的。解u的支集的边界是一条自由边界,通过自由边界

u一般不连续,因此这个方程一般只存在在索伯列夫意义下的广义解,而且由于当u=0时方程蜕化为一阶方程,因此与热传导方程不同,扰动的传播速度是有限的。

反应扩散方程(组)

形如

的半线性抛物型方程组叫做反应扩散方程组。除了研究各种定解问题外,由于(8)的解常具有行波解u(v·x-сt)以及当t→∞时 u(x,t)趋于椭圆型方程组相应的边值问题的解(称为平衡解)这样的性质,因此以研究平衡解的稳定性为核心的各种问题就构成了半线性抛物型方程(组)的定性理论(或叫几何理论)。

抛物型偏微分方程c语言代码,抛物型偏微分方程相关推荐

  1. 解抛物型微分方程c语言程序,偏微分方程的matlab代码_偏微分方程matlab程序,matlab偏微分代码...

    matlab偏微分方程含源代码,大家可以下载运行,里面的都能用 Pie Eais ae h t lat BLep 话k气飞"目基回 古班比比,一生热件身的些懂 1身定 回量1 古典显式格式稳 ...

  2. 深度学习(Deep Ritz,Galerkin,PINN)求解偏微分方程(PDE)实现代码地址

    深度学习求解偏微分方程的实现代码可以在GitHub上查找.下面是几个可能有用的项目地址: Deep Ritz: https://github.com/yseop/DeepRitz Galerkin: ...

  3. 遗传算法c语言程序,遗传算法c语言代码.doc

    遗传算法c语言代码 遗传算法代码 #include #include #include #include #include struct group //染色体的结构 { int city[citie ...

  4. 071_html语言代码

    1. ISO语言代码 1.1. 国际标准化组织(International Organization for Standardization, ISO)简称ISO, 是一个全球性的非政府组织, 是国际 ...

  5. 嵌入式C语言代码规范

    C语言代码规范 参考安富莱C语言编码规范 1.文件与目录 1.文件及目录的命名规定可用的字符集是[A-Z:a-z:0-9:._-]. 2.源文件名后缀用小写字母 .c 和.h. 3.文件的命名要准确清 ...

  6. java 与c 运行效率_Java语言与C语言代码运行效率的比较

    <Java语言与C语言代码运行效率的比较>由会员分享,可在线阅读,更多相关<Java语言与C语言代码运行效率的比较(2页珍藏版)>请在人人文库网上搜索. 1.Java语言与C语 ...

  7. Python语言学习:python语言代码调试—异常处理之详细攻略

    Python语言学习:python语言代码调试-异常处理之详细攻略 目录 python语言代码调试-异常处理 异常捕捉可以使用 try/except 语句 相关文章 Python3 错误和异常 | 菜 ...

  8. 编程笔试(解析及代码实现):猴子吃桃。猴子第一天吃了若干个桃子,当即吃了一半,还不解馋,又多吃了一个…的C++、Java、Python、C#等语言代码实现

    编程笔试(解析及代码实现):猴子吃桃.猴子第一天吃了若干个桃子,当即吃了一半,还不解馋,又多吃了一个. 第二天早上又将剩下的桃子吃了一半,还是不过瘾,又多吃了一个.以后每天都吃前一天剩下的一半再加一个 ...

  9. c语言错误 xef代表什么,单片机C语言代码手册 含100多个经典C程序

    1 单片机单片机 C 语言代码手册语言代码手册 1 LED 灯灯 点亮一个点亮一个 LED include void main while 1 P0 0 x01 P2 0 x7d 流水灯闪烁流水灯闪烁 ...

最新文章

  1. SIM800C SIM卡唯一标识符ICCID
  2. 使用python模拟Simple方式连接ldap
  3. 数值优化(Numerical Optimization)学习系列-文件夹
  4. Re-installation failed due to different application signatures.
  5. DCMTK:DSRListOfItems类的测试程序
  6. C++开发中类的定制
  7. super在python3和python2_Python扩展与 – 使用super()python 3 vs python 2
  8. 数据结构:点之间的最短距离--Floyd算法
  9. git 删除本地和远程分支_如何在本地和远程删除Git分支
  10. Codeforces Beta Round #1 B. Spreadsheets 模拟
  11. 如何解决css样式表在不同浏览器中显示效果不同的问题
  12. uni-app - 最新自定义连接安卓模拟器配置到 MUMU 模拟器运行程序(将uni-app运行到MuMu模拟器上进行调试,配置电脑手机模拟器 MUMU 运行测试)
  13. foxpro导入 mysql_无法导入、导出或链接到 FoxPro 数据库 | Microsoft Docs
  14. 人工智能的主要研究方向都有哪些
  15. 转行的35岁程序员们
  16. 42、PCF8563电子时钟实验
  17. Activiti7工作流引擎-环境搭建及创建数据库
  18. 我的第一个网站——MISDream(四):网站的发布与云服务器的搭建
  19. 【电网规划】基于智能算法的电力系统最优规划matlab仿真
  20. 函数周期表丨值丨数学函数系列

热门文章

  1. max open files mysql_arch-linux – MySQL max_open_files超过1024
  2. 基于51单片机甲醛浓度检测系统仿真设计(仿真+程序+原理图+报告+讲解)
  3. java字符串abcde_关于字符串的实验(abc和abcde的起始地址是否一样?)
  4. MYSQL mediumText大字段存储错误问题排查
  5. 记忆化搜索c语言,2017C语言混搭 标题:10.图形排版 搜索+记忆化
  6. sigslot.h 中文文档
  7. 击败GANs的新生成式模型:score-based model(diffusion model)原理、网络结构、应用、代码、实验、展望
  8. apolloconfig分布式部署
  9. 医院营销 民营医院企划的内涵及形式与内容
  10. c语言字母意义,%C是什么意思? c语言中?和:是什么意思