热传导方程的第一边值问题的计算机模拟大物DTP.docVIP

热传导方程的第一边值问题的计算机模拟 物理学3班 罗杰昌 20082304029 摘要：运用分离变量法求解一维的热传导问题，得出热传导函数，再使用科学计算与模拟平台模拟出一金属杆在一定边值条件和初始条件下温度变化情况，从而体现出计算机在解决实际问题上的巨大作用。 关键词：一维热传导，热传导方程，温度分布，科学计算与模拟平台 1.前言 在热学中,热传导问题是一个基础性的课题,对其进行研究有助于理解热传导的基本规律,澄清基本物理概念。在数学物理方法中,热传导问题是由热传导方程来描述。 热传导方程的求解因边界条件不同而难度各异。有的问题没有严格的解析解,只能求得数值解;有的则存在严格的解析解,其解的形式可能是积分形式，也可能是级数形式。从这些解中很难准确知道温度、热流密度等物理量在空间的分布和随时间的演化规律。本文利用科学计算与模拟平台帮助求解热传导方程的第一边值问题，从而求出动态定解。 2.问题描述 均匀金属杆可看作一维传热系统. 如图1 所示,杆的左端作为坐标原点,则右端的坐标为L ,杆上各点的温度T( x , t) 满足偏微分方程（1）其中,分别为金属杆的热导率、比热容和密度。 金属杆的初始温度为，其两端与0oC的恒温热源接触，杆的侧面与外界环境绝热。现在讨论该金属杆的温度场分布、热流密度分布及其演化情况，其模型为方程 ① 3.分析 对方程①分离变量，并引入参数λ得 ， . 容易看出本征值，对应本证函数为，其中n=1,2,3,….而 ，形式解为, 其中.当具有一、二阶连续导数，三阶导数逐段连续，且时，上述级数一致收敛，形式解即为正规解。 为具体求解，取 计算中，k=1 接着用科学计算与模拟平台模拟出动态定解如下图： n=1 n=3 n=5 4.结论 综上所述，本文对第一类边界条件下的金属杆的一维热传导问题进行了可视化分析，用科学计算与模拟平台分别令n=1，n=3，n=5进行作图，可看出到最终金属杆的温度都趋向0这个定值，使得原来公式化、数值化的结果赋予了生动形象的描述，动态的显示了温度变化的情况。