作业一传热与流体流的数值方法.docVIP

流动与传热的数值方法作业（一） 姓名：徐世杰 学号：120351 题目1： 用Galerkin 方法求以下方程在内部节点的离散化方程。 取线性插值函数，，其中节点间距是均匀的。 题目2：考虑 用控制容积有限差分方法做出内部节点和边界节点的离散化方程； 写出代数方程组的迭代求解程序； 研究空间步长对数值精度和收敛性的影响。 题目一 解：Galerkin 方法就是将对应某个点上的插值函数作为权函数。Galerkin 方法是有限元方法。 可知有： ， i=2，….,n-1 按照习惯，上述积分写成： 可以推得： 由弱解变换可以得 可以得 ，i=2，…,n-1 上式继续推导有： 其中， 化简可以得： i=2，….,n-1 题目二 ①用控制容积有限差分方法做出内部节点和边界节点的离散化方程： 首先进行离散化，先确定节点，再确定控制容积。将0-1划分为N段，共N+1个节点，N个控制容积，其中。 对原方程建立差分方程，内部节点： 边界节点 边界节点 组成代数方程组： 写成矩阵方程组： ②写出代数方程组的迭代求解程序: 用Matlab编写如下