谱元方法课程.ppt

* Tel e-mail：glqin@mail.xjtu.edu.cn * 总体区域类 主要数据及方法包含区域内单元数、总节点数、节点类、标准单元类、总体矩阵等。 总体区域类定义， class CDomain { protected: //整体区域类的属性 区域的属性，单元及网格属性，单元类，标准单元类； 总体矩阵，对角元位置，边界节点序号； public: //整体区域类方法 单元构造，标准单元类的构造，节点类构造； 元素矩阵的形成，属性设置及取值； } * Tel e-mail：glqin@mail.xjtu.edu.cn * 6.2 谱元方法程序设计中的几个问题 插值函数积分的存储 节点的编码及总体矩阵的形成 对称带状矩阵的一维存储 元素的编码 影响元素和影响节点 对角元的位置及一维数组的长度 总刚度矩阵 第一类边界条件的处理 * Tel e-mail：glqin@mail.xjtu.edu.cn * 插值函数积分的存储 在形成元素矩阵时，要计算如下积分， 可以看出， 在(i,j,k)变化时，不必占用(Nx+1)3个存储单元，这样既节约了存储空间又节省了计算时间。Bijk所占的空间可由下式计算， Bijk的检索按如下方式进行，先将i,j,k按由大到小排列，设ijk，则Bijk在数组中的位置为， * Tel e-mail：glqin@mail.xjtu.edu.cn * 6.3 Navier-Stokes方程和能量方程的求解过程 构造标准单元stdElement，求出对各单元不变的量。 构造求解域pGeometry，对求解域进行网格划分，同时构造节点类pDot[NG]，构造各求解单元pElement[Ne]，对节点类中U、V、Φ赋初值，赋零值，Φ在热边界上赋值1，在冷边界上赋0，其余两边为第二类边界条件，在形成单元矩阵时自动满足。 利用Runge-Kutta法求解 、 、 ，将求解结果存于 、 、 解压力Poisson方程，第二类边界条件，解出压力后更新 、 、 求解关于 、 、 的Helmholtz方程，边界条件第一或第二类。解出Helmholtz方程后更新 、 、 ，即完成了一个时间步的求解。 * Tel e-mail：glqin@mail.xjtu.edu.cn * 7 计算实例 常微分方程求解 Helmholtz和Poisson的求解 简单边界Helmholtz方程求解 复杂边界Poisson的求解 移动顶盖方腔驱动非定常流动 封闭方腔中的自然对流 * Tel e-mail：glqin@mail.xjtu.edu.cn * 7.1 常微分方程求解 【例1】用谱方法求解下面方程 该问题的精确解的表达式为： 将求解域[0,-1]等分为10个求解单元,单元内使用5阶Chebyshev多项式进行谱逼近。数值计算结果与精确解的比较如下,比较表明谱元方法具有相当高的精度。 * Tel e-mail：glqin@mail.xjtu.edu.cn * 例1 数值计算结果与精确解的比较如下 序号 x 数值解 解析解 绝对误差 0 1 2 3 4 5 6 7 8 9 10 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 -3.83295107661418e-018 0.155924631208591 0.308292981155188 0.453584301764615 0.588348556299474 0.709240890986803 0.813055055691316 0.896755439706894 0.95750740257533 0.99270559687679 1 0 0.155924631208596 0.308292981155196 0.453584301764625 0.588348556299488 0.709240890986824 0.813055055691339 0.896755439706917 0.957507402575354 0.992705596876806 1 -3.83295107661418e-018 -4.88498130835069e-015 -8.04911692853239e-015 -9.99200722162641e-015 -1.45439