热传导方程

✍ dations ◷ 2025-02-24 10:30:44 #扩散,热传导,抛物型偏微分方程,传热

热传导方程(或称热方程)是一个重要的偏微分方程,它描述一个区域内的温度如何随时间变化。

热传导在三维的各向同性介质里的传播可用以下方程表达:

其中:

热方程是傅里叶冷却律的一个推论(详见条目热传导)。

如果考虑的介质不是整个空间,则为了得到方程的唯一解,必须指定u的边界条件。如果介质是整个空间,为了得到唯一性,必须假定解的增长速度有个指数型的上界,此假定吻合实验结果。

热方程的解具有将初始温度平滑化的特质,这代表热从高温处向低温处传播。一般而言,许多不同的初始状态会趋向同一个稳态(热平衡)。因此我们很难从现存的热分布反解初始状态,即使对极短的时间间隔也一样。

热方程也是抛物线偏微分方程最简单的例子。

利用拉普拉斯算子,热方程可推广为下述形式

其中的 Δ {\displaystyle \Delta } (中译:解析热学)给出。先考虑只有一个空间变数的热方程,这可以当作棍子的热传导之模型。方程如下:

其中 = (, )是和的双变数函数。

假设下述初始条件

其中函数是给定的。再配合下述边界条件

让我们试着找一个非恒等于零的解,使之满足边界条件(3)并具备以下形式:

这套技术称作分离变数法。现在将代回方程(1),

由于等式右边只依赖,而左边只依赖,两边都等于某个常数− λ,于是:

以下将证明(6)没有λ ≤ 0的解:

假设λ < 0,则存在实数、使得

从(3)得到

于是有 = 0 = ,这蕴含恒等于零。

假设λ = 0,则存在实数、使得

仿上述办法可从等式(3)推出恒等于零。

因此必然有λ > 0,此时存在实数、、使得

从等式(3)可知 = 0,因此存在正整数使得

由此得到热方程形如(4)的解。

一般而言,满足(1)与(3)的解相加后仍是满足(1)与(3)的解。事实上可以证明满足(1)、(2)、(3)的解由下述公式给出:

其中

上面采用的方法可以推广到许多不同方程。想法是:在适当的函数空间上,算子 u u x x {\displaystyle u\mapsto u_{xx}} = ,以下函数序列

是Δ的特征矢量。诚然:

此外,任何满足边界条件(0)=()=0的Δ的特征矢量都是某个。令L2(0, )表 上全体平方可积函数的矢量空间。这些函数构成L2(0, )的一组正交归一基。更明白地说:

最后,序列{} ∈ N张出L2(0, )的一个稠密的线性子空间。这就表明我们实际上已将算子Δ 对角化。

一般而言,热传导的研究奠基于以下几个原理。首先注意到热流是能量流的一种形式,因此可以谈论单位时间内流进空间中一块区域的热量。

因此单位时间内进入的热流量也由以下的面积分给出

其中n(x)是在点的向外单位法矢量。

利用格林定理可将之前的面积分转成一个体积分

将以上所有等式合并,便获得支配热流的一般公式。

注记:

在粒子扩散的模型中,我们考虑的方程涉及

不同情况下的方程:

或者

与都是位置与时间的函数。是扩散系数,它控制扩散速度,通常以米/秒为单位。

如果扩散系数依赖于浓度(或第二种情况下的概率密度),则我们得到非线性扩散方程。

单一粒子在粒子扩散方程下的随机轨迹是个布朗运动。

如果一个粒子在时间 t = 0 {\displaystyle t=0} 时,上述 P ( R , t ) {\displaystyle P({\vec {R}},t)}

粒子扩散方程首先由Adolf Fick于1855年导得。

格林函数是扩散方程在粒子位置已知时的解(数学家称之为扩散方程的基本解)。当粒子初始位置在原点 0 {\displaystyle {\vec {0}}} );根据扩散方程对平移的对称性,对一般的已知初始位置 R 0 {\displaystyle {\vec {R}}^{0}} 时有一大群粒子,根据浓度分布的初始值 c ( R , 0 ) {\displaystyle c({\vec {R}},0)} ,浓度分布变为:

在粒子扩散的情形,我们可以将狄拉克δ函数对应的初始条件理解为粒子落在一个已知位置。一般而言,任何扩散过程的解都有这种表法,包括热传导或动量的扩散;后者关系到流体的黏性现象。

以下以简写BC代表边界条件,IC代表初始条件。

(可能的问题:根据上解,u(0)=0)

热方程在许多现象的数学模型中出现,而且常在金融数学中作为期权的模型出现。著名的布莱克-斯科尔斯模型中的差分方程可以转成热方程,并从此导出较简单的解。许多简单期权的延伸模型没有解析解,因此必须以数值方法计算模型给出的定价。热方程可以用Crank-Nicolson法有效地求数值解,此方法也可用于许多无解析解的模型(详见文献Wilmott,1995)。

热方程在流形上的推广是处理阿蒂亚-辛格指标定理的主要工具之一,由此也导向热方程在黎曼几何中有许多应用。

相关

  • 移植排斥移植排斥(英语:transplant rejection)是器官移植后的器官并不被受移植者身体接受的情况。一般来说这是因为免疫系统将移植器官视为异物,如同攻击病毒或细菌一样攻击移植器官所引
  • 木骨架木骨架是一种将木柴利用雌雄榫连结在一起 ,以构成建筑结构的的建筑技术。这一技术的历史可以追溯到几千年前,在世界各地的建筑中都有应用。 欧洲地区主要应用在12世纪至19世纪
  • 沈尹默沈尹默(1883年6月11日-1971年6月1日),原名君默,字中,又字秋明,号君墨,别名匏瓜庵主人,浙江吴兴(今浙江省湖州市吴兴区)人,出生于陕西省汉阴厅(今陕西省汉阴县),中国书法家,早期尝试派诗人。
  • 生命线 (热线)生命线(英文:Lifeline)是一个提供24小时免费的自杀危机防治专线。生命线的志工能透过电话或视讯提供自杀防治服务、心理支持、以及情绪安抚。
  • 国民议会 (奥地利) 政治主题国民议会(德语:Nationalrat)是奥地利议会的下议院,根据宪法,国民议会和联邦议会(上议院)是平等的,但国民议会拥有更大权力,包括可以选举奥地利总理和组成联合政府,联合政府需
  • 雷公连雷公连(学名:)为天南星科雷公连属下的一个种,是中国的特有植物。分布于中国大陆的广西、贵州、四川、湖北西部、湖南西部和南部、云南东南部等地,生长于海拔550米至1,100米的地区
  • 中国野生动物园调查报告《中国野生动物园调查报告》是八位志愿者在调查了北京八达岭野生动物世界、大兴野生动物园、上海野生动物园、广州番禺香江野生动物世界、深圳野生动物园、北京动物园、武汉
  • 斐迪南·菲利普 (奥尔良公爵)奥尔良的斐迪南·菲利普王子(Prince Ferdinand Philippe of Orléans;1810年9月3日-1842年7月13日),全名奥尔良的斐迪南·菲利普·路易·夏尔·埃里克·罗萨利诺(法语:Ferdinand Ph
  • 白雪石白雪石(1915年6月-2011年4月22日),原名白增锐,斋名何须斋,北京人,中华人民共和国画家,第十一届全国政协委员。担任清华大学美术学院咨询委员会委员。2008年,当选第十一届全国政协委员
  • 金彩虹大奖金彩虹大奖(SPECTRUM RADIO GOLDEN RAINBOW AWARD)乃伦敦国际广播电台举办之听众选举活动。每年一月由听众选出多个电影、电视及音乐大奖,并设有金彩虹至尊大奖,表扬过去一年