抛物形扩散方程的有限差分法及数值实例概念.docVIP

偏微分方程数值解 所在学院： 数学与统计学院 课题名称：抛物形方程的有限差分法及数值实例 抛物形方程的有限差分法及数值实例 1.1.1 其中是常数是给定的连续函数按照初边值条件的不同给法可将，满足方程 1.1.1 和初始条件： , 1.1.2 第二，初边值问题（也称混合问题）：求足够光滑的函数，满足方程 1.1.1 和初始条件： , 1.1.3 及边值条件 ， 1.1.4 假定和在相应的区域光滑，并且于，两点满足相容条件,则上述问题有唯一的充分光滑的解。 1.2抛物线扩散方程的求解 下面考虑如下热传导方程 1.2.1 其中，，, 常数 是扩散系数为空间步长，为时间步长，其中，是自然数，用两族平行直线, 和, 将矩形域分割成矩形网格。其中 表示网格节点；表示网格内点（位于开矩形中的网格节点）的集合；表示位于闭矩形中的网格节点的集合；表示-网格边界点的集合。 表示定义在网点处的待求近似解，，。 现在对方程进行差分近似： 向前差分格式 1.2.2 ， 0 1.2.3 计算后得： 1.2.4 其中，，，。 显然这是一个四点显示格式每一层各个节点上的值是通过一个方程组求解到的 1.2.5 若记 ，， 则显格式 1.2.4 可写成向量形式 1.2.6 其中 而对于向前差分格式，当网比时稳定，当时不稳定。这就意味着给定空间步长以后，时间步长必须足够小，才能保证稳定。 1.3抛