1、 2010级数学与应用数学和信息与计算科学专业偏微分方程数值解上机实验实验题目 利用有限元方法和有限差分方法求解偏微分方程 完成日期 2013 年 6 月 17 日 姓 名 班 级 学 号 西北工业大学理学院应用数学系成 绩1偏微分方程数值解上机实验报告实验地点：数学系机房实验时间：第 1415 周，周一、四下午 5、6 节实验分数：占期末考试成绩的 20%一、实验目的及意义掌握有限元方法和有限差分方法的程序实现；学会选择合适的有限差分格式求解一维非线性对流占优的非定常对流扩散问题；学会使用三角线性元、四边形线性元、三角形二次元、四边形二次元的有限元方法求解二维椭圆方程边值问题，并对计算结果进
2、行收敛性分析；尝试采用有限元方法或有限差分方法实现二维初边值抛物型方程的大规模数值求解。通过实验可以提高学生的动手能力，加深学生对算法的理解。二、实验要求在下列给出的三个问题中，最少选择两个问题进行编程实现(问题 2 必选)。要求给出格式的推导过程、算法流程、实现程序、选取的网格参数、以表格或图形的方式给出计算结果、对计算结果进行分析、最后对实验进行总结和讨论。三、实验报告提交实验报告以电子版的形式提交，要求附相应程序。报告最后提交时间为2013 年 6 月 17 日，若在此时间之前不能提交，则实验成绩按零分记。若发现2有抄袭现象，则以作弊论处。问题 1：用差分法求解 Burger 方程： 2
3、,0,(,1)0,(,)sin();01,.uuxttxtt取 =1.0,.,.,问题 2：用三角线性元、四边形线性元、三角形二次元和四边形二次元的有限元方法求解方程： 28cos()in(2),(0,1)(,0),1;(i().uxyxxy取 给出 和 两种范数度量下的计算误差，比较四=1/4, 8/6, 32/4,h2L1H种方法的计算精度，给出相应的数值收敛率，并对所得结果进行分析解释。问题 3：选用合适的数值方法求解方程： 2122(4),()0,1()0,),(0;(sin()cos).yyuxuxtytttuxx取 ，求 时的数值解，请绘制数值解图形，并给出点1/256h0.1,5.t3、 、 处的数值解。1279(,)5631(,)2563(,)2564上机实验（一）一、实验内容二、实验原理三、算法流程四、计算结果及分析五、总结和讨论附录：程序源代码5上机实验（二）一、实验内容二、实验原理三、算法流程四、计算结果及分析五、总结和讨论附录：程序源代码6上机实验（三）一、实验内容二、实验原理三、算法流程四、计算结果及分析五、总结和讨论附录：程序源代码
