偏微分方程在物理学中的完美应用——热方程,推导和示例 2024-08-05 03:00:12 偏微分方程是一个将具有一个以上变量的函数与其偏导数联系起来的方程。为了引入偏微分方程,我们要解决一个简单的问题:模拟薄金属棒内的温度作为位置和时间的函数。在此过程中,我们将从物理原理推导出一维热方程,并求解一些简单的条件:在这个方程中,温度T是位置x和时间T的函数,k、ρ和c分别是金属的热导率、密度和比热容,k/ρc称为扩散系数。 物理过程 我们想要研究,随着时间的增加,热量如何在长为L的金属棒中传导的。金属棒的一端在x=0处,另一端在x=L处。金属棒的长度远大于它的截面半径,所以我们可以把热传导看成是x和t的函数。假设金属棒的比热容是已知的,如果我们能找到温度T(x, t)的函数,我们就能知道热量是如何扩散的。假设棒沿其长度方向是绝热的,因此它只能通过两端吸收或散发热量。这意味着温度分布只取决于以下三个因素: 初始温度分布情况T(x,0)。这叫做初始条件。 金属棒两端的温度,T(0, t)和T(L, t)这些叫做边界条件。 热量在金属棒内由一点传递到另一点的规律。热方程是这种物理定律的数学表示。 对于一组特定的初始和边界条件,求解偏微分方程的问题被称为初始边值问题(IBVP)。在本文中,我们将求解的热方程的初始边值为T(0,t)=T(L,t)=0°C。这些叫做齐次边界条件。 热方程的推导 热方程可以从能量守恒导出:金属杆上某一点储存的热量的时间变化率等于进入该点的净热量流量。这个过程显然符合连续性方程。如果Q是各点处的热量,V是热量流动的矢量场,则:根据热力学第二定律,如果两个相同的物体进行热接触,其中一个比另一个热,那么热量必然以与温度差成比例的速度从较热的物体流向较冷的物体。因此,V与温度的负梯度成正比,所以V=-k∇T,其中k为金属的导热系数。在一维中,它简化为V=(-k∂T/∂x)x,其中x是+x方向的单位向量。Q=ρcT,代入V和Q的表达式,得到热方程: 解热方程 在我们进一步讨论之前,我们需要证明对于任何有物理意义的初始和边界条件,热方程必须存在一个唯一的解。对此的正式证明超出了本文的范围,因此我们将使用一个经验论证。热力学定律告诉我们,无论一开始金属棒的温度分布是怎样的,系统必须经历一个过程,使金属棒达到热平衡,我们在前面讲过这个过程必须服从热方程,因此,对于有物理意义的初始和边界条件,热方程的解是存在的。此外,经典物理学的基本假设之一是,相同的实验条件必然会导致相同的结果,因此,金属棒进入热平衡的特定方式,由初始条件和边界条件所唯一规定。这意味着,对于热方程,如果f(x,t)和g(x,t)是两个不同的函数且满足相同的IBVP,那么f和g有相同的形式。此外,热方程是线性的,因此如果f和g是解,α和β是任何实数,那么αf+βg也是一个解。所以我们可以得出结论,解是相同形式的函数的线性组合。考虑下面的函数,我们可以通过试错来推测:其中n是大于0的正整数。该函数满足热方程:这个函数也满足边界条件,因为sin(0)=sin(nπ)=0。因此通解为:如果我们能找到系数A_n,使这个通解满足初始条件,问题就解决了。也就是说,我们需要找到一个A_n,这样:这叫做初始条件下的傅里叶正弦级数展开式。系数A_n叫做傅里叶系数。 计算傅里叶系数 初始条件T(x,0)是区间[0,L]上的分段连续函数,且在边界处为零。结果证明具有这些性质的函数集合是加法和标量乘法下的向量空间。我们称这个向量空间为:这个向量空间有一个内积。对于f,g∈𝕊ᴸ(加法和标量乘法下的向量空间),一个可能的内积为:我们可以通过使用单位向量的点积将其投影到轴上来找到几何向量的分量,单位向量构成了ℝⁿ的基。同样,如果我们能为𝕊ᴸ找到一个基,那么我们可以将任何f∈𝕊ᴸ投影到基函数上,以便将f表示为基函数的线性组合。 求v的分量 对于整数m,n>0,函数sin(nπx/L)是标准正交的:因此,我们可以将任意函数f∈𝕊ᴸ表示为基集中函数的线性组合:线性组合的系数由欧拉积分给出:作为演示,让我们找到一个单位锯齿脉冲的傅立叶系数:显然,saw(x)∈𝕊¹,因此:傅里叶系数为:因此,锯齿波的傅里叶级数展开为:下面的动画展示了,随着sin项的增加,傅里叶级数如何接近锯齿状(波形)的。x=1附近的误差称为吉布斯现象。吉布斯现象是一种不可避免的误差,它使不连续函数的傅里叶级数将不连续时的函数值高估约9%。吉布斯现象永远不能完全消除,但当傅里叶级数中的项数接近无穷时,误差收敛到完全局限于不连续点。例如,如果在锯齿形的傅里叶级数展开式中包含无限项,我们会发现当 0≤x<1 时,级数将完全等于x,而在x=1时,级数的值约为1.09。这告诉我们,求解热方程的齐次IBVP等于使用欧拉积分来求傅里叶系数:例:金属棒初始温度是均匀的假设一个绝热的,一米长的金属棒扩散系数为k/ρc=0.1m²/s(不现实,为了方便作图),最初温度为100°C,在温度为0℃、t=0时夹紧冷却元件。初始条件和边界条件为:傅里叶系数是:让我们来验证这个匹配初始和边界条件的傅里叶级数:所以,解是:现在让我们画出解:作为一个三维图:例:温度“尖峰”现在假设金属棒每个地方的初始温度都是0°C,除了中间10厘米的温度是100°C。这次的扩散系数是0.0075m²/s。计算傅里叶系数最简单的方法是把它们转换成更一般的形式:所以解是:虽然金属棒最终会达到热平衡状态,但是5秒后温度下降的非常缓慢,所以动画值展示5秒。下面是3D图:例:随机热分布这一次,这根金属棒,我们假设是铜做的,扩散系数为1.11×10⁻⁴m²/s,它的热量分布是随机的(以10厘米为一段),如下表:在这种情况下,最好是通过数值积分来求傅里叶系数而不是试图找到一个封闭式的表达式。下面是一个解的动画:240秒后,温度的变化非常缓慢。最有趣的行为发生在最初的60秒: 结束语 这就是第一部分的内容。现在你知道了如何解最简单的情况下的热方程,你可以使用热方程来分析更有趣的问题。 赞 (0) 相关推荐 如何解常系数的线性常微分方程|高数札记(1) 可乐数学按:我们将陆续刊登张博士的高数札记系列文章,其中记录了他教高等数学的一些感悟.第一篇中他介绍了求解一.二阶常系数线性常微分方程的方法. 高数札记1 @21世纪学习数学的人 这学期第一次教高等数 ... 方程的常点与奇点 方程的常点与奇点 微分方程 "线性"是指函数y及其n阶导数的幂都为1: "常系数"是指函数y及其n阶导数前的系数都为常数: "微分方程"即以自变量x,函数y及其n阶导数 ... 傅里叶变换的应用Ⅰ——热传导与扩散 事实上,傅里叶变换本来就是傅里叶在解决热传导问题时提出的理论.我们也从这个问题入手来讲解傅里叶变换的应用. 1 傅里叶热传导定律 热传导是一种非常常见的现象.当系统内温度分布不均时,热量就会自发地从高 ... 原来鹿晗心目中的完美女神是她?热巴知道后要吃醋啊 本文由明星粉丝团作者柴扉原创,未经允许不得转载 自从迪丽热巴参加完<奔跑吧兄弟>之后,如今每次提到鹿晗就不得不说迪丽热巴了.不少网友看过这一季的跑男之后,大家就喜欢上了陆地夫妇.总觉得他们 ... 物理学中最深刻的联系之一——对称性与守恒定律之间的关系|牛顿|动量_网易订阅 2020-10-15 18:45:02 来源: 老胡说科学 守恒定律说,孤立系统的某些性质不会随着系统的演化而改变.作为诺特定理的结果,守恒定律"与底层物理中的对称性有关".换句话 ... 不完美中的完美 杨强劲 说到课改,我们总会想到诸多的理论概念和创新词汇,也容易将之放在一个宏观的背景下,给出学术味十足的观照.实际上对于基础教育学校一线的管理者和教师来说,如何将课改理念有效地落实到课堂教学中,从微观 ... 译研动态|“大使荐书”创新合作模式,中尼文学作品引热议 2017年12月21日,由中国文化译研网(CCTSS)与尼泊尔驻华使馆联合主办的"文明的交响·中尼文化沙龙"在京成功举行.文化部.中华文化发展委员会.统战部.国家宗教局等部委,阿富 ... 明暗双色修容棒,高光加隐影,让你即刻打造轮廓妆容,立体V字妆容,轻松搞定,打造出理想中的完美脸部线条 修容棒是近期最火热的修容产品,质地细腻,不容易晕染,而且上妆持久,是很多专业化妆师的不二选择.越来越多的美女们偏爱修容棒,但是很多化妆小白还没有学会怎么使用修容棒,那么跟着小编带来的修容棒教程,一起与 ... 《伤寒论》中用于治疗血热互结的蓄血证,是... <伤寒论>中用于治疗血热互结的蓄血证,是泄热逐瘀代表方. 糖尿病常出现多饮.多食.多尿,大便干燥,便秘,为胃肠燥热.灼伤阴津所致. 又多见口唇紫暗,舌质暗红,边有瘀斑.舌下静脉青紫,脉沉涩 ... 词源趣谈:focus——从家中的“火炉”到物理学中的“焦点”2 英语单词focus直接来自拉丁语,由词根foc-加名词后缀-us组成,本意是"火炉.壁炉". 17世纪初,著名物理学家开普勒最早用focus一词来表示光线的汇聚点,即所谓的&quo ... 【左岸风文学】翠儿|依然是完美中的完美 外一首 依然是完美中的完美 外一首 作者:翠儿 盏无水痕为妙,谷雨茶滋味鲜活 一芽一嫩叶,一芽两嫩叶,雀舌,旗枪 心中有佛见佛,有魔见魔 浅浅而尝,不似液体,而像饮下一朵云 暖而轻柔,一路袅娜 温润调畅.清 ... 《病机十九条》中属“火、热”的病机和病症该如何理解? 在<素问·至真要大论>中,论述属热的有四条,分别为"诸病有声,鼓之如鼓,皆属于热","诸腹胀大,皆属于热","诸转反戾,水液混浊,皆属于热 ...